# Run TEBD
pepsl = PEPS(Nx,
Ny,
d,
D[0],
chi[0],
fname=fnamel,
fdir=savedir,
norm_tol=0.5,
norm_bs_upper=3.,
norm_bs_lower=0.)
# Loop over all optimizaton parameters
for ind in range(len(D)):
# Update PEPS Parameters
pepsl.D = D[ind]
pepsl.chi = chi[ind]
pepsl.fname = prepend + "Nx{}_Ny{}_sx{}_sy{}_D{}_chi{}_run_left".format(
Nx, Ny, sxind, syind, D[ind], chi[ind])
# --------------------------------------------------------------------
# Calculate left eigenstate
for i in range(5):
try:
Efl, pepsl = run_tebd(Nx,
Ny,
d,
opsl,
peps=pepsl,
D=D[ind],
chi=chi[ind],
norm_bs_lower=0.)
pepsl = PEPS(Nx,
Ny,
d,
D[0],
chi[0],
fname=fnamel,
fdir=savedir,
norm_tol=0.5,
norm_bs_upper=3.,
norm_bs_lower=0.)
# Loop over all optimizaton parameters
for ind in range(len(D)):
# Update PEPS Parameters
peps.D = D[ind]
peps.chi = chi[ind]
peps.fname = "Nx{}_Ny{}_sx{}_sy{}_D{}_chi{}_run_right".format(
Nx, Ny, sxind, syind, D[ind], chi[ind])
pepsl.fname = "Nx{}_Ny{}_sx{}_sy{}_D{}_chi{}_run_left".format(
Nx, Ny, sxind, syind, D[ind], chi[ind])
# --------------------------------------------------------------------
# Calculate right eigenstate
for i in range(5):
try:
Ef, peps = run_tebd(Nx,
Ny,
d,
ops,
peps=peps,