Skip to content

Add tests #12

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
merged 8 commits into from
Apr 10, 2025
Merged

Add tests #12

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
8 commits
Select commit Hold shift + click to select a range
1784a62
add tests on exact diagonalisation
sanderdemeyer Apr 1, 2025
092b7be
error handling Hubbard model
sanderdemeyer Apr 1, 2025
dc9c39c
add tJ test
sanderdemeyer Apr 1, 2025
4926222
fix formatting
sanderdemeyer Apr 1, 2025
39505ef
increase patch coverage
sanderdemeyer Apr 1, 2025
25a4075
increase patch coverage some more
sanderdemeyer Apr 1, 2025
8cdb0f4
Update src/hubbardoperators.jl
sanderdemeyer Apr 10, 2025
ef3e7c4
change to @test_throws
sanderdemeyer Apr 10, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 6 additions & 6 deletions src/hubbardoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -245,7 +245,7 @@ function u_plus_u_min(T, ::Type{Trivial}, ::Type{U1Irrep})
return t
end
function u_plus_u_min(T, ::Type{Trivial}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`u_plus_u_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
function u_plus_u_min(T, ::Type{U1Irrep}, ::Type{Trivial})
t = two_site_operator(T, U1Irrep, Trivial)
Expand All @@ -266,7 +266,7 @@ function u_plus_u_min(T, ::Type{U1Irrep}, ::Type{U1Irrep})
return t
end
function u_plus_u_min(T, ::Type{U1Irrep}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`u_plus_u_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
function u_plus_u_min(T, ::Type{SU2Irrep}, ::Type{Trivial})
return error("Not implemented")
Expand All @@ -275,7 +275,7 @@ function u_plus_u_min(T, ::Type{SU2Irrep}, ::Type{U1Irrep})
return error("Not implemented")
end
function u_plus_u_min(T, ::Type{SU2Irrep}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`u_plus_u_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
const u⁺u⁻ = u_plus_u_min

Expand Down Expand Up @@ -305,7 +305,7 @@ function d_plus_d_min(T, ::Type{Trivial}, ::Type{U1Irrep})
return t
end
function d_plus_d_min(T, ::Type{Trivial}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`d_plus_d_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
function d_plus_d_min(T, ::Type{U1Irrep}, ::Type{Trivial})
t = two_site_operator(T, U1Irrep, Trivial)
Expand All @@ -326,7 +326,7 @@ function d_plus_d_min(T, ::Type{U1Irrep}, ::Type{U1Irrep})
return t
end
function d_plus_d_min(T, ::Type{U1Irrep}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`d_plus_d_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
function d_plus_d_min(T, ::Type{SU2Irrep}, ::Type{Trivial})
return error("Not implemented")
Expand All @@ -335,7 +335,7 @@ function d_plus_d_min(T, ::Type{SU2Irrep}, ::Type{U1Irrep})
return error("Not implemented")
end
function d_plus_d_min(T, ::Type{SU2Irrep}, ::Type{SU2Irrep})
return error("Not implemented")
throw(ArgumentError("`d_plus_d_min` is not symmetric under `SU2Irrep` spin symmetry"))
end
const d⁺d⁻ = d_plus_d_min

Expand Down
27 changes: 25 additions & 2 deletions test/bosonoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6,9 +6,9 @@ using .TensorKitTensorsTestSetup
using TensorKitTensors.BosonOperators
using StableRNGs

const cutoff = 4

@testset "Non-symmetric bosonic operators" begin
cutoff = 4

# inferrability
A = @inferred a(; cutoff)
A⁺ = @inferred a⁺(; cutoff)
Expand Down Expand Up @@ -42,6 +42,8 @@ const cutoff = 4
end

@testset "U1-symmetric bosonic operators" begin
cutoff = 4

rng = StableRNG(123)
# inferrability
N = @inferred n(U1Irrep; cutoff)
Expand All @@ -66,3 +68,24 @@ end

test_operator(O_u1, O_triv; L)
end

@testset "Exact Diagonalization" begin
cutoff = 1
L = 2
for symmetry in [Trivial U1Irrep]
rng = StableRNG(123)
# inferrability
N = @inferred n(U1Irrep; cutoff)
A⁺A = @inferred a⁺a(U1Irrep; cutoff)
AA⁺ = @inferred aa⁺(U1Irrep; cutoff)
V = @inferred boson_space(U1Irrep; cutoff)

a_pm, a_mp, a_n = rand(rng, 3)
O = (N ⊗ id(V) + id(V) ⊗ N) * a_n + AA⁺ * a_mp + A⁺A * a_pm

true_eigenvals = sort([0, 2 * a_n, a_n + sqrt(a_mp * a_pm),
a_n - sqrt(a_mp * a_pm)])
eigenvals = expanded_eigenvalues(O; L)
@test eigenvals ≈ true_eigenvals
end
end
6 changes: 3 additions & 3 deletions test/fermionoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,9 +37,9 @@ end
H = -t * (c⁻c⁺() + c⁺c⁻()) +
V * ((n() - 0.5 * id(pspace)) ⊗ (n() - 0.5 * id(pspace))) -
0.5 * mu * (n() ⊗ id(pspace) + id(pspace) ⊗ n())
# Values based on https://arxiv.org/abs/1610.05003v1. Half-Chain Entanglement Entropy in the One-Dimensional Spinless Fermion Model
true_eigenvalues = sort([V / 4, V / 4 - mu, -V / 4 - mu / 2 + t, -V / 4 - mu / 2 - t])

# Values based on https://arxiv.org/abs/1610.05003v1. Half-Chain Entanglement Entropy in the One-Dimensional Spinless Fermion Model
true_eigenvals = sort([V / 4, V / 4 - mu, -V / 4 - mu / 2 + t, -V / 4 - mu / 2 - t])
eigenvals = expanded_eigenvalues(H; L)
@test eigenvals ≈ true_eigenvalues
@test eigenvals ≈ true_eigenvals
end
30 changes: 30 additions & 0 deletions test/hubbardoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,10 +64,15 @@ end
d_num(particle_symmetry, spin_symmetry) ≈
d_num(particle_symmetry, spin_symmetry) *
u_num(particle_symmetry, spin_symmetry)
else
@test_throws ArgumentError u_plus_u_min(particle_symmetry, spin_symmetry)
@test_throws ArgumentError d_plus_d_min(particle_symmetry, spin_symmetry)
end
else
@test_broken c_plus_c_min(particle_symmetry, spin_symmetry)
@test_broken c_min_c_plus(particle_symmetry, spin_symmetry)
@test_broken d_plus_d_min(particle_symmetry, spin_symmetry)
@test_broken u_plus_u_min(particle_symmetry, spin_symmetry)
end
end
end
Expand Down Expand Up @@ -114,3 +119,28 @@ end
@test vals_triv ≈ vals_symm
end
end

@testset "Exact diagonalisation" begin
for particle_symmetry in [Trivial, U1Irrep, SU2Irrep],
spin_symmetry in [Trivial, U1Irrep, SU2Irrep]

if (particle_symmetry, spin_symmetry) in implemented_symmetries
rng = StableRNG(123)

L = 2
t, U = rand(rng, 5)
mu = 0.0
E⁻ = U / 2 - sqrt((U / 2)^2 + 4 * t^2)
E⁺ = U / 2 + sqrt((U / 2)^2 + 4 * t^2)
H_triv = hubbard_hamiltonian(particle_symmetry, spin_symmetry; t, U, mu, L)

# Values based on https://arxiv.org/pdf/0807.4878. Introduction to Hubbard Model and Exact Diagonalization
true_eigenvals = sort(vcat(repeat([-t], 2), [E⁻], repeat([0], 4),
repeat([t], 2),
repeat([U - t], 2), [U], [E⁺], repeat([U + t], 2),
[2 * U]))
eigenvals = expanded_eigenvalues(H_triv; L)
@test eigenvals ≈ true_eigenvals
end
end
end
17 changes: 17 additions & 0 deletions test/spinoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -125,3 +125,20 @@ end
S_exchange(; spin) * a_exchange
test_operator(O_u1, O_triv; L)
end

@testset "Exact diagonalisation for $sector symmetry" for sector in [Trivial U1Irrep]
spin = 1

ZZ = @inferred S_zz(sector; spin)
plusmin = @inferred S_plusmin(sector; spin)
minplus = @inferred S_minplus(sector; spin)
O = ZZ + 0.5 * (plusmin + minplus)

true_eigenvals = vcat([-2.0], repeat([-1.0], 3), repeat([1.0], 5))
eigenvals = expanded_eigenvalues(O; L=2)
@test eigenvals ≈ true_eigenvals

# Value based on https://doi.org/10.1088/0953-8984/2/26/010. Exact diagonalisations of open spin-1 chains
eigenvals = expanded_eigenvalues(O; L=4)
@test eigenvals[2] - eigenvals[1] ≈ 0.509170 atol = 1e-6
end
69 changes: 51 additions & 18 deletions test/tjoperators.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -53,14 +53,18 @@ end
@test u_plus_u_min(particle_symmetry, spin_symmetry; slave_fermion)' ≈
u_min_u_plus(particle_symmetry, spin_symmetry; slave_fermion)
else
@test_broken d_plus_d_min(particle_symmetry, spin_symmetry;
slave_fermion)
@test_broken d_min_d_plus(particle_symmetry, spin_symmetry;
slave_fermion)
@test_broken u_plus_u_min(particle_symmetry, spin_symmetry;
slave_fermion)
@test_broken u_min_u_plus(particle_symmetry, spin_symmetry;
slave_fermion)
@test_throws ArgumentError d_plus_d_min(particle_symmetry,
spin_symmetry;
slave_fermion)
@test_throws ArgumentError d_min_d_plus(particle_symmetry,
spin_symmetry;
slave_fermion)
@test_throws ArgumentError u_plus_u_min(particle_symmetry,
spin_symmetry;
slave_fermion)
@test_throws ArgumentError u_min_u_plus(particle_symmetry,
spin_symmetry;
slave_fermion)
end

# test number operator
Expand All @@ -77,8 +81,10 @@ end
c_num_hole(particle_symmetry, spin_symmetry; slave_fermion) +
c_num(particle_symmetry, spin_symmetry; slave_fermion)
else
@test_broken u_num(particle_symmetry, spin_symmetry; slave_fermion)
@test_broken d_num(particle_symmetry, spin_symmetry; slave_fermion)
@test_throws ArgumentError u_num(particle_symmetry, spin_symmetry;
slave_fermion)
@test_throws ArgumentError d_num(particle_symmetry, spin_symmetry;
slave_fermion)
end

# test spin operator
Expand All @@ -88,11 +94,12 @@ end
d_min_u_min(particle_symmetry, spin_symmetry; slave_fermion)) /
sqrt(2)
else
@test_broken c_singlet(particle_symmetry, spin_symmetry; slave_fermion)
@test_broken u_min_d_min(particle_symmetry, spin_symmetry;
slave_fermion)
@test_broken d_min_u_min(particle_symmetry, spin_symmetry;
slave_fermion)
@test_throws ArgumentError c_singlet(particle_symmetry, spin_symmetry;
slave_fermion)
@test_throws ArgumentError u_min_d_min(particle_symmetry, spin_symmetry;
slave_fermion)
@test_throws ArgumentError d_min_u_min(particle_symmetry, spin_symmetry;
slave_fermion)
end

if spin_symmetry == Trivial
Expand Down Expand Up @@ -166,12 +173,11 @@ function tjhamiltonian(particle_symmetry, spin_symmetry; t, J, mu, L, slave_ferm
end

@testset "spectrum" begin
rng = StableRNG(123)
L = 4
t = randn()
J = randn()
mu = randn()

for slave_fermion in (false, true)
t, J, mu = rand(rng, 3)
H_triv = tjhamiltonian(Trivial, Trivial; t, J, mu, L, slave_fermion)
vals_triv = mapreduce(vcat, eigvals(H_triv)) do (c, v)
return repeat(real.(v), dim(c))
Expand Down Expand Up @@ -199,3 +205,30 @@ end
end
end
end

@testset "Exact Diagonalisation" begin
rng = StableRNG(123)
L = 2

for particle_symmetry in [Trivial, U1Irrep],
spin_symmetry in [Trivial, U1Irrep, SU2Irrep]

for slave_fermion in (false, true)
if (particle_symmetry, spin_symmetry) in implemented_symmetries
t, J = rand(rng, 2)
num = c_num(particle_symmetry, spin_symmetry; slave_fermion)
H = (-t) *
(c_plus_c_min(particle_symmetry, spin_symmetry; slave_fermion) +
c_min_c_plus(particle_symmetry, spin_symmetry; slave_fermion)) +
J *
(S_exchange(particle_symmetry, spin_symmetry; slave_fermion) -
(1 / 4) * (num ⊗ num))

true_eigenvals = sort(vcat([-J], repeat([-t], 2), repeat([t], 2),
repeat([0.0], 4)))
eigenvals = expanded_eigenvalues(H; L)
@test eigenvals ≈ true_eigenvals
end
end
end
end