Merge branch 'main' into bd/printing

Author: borisdevos

.gitignore

@@ -2,5 +2,6 @@
 *.jl.cov
 *.jl.mem
 Manifest.toml
+.vscode/
 .DS_Store
 build

src/fermions.jl

@@ -55,7 +55,8 @@ twist(a::FermionParity) = a.isodd ? -1 : +1
 
 function fusiontensor(a::I, b::I, c::I) where {I <: FermionParity}
     @warn "FermionParity Arrays do not preserve categorical properties." maxlog = 1
-    return fill(Int(Nsymbol(a, b, c)), (1, 1, 1, 1))
+    Nabc = Nsymbol(a, b, c)
+    return fill(Int(Nabc), (1, 1, 1, Nabc))
 end
 
 function Base.show(io::IO, a::FermionParity)

src/irreps/cu1irrep.jl

@@ -218,8 +218,9 @@ function Rsymbol(a::CU1Irrep, b::CU1Irrep, c::CU1Irrep)
 end
 
 function fusiontensor(a::CU1Irrep, b::CU1Irrep, c::CU1Irrep)
-    C = fill(zero(sectorscalartype(CU1Irrep)), dim(a), dim(b), dim(c), 1)
-    !Nsymbol(a, b, c) && return C
+    Nabc = Nsymbol(a, b, c)
+    C = fill(zero(sectorscalartype(CU1Irrep)), dim(a), dim(b), dim(c), Nabc)
+    Nabc || return C
     if c.j == 0
         if a.j == b.j == 0
             C[1, 1, 1, 1] = 1.0

src/irreps/dnirrep.jl

@@ -198,8 +198,9 @@ end
 
 function fusiontensor(a::I, b::I, c::I) where {N, I <: DNIrrep{N}}
     T = sectorscalartype(I)
-    C = zeros(T, dim(a), dim(b), dim(c), 1)
-    Nsymbol(a, b, c) || return C
+    Nabc = Nsymbol(a, b, c)
+    C = zeros(T, dim(a), dim(b), dim(c), Nabc)
+    Nabc || return C
 
     if c.j == 0
         if a.j == b.j == 0 || (2 * a.j == 2 * b.j == N)

src/irreps/irreps.jl

@@ -66,7 +66,8 @@ Rsymbol(a::I, b::I, c::I) where {I <: AbelianIrrep} = braidingscalartype(I)(Nsym
 
 function fusiontensor(a::I, b::I, c::I) where {I <: AbelianIrrep}
     T = sectorscalartype(I)
-    return fill(T(Nsymbol(a, b, c)), (1, 1, 1, 1))
+    Nabc = Nsymbol(a, b, c)
+    return fill(T(Nabc), (1, 1, 1, Nabc))
 end
 
 include("znirrep.jl")

src/irreps/su2irrep.jl

@@ -72,7 +72,9 @@ end
 
 function fusiontensor(a::SU2Irrep, b::SU2Irrep, c::SU2Irrep)
     T = sectorscalartype(SU2Irrep)
-    C = Array{T}(undef, dim(a), dim(b), dim(c), 1)
+    Nabc = Nsymbol(a, b, c)
+    C = Array{T}(undef, dim(a), dim(b), dim(c), Nabc)
+    Nabc || return C
     ja, jb, jc = a.j, b.j, c.j
 
     for kc in 1:dim(c), kb in 1:dim(b), ka in 1:dim(a)

src/sectors.jl

@@ -691,15 +691,15 @@ hopflink(a::I, b::I) where {I <: Sector} = sum(dim(c) * tr(Rsymbol(a, b, c) * Rs
 """
      Smatrix(::Type{I}) where {I <: Sector}
 
-Return the S-matrix of the sector type `I`, which is a matrix containing the hopflinks of all pairs of sectors of type `I`.
+Return the S-matrix of the sector type `I`, which is a matrix containing the hopflinks of all pairs of sectors of type `I`, with the second sector being taken dual.
 The S-matrix is not normalized by the total quantum dimension here.
 """
 function Smatrix(::Type{I}) where {I <: Sector}
     Base.IteratorSize(values(I)) isa Base.IsInfinite &&
         throw(ArgumentError("Only defined for sectors with a finite number of simple objects"))
     vals = values(I)
     l = length(vals)
-    return reshape([hopflink(a, b) for a in vals, b in vals], (l, l))
+    return reshape([hopflink(a, dual(b)) for a in vals, b in vals], (l, l))
 end
 
 """
@@ -715,13 +715,14 @@ end
 """
     topological_central_charge(::Type{I}) where {I <: Sector}
 
-Return the topological central charge c of the modular sector type `I`, where c is determined mod 8.
+Return the topological central charge c of the braided sector type `I`, where c is determined mod 8.
 We choose convention by restrict the returning value as rational numbers in (-4, 4].
 """
 function topological_central_charge(::Type{I}) where {I <: Sector}
-    ξ = sum(dim(a)^2 * twist(a) for a in values(I)) / dim(I)
-    @assert isapprox(abs(ξ), 1) "Sector $I is not modular"
-    c_float = angle(ξ) * 8 / (2π)
+    gauss_sum = sum(dim(a)^2 * twist(a) for a in values(I))
+    Θ = abs(gauss_sum)
+    @assert Θ > sqrt(eps(float(Θ))) "Topological central charge is not defined for sector type $I" # For non-modular categories, central charge is also meaningful. See https://arxiv.org/pdf/1602.05946. For super modular category, Gauss sum vanishes, and its central charge needs to be defined in another manner: https://arxiv.org/pdf/1603.09294.
+    c_float = angle(gauss_sum) * 8 / (2π)
 
     isapprox(c_float, -4) && return 4 // 1
 

test/newsectors.jl

@@ -65,7 +65,9 @@ end
 dim(s::NewSU2Irrep) = twice(s.j) + 1
 
 function fusiontensor(a::NewSU2Irrep, b::NewSU2Irrep, c::NewSU2Irrep)
-    C = Array{Float64}(undef, dim(a), dim(b), dim(c), 1)
+    Nabc = Nsymbol(a, b, c)
+    C = Array{Float64}(undef, dim(a), dim(b), dim(c), Nabc)
+    Nabc == 0 && return C
     ja, jb, jc = a.j, b.j, c.j
 
     for kc in 1:dim(c), kb in 1:dim(b), ka in 1:dim(a)

test/runtests.jl

@@ -266,15 +266,22 @@ end
 end
 
 @testset "Ismodular" begin
-    @test !ismodular(Z2Irrep)
-    @test !ismodular(Z3Irrep)
-    @test !ismodular(FermionParity)
-    @test !ismodular(A4Irrep)
-    @test !ismodular(IsingAnyon ⊠ Z2Irrep)
-    @test ismodular(IsingAnyon)
-    @test ismodular(FibonacciAnyon)
-    @test ismodular(TimeReversed{IsingAnyon})
-    @test ismodular(IsingAnyon ⊠ TimeReversed{IsingAnyon})
+    Tannakian_list = [Z2Irrep, Z3Irrep, A4Irrep, D3Irrep, D4Irrep, Z2Irrep ⊠ D4Irrep, D3Irrep ⊠ A4Irrep]
+    Super_Tannakian_list = [FermionParity, FermionParity ⊠ A4Irrep, FermionParity ⊠ Z3Irrep, D4Irrep ⊠ FermionParity]
+    UMTC_list = [
+        IsingAnyon, FibonacciAnyon, TimeReversed{IsingAnyon}, TimeReversed{FibonacciAnyon},
+        FibonacciAnyon ⊠ FibonacciAnyon, FibonacciAnyon ⊠ IsingAnyon, IsingAnyon ⊠ TimeReversed{IsingAnyon},
+        TimeReversed{FibonacciAnyon} ⊠ IsingAnyon, IsingAnyon ⊠ IsingAnyon ⊠ IsingAnyon,
+        IsingAnyon ⊠ FibonacciAnyon ⊠ IsingAnyon ⊠ TimeReversed{IsingAnyon} ⊠ TimeReversed{FibonacciAnyon},
+    ]
+    UMTC_over_RepG_list = [Z2Irrep ⊠ IsingAnyon, Z3Irrep ⊠ FibonacciAnyon, D3Irrep ⊠ TimeReversed{IsingAnyon}, A4Irrep ⊠ FibonacciAnyon ⊠ TimeReversed{IsingAnyon}]
+    UMTC_over_sRepG_list = [Z2Irrep ⊠ FermionParity ⊠ IsingAnyon, FermionParity ⊠ Z3Irrep ⊠ FibonacciAnyon, D3Irrep ⊠ TimeReversed{IsingAnyon} ⊠ FermionParity, A4Irrep ⊠ FibonacciAnyon ⊠ FermionParity ⊠ TimeReversed{IsingAnyon}]
+    for sect in [Tannakian_list..., Super_Tannakian_list..., UMTC_over_RepG_list..., UMTC_over_sRepG_list...]
+        @test !ismodular(sect)
+    end
+    for sect in UMTC_list
+        @test ismodular(sect)
+    end
 end
 
 @testset "Total quantum dimension" begin
@@ -290,6 +297,10 @@ end
 end
 
 @testset "Topological central charge" begin
+    @test topological_central_charge(Z2Irrep) == 0 // 1
+    @test topological_central_charge(A4Irrep) == 0 // 1
+    @test topological_central_charge(Z2Irrep ⊠ IsingAnyon) == 1 // 2
+    @test topological_central_charge(FibonacciAnyon ⊠ D4Irrep) == - 14 // 5
     @test topological_central_charge(IsingAnyon) == 1 // 2
     @test topological_central_charge(TimeReversed{IsingAnyon}) == - 1 // 2
     @test topological_central_charge(IsingAnyon ⊠ TimeReversed{IsingAnyon}) == 0 // 1