# 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,