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

Moshinsky brackets comparison with Buck et al.

This commit is contained in:
ysyapa 2024-04-16 20:07:05 +00:00
parent d987b02290
commit f956f92546
1 changed files with 25 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
test/buck_et_al.jl Normal file
View File

@ -0,0 +1,25 @@
using NuclearToolkit
# Table 1 from Buck et al.
n1 = [0,0,0,0,0,0,2,2,2,2];
l1 = [0,1,1,2,2,2,2,2,2,2];
n2 = [0,0,0,0,0,0,1,1,1,1];
l2 = [0,3,5,2,4,5,3,3,4,4];
N = [0,0,0,0,1,0,0,1,0,3];
L = [0,2,1,1,3,5,3,0,2,2];
n = [0,1,0,0,0,0,1,2,4,0];
l = [0,0,5,3,1,2,6,5,2,4];
Λ = [0,2,6,4,3,4,4,5,2,4];
Elhs = l1 .+ l2 .+ 2 .* (n1 .+ n2)
Erhs = l .+ L .+ 2 .* (n .+ N)
E_max = max(maximum(Elhs), maximum(Erhs))
j_max = 2 * E_max + 1
l_max = j_max
dtri = NuclearToolkit.prep_dtri(l_max + 1);
dcgm0 = NuclearToolkit.prep_dcgm0(l_max);
d6j = nothing # NuclearToolkit.prep_d6j_int(E_max, j_max, to);
bracs = NuclearToolkit.HObracket_d6j.(N, L, n, l, n1, l1, n2, l2, Λ, 1.0, Ref(dtri), Ref(dcgm0), Ref(d6j))
display(bracs)