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

Moshinsky partially optimized

This commit is contained in:
Nuwan Yapa 2024-03-06 14:27:02 -05:00
parent a069bef7f6
commit 0ed7a1f25d
1 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
ho_basis.jl
View File

@ -76,15 +76,15 @@ end
get_H_matrix(V_l, ns, ls) = get_T_matrix(ns, ls) + get_V_matrix(V_l, ns, ls) get_H_matrix(V_l, ns, ls) = get_T_matrix(ns, ls) + get_V_matrix(V_l, ns, ls)
function Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ) function Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)
l_max = 2*max(maximum(l1s), maximum(l2s)) # OPTIMIZE E_max = maximum(Es)
E_max = maximum(Es) # OPTIMIZE j_max = 2 * E_max + 1
j_max = l_max # OPTIMIZE l_max = j_max
to = 0 # unused to = 0 # unused
dtri = NuclearToolkit.prep_dtri(l_max) ; println("dtri prepared") dtri = NuclearToolkit.prep_dtri(l_max + 1);
dcgm0 = NuclearToolkit.prep_dcgm0(l_max) ; println("dcgm0 prepared") dcgm0 = NuclearToolkit.prep_dcgm0(l_max);
d6j = NuclearToolkit.prep_d6j_int(E_max, j_max, to) ; println("d6j prepared") d6j = NuclearToolkit.prep_d6j_int(E_max, j_max, to);
mat = zeros(length(Es), length(Es)) mat = zeros(length(Es), length(Es))
s = hcat(Es, n1s, l1s, n2s, l2s) s = hcat(Es, n1s, l1s, n2s, l2s)
for idx in CartesianIndices(mat) for idx in CartesianIndices(mat)