diff --git a/src/caches/firk_caches.jl b/src/caches/firk_caches.jl index 536f5bffde..a0b8b46b08 100644 --- 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)