On an asymptotic method for computing the modified energy for symplectic methods
Numerical Analysis
2014-01-06 v1
Abstract
We revisit an algorithm by Skeel et al. for computing the modified, or shadow, energy associated with the symplectic discretization of Hamiltonian systems. By rephrasing the algorithm as a Richardson extrapolation scheme arbitrary high order of accuracy is obtained, and provided error estimates show that it does capture the theoretical exponentially small drift associated with such discretizations. Several numerical examples illustrate the theory.
Cite
@article{arxiv.1304.0673,
title = {On an asymptotic method for computing the modified energy for symplectic methods},
author = {Per Christian Moan and Jitse Niesen},
journal= {arXiv preprint arXiv:1304.0673},
year = {2014}
}
Comments
16 pages, 7 figures