Krylov integrators for Hamiltonian systems (Q1731604)

The authors consider Arnoldi-like processes to obtain symplectic subspaces for Hamiltonian systems. They use local low-dimensional approximations of large systems, with special emphasis on Krylov subspaces and some of their extensions. These subspaces can be used in two ways: either by numerically solving low-dimensional systems and then mapping them back to higher dimension, or for approximating the necessary functions in exponential integrators applied to higher dimensional systems. The latter case is used for linear systems while, in the former case, one can expect excellent energy preservation. The authors consider second-order exponential integrators which solve linear systems exactly, and for which these two approaches are in a certain sense equivalent. The authors also discuss the time symmetry-preservation properties of the integrators. They further present how numerical tests of these methods combined with symplectic subspaces also show promising behavior when applied to nonlinear Hamiltonian problems.