From 870629e29ecd9af5cc5361db9b0fe1fae0225d05 Mon Sep 17 00:00:00 2001
From: Christopher Rackauckas <accounts@chrisrackauckas.com>
Date: Thu, 19 Oct 2023 06:33:48 -0400
Subject: [PATCH] Switch to robust multi-newton to improve initialization
 robustness

---
 src/initialize_dae.jl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/initialize_dae.jl b/src/initialize_dae.jl
index af4f96ceae..37df0a6ebf 100644
--- 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)