Abstract
We present a GPU accelerated N-body integrator using the Bulirsch–Stoer method, called
GANBISS (GPU accelerated n-body code for binary star systems). It is designed to simulate
the dynamical evolution of planetesimal disks in binary star systems which contain some
thousand disk objects. However, it can also be used for studies of non-interacting massless
bodies where up to 50 million objects can be studied in a simulation. GANBISS shows the
energy and angular momentum conservation behavior of non-symplectic integration methods.
The code is written in CUDA C and can be run on NVIDIA GPUs of compute capability of
at least 3.5. A comparison of GPU and CPU computations indicates a speed-up of the GPU
performance of up to 100 times—depending on the number of disk objects.
GANBISS (GPU accelerated n-body code for binary star systems). It is designed to simulate
the dynamical evolution of planetesimal disks in binary star systems which contain some
thousand disk objects. However, it can also be used for studies of non-interacting massless
bodies where up to 50 million objects can be studied in a simulation. GANBISS shows the
energy and angular momentum conservation behavior of non-symplectic integration methods.
The code is written in CUDA C and can be run on NVIDIA GPUs of compute capability of
at least 3.5. A comparison of GPU and CPU computations indicates a speed-up of the GPU
performance of up to 100 times—depending on the number of disk objects.
| Original language | English |
|---|---|
| Article number | 33 |
| Number of pages | 14 |
| Journal | Celestial Mechanics and Dynamical Astronomy |
| Volume | 135 |
| Issue number | 3 |
| Publication status | Published - May 2023 |
Austrian Fields of Science 2012
- 103004 Astrophysics
- 102009 Computer simulation