enckehh : an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour

We present enckehh, a new, highly accurate code for orbital dynamics of perturbed Keplerian systems such as planetary systems or galactic centre systems. It solves Encke’s equations of motion, which assume perturbed Keplerian orbits. By incorporating numerical techniques, we have made the code follow optimal roundoff error growth. In a 1012 d integration of the outer Solar system, enckehh was 3.5 orders of magnitude more accurate than ias15 in a fixed timestep test. Adaptive steps are recommended for ias15. Through study of efficiency plots, we show that enckehh reaches significantly higher accuracy than the rebound integrators ias15 and whckl for fixed step size.