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

Prettier monitoring

This commit is contained in:
Nuwan Yapa 2025-02-10 18:58:36 -05:00
parent 664e38d624
commit aeb3225254
1 changed files with 11 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
system.jl
View File

@ -39,6 +39,10 @@ function solve_system(s::system, initial_dens=nothing, initial_flds=(zeros(1 + s
(Φ0s, W0s, B0s, A0s) = initial_flds (Φ0s, W0s, B0s, A0s) = initial_flds
if monitor_plot
p = plot(legends=false, layout=(2, 2), title=["ρ_sp" "ρ_vp" "ρ_sn" "ρ_vn"])
end
E_total_previous = NaN E_total_previous = NaN
while true while true
@ -59,7 +63,13 @@ function solve_system(s::system, initial_dens=nothing, initial_flds=(zeros(1 + s
occs_n = fillNucleons(s.N, κs_n, Es_n) occs_n = fillNucleons(s.N, κs_n, Es_n)
(ρ_sp, ρ_vp) = calculateNucleonDensity(κs_n, Es_n, occs_n, false, S_interp, V_interp, R_interp, A_interp, s.r_max, s.divs) (ρ_sp, ρ_vp) = calculateNucleonDensity(κs_n, Es_n, occs_n, false, S_interp, V_interp, R_interp, A_interp, s.r_max, s.divs)
monitor_plot && plot(rs(s), hcat(ρ_sp, ρ_vp, ρ_sn, ρ_vn), layout=(2, 2), title=["ρ_sp" "ρ_vp" "ρ_sn" "ρ_vn"]) |> display if monitor_plot
for s in p.series_list
s.plotattributes[:linecolor] = :gray
end
plot!(p, rs(s), hcat(ρ_sp, ρ_vp, ρ_sn, ρ_vn), linecolor=:red)
display(p)
end
# total binding energy of nucleons # total binding energy of nucleons
E_total = sum(M_p .- Es_p) + sum(M_n .- Es_n) E_total = sum(M_p .- Es_p) + sum(M_n .- Es_n)