Skip to content

[Fix] Transformers for non-abelian non-selfadjoint sectors #256

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
lkdvos merged 4 commits into from
Jul 13, 2025
Merged

[Fix] Transformers for non-abelian non-selfadjoint sectors #256

Show file tree
Hide file tree
Changes from 3 commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
150f4d6
Correctly handle dual of sectors
lkdvos Jul 12, 2025
0242b64
Fix debug message for transformers of length 1
lkdvos Jul 12, 2025
6f52ba0
Add non-abelian non-selfadjoint test
lkdvos Jul 12, 2025
c33606d
Bump v0.14.8
lkdvos Jul 13, 2025
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
6 changes: 3 additions & 3 deletions src/tensors/treetransformers.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -65,12 +65,12 @@ function GenericTreeTransformer(transform, p, Vdst, Vsrc)
I = sectortype(Vsrc)

uncoupleds_src = map(structure_src.fusiontreelist) do (f₁, f₂)
return TupleTools.vcat(f₁.uncoupled, f₂.uncoupled)
return TupleTools.vcat(f₁.uncoupled, dual.(f₂.uncoupled))
end
uncoupleds_src_unique = unique(uncoupleds_src)

uncoupleds_dst = map(structure_dst.fusiontreelist) do (f₁, f₂)
return TupleTools.vcat(f₁.uncoupled, f₂.uncoupled)
return TupleTools.vcat(f₁.uncoupled, dual.(f₂.uncoupled))
end

T = sectorscalartype(I)
Expand Down Expand Up @@ -116,7 +116,7 @@ function GenericTreeTransformer(transform, p, Vdst, Vsrc)

@debug("TreeTransformer for $Vsrc to $Vdst via $p",
nblocks = length(data),
sz_median = size(data[end ÷ 2][1], 1),
sz_median = size(data[cld(end, 2)][1], 1),
sz_max = size(data[1][1], 1),
Δt)

Expand Down
6 changes: 6 additions & 0 deletions test/runtests.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -100,6 +100,12 @@ VfSU₂ = (ℂ[FermionSpin](0 => 3, 1 // 2 => 1),
ℂ[FermionSpin](1 // 2 => 1, 1 => 1)',
ℂ[FermionSpin](0 => 2, 1 // 2 => 2),
ℂ[FermionSpin](0 => 1, 1 // 2 => 1, 3 // 2 => 1)')
VSU₂U₁ = (Vect[SU2Irrep ⊠ U1Irrep]((0, 0) => 1, (1 // 2, -1) => 1),
Vect[SU2Irrep ⊠ U1Irrep]((0, 0) => 2, (0, 2) => 1, (1, 0) => 1, (1, -2) => 1,
(1 // 2, -1) => 1),
Vect[SU2Irrep ⊠ U1Irrep]((1 // 2, 1) => 1, (1, -2) => 1)',
Vect[SU2Irrep ⊠ U1Irrep]((0, 0) => 2, (0, 2) => 1, (1 // 2, 1) => 1),
Vect[SU2Irrep ⊠ U1Irrep]((0, 0) => 1, (1 // 2, 1) => 1)')
# VSU₃ = (ℂ[SU3Irrep]((0, 0, 0) => 3, (1, 0, 0) => 1),
# ℂ[SU3Irrep]((0, 0, 0) => 3, (2, 0, 0) => 1)',
# ℂ[SU3Irrep]((1, 1, 0) => 1, (2, 1, 0) => 1),
Expand Down
10 changes: 5 additions & 5 deletions test/tensors.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
for V in (Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂)#, VSU₃)
for V in (Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂, VSU₂U₁)#, VSU₃)
V1, V2, V3, V4, V5 = V
@assert V3 * V4 * V2 ≿ V1' * V5' # necessary for leftorth tests
@assert V3 * V4 ≾ V1' * V2' * V5' # necessary for rightorth tests
Expand All @@ -10,15 +10,15 @@ spacelist = try
if Sys.iswindows()
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂)
elseif Sys.isapple()
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VfU₁, VfSU₂)#, VSU₃)
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VfU₁, VfSU₂, VSU₂U₁)#, VSU₃)
else
(Vtr, Vℤ₂, Vfℤ₂, VU₁, VCU₁, VSU₂, VfSU₂)#, VSU₃)
(Vtr, Vℤ₂, Vfℤ₂, VU₁, VCU₁, VSU₂, VfSU₂, VSU₂U₁)#, VSU₃)
end
else
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂)#, VSU₃)
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂, VSU₂U₁)#, VSU₃)
end
catch
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂)#, VSU₃)
(Vtr, Vℤ₂, Vfℤ₂, Vℤ₃, VU₁, VfU₁, VCU₁, VSU₂, VfSU₂, VSU₂U₁)#, VSU₃)
end

for V in spacelist
Expand Down
Loading