From 3e6c28ef17826371b214b4ce9652876ae4373f02 Mon Sep 17 00:00:00 2001 From: Nuwan Yapa Date: Mon, 10 Feb 2025 20:30:33 -0500 Subject: [PATCH] Dirty fix for stability --- system.jl | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/system.jl b/system.jl index ac61ab8..89d23ec 100644 --- a/system.jl +++ b/system.jl @@ -46,7 +46,7 @@ function solve_system(s::system, initial_dens=nothing, initial_flds=(zeros(1 + s E_total_previous = NaN while true - (Φ0s, W0s, B0s, A0s) = solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, s.r_max, s.divs, 50; initial_sol = (Φ0s, W0s, B0s, A0s)) + (Φ0s, W0s, B0s, A0s) = solveMesonWfs(ρ_sp, ρ_vp, ρ_sn, ρ_vn, s.r_max, s.divs, isnan(E_total_previous) ? 50 : 5; initial_sol = (Φ0s, W0s, B0s, A0s)) S_interp = linear_interpolation(rs(s), Φ0s) V_interp = linear_interpolation(rs(s), W0s)