Hi all, The DifferentialEquations’ suggested integration methods when trying to solve this stiff problem all yield the same error. The native Julia ones, that is. Sundial’s CVODE_BDF as well as lsoda are able to solve this with no problems. My code: using DifferentialEquations using Plots gr() us…

Your issue is that you’re not handling autodiff’d values correctly. [image] Andre_Mello: y = zeros(102) That assumes that all of your values will be Float64, which isn’t true during autodiff. So either turn autodiff off like Rosenbrock23(autodiff=false), or make sure your types are correct, …

Error while using DifferentialEquations.jl stiff methods

Specific Domains Numerics

Andre_Mello March 18, 2019, 11:10pm 3

I see!
Thanks @ChrisRackauckas, I’ll be sure to keep my typing in mind from now on.

Topic		Replies	Views
Strange Behaviour of Stiff Solvers in DifferentialEquations.jl Numerics diffeq	2	676	November 23, 2019
Solver for stiff differential equations New to Julia question	13	1367	May 10, 2022
(Stiffness?) Problem with Differential Equation System New to Julia diffeq	11	2586	November 25, 2018
ODE solving is fast on DifferentialEquation@7.9.0 but significantly slower on later versions Numerics package , differentialequation	10	437	January 6, 2025
Domain Troubles with Differential Equations Solver (with non-auto-differentiable function inside the ODE def) Numerics question , diffeq	33	3635	August 19, 2018

Error while using DifferentialEquations.jl stiff methods

Related topics