Switch to robust multi-newton to improve initialization robustness

SciML · Oct 19, 2023 · 870629e · 870629e
1 parent d763352
commit 870629e
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/initialize_dae.jl b/src/initialize_dae.jl
@@ -23,7 +23,7 @@ BrownFullBasicInit(abstol) = BrownFullBasicInit(; abstol = abstol, nlsolve = not
 using SciMLNLSolve
 default_nlsolve(alg, isinplace, u, autodiff = false) = alg
 function default_nlsolve(::Nothing, isinplace, u, autodiff = false)
-    NLSolveJL(autodiff = autodiff ? :forward : :central)
+    RobustMultiNewton(autodiff = autodiff)
 end
 function default_nlsolve(::Nothing, isinplace::Val{false}, u::StaticArray, autodiff = false)
     SimpleNewtonRaphson(autodiff = autodiff)