Fix sectorscalartype fallback definition for UniqueFusion() sector types

Project.toml

@@ -1,7 +1,7 @@
 name = "TensorKitSectors"
 uuid = "13a9c161-d5da-41f0-bcbd-e1a08ae0647f"
 authors = ["Lukas Devos", "Jutho Haegeman"]
-version = "0.1.2"
+version = "0.1.3"
 
 [deps]
 HalfIntegers = "f0d1745a-41c9-11e9-1dd9-e5d34d218721"

src/sectors.jl

@@ -105,14 +105,21 @@
 Return the scalar type of the topological data (Fsymbol, Rsymbol) of the sector `I`.
 """
 function sectorscalartype(::Type{I}) where {I<:Sector}
-    if BraidingStyle(I) isa NoBraiding
-        return eltype(Core.Compiler.return_type(Fsymbol, NTuple{6,I}))
+    if BraidingStyle(I) === NoBraiding()
+        return _Fscalartype(I)
     else
-        Feltype = eltype(Core.Compiler.return_type(Fsymbol, NTuple{6,I}))
-        Reltype = eltype(Core.Compiler.return_type(Rsymbol, NTuple{3,I}))
-        return Base.promote_op(*, Feltype, Reltype)
+        return Base.promote_op(*, _Fscalartype(I), _Rscalartype(I))
     end
 end
+function _Fscalartype(::Type{I}) where {I<:Sector}
+    Ftype = Core.Compiler.return_type(Fsymbol, NTuple{6,I})
+    return FusionStyle(I) === UniqueFusion() ? Ftype : eltype(Ftype)
+end
+function _Rscalartype(::Type{I}) where {I<:Sector}
+    BraidingStyle(I) === NoBraiding() && throw(ArgumentError("No braiding for sector $I"))
+    Rtype = Core.Compiler.return_type(Rsymbol, NTuple{3,I})
+    return FusionStyle(I) === UniqueFusion() ? Rtype : eltype(Rtype)
+end
 
 """
     isreal(::Type{<:Sector}) -> Bool

test/sectors.jl

@@ -10,9 +10,14 @@ Istr = TKS.type_repr(I)
     @constinferred dim(s[1])
     @constinferred frobeniusschur(s[1])
     @constinferred Nsymbol(s...)
-    @constinferred Rsymbol(s...)
-    @constinferred Bsymbol(s...)
-    @constinferred Fsymbol(s..., s...)
+    R = @constinferred Rsymbol(s...)
+    B = @constinferred Bsymbol(s...)
+    F = @constinferred Fsymbol(s..., s...)
+    if FusionStyle(I) === SimpleFusion()
+        @test typeof(R * F) <: @constinferred sectorscalartype(I)
+    else
+        @test Base.promote_op(*, eltype(R), eltype(F)) <: @constinferred sectorscalartype(I)
+    end
     it = @constinferred s[1] ⊗ s[2]
     @constinferred ⊗(s..., s...)
 end