Added TimeReversed and its related functions (#11)

Author: Chenqitrg

.gitignore

@@ -2,3 +2,4 @@
 *.jl.cov
 *.jl.mem
 /Manifest.toml
+.DS_Store

src/TensorKitSectors.jl

@@ -18,7 +18,7 @@ export rightone, leftone
 export pentagon_equation, hexagon_equation
 
 export Trivial, Z2Irrep, Z3Irrep, Z4Irrep, ZNIrrep, U1Irrep, SU2Irrep, CU1Irrep
-export ProductSector
+export ProductSector, TimeReversed
 export FermionParity, FermionNumber, FermionSpin
 export PlanarTrivial, FibonacciAnyon, IsingAnyon
 
@@ -58,7 +58,7 @@ include("precompile.jl")
 function __precompile__()
     for I in (Trivial, Z2Irrep, Z3Irrep, Z4Irrep, ZNIrrep, U1Irrep, SU2Irrep, CU1Irrep,
               FermionParity, FermionNumber, FermionSpin, PlanarTrivial, FibonacciAnyon,
-              IsingAnyon)
+              IsingAnyon, TimeReversed{IsingAnyon}, TimeReversed{FibonacciAnyon})
         precompile_sector(I)
     end
 end

src/sectors.jl

@@ -459,3 +459,62 @@ function Base.iterate(s::SectorSet{I}, args...) where {I<:Sector}
     val, state = next
     return convert(I, s.f(val)), state
 end
+
+# Time reversed sector
+struct TimeReversed{I<:Sector} <: Sector
+    a::I
+    function TimeReversed{I}(a::I) where {I<:Sector}
+        if BraidingStyle(I) isa NoBraiding
+            throw(ArgumentError("TimeReversed is not defined for sectors $I with no braiding"))
+        end
+        return new{I}(a)
+    end
+end
+FusionStyle(::Type{TimeReversed{I}}) where {I<:Sector} = FusionStyle(I)
+BraidingStyle(::Type{TimeReversed{I}}) where {I<:Sector} = BraidingStyle(I)
+function Nsymbol(a::TimeReversed{I}, b::TimeReversed{I},
+                 c::TimeReversed{I}) where {I<:Sector}
+    return Nsymbol(a.a, b.a, c.a)
+end
+
+function Fsymbol(a::TimeReversed{I}, b::TimeReversed{I}, c::TimeReversed{I},
+                 d::TimeReversed{I}, e::TimeReversed{I},
+                 f::TimeReversed{I}) where {I<:Sector}
+    return Fsymbol(a.a, b.a, c.a, d.a, e.a, f.a)
+end
+function Rsymbol(a::TimeReversed{I}, b::TimeReversed{I},
+                 c::TimeReversed{I}) where {I<:Sector}
+    return adjoint(Rsymbol(a.a, b.a, c.a))
+end
+
+Base.one(::Type{TimeReversed{I}}) where {I<:Sector} = TimeReversed{I}(one(I))
+Base.conj(c::TimeReversed{I}) where {I<:Sector} = TimeReversed{I}(conj(c.a))
+function ⊗(c1::TimeReversed{I}, c2::TimeReversed{I}) where {I<:Sector}
+    return Iterators.map(TimeReversed{I}, c1.a ⊗ c2.a)
+end
+function Base.IteratorSize(::Type{SectorValues{TimeReversed{I}}}) where {I<:Sector}
+    return Base.IteratorSize(values(I))
+end
+function Base.length(::SectorValues{TimeReversed{I}}) where {I<:Sector}
+    return length(values(I))
+end
+function Base.getindex(::SectorValues{TimeReversed{I}}, i::Int) where {I<:Sector}
+    return TimeReversed{I}(getindex(values(I), i))
+end
+function Base.iterate(::SectorValues{TimeReversed{I}}, state...) where {I<:Sector}
+    next = iterate(values(I), state...)
+    if isnothing(next)
+        return nothing
+    else
+        obj, nextstate = next
+        return TimeReversed{I}(obj), nextstate
+    end
+end
+function findindex(::SectorValues{TimeReversed{I}},
+                   a::TimeReversed{I}) where {I<:Sector}
+    return findindex(values(I), a.a)
+end
+
+function Base.isless(c1::TimeReversed{I}, c2::TimeReversed{I}) where {I<:Sector}
+    return isless(c1.a, c2.a)
+end

test/runtests.jl

@@ -20,7 +20,21 @@ const sectorlist = (Z2Irrep, Z3Irrep, Z4Irrep, U1Irrep, CU1Irrep, SU2Irrep, NewS
                     Z3Irrep ⊠ Z4Irrep, FermionParity ⊠ U1Irrep ⊠ SU2Irrep,
                     FermionParity ⊠ SU2Irrep ⊠ SU2Irrep, NewSU2Irrep ⊠ NewSU2Irrep,
                     NewSU2Irrep ⊠ SU2Irrep, FermionParity ⊠ SU2Irrep ⊠ NewSU2Irrep,
-                    FibonacciAnyon ⊠ FibonacciAnyon ⊠ Z2Irrep)
+                    FibonacciAnyon ⊠ FibonacciAnyon ⊠ Z2Irrep,
+                    TimeReversed{Z2Irrep},
+                    TimeReversed{Z3Irrep}, TimeReversed{Z4Irrep},
+                    TimeReversed{U1Irrep}, TimeReversed{CU1Irrep}, TimeReversed{SU2Irrep},
+                    TimeReversed{FibonacciAnyon}, TimeReversed{IsingAnyon},
+                    TimeReversed{FermionParity},
+                    TimeReversed{FermionParity ⊠ FermionParity},
+                    TimeReversed{Z2Irrep ⊠ Z3Irrep ⊠ Z4Irrep},
+                    TimeReversed{Z2Irrep} ⊠ TimeReversed{Z3Irrep} ⊠ TimeReversed{Z4Irrep},
+                    TimeReversed{NewSU2Irrep ⊠ NewSU2Irrep},
+                    TimeReversed{Z2Irrep ⊠ FibonacciAnyon ⊠ FibonacciAnyon},
+                    TimeReversed{NewSU2Irrep ⊠ SU2Irrep},
+                    TimeReversed{FermionParity ⊠ U1Irrep ⊠ SU2Irrep},
+                    TimeReversed{FermionParity ⊠ SU2Irrep ⊠ SU2Irrep},
+                    TimeReversed{FermionParity ⊠ SU2Irrep ⊠ NewSU2Irrep})
 
 @testset "$(TensorKitSectors.type_repr(I))" for I in sectorlist
     @include("sectors.jl")