WHFast512: A symplectic N-body integrator for planetary systems optimized with AVX512 instructions

Pejvak Javaheri,Hanno Rein,Dan Tamayo
DOI: https://doi.org/10.21105/astro.2307.05683
2023-07-23
Abstract:We describe the implementation of the symplectic N-body integrator WHFast512 using Single Instruction Multiple Data (SIMD) parallelism and 512-bit Advanced Vector Extensions (AVX512). We are able to speed up integrations of planetary systems by up to 4.7x compared to the non-vectorized version of WHFast. WHFast512 can integrate the Solar System with 8 planets for 5 billion years in less than 1.4 days. To our knowledge, this makes WHFast512 the fastest direct N-body integrator for systems of this kind. As an example, we present an ensemble of 40-Gyr integrations of the Solar System. Ignoring the Sun's post-main sequence evolution, we show that the instability rate is well captured by a diffusion model. WHFast512 is freely available within the REBOUND package.
Earth and Planetary Astrophysics,Instrumentation and Methods for Astrophysics,Computational Physics
What problem does this paper attempt to address?