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 Direct N-Body Algorithms on Special-Purpose Supercomputers

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

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

Large Scale Numerical Simulation Via Parallelization and Reconfigurable Computing Hardware

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

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

N-Body Simulations on GPUs

GRAPE-6A: A single-card GRAPE-6 for parallel PC-GRAPE cluster system

GRAPE-5: A Special-Purpose Computer for N-body Simulation

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

Astrophysical N-Body Simulations: Algorithms and Challenges

From Newton to Einstein –N-Body Dynamics in Galactic Nuclei and SPH Using New Special Hardware and Astrogrid-D

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

How to Build and Use Special Purpose PC Clusters in Stellar Dynamics

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

Monitoring the Postural Stability of Planar Bipedal Robots using the Moment-Height Stability Measure

Direct N-body application on low-power and energy-efficient parallel architectures

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

The Formation of Galaxies: A Challenge for Supercomputers -- a Simple Task for GRAPE ?

High Performance Gravitational N-body Simulations on a Planet-wide Distributed Supercomputer

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