print_str = ''
for j in range(Ny-1):
print_str += '{} '.format(density_top[0][i][j])
print(print_str)
for i in range(Nx):
print_str = ''
for j in range(Ny-1):
print_str += '{} '.format(density_bot[0][i][j])
print(print_str)
print('Horizontal Density')
for i in range(Ny):
print_str = ''
for j in range(Nx-1):
print_str += '{} '.format(density_top[1][i][j])
print(print_str)
for i in range(Ny):
print_str = ''
for j in range(Nx-1):
print_str += '{} '.format(density_bot[1][i][j])
print(print_str)
# Increase PEPS D if wanted
if (len(D)-1 > ind):
if (D[ind+1] > D[ind]):
peps.increase_mbd(D[ind+1],chi=chi[ind+1])
peps.chi = chi[ind+1]
if (D[ind+1] > D[ind]):
pepsl.increase_mbd(D[ind+1],chi=chi[ind+1])
pepsl.chi = chi[ind+1]
peps.normalize()
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],
n_step=n_step[ind],
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,
D=D[ind],