Skip to content

Make SUWeight axis order definite #315

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
Yue-Zhengyuan merged 13 commits into from
Jan 6, 2026
Merged

Make SUWeight axis order definite #315

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
79297a0
Make SUWeight axis order definite
Yue-Zhengyuan Dec 12, 2025
638e9d5
Update related tests
Yue-Zhengyuan Dec 12, 2025
259c97f
Update docs on CTMRGEnv(wts)
Yue-Zhengyuan Dec 12, 2025
2b3a532
Update SU and bond truncation
Yue-Zhengyuan Dec 12, 2025
6ec6239
Test SU with non-standard virtual dualness
Yue-Zhengyuan Dec 12, 2025
f230f2f
Undo metadata changes
Yue-Zhengyuan Dec 12, 2025
d642197
Update docstring
Yue-Zhengyuan Dec 12, 2025
dcbe1ac
Refactor SUWeight rotation
Yue-Zhengyuan Dec 15, 2025
07f3fc6
Apply suggestions and add SUWeight random init
Yue-Zhengyuan Dec 16, 2025
b8f6686
Streamline CTMRGEnv(wts)
Yue-Zhengyuan Dec 16, 2025
1a3e816
Correct a comment
Yue-Zhengyuan Dec 16, 2025
12f9dc3
Apply suggestions
Yue-Zhengyuan Jan 2, 2026
a8aae2b
Improve `rand!(wts)`
Yue-Zhengyuan Jan 3, 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
2 changes: 1 addition & 1 deletion docs/src/examples/hubbard_su/index.md
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,7 +130,7 @@ is achieved by using `trunc=truncrank(χ)` with different `χ`s in the CTMRG run
χenv₀, χenv = 6, 16
env_space = Vect[fℤ₂](0 => χenv₀ / 2, 1 => χenv₀ / 2)
normalize!.(peps.A, Inf)
env = CTMRGEnv(wts, peps)
env = CTMRGEnv(wts)
for χ in [χenv₀, χenv]
global env, = leading_boundary(
env, peps; alg = :sequential, tol = 1.0e-8, maxiter = 50, trunc = truncrank(χ)
Expand Down
2 changes: 1 addition & 1 deletion docs/src/examples/hubbard_su/main.ipynb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -140,7 +140,7 @@
"χenv₀, χenv = 6, 16\n",
"env_space = Vect[fℤ₂](0 => χenv₀ / 2, 1 => χenv₀ / 2)\n",
"normalize!.(peps.A, Inf)\n",
"env = CTMRGEnv(wts, peps)\n",
"env = CTMRGEnv(wts)\n",
"for χ in [χenv₀, χenv]\n",
" global env, = leading_boundary(\n",
" env, peps; alg = :sequential, tol = 1.0e-8, maxiter = 50, trunc = truncrank(χ)\n",
Expand Down
2 changes: 1 addition & 1 deletion examples/hubbard_su/main.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -83,7 +83,7 @@ is achieved by using `trunc=truncrank(χ)` with different `χ`s in the CTMRG run
χenv₀, χenv = 6, 16
env_space = Vect[fℤ₂](0 => χenv₀ / 2, 1 => χenv₀ / 2)
normalize!.(peps.A, Inf)
env = CTMRGEnv(wts, peps)
env = CTMRGEnv(wts)
for χ in [χenv₀, χenv]
global env, = leading_boundary(
env, peps; alg = :sequential, tol = 1.0e-8, maxiter = 50, trunc = truncrank(χ)
Expand Down
1 change: 1 addition & 0 deletions src/PEPSKit.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,7 @@
module PEPSKit

using LinearAlgebra, Statistics, Base.Threads, Base.Iterators, Printf
using Random
using Compat
using Accessors: @set, @reset
using VectorInterface
Expand Down
2 changes: 1 addition & 1 deletion src/algorithms/time_evolution/evoltools.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -233,7 +233,7 @@ function _apply_gate(

a, b = absorb_s(a, s, b)
if need_flip
a, s, b = flip_svd(a, s, b)
a, s, b = flip(a, numind(a)), _fliptwist_s(s), flip(b, 1)
end
return a, s, b, ϵ
end
Expand Down
2 changes: 1 addition & 1 deletion src/algorithms/time_evolution/simpleupdate3site.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -479,7 +479,7 @@ function _su3site_se!(
Ms = [first(_unfuse_physicalspace(M, Vphy)) for (M, Vphy) in zip(Ms, Vphys)]
end
for (wt, wt_idx, flip) in zip(wts, wt_idxs, flips)
env[CartesianIndex(wt_idx)] = normalize(flip ? _flip_s(wt) : wt, Inf)
env[CartesianIndex(wt_idx)] = normalize(flip ? _fliptwist_s(wt) : wt, Inf)
end
end
# restore virtual arrows in `Ms`
Expand Down
4 changes: 2 additions & 2 deletions src/algorithms/truncation/bond_truncation.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -139,7 +139,7 @@ function bond_truncate(
a, s, b = svd_trunc!(permute(_combine_ab(a, b), perm_ab); trunc = alg.trunc)
a, b = absorb_s(a, s, b)
if need_flip
a, s, b = flip_svd(a, s, b)
a, s, b = flip(a, numind(a)), _fliptwist_s(s), flip(b, 1)
end
return a, s, b, (; fid, Δfid, Δs)
end
Expand Down Expand Up @@ -186,7 +186,7 @@ function bond_truncate(
@tensor a[-1 -2; -3] := Qa[-1 -2 3] * u[3 -3]
@tensor b[-1; -2 -3] := vh[-1 1] * Qb[1 -2 -3]
if need_flip
a, s, b = flip_svd(a, s, b)
a, s, b = flip(a, numind(a)), _fliptwist_s(s), flip(b, 1)
end
return a, s, b, info
end
Loading