Interested in numerically solving non-separable Hamiltonians. The Wolfram documentation states that the symplectic schemes can only be implemented for separable Hamiltonians of the form

$$ H(p,q) = T(p) + V(q).$$

Even this, I believe has some bugs which are discussed here. Now my thoughts are to simply let Mathematica do the leg work for e.g. generating the coefficients of the butcher tableau etc. etc. and then use this in my "own" designed integrator. Basically sidestepping the limitation of checking whether a Hamiltonian is separable or not...


  • Is there an alternative method to what I have suggested?

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.