English

Accelerating NBODY6 with Graphics Processing Units

Instrumentation and Methods for Astrophysics 2015-06-05 v1 Computational Physics

Abstract

We describe the use of Graphics Processing Units (GPUs) for speeding up the code NBODY6 which is widely used for direct NN-body simulations. Over the years, the N2N^2 nature of the direct force calculation has proved a barrier for extending the particle number. Following an early introduction of force polynomials and individual time-steps, the calculation cost was first reduced by the introduction of a neighbour scheme. After a decade of GRAPE computers which speeded up the force calculation further, we are now in the era of GPUs where relatively small hardware systems are highly cost-effective. A significant gain in efficiency is achieved by employing the GPU to obtain the so-called regular force which typically involves some 99 percent of the particles, while the remaining local forces are evaluated on the host. However, the latter operation is performed up to 20 times more frequently and may still account for a significant cost. This effort is reduced by parallel SSE/AVX procedures where each interaction term is calculated using mainly single precision. We also discuss further strategies connected with coordinate and velocity prediction required by the integration scheme. This leaves hard binaries and multiple close encounters which are treated by several regularization methods. The present nbody6-GPU code is well balanced for simulations in the particle range 1042×10510^4-2 \times 10^5 for a dual GPU system attached to a standard PC.

Keywords

Cite

@article{arxiv.1205.1222,
  title  = {Accelerating NBODY6 with Graphics Processing Units},
  author = {Keigo Nitadori and Sverre J. Aarseth},
  journal= {arXiv preprint arXiv:1205.1222},
  year   = {2015}
}

Comments

8 pages, 3 figures, 2 tables, MNRAS accepted

R2 v1 2026-06-21T20:59:14.889Z