Global modified Hamiltonian for constrained symplectic integrators
Hairer, Ernst
In: Numerische Mathematik, 2003, vol. 95, no. 2, p. 325-336
Ajouter à la liste personnelle- Summary
- Summary.: We prove that the numerical solution of partitioned Runge-Kutta methods applied to constrained Hamiltonian systems (e.g., the Rattle algorithm or the Lobatto IIIA-IIIB pair) is formally equal to the exact solution of a constrained Hamiltonian system with a globally defined modified Hamiltonian. This property is essential for a better understanding of their longtime behaviour. As an illustration, the equations of motion of an unsymmetric top are solved using a parameterization with Euler parameters