From 7822c0f168da632ea12cd9842833af1b8f707545 Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Tue, 28 Jan 2025 12:20:30 -0500 Subject: [PATCH] Compare calculated fields against benchmark --- test/Pb208_flds.jl | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/test/Pb208_flds.jl b/test/Pb208_flds.jl index caa026b..b94e593 100644 --- a/test/Pb208_flds.jl +++ b/test/Pb208_flds.jl @@ -17,3 +17,14 @@ divs = length(xs) - 1 plot(xs, hcat(Φ0, W0, B0, A0), layout=4, label=["Φ0" "W0" "B0" "A0"]) 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"])