Add safety zeroing in RadauIIA5 non fully defined Jacobians

More generally fixes SciML/DifferentialEquations.jl#1001
SciML · Nov 26, 2023 · 680301f · 680301f
1 parent e11d38a
commit 680301f
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/caches/firk_caches.jl b/src/caches/firk_caches.jl
@@ -239,6 +239,7 @@ function alg_cache(alg::RadauIIA5, u, rate_prototype, ::Type{uEltypeNoUnits},
         error("Non-concrete Jacobian not yet supported by RadauIIA5.")
     end
     W2 = similar(J, Complex{eltype(W1)})
+    recursivefill!(W2, false)
 
     du1 = zero(rate_prototype)