Dense output documentation
Opened this issue · 3 comments
Can you add some more documentation for how to use dense output? Suppose I want to integrate and ODE, and determine the solutions at times [1.0, 2.0, 3.0, 4.0, 5.0, 6.0]. How could this be achieved using dense output interpolation?
Right now I am calling the integrate function, 6 different times to get the solutions at my six desired time points. But I think this is pretty inefficient because, the integrator step size appears to restart every time I call integrate.
There's a better way to do it. You have to use the solout function. I can try to make a simple example...
BTW, I've made a speedy python wrapper to this package: https://github.com/Nicholaswogan/numblsoda
My hope is to make a Python interface to a explicit ODE solver which compares similarly to Julia's explicit methods.
Currently the wrapper is repeated calling integrate, instead of using dense output.
I figured out how to do the dense output stuff: https://github.com/Nicholaswogan/numblsoda/blob/main/src/dop853_c_interface.f90
An example would still be nice for documentation purposes