Mario Mulansky

Do you mean the exception raised when the step size adjustment fails? This has been modularized in f5ee194e7c45ecd5fa503372d7e3a74f620ce16e, but is not yet customizable. However, I'm not sure this is necessary...

is there a difficulty? you just give the time as a member of the device functor in the system function, or am i missing something?

Differential equations of the form M x' = F(x) are, in general, not ODEs - so odeint does not provide solvers for this. If M is constant and non-singular, the...

Well the only way I see is to invert M and treat it as an ODE. You can use Eigen to invert M, for example.

> One could also solve for `y` in `My = F(x)` using an iterative solver. Depending on `M` this could be more memory efficient. Right, but then you would need...

I agree that providing a numerical jacobian is definitely doable and should be added to odeint. We will try to include that at some point. Thanks for pointing that out...

For some reason, we called the type ```result_type```, not ```return_type```. Try using this name in the typedef and see if it works. Otherwise, your approach is correct.

Did you try defining it as a const ref? This should be the right way.

Hi everyone, After some discussion I've revisited this problem and added a short example based on dopri5, find_if and a final bisection: 766d2cb0631b89f8f7f5c8187b792a216e691bd9 Quite straight-forward and shouldnt be too difficult...

Can you build basic CUDA examples (i.e. that dont use odeint)? What version of odeint are you using? The latest from github or some that shipped with boost, which boost...