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.

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

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

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

From Thread to Transcontinental Computer: Disturbing Lessons in Distributed Supercomputing

4.45 Pflops Astrophysical N-Body Simulation on K computer -- The Gravitational Trillion-Body Problem

Astrophysical particle simulations with large custom GPU clusters on three continents

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

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

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

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

Swphotons: Toward Trillion‐body‐scale Cosmological N‐body Simulations on Sunway TaihuLight Supercomputer

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

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

Astrophysical Particle Simulations with Custom GPU Clusters

Fifteen Minutes of Fame: The Dynamics of Information Access on the Web

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

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

N-Body Simulations on GPUs

The New Worlds Simulations: Large-scale Simulations across Three Cosmologies

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

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