Session CD - Chaos and Statistical Physics.
MIXED session, Saturday morning, November 02
Room 219, Morton Hall
In this talk, we compare the performance of symplectic integrators and ordinary-differential-equation (ODE) solvers for some low dimensional Hamiltonian systems. Considered here are the symplectic integrators of second-order (Feng, Ruth, Friedman amp; Auerbach), fourth-order (Forest amp; Ruth, Candy amp; Rozmus) and sixth-order (Yoshida), and the ODE solvers of DEABM and DERKF (Shampine) of the SLATEK library and LSODE (Hindmarsh) of ODEPACK. These integrators and solvers are used for trajectory computation of the Holmes oscillator, two-body central force field motion, and Henon-Heiles system. The main observation is that ODE solvers outperform the symplectic integrators in accuracy, but fall short in efficiency.
