From 1c195cd5da5dcd2e071c44e906aec4b33f5e518f Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Fri, 21 Feb 2025 10:24:43 -0500 Subject: [PATCH] Minor change in a test --- test/Pb208_nucleon_dens.jl | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/test/Pb208_nucleon_dens.jl b/test/Pb208_nucleon_dens.jl index 9c1c277..63575fe 100644 --- a/test/Pb208_nucleon_dens.jl +++ b/test/Pb208_nucleon_dens.jl @@ -21,15 +21,15 @@ s.W0 = Vs s.B0 = Rs s.A0 = As -(ρ_sp, ρ_vp) = solveNucleonDensity(true, s) +(s.ρ_sp, s.ρ_vp) = solveNucleonDensity(true, s) -p_sp = plot(rs(s), ρ_sp, xlabel="r (fm)", label="ρₛₚ(r) calculated") -p_vp = plot(rs(s), ρ_vp, xlabel="r (fm)", label="ρᵥₚ(r) calculated") +p_sp = plot(rs(s), s.ρ_sp, xlabel="r (fm)", label="ρₛₚ(r) calculated") +p_vp = plot(rs(s), s.ρ_vp, xlabel="r (fm)", label="ρᵥₚ(r) calculated") -(ρ_sn, ρ_vn) = solveNucleonDensity(false, s) +(s.ρ_sn, s.ρ_vn) = solveNucleonDensity(false, s) -p_sn = plot(rs(s), ρ_sn, xlabel="r (fm)", label="ρₛₙ(r) calculated") -p_vn = plot(rs(s), ρ_vn, xlabel="r (fm)", label="ρᵥₙ(r) calculated") +p_sn = plot(rs(s), s.ρ_sn, xlabel="r (fm)", label="ρₛₙ(r) calculated") +p_vn = plot(rs(s), s.ρ_vn, xlabel="r (fm)", label="ρᵥₙ(r) calculated") # benchmark data generated from Hartree.f # format: x Rhos(n) Rhov(n) Rhot(n) Rhos(p) Rhov(p) Rhot(p)