fix: use default solver tolerance if not provided in OverrideInit

SciML · Nov 12, 2024 · efaa9f2 · efaa9f2
1 parent 705bf97
commit efaa9f2
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/initialization.jl b/src/initialization.jl
@@ -182,7 +182,11 @@ function get_initial_values(prob, valp, f, alg::OverrideInit,
         initdata.update_initializeprob!(initprob, valp)
     end
 
-    nlsol = solve(initprob, nlsolve_alg; abstol = alg.abstol)
+    if alg.abstol !== nothing
+        nlsol = solve(initprob, nlsolve_alg; abstol = alg.abstol)
+    else
+        nlsol = solve(initprob, nlsolve_alg)
+    end
 
     u0 = initdata.initializeprobmap(nlsol)
     if initdata.initializeprobpmap !== nothing