Discrete gradient algorithms of high order for one-dimensional systems
Computational Physics
2010-08-24 v1 Numerical Analysis
Numerical Analysis
Abstract
We show how to increase the order of one-dimensional discrete gradient numerical integrator without losing its advantages, such as exceptional stability, exact conservation of the energy integral and exact preservation of the trajectories in the phase space. The accuracy of our integrators is higher by several orders of magnitude as compared with the standard discrete gradient scheme (modified midpoint rule) and, what is more, our schemes have very high accuracy even for large time steps.
Cite
@article{arxiv.1008.3895,
title = {Discrete gradient algorithms of high order for one-dimensional systems},
author = {Jan L. Cieśliński and Bogusław Ratkiewicz},
journal= {arXiv preprint arXiv:1008.3895},
year = {2010}
}
Comments
22 pages, 9 figures