yapanuwan
/
NuclearRMF
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Missed change

This commit is contained in:
Nuwan Yapa 2025-02-18 17:11:02 -05:00
parent b7dd47eab8
commit b04b41373b
1 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
test/Linear_nucleon_spectrum.jl
View File

@ -10,20 +10,20 @@ Vs = test_data[:, 3]
Rs = test_data[:, 4] Rs = test_data[:, 4]
As = test_data[:, 5] As = test_data[:, 5]
S_interp = linear_interpolation(xs, Ss)
V_interp = linear_interpolation(xs, Vs)
R_interp = linear_interpolation(xs, Rs)
A_interp = linear_interpolation(xs, As)
p = false p = false
r_max = maximum(xs) r_max = maximum(xs)
divs = length(xs) - 1 divs = length(xs) - 1
s = system(8, 8, r_max, divs) s = system(8, 8, r_max, divs)
s.Φ0 = Ss
s.W0 = Vs
s.B0 = Rs
s.A0 = As
E_min = 860 E_min = 860
E_max = 939 E_max = 939
κs, Es = findAllOrbitals(p, S_interp, V_interp, R_interp, A_interp, s, E_min, E_max) κs, Es = findAllOrbitals(p, s, E_min, E_max)
Ebinds = (p ? M_p : M_n) .- Es Ebinds = (p ? M_p : M_n) .- Es
occ = fillNucleons(Z_or_N(s, p), κs, Es) occ = fillNucleons(Z_or_N(s, p), κs, Es)