Skip to content

Fix bug impure tensors complex phi4 + added test for impure tensors for real and complex phi4 #168

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
VictorVanthilt merged 4 commits into from
Apr 16, 2026
Merged

Fix bug impure tensors complex phi4 + added test for impure tensors for real and complex phi4 #168

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
e007891
Fix bug impure tensors
JaridPiceu Apr 15, 2026
c1c9c4e
Add test for impure tensor real and complex phi4
JaridPiceu Apr 15, 2026
f80a5c2
Make compatible with Runic
JaridPiceu Apr 15, 2026
9c763a0
Lowering bond dimension in test
JaridPiceu Apr 15, 2026
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
2 changes: 1 addition & 1 deletion src/models/phi4_complex.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -105,7 +105,7 @@ function moment_matrix(N, μ0, λ; rtol = 1.0e-8)
return M
end

function phi4_complex_tensor(f::TensorMap{TT, SS, NN, NN}, weights::Matrix{Float64}; T::Type{<:Number} = ComplexF64) where {TT, SS, NN} # for Trivial symmetry
function phi4_complex_tensor(f::TensorMap{TT, SS, NN, NN}, weights::Matrix{<:Number}; T::Type{<:Number} = ComplexF64) where {TT, SS, NN} # for Trivial symmetry
Comment thread
VictorVanthilt marked this conversation as resolved.
K = size(weights, 1)
N = K^2
perms = collect(permutations(1:4)) # 24 total
Expand Down
53 changes: 53 additions & 0 deletions test/models.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -84,3 +84,56 @@ for (model, temp, answer, description) in model_temp_answer_string_3d
@test free_energy(data, temp; scalefactor = 8.0) ≈ answer rtol = 1.0e-3
end
end


# Test for impurity tensors
# Real phi^4
@testset "Real φ⁴ model - Impure" begin
# Disordered phase
λ = 1.0
μ0 = 0.0

Tpure = phi4_real(Trivial, 10, μ0, λ)
T_imp1 = phi4_real_imp1(Trivial, 10, μ0, λ)
T_imp2 = phi4_real_imp2(Trivial, 10, μ0, λ)

scheme = ImpurityHOTRG(Tpure, T_imp1, T_imp1, T_imp2)

data = run!(scheme, truncrank(16), maxiter(25))

order_para = data[end][4] / data[end][1]
@test order_para ≈ 0.0 atol = 1.0e-3

# Ordered phase
μ0 = -2.0

Tpure = phi4_real(Trivial, 10, μ0, λ)
T_imp1 = phi4_real_imp1(Trivial, 10, μ0, λ)
T_imp2 = phi4_real_imp2(Trivial, 10, μ0, λ)

scheme = ImpurityHOTRG(Tpure, T_imp1, T_imp1, T_imp2)

data = run!(scheme, truncrank(16), maxiter(25))

order_para = data[end][4] / data[end][1]
@test order_para ≈ 1.5317112652447245 rtol = 1.0e-3
end

# Complex phi^4
@testset "Complex φ⁴ model - Impure" begin
# Disordered phase
λ = 1.0
μ0 = 0.0

Tpure = phi4_complex(Trivial, 6, μ0, λ; T = ComplexF64)
T_imp11 = phi4_complex_impϕ(6, μ0, λ)
T_imp12 = phi4_complex_impϕdag(6, μ0, λ)
T_imp2 = phi4_complex_impϕ2(6, μ0, λ)

scheme = ImpurityHOTRG(Tpure, T_imp11, T_imp12, T_imp2)

data = run!(scheme, truncrank(16), maxiter(25))

order_para = data[end][4] / data[end][1]
@test order_para ≈ 0.0 atol = 1.0e-3
end
Loading