fix sectorscalartype for product sectors (#72)

Author: borisdevos

src/product.jl

@@ -228,7 +228,11 @@ function braidingscalartype(::Type{<:ProductSector{T}}) where {T <: SectorTuple}
     return typeof(prod(zero ∘ braidingscalartype, _sectors(T)))
 end
 function sectorscalartype(::Type{<:ProductSector{T}}) where {T <: SectorTuple}
-    return typeof(prod(zero ∘ sectorscalartype, _sectors(T)))
+    return if BraidingStyle(ProductSector{T}) == NoBraiding()
+        typeof(prod(zero ∘ fusionscalartype, _sectors(T)))
+    else
+        typeof(prod(zero ∘ sectorscalartype, _sectors(T)))
+    end
 end
 function dimscalartype(::Type{<:ProductSector{T}}) where {T <: SectorTuple}
     return typeof(prod(zero ∘ dimscalartype, _sectors(T)))

test/runtests.jl

@@ -11,7 +11,7 @@ const sectorlist = (
     DNIrrep{3}, DNIrrep{4}, DNIrrep{5}, CU1Irrep,
     A4Irrep, SU2Irrep, NewSU2Irrep,
     FibonacciAnyon, IsingAnyon, FermionParity,
-    FermionParity ⊠ FermionParity,
+    FermionParity ⊠ FermionParity, FibonacciAnyon ⊠ PlanarTrivial,
     Z3Irrep ⊠ Z4Irrep, FermionParity ⊠ U1Irrep ⊠ SU2Irrep,
     FermionParity ⊠ SU2Irrep ⊠ SU2Irrep, NewSU2Irrep ⊠ NewSU2Irrep,
     NewSU2Irrep ⊠ SU2Irrep, FermionParity ⊠ SU2Irrep ⊠ NewSU2Irrep,