Authors: Paola Console Ernst Hairer Christian Lubich
Publish Date: 2013/02/15
Volume: 124, Issue: 3, Pages: 517-539
Abstract
A method of choice for the longtime integration of constrained Hamiltonian systems is the Rattle algorithm It is symmetric symplectic and nearly preserves the Hamiltonian but it is only of order two and thus not efficient for high accuracy requirements In this article we prove that certain symmetric linear multistep methods have the same qualitative behavior and can achieve an arbitrarily high order with a computational cost comparable to that of the Rattle algorithm
Keywords: