@@ -551,9 +551,6 @@ const TransposeKey{I<:Sector,N₁,N₂} = Tuple{<:FusionTree{I},<:FusionTree{I},
551551function _transpose ((f₁, f₂, p1, p2):: TransposeKey{I,N₁,N₂} ) where {I<: Sector ,N₁,N₂}
552552 N = N₁ + N₂
553553 p = linearizepermutation (p1, p2, length (f₁), length (f₂))
554- @info " _transpose"
555- @show f₁
556- @show f₂
557554 newtrees = repartition (f₁, f₂, N₁)
558555 length (p) == 0 && return newtrees
559556 i1 = findfirst (== (1 ), p)
@@ -623,12 +620,6 @@ function planar_trace(f₁::FusionTree{I}, f₂::FusionTree{I},
623620 newtrees = FusionTreeDict {Tuple{F₁,F₂},T} ()
624621 for ((f₁′, f₂′), coeff′) in repartition (f₁, f₂, N)
625622 for (f₁′′, coeff′′) in planar_trace (f₁′, q1′, q2′) # errors in this planar_trace first
626- @info " planar_trace"
627- @show f₁
628- @show f₂
629- @show f₁′
630- @show f₁′′
631- @show f₂′
632623 for (f12′′′, coeff′′′) in transpose (f₁′′, f₂′, p1′, p2′) # for a different unit errors here
633624 coeff = coeff′ * coeff′′ * coeff′′′
634625 if ! iszero (coeff)
@@ -711,16 +702,13 @@ function elementary_trace(f::FusionTree{I,N}, i) where {I<:Sector,N}
711702 F = fusiontreetype (I, N - 2 )
712703 newtrees = FusionTreeDict {F,T} ()
713704 _one = f. coupled # otherwise ArgumentError above thrown
714- @info " elementary_trace"
715- @show f
716705
717706 j = mod1 (i + 1 , N)
718707 b = f. uncoupled[i]
719708 b′ = f. uncoupled[j]
720709 # if trace is zero, return empty dict
721710 (b == dual (b′) && f. isdual[i] != f. isdual[j]) || return newtrees
722711 if i < N
723- @info " elementary_trace: i < N"
724712 inner_extended = (leftone (f. uncoupled[1 ]), f. uncoupled[1 ], f. innerlines... , f. coupled)
725713 a = inner_extended[i]
726714 d = inner_extended[i + 2 ]
@@ -755,11 +743,9 @@ function elementary_trace(f::FusionTree{I,N}, i) where {I<:Sector,N}
755743 if f. isdual[i]
756744 coeff *= frobeniusschur (b)
757745 end
758- @show f′
759746 push! (newtrees, f′ => coeff)
760747 return newtrees
761748 else # i == N
762- @info " elementary_trace: i == N"
763749 if N == 2
764750 f′ = FusionTree {I} ((), _one, (), (), ()) # or leftone(f.uncoupled[1]) == rightone(f.uncoupled[2])
765751 coeff = sqrtdim (b)
@@ -777,8 +763,6 @@ function elementary_trace(f::FusionTree{I,N}, i) where {I<:Sector,N}
777763 vertices_ = TupleTools. front (f. vertices)
778764 f_ = FusionTree (uncoupled_, coupled_, isdual_, inner_, vertices_)
779765 fs = FusionTree ((b,), b, (! f. isdual[1 ],), (), ())
780- @show f_
781- @show fs
782766 unit = leftone (fs. coupled)
783767 for (f_′, coeff) in merge (fs, f_, unit, 1 ) # coloring gets reversed here, should be the other unit
784768 f_′. innerlines[1 ] == unit || continue # is this one valid?
@@ -787,7 +771,6 @@ function elementary_trace(f::FusionTree{I,N}, i) where {I<:Sector,N}
787771 inner′ = N <= 4 ? () : Base. tail (Base. tail (f_′. innerlines))
788772 vertices′ = N <= 3 ? () : Base. tail (Base. tail (f_′. vertices))
789773 f′ = FusionTree (uncoupled′, unit, isdual′, inner′, vertices′) # and this one?
790- @show f′
791774 coeff *= sqrtdim (b)
792775 if ! (f. isdual[N])
793776 coeff *= conj (frobeniusschur (b))
0 commit comments