Polynomial map factorization of symplectic maps

Long-term stability studies of nonlinear Hamiltonian systems require symplectic integration algorithms which are both fast and accurate. In this paper, we study a symplectic integration method wherein the symplectic map representing the Hamiltonian system is refactorized using polynomial symplectic maps. This method is analyzed for the three degree of freedom case. Finally, we apply this algorithm to study a large particle storage ring.