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.

Graphic-Card Cluster for Astrophysics (GraCCA) -- Performance Tests

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

Astrophysical Particle Simulations with Custom GPU Clusters

AHA-GRAPE: Adaptive Hydrodynamic Architecture - GRAvity PipE

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

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

Astrophysical particle simulations with large custom GPU clusters on three continents

GAMER: a GPU-Accelerated Adaptive Mesh Refinement Code for Astrophysics

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

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

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

Scalable Clustering Using Graphics Processors

GRAPESPH: Cosmological SPH simulations with the special purpose hardware GRAPE

GPU Accelerated Non-Illuminated Graphical Electromagnetic Computing Method with High Accuracy

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

N-Body Simulations on GPUs

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

CPU/GPU Computing for Long-Wave Radiation Physics on Large GPU Clusters.

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

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