update CGC orthogonality test for GenericFusion

Author: lkdvos

test/sectors.jl

@@ -92,11 +92,18 @@ end
     hasfusiontensor(I) || return nothing
     for a in smallset(I), b in smallset(I)
         cs = vec(collect(a ⊗ b))
-        CGCs = map(c -> reshape(fusiontensor(a, b, c), :, dim(c)), cs)
-        M = map(Iterators.product(CGCs, CGCs)) do (cgc1, cgc2)
-            return LinearAlgebra.norm(cgc1' * cgc2)
+        cgcs = map(c -> fusiontensor(a, b, c), cs)
+        for (c, cgc) in zip(cs, cgcs), (c′, cgc′) in zip(cs, cgcs)
+            for μ in 1:Nsymbol(a, b, c), ν in 1:Nsymbol(a, b, c′)
+                @tensor overlap[mc mc'] := conj(view(cgc, :, :, :, μ)[ma mb mc]) *
+                    view(cgc′, :, :, :, ν)[ma mb mc']
+                if μ == ν && c == c′
+                    @test isapprox(overlap, LinearAlgebra.I; atol = 1.0e-12)
+                else
+                    @test isapprox(LinearAlgebra.norm(overlap), 0; atol = 1.0e-12)
+                end
+            end
         end
-        @test isapprox(M' * M, LinearAlgebra.Diagonal(dim.(cs)); atol = 1.0e-12)
     end
 end