Closed henry2004y closed 7 months ago
For example, in demo_proton_electron, if we use
sol_e = solve(prob_e, Vern9(); save_idxs=[1,2,3])
instead of
sol_e = solve(prob_e, Vern9())
The display result would be strange.
Related to https://github.com/SciML/OrdinaryDiffEq.jl/issues/2140
Chris explained that this bug originates from doing lazy interpolations some of the schemes in the above link.
For example, in demo_proton_electron, if we use
instead of
The display result would be strange.
Related to https://github.com/SciML/OrdinaryDiffEq.jl/issues/2140