main.py

# %%
import jackknife
import adga
import numpy as np
import sys


# %%
if len(sys.argv) == 4:
    adga_root = sys.argv[1]
    adga_config = sys.argv[2]
    two_particle_file = sys.argv[3]
elif len(sys.argv) < 4:
    print("Command line argument missing.")
    print("""Please specify the ADGA root directory, the configuration file and
          the 2-particle Green's function file""")
    sys.exit(1)
else:
    print("Too many command line arguments.")
    print("""Please specify only the ADGA root directory, the configuration
          file and the 2-particle Green's function file""")
    sys.exit(2)

abinitio_dga = adga.Adga(adga_root, adga_config, two_particle_file)
x_generator = abinitio_dga.get_worm_sample_generator()
n = abinitio_dga.n_worm_samples
f = abinitio_dga.__call__

# %%
# Do the Jackknife estimation
jackknife.do_jackknife_estimation(x_generator, n, f)