the rise of computing power has yielded many important insights. There is therefore considerable interest in developing fast and accurate numerical integrators.

A large number of such integrators have been developed over the years to perform this task. For many long term integrations, symplectic integrators have proven to be a favourable choice. Symplectic schemes incorporate the symmetries of Hamiltonian systems, and therefore typically conserve quantities like the energy and angular momentum better than non-symplectic integrators.

For integrations of planetary systems, [@WisdomHolman1991], and independently [@Kinoshita1991], developed a widely used class of symplectic integrators. The ideas of [@WisdomHolman1991] developed from the original ideas of the mapping method of [@Wisdom1981]. We refer to these as a Wisdom-Holman mapping or a Wisdom-Holman integrator. Since then, many authors have modified and built upon this method, and several have made their integrators publicly available to the astrophysics community [e.g. @Chambers1997; @Duncan1998].

The Wisdom-Holman integrator exploits the intuition from perturbation theory that one can separate the problem into a system of Keplerian orbits about the Sun, modified by small perturbations among the planets. The nuisance is that while Newton provided us the solution to the two-body problem, Poincar[é]{} showed that the remaining superimposed perturbations are not integrable. 