Skip to content

Cover force_mpo in simple update test #350

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
leburgel merged 3 commits into from
Apr 10, 2026
Merged

Cover force_mpo in simple update test #350

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
ca5f3af
Remove bipartite state restriction on 2-site MPO gate
Yue-Zhengyuan Apr 8, 2026
4789be7
Update test to catch force_mpo error for finite-T Ising
Yue-Zhengyuan Apr 9, 2026
2a0f3c5
Merge remote-tracking branch 'upstream/master' into su-bipartite
Yue-Zhengyuan Apr 10, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
5 changes: 1 addition & 4 deletions src/algorithms/time_evolution/simpleupdate.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down
7 changes: 4 additions & 3 deletions test/timeevol/tf_ising_finiteT.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -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
Comment thread
leburgel marked this conversation as resolved.

# 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
Expand All @@ -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)
Expand Down
Loading