Merge remote-tracking branch 'upstream/master' into lb/initialize_env

Author: Yue-Zhengyuan

src/algorithms/time_evolution/simpleupdate.jl

@@ -142,10 +142,7 @@ function su_iter(
             state2[r, c] = _apply_sitegate(state2[r, c], gate; alg.purified)
         elseif length(sites) == 2
             (d, r, c), = _nn_bondrev(sites..., (Nr, Nc))
-            if alg.bipartite
-                length(sites) > 2 && error("Multi-site MPO gates are not compatible with bipartite states.")
-                r > 1 && continue
-            end
+            alg.bipartite && r > 1 && continue
             ϵ′ = _su_iter!(state2, gate, env2, sites, alg)
             ϵ = max(ϵ, ϵ′)
             (!alg.bipartite) && continue

test/examples/heisenberg.jl

@@ -55,7 +55,7 @@ end
 @testset "C4v AD optimization with scalartype T=$T and projector_alg=$projector_alg" for (T, projector_alg) in
     Iterators.product([Float64, ComplexF64], [:c4v_eigh, :c4v_qr])
     # initialize symmetric states
-    Random.seed!(1234567)
+    Random.seed!(123456789)
     symm = RotateReflect()
     H′ = heisenberg_XYZ_c4v(InfiniteSquare())
     H = T <: Real ? real(H′) : H′

test/timeevol/tf_ising_finiteT.jl

@@ -44,12 +44,13 @@ dt, nstep = 1.0e-3, 400
 β = dt * nstep
 
 # when g = 2, β = 0.4 and 2β = 0.8 belong to two phases (without and with nonzero σᶻ)
-@testset "Finite-T SU (bipartite = $(bipartite))" for bipartite in (true, false)
+@testset "Finite-T SU (force_mpo = $(force_mpo))" for force_mpo in (false, true)
     # use second order Trotter decomposition
     symmetrize_gates = true
+    bipartite = true
 
     # PEPO approach: results at β, or T = 2.5
-    alg = SimpleUpdate(; trunc = trunc_pepo, purified = false, bipartite)
+    alg = SimpleUpdate(; trunc = trunc_pepo, purified = false, bipartite, force_mpo)
     pepo, wts, info = time_evolve(pepo0, ham, dt, nstep, alg, wts0; symmetrize_gates)
 
     ## BP gauge fixing
@@ -72,7 +73,7 @@ dt, nstep = 1.0e-3, 400
     @test isapprox(abs.(result_2β), bm_2β, rtol = 1.0e-4)
 
     # Purification approach: results at 2β, or T = 1.25
-    alg = SimpleUpdate(; trunc = trunc_pepo, purified = true, bipartite)
+    alg = SimpleUpdate(; trunc = trunc_pepo, purified = true, bipartite, force_mpo)
     pepo, wts, info = time_evolve(pepo0, ham, dt, 2 * nstep, alg, wts0; symmetrize_gates)
     env = converge_env(InfinitePEPS(pepo), 8)
     result_2β′ = measure_mag(pepo, env; purified = true)