A Very Fast Symplectic Integrator for Planetary Systems: WHFast512 in x86 Assembly

Abstract

Long-term integrations of planetary systems are used in a variety of astrophysical applications, for example to determine a system's stability. In this paper, we describe significant speed and accuracy improvements to the WHFast512 integrator. By completely rewriting WHFast512 in x86 assembly, we gain fine-grained control over CPU instructions and are able to keep the entire simulation state in CPU registers, only writing to memory when an output is required. We also reduce the number of branche…