Dense output documentation

[Nicholaswogan]

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.

[jacobwilliams]

There's a better way to do it. You have to use the solout function. I can try to make a simple example...

[Nicholaswogan]

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.

[Nicholaswogan]

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