Numerical integration of Hamiltonian problems by G-symplectic methods

Abstract

It is the purpose of this paper to consider the employ of General Linear Methods (GLMs) as geometric numerical solvers for the treatment of Hamiltonian problems. Indeed, even if the numerical flow generated by a GLM cannot be symplectic, we exploit here a concept of near conservation for such methods which, properly combined with other desirable features (such as symmetry and boundedness of parasitic components), allows to achieve an accurate conservation of the Hamiltonian. In this paper we focus our attention on the connection between order of convergence and Hamiltonian deviation by multivalue methods. Moreover, we derive a semi-implicit GLM which results competitive to symplectic Runge-Kutta methods, which are notoriously implicit.