Sapporo2: A versatile direct $N$-body library

Abstract

Astrophysical direct $N$-body methods have been one of the first production algorithms to be implemented using NVIDIA's CUDA architecture. Now, almost seven years later, the GPU is the most used accelerator device in astronomy for simulating stellar systems. In this paper we present the implementation of the Sapporo2 $N$-body library, which allows researchers to use the GPU for $N$-body simulations with little to no effort. The first version, released five years ago, is actively used, but lacks advanced features and versatility in numerical precision and support for higher order integrators. In this updated version we have rebuilt the code from scratch and added support for OpenCL, multi-precision and higher order integrators. We show how to tune these codes for different GPU architectures and present how to continue utilizing the GPU optimal even when only a small number of particles ($N<100$) is integrated. This careful tuning allows Sapporo2 to be faster than Sapporo1 even with the added options and double precision data loads. The code runs on a range of NVIDIA and AMD GPUs in single and double precision accuracy. With the addition of OpenCL support the library is also able to run on CPUs and other accelerators that support OpenCL.

15 pages, 7 figures. Accepted for publication in Computational Astrophysics and Cosmology

References in corpus (9)

• Accelerating NBODY6 with Graphics Processing Units
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
• Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
• Harvesting graphics power for MD simulations
• N-body simulations of gravitational dynamics
• 6th and 8th Order Hermite Integrator for N-body Simulations
• A fully parallel, high precision, N-body code running on hybrid computing platforms
• The formation of young dense star clusters through mergers
• A Performance Comparison of Different Graphics Processing Units Running Direct N-Body Simulations