Backward error analysis for conjugate symplectic methods
Loading...
Date
2023-12-08
Open Access Location
Authors
Journal Title
Journal ISSN
Volume Title
Publisher
American Institute of Mathematical Sciences (AIMS)
Rights
(c) 2023 The Author/s
CC BY 4.0
CC BY 4.0
Abstract
The numerical solution of an ordinary differential equation can be interpreted as the exact solution of a nearby modified equation. Investigating the behaviour of numerical solutions by analysing the modified equation is known as backward error analysis. If the original and modified equation share structural properties, then the exact and approximate solution share geometric features such as the existence of conserved quantities. Conjugate symplectic methods preserve a modified symplectic form and a modified Hamiltonian when applied to a Hamiltonian system. We show how a blended version of variational and symplectic techniques can be used to compute modified symplectic and Hamiltonian structures. In contrast to other approaches, our backward error analysis method does not rely on an ansatz but computes the structures systematically, provided that a variational formulation of the method is known. The technique is illustrated on the example of symmetric linear multistep methods with matrix coefficients.
Description
Keywords
variational integrators, backward error analysis, Euler–Lagrange equations, multistep methods, conjugate symplectic methods
Citation
McLachlan RI, Offen C. (2023). Backward error analysis for conjugate symplectic methods. Journal of Geometric Mechanics. 15. 1. (pp. 98-115).