How to regularize a symplectic-energy-momentum integrator
arXiv:math/0507483
Abstract
We identify ghost trajectories of symplectic-energy-momentum (SEM) integration and show that the ghost trajectories are not time reversible. We explain how SEM integration can be regularized, in a SEM preserving manner, so that it is time reversible. We describe an algorithm for implementing the regularized SEM integrator. Simulation results for the pendulum are given. Coordinate invariance of the regularized SEM integrator is briefly discussed.
17 pages, see http://rose-hulman.edu/~shibberu/DTH_Dynamics/DTH_Dynamics.htm