update frobeniusschur

Author: lkdvos

src/TensorKit.jl

@@ -41,7 +41,8 @@ export infimum, supremum, isisomorphic, ismonomorphic, isepimorphic
 
 # methods for sectors and properties thereof
 export sectortype, sectors, hassector, Nsymbol, Fsymbol, Rsymbol, Bsymbol,
-    frobeniusschur, twist, otimes, sectorscalartype, deligneproduct
+    frobeniusschur, frobenius_schur_phase, frobenius_schur_indicator,
+    twist, otimes, sectorscalartype, deligneproduct
 export fusiontrees, braid, permute, transpose
 export ZNSpace, SU2Irrep, U1Irrep, CU1Irrep
 # other fusion tree manipulations, should not be exported:

src/fusiontrees/manipulations.jl

@@ -246,7 +246,7 @@ function flip(f₁::FusionTree{I, N₁}, f₂::FusionTree{I, N₂}, i::Int; inv:
     @assert 0 < i ≤ N₁ + N₂
     if i ≤ N₁
         a = f₁.uncoupled[i]
-        χₐ = frobeniusschur(a)
+        χₐ = frobenius_schur_phase(a)
         θₐ = twist(a)
         if !inv
             factor = f₁.isdual[i] ? χₐ * θₐ : one(θₐ)
@@ -259,7 +259,7 @@ function flip(f₁::FusionTree{I, N₁}, f₂::FusionTree{I, N₂}, i::Int; inv:
     else
         i -= N₁
         a = f₂.uncoupled[i]
-        χₐ = frobeniusschur(a)
+        χₐ = frobenius_schur_phase(a)
         θₐ = twist(a)
         if !inv
             factor = f₂.isdual[i] ? χₐ * one(θₐ) : θₐ
@@ -302,7 +302,7 @@ function bendright(f₁::FusionTree{I, N₁}, f₂::FusionTree{I, N₂}) where {
 
     coeff₀ = sqrtdim(c) * invsqrtdim(a)
     if f₁.isdual[N₁]
-        coeff₀ *= conj(frobeniusschur(dual(b)))
+        coeff₀ *= conj(frobenius_schur_phase(dual(b)))
     end
     if FusionStyle(I) isa MultiplicityFreeFusion
         coeff = coeff₀ * Bsymbol(a, b, c)
@@ -343,7 +343,7 @@ function foldright(f₁::FusionTree{I, N₁}, f₂::FusionTree{I, N₂}) where {
     isduala = f₁.isdual[1]
     factor = sqrtdim(a)
     if !isduala
-        factor *= conj(frobeniusschur(a))
+        factor *= conj(frobenius_schur_phase(a))
     end
     c1 = dual(a)
     c2 = f₁.coupled
@@ -767,7 +767,7 @@ function elementary_trace(f::FusionTree{I, N}, i) where {I <: Sector, N}
             end
         end
         if f.isdual[i]
-            coeff *= frobeniusschur(b)
+            coeff *= frobenius_schur_phase(b)
         end
         push!(newtrees, f′ => coeff)
         return newtrees
@@ -777,7 +777,7 @@ function elementary_trace(f::FusionTree{I, N}, i) where {I <: Sector, N}
             f′ = FusionTree{I}((), unit, (), (), ())
             coeff = sqrtdim(b)
             if !(f.isdual[N])
-                coeff *= conj(frobeniusschur(b))
+                coeff *= conj(frobenius_schur_phase(b))
             end
             push!(newtrees, f′ => coeff)
             return newtrees
@@ -798,7 +798,7 @@ function elementary_trace(f::FusionTree{I, N}, i) where {I <: Sector, N}
             f′ = FusionTree(uncoupled′, unit, isdual′, inner′, vertices′)
             coeff *= sqrtdim(b)
             if !(f.isdual[N])
-                coeff *= conj(frobeniusschur(b))
+                coeff *= conj(frobenius_schur_phase(b))
             end
             newtrees[f′] = get(newtrees, f′, zero(coeff)) + coeff
         end

src/tensors/braidingtensor.jl

@@ -385,10 +385,10 @@ end
 #                 f₁′′.innerlines[1] == u || continue
 #                 coeff = coeff′ * coeff′′ * sqrtdim(a) * sqrtdim(b)
 #                 if f₁′′.isdual[1]
-#                     coeff *= frobeniusschur(a)
+#                     coeff *= frobenius_schur_phase(a)
 #                 end
 #                 if f₁′′.isdual[3]
-#                     coeff *= frobeniusschur(b)
+#                     coeff *= frobenius_schur_phase(b)
 #                 end
 #                 f12 = (f₀, f₂′)
 #                 if @isdefined newtrees
@@ -453,10 +453,10 @@ end
 #                 f₂′′.innerlines[1] == u || continue
 #                 coeff = coeff′ * conj(coeff′′ * sqrtdim(a) * sqrtdim(b))
 #                 if f₂′′.isdual[1]
-#                     coeff *= conj(frobeniusschur(a))
+#                     coeff *= conj(frobenius_schur_phase(a))
 #                 end
 #                 if f₂′′.isdual[3]
-#                     coeff *= conj(frobeniusschur(b))
+#                     coeff *= conj(frobenius_schur_phase(b))
 #                 end
 #                 f12 = (f₁′, f₀)
 #                 if @isdefined newtrees
@@ -517,7 +517,7 @@ end
 #                 f₁′′.innerlines[1] == u || continue
 #                 coeff = coeff′ * coeff′′ * sqrtdim(a)
 #                 if f₁′′.isdual[1]
-#                     coeff *= frobeniusschur(a)
+#                     coeff *= frobenius_schur_phase(a)
 #                 end
 #                 f₁′′′ = FusionTree{I}((b,), b, (f₁′′.isdual[3],), (), ())
 #                 f12 = (f₁′′′, f₂′)
@@ -579,7 +579,7 @@ end
 #                 f₂′′.innerlines[1] == u || continue
 #                 coeff = coeff′ * conj(coeff′′ * sqrtdim(a))
 #                 if f₂′′.isdual[1]
-#                     coeff *= conj(frobeniusschur(a))
+#                     coeff *= conj(frobenius_schur_phase(a))
 #                 end
 #                 f₂′′′ = FusionTree{I}((b,), b, (f₂′′.isdual[3],), (), ())
 #                 f12 = (f₁′, f₂′′′)