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	672	November 23, 2019
Domain Troubles with Differential Equations Solver (with non-auto-differentiable function inside the ODE def) Numerics question , diffeq	33	3576	August 19, 2018
(Stiffness?) Problem with Differential Equation System New to Julia diffeq	11	2529	November 25, 2018
DifferentialEquations -- what goes wrong? General Usage	14	2435	November 13, 2019
ERROR: InexactError: Int64(Int64, 1//1000000) New to Julia diffeq	9	6592	November 8, 2018

Error while using DifferentialEquations.jl stiff methods

Related topics