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

Compare calculated fields against benchmark

This commit is contained in:
Nuwan Yapa 2025-01-28 12:20:30 -05:00
parent de1b7b8092
commit 7822c0f168
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
test/Pb208_flds.jl
View File

@ -17,3 +17,14 @@ divs = length(xs) - 1
plot(xs, hcat(Φ0, W0, B0, A0), layout=4, label=["Φ0" "W0" "B0" "A0"]) plot(xs, hcat(Φ0, W0, B0, A0), layout=4, label=["Φ0" "W0" "B0" "A0"])
xlabel!("r (fm)") xlabel!("r (fm)")
# test data generated from Hartree.f
# format: x S(x) V(x) R(x) A(x)
test_data = readdlm("test/Pb208FldsFSUGarnet.csv")
xs_bench = test_data[:, 1]
Φ0_bench = test_data[:, 2]
W0_bench = test_data[:, 3]
B0_bench = test_data[:, 4]
A0_bench = test_data[:, 5]
plot!(xs_bench, hcat(Φ0_bench, W0_bench, B0_bench, A0_bench), layout=4, label=["Φ0 benchmark" "W0 benchmark" "B0 benchmark" "A0 benchmark"])