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

* Correctly handle dual of sectors

* Fix debug message for transformers of length 1

* Add non-abelian non-selfadjoint test

* Bump v0.14.8

Author: lkdvos

Project.toml

@@ -1,7 +1,7 @@
 name = "TensorKit"
 uuid = "07d1fe3e-3e46-537d-9eac-e9e13d0d4cec"
 authors = ["Jutho Haegeman"]
-version = "0.14.7"
+version = "0.14.8"
 
 [deps]
 LRUCache = "8ac3fa9e-de4c-5943-b1dc-09c6b5f20637"

src/tensors/treetransformers.jl

@@ -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)
@@ -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)
 

test/runtests.jl

@@ -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),

test/tensors.jl

@@ -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
@@ -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