diff --git a/src/basis_embedding.jl b/src/basis_embedding.jl
index b9b28c8..ff7b157 100644
--- 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
index 13744f9..4862bbe 100644
--- 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()