Conjugate-symplecticity properties of Euler-Maclaurin methods and their implementation on the infinity computer (Q1989370)

This paper studies the conservation properties of Euler-Maclaurin formulae and their implementation on the Infinity computer [\textit{Y.D. Sergeyev}, Optim. Lett. 5, No. 4, 575--585 (2011; Zbl 1230.65028)] to solve canonical Hamiltonian dynamical systems. The authors show that the Euler-MacLaurin method of order \(p\) is conjugate-symplectic up to order \(p + 2\). This property makes them suitable for integrating canonical Hamiltonian systems over long times. The new approach is able to work with black-box representations of the functions to compute the exact higher order derivatives using numerical infinities and infinitesimals. The new method belongs to a family of multi-derivative one-step methods containing the classical trapezoidal method as seed technique. The Euler-Maclaurin formulae of order four and six are applied to the nonlinear pendulum, the Kepler problem, the Fermi-Pasta-Ulam problem and to a non-separable Hamiltonian problem. The long-time behavior of their numerical solutions is compared with that of the numerical solutions computed by the symplectic Gauss methods of order four and six.