From 220dd0eae36b64a9002b7226051dd640222e52db Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Wed, 29 Jan 2025 18:28:59 -0500 Subject: [PATCH] Unit fix (need to understand) --- mesons.jl | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/mesons.jl b/mesons.jl index 7ac95af..217be65 100644 --- a/mesons.jl +++ b/mesons.jl @@ -48,20 +48,20 @@ end 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) # A0 doesn't need iterations - src_A0 = -sqrt(g2_γ) .* ρ_vp + src_A0 = -g2_γ .* ρ_vp .* ħc A0 = solveKG(0, src_A0, r_max) (Φ0, W0, 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 - @. 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) - @. 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) - @. 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) end