remove warning, add error in lua for safety

fastpm · Jan 8, 2024 · c90d61a · c90d61a
1 parent a3273eb
commit c90d61a
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 1 deletion.
diff --git a/libfastpm/cosmology.c b/libfastpm/cosmology.c
@@ -356,7 +356,6 @@ static ode_soln growth_ode_solve(double a, FastPMCosmology * c)
         if (a >= aini) {
             fastpm_raise(-1, "Growth ODE unsuccesful at a=%g.", a);
         } else {
-            fastpm_info("WARNING: Growth ODE unsuccesful at a=%g. Returning growth of 0.", a);
             soln.y0 = 0;
             soln.y1 = 0;
             soln.y2 = 0;

diff --git a/src/lua-runtime-fastpm.lua b/src/lua-runtime-fastpm.lua
@@ -70,6 +70,13 @@ schema.growth_mode.choices = {
     ODE = 'FASTPM_GROWTH_MODE_ODE',
 }
 
+-- enforce a<0.00625 (hard coded) when using ODE growth mode
+function schema.time_step.action (time_step)
+    if time_step[1] < 0.00625 then
+        error("Cannot start the simulation at a<0.00625 when growth_mode=='ODE'.")
+    end
+end
+
 -- enforce Omega_m
 function schema.omega_m.action (value)
     if value ~= nil then