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

Unit fix (need to understand)

This commit is contained in:
Nuwan Yapa 2025-01-29 18:28:59 -05:00
parent 4f61e8a51f
commit 220dd0eae3
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
mesons.jl
View File

@ -48,20 +48,20 @@ end
Reference: P. Giuliani, K. Godbey, E. Bonilla, F. Viens, and J. Piekarewicz, Frontiers in Physics 10, (2023)" Reference: P. Giuliani, K. Godbey, E. Bonilla, F. Viens, and J. Piekarewicz, Frontiers in Physics 10, (2023)"
function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3) function solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, r_max, divs, iterations=3)
# A0 doesn't need iterations # A0 doesn't need iterations
src_A0 = -sqrt(g2_γ) .* ρ_vp src_A0 = -g2_γ .* ρ_vp .* ħc
A0 = solveKG(0, src_A0, r_max) A0 = solveKG(0, src_A0, r_max)
(Φ0, W0, B0) = (zeros(1 + divs) for _ in 1:3) (Φ0, W0, B0) = (zeros(1 + divs) for _ in 1:3)
(src_Φ0, src_W0, src_B0) = (zeros(1 + divs) for _ in 1:3) (src_Φ0, src_W0, src_B0) = (zeros(1 + divs) for _ in 1:3)
for _ in 1:iterations for _ in 1:iterations
@. src_Φ0 = g2_s * ((κ/2) * Φ0^2 + (λ/6) * Φ0^3 - (ρ_sp + ρ_sn)) @. src_Φ0 = g2_s * ((κ/2) * Φ0^2 + (λ/6) * Φ0^3 - (ρ_sp + ρ_sn)) * ħc
Φ0 .= solveKG(m_s, src_Φ0, r_max) Φ0 .= solveKG(m_s, src_Φ0, r_max)
@. src_W0 = g2_v * ((ζ/6) * W0^3 + 2 * Λv * B0^2 * W0 - (ρ_vp + ρ_vn)) @. src_W0 = g2_v * ((ζ/6) * W0^3 + 2 * Λv * B0^2 * W0 - (ρ_vp + ρ_vn)) * ħc
W0 .= solveKG(m_ω, src_W0, r_max) W0 .= solveKG(m_ω, src_W0, r_max)
@. src_B0 = g2_ρ * (2 * Λv * W0^2 * B0 - (ρ_vp - ρ_vn) / 2) @. src_B0 = g2_ρ * (2 * Λv * W0^2 * B0 - (ρ_vp - ρ_vn) / 2) * ħc
B0 .= solveKG(m_ρ, src_B0, r_max) B0 .= solveKG(m_ρ, src_B0, r_max)
end end