Now Emax = max(l1 + l2 + 2*(n1 + n2))

2024-04-19 15:31:38 -04:00 · 2024-04-19 15:31:38 -04:00 · 4fae06f8cd
parent 43c4a5941c
commit 4fae06f8cd
3 changed files with 3 additions and 3 deletions
--- a/ho_basis.jl
+++ b/ho_basis.jl
@ -35,7 +35,7 @@ function get_2p_basis(E_max)
    l2s = Int[]

    # E = 2*n1 + l1 + 2*n2 + l2
-    for E in 0 : 2 * E_max
+    for E in 0 : E_max
        for n1 in 0 : E ÷ 2
            for n2 in 0 : (E - 2*n1) ÷ 2
                for l1 in 0 : (E - 2*n1 - 2*n2)
--- a/ho_basis_3body.jl
+++ b/ho_basis_3body.jl
@ -2,7 +2,7 @@ using Arpack, SparseArrays
 include("ho_basis.jl")
 include("p_space.jl")

-E_max = 20
+E_max = 40
 μ1ω1_gen = 1/2
 μ2ω2_gen = 2
 Λ = 0
--- a/test/moshinsky.jl
+++ b/test/moshinsky.jl
@ -4,7 +4,7 @@ using LinearAlgebra

 include("../ho_basis.jl")

-E_max = 15
+E_max = 30
 Λ = 0

 println("No of threads = ", Threads.nthreads())