some units, leftunits and isunits

Author: borisdevos

src/algorithms/excitation/chepigaansatz.jl

@@ -35,10 +35,10 @@ end
 
 function excitations(
         H, alg::ChepigaAnsatz, ψ::FiniteMPS, envs = environments(ψ, H);
-        sector = one(sectortype(ψ)), num::Int = 1, pos::Int = length(ψ) ÷ 2
+        sector = leftunit(ψ), num::Int = 1, pos::Int = length(ψ) ÷ 2
     )
     1 ≤ pos ≤ length(ψ) || throw(ArgumentError("invalid position $pos"))
-    sector == one(sector) || error("not yet implemented for charged excitations")
+    isunit(sector) || error("not yet implemented for charged excitations")
 
     # add random offset to kickstart Krylov process:
     AC = ψ.AC[pos]
@@ -100,10 +100,10 @@ end
 
 function excitations(
         H, alg::ChepigaAnsatz2, ψ::FiniteMPS, envs = environments(ψ, H);
-        sector = one(sectortype(ψ)), num::Int = 1, pos::Int = length(ψ) ÷ 2
+        sector = leftunit(ψ), num::Int = 1, pos::Int = length(ψ) ÷ 2
     )
     1 ≤ pos ≤ length(ψ) - 1 || throw(ArgumentError("invalid position $pos"))
-    sector == one(sector) || error("not yet implemented for charged excitations")
+    isunit(sector) || error("not yet implemented for charged excitations")
 
     # add random offset to kickstart Krylov process:
     @plansor AC2[-1 -2; -3 -4] := ψ.AC[pos][-1 -2; 1] * ψ.AR[pos + 1][1 -4; -3]

src/states/quasiparticle_state.jl

@@ -35,7 +35,7 @@ end
 #constructors
 function LeftGaugedQP(
         datfun, left_gs, right_gs = left_gs;
-        sector = one(sectortype(left_gs)), momentum = 0.0
+        sector = leftunit(left_gs), momentum = 0.0
     )
     # find the left null spaces for the TNS
     excitation_space = Vect[typeof(sector)](sector => 1)
@@ -56,7 +56,7 @@ function LeftGaugedQP(
 end
 function LeftGaugedQP(
         datfun, left_gs::MultilineMPS, right_gs::MultilineMPS = left_gs;
-        sector = one(sectortype(left_gs)), momentum = 0.0
+        sector = leftunit(left_gs), momentum = 0.0
     )
     # not sure why this is needed for type stability
     Tresult = leftgaugedqptype(eltype(parent(left_gs)), typeof(momentum))
@@ -69,7 +69,7 @@ end
 
 function RightGaugedQP(
         datfun, left_gs, right_gs = left_gs;
-        sector = one(sectortype(left_gs)), momentum = 0.0
+        sector = leftunit(left_gs), momentum = 0.0
     )
     # find the left null spaces for the TNS
     excitation_space = Vect[typeof(sector)](sector => 1)
@@ -226,7 +226,7 @@ auxiliarysector(state::QP) = only(sectors(auxiliaryspace(state)))
 eachsite(state::QP) = eachsite(state.left_gs)
 
 istopological(qp::QP) = qp.left_gs !== qp.right_gs
-istrivial(qp::QP) = !istopological(qp) && isone(auxiliarysector(qp))
+istrivial(qp::QP) = !istopological(qp) && isunit(auxiliarysector(qp))
 
 Base.copy(a::QP) = copy!(similar(a), a)
 Base.copyto!(a::QP, b::QP) = copy!(a, b)

src/utility/plotting.jl

@@ -119,7 +119,7 @@ function transferplot end
         sector_formatter = string
     )
     if sectors === nothing
-        sectors = [one(sectortype(h.args[1]))]
+        sectors = [unit(sectortype(h.args[1]))]
     end
 
     for sector in sectors

test/algorithms.jl

@@ -1057,7 +1057,7 @@ module TestAlgorithms
             FiniteMPS(physicalspace(H), maxVspaces[2:(end - 1)]), H; verbosity = 0
         )
         E₀ = expectation_value(gs, H)
-        @test E₀ ≈ first(vals_dense[one(U1Irrep)])
+        @test E₀ ≈ first(vals_dense[unit(U1Irrep)])
 
         for (sector, vals) in vals_dense
             # ED tests
@@ -1074,7 +1074,7 @@ module TestAlgorithms
             Es, Bs = excitations(H, QuasiparticleAnsatz(; tol), gs; sector, num = 1)
             Es = Es .+ E₀
             # first excited state is second eigenvalue if sector is trivial
-            @test Es[1] ≈ vals[isone(sector) ? 2 : 1] atol = 1.0e-8
+            @test Es[1] ≈ vals[isunit(sector) ? 2 : 1] atol = 1.0e-8
         end
 
         # shifted charges tests

test/multifusion.jl

@@ -32,7 +32,7 @@ module TestMultifusion
         P = Vect[I](D0 => 1, D1 => 1)
         X = zeros(T, P ← P)
         for (s, f) in fusiontrees(X)
-            isone(only(f.uncoupled)) ? X[s, f] .= g : X[s, f] .= -g
+            isunit(only(f.uncoupled)) ? X[s, f] .= g : X[s, f] .= -g
         end
         ZZ = zeros(T, P^2 ← P^2)
         for (s, f) in fusiontrees(ZZ)