Exemplo n.º 1
0
    backup_interval = 100

if "-loginterval" in sys.argv:
    log_interval = int(sys.argv[sys.argv.index("-loginterval") + 1])
else:
    log_interval = 100
    
mode = "fu"
if "-mode" in sys.argv:
    mode = sys.argv[sys.argv.index("-mode") + 1]

name_suffix = ""
if "-namesuffix" in sys.argv:
    name_suffix = sys.argv[sys.argv.index("-namesuffix") + 1]

globallog.write("tfi_gs_2d.py, D={:d}, chi={:d}, h={:f}, tau={:.0e}, iterations={:d}, trotter order {:d}{:s}{:s}\n".format(D, chi, h, tau, maxiterations, 2 if trotter_second_order else 1, ", ffu" if fast_full_update else "", "" if name_suffix == "" else ", name_suffix=" + name_suffix))

basepath = "output_tfi/"
if mode == "su":
    basepath = "output_tfi_su/"

if name_suffix != "":
    name_suffix = "_" + name_suffix

if not output_to_terminal:
    f = open(basepath + "log_tfi_gs_2d_D={:d}_chi={:d}_h={:f}_tau={:.0e}{:s}.txt".format(D, chi, h, tau, name_suffix), "a")
    sys.stdout = f
    sys.stderr = f

if os.path.isfile(basepath + statefile):
    a, nns = peps.load(basepath + statefile)
Exemplo n.º 2
0
import globallog
from time import time
from polish import polish

t0 = time()

chi = int(sys.argv[1])
h = float(sys.argv[2])
statefile = sys.argv[3]
output_to_terminal = "-writehere" in sys.argv

num_workers = 1
if "-workers" in sys.argv:
    num_workers = int(sys.argv[sys.argv.index("-workers") + 1])

globallog.write('tfi_gs_polish.py, chi={0:d}, h={1:f}, statefile="{2:s}"\n'.format(chi, h, statefile))

basepath_in = "output_tfi/"
basepath_out = "output_tfi_polish/"

if not output_to_terminal:
    logfile = open(basepath_out + statefile[:-5] + ".log", "a")
    sys.stdout = logfile
    sys.stderr = logfile


def test_fct(a, A):
    n = len(a)
    X = map(lambda b: peps.make_double_layer(b, o=gates.sigmax), a)
    Z = map(lambda b: peps.make_double_layer(b, o=gates.sigmaz), a)
    AT = map(lambda B: B.transpose([1, 2, 3, 0]), A)
Exemplo n.º 3
0
# from polish import polish
import polish

t0 = time()

chi = int(sys.argv[1])
J = int(sys.argv[2])
h = float(sys.argv[3])
statefile = sys.argv[4]

num_workers = 1
if "-workers" in sys.argv:
    num_workers = int(sys.argv[sys.argv.index("-workers") + 1])

globallog.write(
    'heisenberg_gs_polish.py, chi={0:d}, J={1:d}, h={2:f}, statefile="{3:s}"\n'.format(chi, J, h, statefile)
)

basepath = "output_heisenberg_polish/"

if not ("-writehere" in sys.argv):
    logfile = open(basepath + statefile[statefile.rfind("/") + 1 : -5] + ".log", "a")
    sys.stdout = logfile
    sys.stderr = logfile


def test_fct(a, A):
    n = len(a)
    Z = map(lambda b: peps.make_double_layer(b, o=gates.sigmaz), a)
    AT = map(lambda B: B.transpose([1, 2, 3, 0]), A)
    ZT = map(lambda B: B.transpose([1, 2, 3, 0]), Z)