Closed vancleve closed 3 weeks ago
yeah, this looks like a simple typo.
I also get the error with alg=Rodas5()
fwiw.
In line 1173 of stiff_addsteps instead of
@.. broadcast=false tmp=uprev + a71 * k1 + a72 * k2 + a73 * k3 + a74 * k4 +
a75 * k5 + a76 * k6 + k7
it must be
@.. broadcast=false tmp .+= k7
I believe that is correct. Fix incoming.
Stacktrace starts here: https://github.com/SciML/OrdinaryDiffEq.jl/blob/7f0eb18c635fdbe7972f1ddb1ce4dab3bd28d79e/src/dense/stiff_addsteps.jl#L1173
Minimish working example
macOS 14.5 Julia 1.10.4 DifferentialEquations v7.13.0