Add fusionproduct

src/TensorKit.jl

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

src/sectors/sectors.jl

@@ -120,6 +120,7 @@ Base.isless(::Trivial, ::Trivial) = false
 #---------------------------------------------
 """
     ⊗(a::I, b::I) where {I<:Sector}
+    fusionproduct(a::I, b::I) where {I<:Sector}
 
 Return an iterable of elements of `c::I` that appear in the fusion product `a ⊗ b`.
 
@@ -128,6 +129,7 @@ Note that every element `c` should appear at most once, fusion degeneracies (if
 """
 ⊗(::Trivial, ::Trivial) = (Trivial(),)
 ⊗(I::Sector) = (I,)
+fusionproduct(a::I, b::Vararg{I}) where {I<:Sector} = ⊗(a, b...)
 
 """
     Nsymbol(a::I, b::I, c::I) where {I<:Sector} -> Integer