GANBISS: a new GPU accelerated N-body code for binary star systems

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.