help="return max eigenvalue")
    parser.add_argument("-o", "--output", type=str, required=False,
                        help="file to write output to")
    parser.add_argument("-c", "--conditionnumber", type=float, required=False,
                        help="condition number to determine cutoff")
    args = parser.parse_args()

    if args.verbose:
        logging.basicConfig(format='%(levelname)s: %(message)s', level=logging.DEBUG)
        logging.debug("Verbose debuging mode activated")
    else:
        logging.basicConfig(format='%(levelname)s: %(message)s', level=logging.INFO)

    if not args.operators:
        logging.info("Operators not specified, attempting to automagically determine")
        ops = determine_operators.matching_operators(args.input_dir, args.filewild)
        logging.info("Operators automagically found to be {}".format(",".join(ops)))
        if not ops:
            logging.error("Error: no operators found")
            parser.print_help()
            parser.exit()
        args.operators = ops
    else:
        ops = readops(args.operators)
        logging.debug("found operators: {}".format(','.join(ops)))
        args.operators = ops

    cor_matrix = {}
    for snk in args.operators:
        for src in args.operators:
            filename = args.input_dir + args.filewild.format(snk, src)
Exemple #2
0
                    help="Period in time direction (not required for all functions)")
parser.add_argument("--periodic", action="store_true", help="use periodic effective mass functions")
parser.add_argument("--prune", action="store_true", help="don't plot stuff where the correlator value is too close to zero (with error)")
parser.add_argument("-c", "--configs", type=int, required=False, help="specify the configs to be used\n")
parser.add_argument("-t", "--times", required=False, type=int, help="specify the times to be used\n")


args = parser.parse_args()

if args.fit:
    raise DeprecationWarning("Fit no long works (plotting fits was reworked)")
    funct = functions[args.function](Nt=args.period)

if not args.operators:
    print "Operators not specified, attempting to automagically determine"
    ops = determine_operators.matching_operators(args.input_dir, args.format)
    print ops
    if not ops:
        print "Error: no operators found"
        parser.print_help()
        parser.exit()
    args.operators = ops

cor_template = args.format
if (not args.make_from_operators) and (not args.format_vev and not args.no_vev):
    print "Error: must specify vev format or no-vev (-nv)"
    parser.print_help()
    print "\nError: must specify vev format or no-vev (-nv)"
    parser.exit()
vev_template = args.format_vev