I am using OrdinaryDiffEq with Sundials.
While other options such as
- stability_limit_detect
- max_hnil_warns
- max_error_test_failures
seem to be working, the option maxiters doesn’t seem to stop iteration at all; The iteration even exceeds the default limits which are 500 for sundials and 1e5 for the Julia interface.
I have made an example:
using OrdinaryDiffEq, Sundials
function osc(du,u,p,t)
du[1] = u[2]
du[2] = -u[1]
end
u0 = [1.0;0.0]
tspan = (0.0,1e8)
prob = ODEProblem(osc,u0,tspan)
sol = solve(prob,CVODE_BDF(), maxiters=10)
sol.t|>length
this returns
2086564 (>1e5)
on my PC.
the version of packages are
Sundials v4.2.3
OrdinaryDiffEq v5.40.0
on
Julia v"1.4.2"
I know there is a similar issue without a solution yet. I wonder if there is a work-around for this problem because my original problem takes too much time using other solvers.
https://github.com/SciML/Sundials.jl/issues/269