fix construct and remove braidingtensor preprocessors

Jutho · Jan 3, 2024 · fff1008 · fff1008
1 parent 8a2df70
commit fff1008
Show file tree

Hide file tree

Showing 5 changed files with 167 additions and 157 deletions.
diff --git a/src/planar/macros.jl b/src/planar/macros.jl
@@ -102,7 +102,7 @@ function _plansor(expr, kwargs...)
 
     tparser = TO.tensorparser(expr, kwargs...)
     pparser = planarparser(expr, kwargs...)
-    insert!(tparser.preprocessors, 5, _remove_braidingtensors)
+    insert!(tparser.preprocessors, 4, _remove_braidingtensors)
     tensorex = tparser(expr)
     planarex = pparser(expr)
 

diff --git a/src/planar/planaroperations.jl b/src/planar/planaroperations.jl
@@ -61,6 +61,7 @@ function planarcontract!(C::AbstractTensorMap{S,N₁,N₂},
     codB, domB = codomainind(B), domainind(B)
     oindA, cindA = pA
     cindB, oindB = pB
+    # @show codA, domA, codB, domB, oindA, cindA, oindB, cindB, pAB
     oindA, cindA, oindB, cindB = reorder_indices(codA, domA, codB, domB, oindA, cindA,
                                                  oindB, cindB, pAB...)
 
@@ -91,6 +92,7 @@ _cyclicpermute(t::Tuple{}) = ()
 function reorder_indices(codA, domA, codB, domB, oindA, oindB, p1, p2)
     N₁ = length(oindA)
     N₂ = length(oindB)
+    # @show codA, domA, codB, domB, oindA, oindB, p1, p2
     @assert length(p1) == N₁ && all(in(p1), 1:N₁)
     @assert length(p2) == N₂ && all(in(p2), N₁ .+ (1:N₂))
     oindA2 = TupleTools.getindices(oindA, p1)