headmyshoulder

No, there is absolutely no difficulty. This question arised on code project and this issue is a reminder to add an example to our docs and tutorials.

Why an iterative solver? If M is small enough, LU might be sufficient and the decomposition can be computed only once. I also think that the precision is better you...

Unfortunately odeint does not provide solvers for non-separable Hamiltonian systems. I think the book about geometrical integration from E. Hairer could help you with more advances methods for non-separable Hamiltonians....

Ok, I see. Btw. maybe the Velocity-Verlet integrator fits your needs. I think it was developed for particle methods and it should be energy conserving (at least in the average).

Moved into feature wishlist in the wiki!

One could use a tuple of tuples. It is not very elegant, but I don't see any other way.

Hi, thank you for your suggestion. We will look on it. But be warned, it will take some time.

Hi Janick, yesterday we discussed a similar problem. I think this kind of application needs a special integrate function which handles the boundaries. This function might be implemented be the...

Hmm, after reading about your problem I think that other solutions might also fit well for you. You might implement a step-size controller proxy which uses a controlled_runge_kutta stepper with...

Hi, sorry there are no news about conditional integrators. It is still on our list. Actually I first want to finalize the iterators before implementing the conditional integrate. The Iterators...