DifferentialEquations: Sundials and LSODA report convergence failure, Julia solvers silently fail?

[ChrisRackauckas]

@qtoloza

The real problem is that solve() does not throw an error and halt integration when I use a native Julia solver (I’ve tried AutoDP5(Rosenbrock23), Tsit5, Rosenbrock23, Vern6, etc.) - instead, it seemingly continues to integrate forever while accumulating allocations. CVODE_BDF and lsoda both report convergence failure and stop early.

Every solver cuts off slightly differently. Runge-Kutta methods can keep iterating until maxiters, while more stable methods like Rosenbrock23 can very accurately try to get a diverging trajectories. You can define things like unstable_check to customize the divergence checking.