yapanuwan
/
BergEC-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Merge branch 'main' into berggren

This commit is contained in:
ysyapa 2024-05-31 20:13:38 +00:00
parent 29ed21f035 58e32f4950
commit 59f74b0263
4 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
3body_resonance.jl
View File

@ -17,7 +17,7 @@ E_max = 30
println("No of threads = ", Threads.nthreads())
@time "Basis" begin
Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max)
Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max, Λ)
l_max = max(maximum(l1s), maximum(l2s))
n_max = max(maximum(n1s), maximum(n2s))
mask1 = (n2s .== n2s') .&& (l2s .== l2s')

3
ho_basis.jl
View File

@ -29,7 +29,7 @@ function get_sp_basis(E_max)
return (Es, ns, ls)
end
function get_2p_basis(E_max)
function get_2p_basis(E_max, Λ=-1)
Es = Int[]
n1s = Int[]
l1s = Int[]
@ -42,6 +42,7 @@ function get_2p_basis(E_max)
for n2 in 0 : (E - 2*n1) ÷ 2
for l1 in 0 : (E - 2*n1 - 2*n2)
l2 = E - 2*n1 - 2*n2 - l1
if Λ≥0 && !triangle_ineq(l1, l2, Λ); continue; end
push!(Es, E)
push!(n1s, n1)
push!(l1s, l1)

2
ho_basis_3body.jl
View File

@ -18,7 +18,7 @@ E_max = 40
println("No of threads = ", Threads.nthreads())
@time "Basis" begin
Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max)
Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max, Λ)
mask1 = (n2s .== n2s') .&& (l2s .== l2s')
mask2 = (n1s .== n1s') .&& (l1s .== l1s')
end

2
test/moshinsky.jl
View File

@ -9,7 +9,7 @@ E_max = 30
println("No of threads = ", Threads.nthreads())
@time "Basis" Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max)
@time "Basis" Es, n1s, l1s, n2s, l2s = get_2p_basis(E_max, Λ)
println("Basis size = ", length(Es))