defined set_sector_label

Author: ogauthe

NDTensors/src/lib/SymmetrySectors/src/sector_definitions/u1.jl

@@ -11,11 +11,11 @@ struct U1{T} <: AbstractSector
 end
 
 SymmetryStyle(::Type{<:U1}) = AbelianStyle()
-
-GradedAxes.dual(u::U1) = U1(-u.n)
-
 sector_label(u::U1) = u.n
 
+set_sector_label(s::U1, sector_label) = typeof(s)(sector_label)
+GradedAxes.dual(s::U1) = set_sector_label(s, -sector_label(s))
+
 trivial(::Type{U1}) = trivial(U1{Int})
 trivial(::Type{U1{T}}) where {T} = U1(zero(T))
 

NDTensors/src/lib/SymmetrySectors/src/sector_definitions/zn.jl

@@ -6,20 +6,20 @@ using ...GradedAxes: GradedAxes
 
 struct Z{N} <: AbstractSector
   m::Int
-  Z{N}(m) where {N} = new{N}(m % N)
+  Z{N}(m) where {N} = new{N}(mod(m, N))
 end
 
-SymmetryStyle(::Type{<:Z}) = AbelianStyle()
+modulus(::Type{Z{N}}) where {N} = N
+modulus(c::Z) = modulus(typeof(c))
 
+SymmetryStyle(::Type{<:Z}) = AbelianStyle()
 sector_label(c::Z) = c.m
-modulus(::Type{Z{N}}) where {N} = N
 
-modulus(c::Z) = modulus(typeof(c))
+set_sector_label(s::Z, sector_label) = typeof(s)(sector_label)
+GradedAxes.dual(s::Z) = set_sector_label(s, -sector_label(s))
 
 trivial(sector_type::Type{<:Z}) = sector_type(0)
 
 function abelian_label_fusion_rule(sector_type::Type{<:Z}, n1, n2)
-  return sector_type((n1 + n2) % modulus(sector_type))
+  return sector_type(n1 + n2)
 end
-
-GradedAxes.dual(c::Z) = typeof(c)(mod(-sector_label(c), modulus(c)))

NDTensors/src/lib/SymmetrySectors/test/test_simple_sectors.jl

@@ -70,6 +70,8 @@ using Test: @inferred, @test, @testset, @test_throws
     @test dual(Z{2}(1)) == Z{2}(1)
     @test isless(Z{2}(0), Z{2}(1))
     @test !isless(Z{2}(1), Z{2}(0))
+    @test Z{2}(0) == z0
+    @test Z{2}(-3) == z1
 
     @test Z{2}(0) == TrivialSector()
     @test TrivialSector() < Z{2}(1)