Commutators now in PConv and fix faulty automatic indexing

lukasgrunwald · Mar 20, 2024 · 97dfa6d · 97dfa6d
1 parent 9b311f8
commit 97dfa6d
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 7 deletions.
diff --git a/src/basis_embedding.jl b/src/basis_embedding.jl
@@ -74,8 +74,6 @@ Map operator `op` in Full-Hilbert space into half-filled U(1) subspace.
 function map_u1(op)
     if size(op, 1) == 4^2 # Fermions
         idx = idx_fermi
-    elseif size(op, 1) == 16^2 # 2-band Fermions
-        idx = idx_multi
     else
         error("No automatic mapping implemented! Call function with idx instead.")
     end

diff --git a/test/embedding.jl b/test/embedding.jl
@@ -6,12 +6,8 @@ using LadderOperators
 using LinearAlgebra: I
 using Test
 
-# Convenience functions for checking the functions
-commutator(A, B) = A*B - B * A
-anti_commutator(A, B) = A*B + B * A
-
 # —————————————————————————————————— Spinless Fermions ————————————————————————————————— #
-Nsites = 3
+Nsites = 4
 Cdag, C, Z, N = local_spinless_fermions()
 Cdag_, C_ = ladder_operators(:spinless, Nsites)
 
@@ -32,6 +28,7 @@ N = embed_fock(Nsites, N) # directly translate to fock space
     @test C == C
 end
 
+##
 # —————————————————————————————————— Spinful Fermions —————————————————————————————————— #
 Nsites = 2
 Cdag, C, Z, N = local_fermions()