use SymmetrySectors.nsymbol

Author: ogauthe

Project.toml

@@ -29,7 +29,7 @@ LRUCache = "1.6.1"
 LabelledNumbers = "0.1.0"
 LinearAlgebra = "1.10.0"
 Strided = "2.2.0"
-SymmetrySectors = "0.1.7"
+SymmetrySectors = "0.1.9"
 TensorAlgebra = "0.2.4"
 TensorProducts = "0.1.3"
 TypeParameterAccessors = "0.2.0"

src/fusion_trees/fusiontree.jl

@@ -5,9 +5,16 @@
 
 using GradedUnitRanges: AbstractGradedUnitRange, GradedUnitRanges, flip, isdual, sector_type
 using SymmetrySectors:
-  ×, AbstractSector, SectorProduct, SymmetrySectors, arguments, trivial, to_gradedrange
+  ×,
+  AbstractSector,
+  SectorProduct,
+  SymmetrySectors,
+  arguments,
+  nsymbol,
+  to_gradedrange,
+  trivial
 using TensorAlgebra: flatten_tuples
-using TensorProducts: ⊗, tensor_product
+using TensorProducts: ⊗
 
 #
 # A fusion tree fuses N sectors sec1, secN  onto one sector fused_sec. A given set of
@@ -185,14 +192,6 @@ function outer_multiplicity_kron(
   return linear_inds[outer_multiplicity1, outer_multiplicity2]
 end
 
-# TODO move to GradedUnitRanges
-function nsymbol(s1::AbstractSector, s2::AbstractSector, s3::AbstractSector)
-  full_space = to_gradedrange(s1 ⊗ s2)
-  x = findfirst(==(s3), blocklabels(full_space))
-  isnothing(x) && return 0  # OR labelled(0, s3)?
-  return Int(blocklengths(full_space)[x])
-end
-
 function outer_multiplicity_split(
   sec1::S, sec2::S, fused::S, outer_mult_index::Integer
 ) where {S<:SectorProduct}

test/Project.toml

@@ -16,4 +16,5 @@ Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 Aqua = "0.8.9"
 SafeTestsets = "0.1"
 Suppressor = "0.2"
+SymmetrySectors = "0.1.9"
 Test = "1.10"