Abstract:We present the results of gravitational direct $N$-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the $N$-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor-corrector Hermite integrator with block time steps. We find that for a large number of particles ($N \apgt 10^4$) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in $10^6$ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around $10^{-11}$, whereas for the GPU machine precision saturates around $10^{-6}$ For $N\apgt 10^6$ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.

Performance Analysis of Parallel Gravitationaln-Body Codes on Large GPU Clusters

High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA

Heterogeneous Programming and Optimization of Gyrokinetic Toroidal Code and Large-Scale Performance Test on TH-1A.

Nbody6++gpu: Ready for the Gravitational Million-Body Problem

High Performance Direct Gravitational N-body Simulations on Graphics Processing Units

High performance massively parallel direct N-body simulations on large GPU clusters

High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg

N-Body Simulations on GPUs

Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers

Supermassive Black Hole Binaries in High Performance Massively Parallel Direct N-body Simulations on Large GPU Clusters

Acceleration of Hybrid MPI Parallel NBODY6++ for Largen-Body Globular Cluster Simulations

Efficiency of parallel computations of gravitational forces by TreeCode method in N-body models

Astrophysical particle simulations with large custom GPU clusters on three continents

Parallel Time-Space Processing Model Based Fast N -Body Simulation on GPUs

SAPPORO: A way to turn your graphics cards into a GRAPE-6

Up to 700k GPU cores, Kepler, and the Exascale future for simulations of star clusters around black holes

Accelerating Astrophysical Particle Simulations with Programmable Hardware (FPGA and GPU)

Astrophysical Particle Simulations with Custom GPU Clusters

The Chamomile Scheme: An Optimized Algorithm for N-body simulations on Programmable Graphics Processing Units

Up to 700k GPU Cores, Kepler, and the Exascale Future for Simulations of Star Clusters Around Black Holes.

Parallelization of the Symplectic Massive Body Algorithm (SyMBA) $N$-body Code