Fourth-order symplectic integration
Physica D
A symplectic integration algorithm for separable Hamiltonian functions
Journal of Computational Physics
Almost Poisson integration of rigid body systems
Journal of Computational Physics
Local path fitting: A new approach to variational integrators
Journal of Computational and Applied Mathematics
Hi-index | 31.45 |
Symplectic Runge-Kutta schemes for the integration of general Hamiltonian systems are implicit. In practice, one has to solve the implicit algebraic equations using some iterative approximation method, in which case the resulting integration scheme is no longer symplectic. In this paper, the preservation of the symplectic structure is analyzed under two popular approximation schemes, fixed-point iteration and Newton's method, respectively. Error bounds for the symplectic structure are established when N fixed-point iterations or N iterations of Newton's method are used. The implications of these results for the implementation of symplectic methods are discussed and then explored through numerical examples. Numerical comparisons with non-symplectic Runge-Kutta methods and pseudo-symplectic methods are also presented.