GANBISS: A new GPU accelerated N-body code for Binary Star Systems

Abstract

We present a GPU acclerated N-body integrator using the Bulirsch-Stoer method, called GANBISS(GPU Accelerated N-body code for Binary Star Systems). It is designed to (i) simulate the dynamics and evolution of planetesimal disks in binary star systems which contains of some thousand disk objects and (ii) supports the integration of up to 40 million non-interacting and massless bodies. GANBISS shows the energy and angular momentum conservation behaviour of non-symplectic integration methods. We compare the performance to a CPU implementation and find a speed up of up to 100 times faster for the GPU implementation, depending on the number of integrated bodies. The code is written in CUDA C and can be run on NVIDIA GPUs of compute capability of at least 3.5.