Christopher Rackauckas

https://github.com/wildart/Evolutionary.jl/issues/16

No, it hasn't been updated still?

Yes, this makes sense. I am a little worried about compile times, but maybe it all just quickly compiles away.

We can also hack it with getproperty overloading

@mschauer has a goldmine of papers from which we can pull from: http://pub.math.leidenuniv.nl/~schauermr/

https://www.tandfonline.com/doi/full/10.1080/17442508.2017.1381097

Have you tried a direct check on the gradient defined by autodiff? We should add a test that it's computing the right thing in that case since it's quite nested....

The biggest issue, the reason for the two integrators, is because during the iterations it will need to be disconnected so that way you're building the new interval solution while...

>would still be interesting to restrict the dense solution to certain indices That could be done by allowing a choice of indices for which the `k` is then copied over...

If we wanted, we can make `(du,u,integrator,t) -> f(du,u,integrator,p,t)` or even build the history function there, so the interface is similar but using the OrdinaryDiffEq integrator passing instead of enclosing.