WHFast512: A symplectic N-body integrator for planetary systems optimized with AVX512 instructions
arXiv:2307.05683 · doi:10.21105/astro.2307.05683
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.
arXiv comment8 pages, 4 figures, 4.7x faster, accepted by the Open Journal of Astrophysics, WHFast512 source code available at https://github.com/hannorein/rebound/
References in corpus (6)
- A study of symplectic integrators for planetary system problems: error analysis and comparisons
- Long-term dynamics of the solar system inner planets
- Stepsize errors in the -body problem: discerning Mercury's true possible long-term orbits
- Simple physics and integrators accurately reproduce Mercury instability statistics
- orbitN: A symplectic integrator for planetary systems dominated by a central mass -- Insight into long-term solar system chaos
- Mercury's chaotic secular evolution as a subdiffusive process