Reconsidering trigonometric integrators

Robert I McLachlan, Dion R O'Neale


In this paper we look at the performance of trigonometric integrators applied to highly oscillatory differential equations. It is widely known that some of the trigonometric integrators suffer from low order resonances for particular step sizes. We show here that, in general, trigonometric integrators also suffer from higher order resonances which can lead to loss of nonlinear stability. We illustrate this with the Fermi-Pasta-Ulam problem, a highly oscillatory Hamiltonian system. We also show that in some cases trigonometric integrators preserve invariant or adiabatic quantities but at the wrong values. We use statistical properties such as time averages to further evaluate the performance of the trigonometric methods and compare the performance with that of the mid-point rule.



numerical integration, initial value problems, highly oscillatory Hamiltonian systems, trigonometric integrators, symplectic algorithms.


Remember, for most actions you have to record/upload into this online system
and then inform the editor/author via clicking on an email icon or Completion button.
ANZIAM Journal, ISSN 1446-8735, copyright Australian Mathematical Society.