A fully parallel, high precision, N-body code running on hybrid computing platforms

Abstract

We present a new implementation of the numerical integration of the classical, gravitational, N-body problem based on a high order Hermite's integration scheme with block time steps, with a direct evaluation of the particle-particle forces. The main innovation of this code (called HiGPUs) is its full parallelization, exploiting both OpenMP and MPI in the use of the multicore Central Processing Units as well as either Compute Unified Device Architecture (CUDA) or OpenCL for the hosted Graphic Processing Units. We tested both performance and accuracy of the code using up to 256 GPUs in the supercomputer IBM iDataPlex DX360M3 Linux Infiniband Cluster provided by the italian supercomputing consortium CINECA, for values of N up to 8 millions. We were able to follow the evolution of a system of 8 million bodies for few crossing times, task previously unreached by direct summation codes. The code is freely available to the scientific community.

Paper submitted to Journal of Computational Physics consisting in 28 pages, 9 figures.The previous submitted version was lacking of the bibliography, for a Tex problem

References in corpus (6)

• Accelerating NBODY6 with Graphics Processing Units
• N-body simulations of gravitational dynamics
• 6th and 8th Order Hermite Integrator for N-body Simulations
• Merging of globular clusters within inner galactic regions. II. The Nuclear Star Cluster formation
• N-Body Simulations on GPUs
• A Performance Comparison of Different Graphics Processing Units Running Direct N-Body Simulations

Cited by in corpus (44)

• Very massive stars, pair-instability supernovae and intermediate-mass black holes with the SEVN code
• The Astrophysical Multipurpose Software Environment
• Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code
• Gravitational wave sources from inspiralling globular clusters in the Galactic Centre and similar environments
• Collisionless loss-cone refilling: there is no final parsec problem
• The MEGaN project II. Gravitational waves from intermediate mass- and binary black holes around a supermassive black hole
• Hierarchical black hole triples in young star clusters: impact of Kozai-Lidov resonance on mergers
• A Semi-Analytic dynamical friction model that reproduces core stalling
• Brownian motion of massive black hole binaries and the final parsec problem
• On the formation of compact, massive sub-systems in stellar clusters and its relation with intermediate mass black holes
• Compact Binary Coalescences: Astrophysical Processes and Lessons Learned
• Gamma-ray and X-ray emission from the Galactic Centre: hints on the nuclear star cluster formation history
• The influence of Massive Black Hole Binaries on the Morphology of Merger Remnants
• Lack of nuclear clusters in dwarf spheroidal galaxies: implications for massive black holes formation and the cusp/core problem
• Supermassive black holes coalescence mediated by massive perturbers: implications for gravitational waves emission and nuclear cluster formation
• The dearth of nuclear star clusters in bright galaxies
• Supernova Kicks and Dynamics of Compact Remnants in the Galactic Centre
• The MEGaN project I. Missing formation of massive nuclear clusters and tidal disruption events by star clusters - massive black hole interactions
• Star cluster disruption by a massive black hole binary
• A Performance Comparison of Different Graphics Processing Units Running Direct N-Body Simulations
• From hydrodynamics to N-body simulations of star clusters: mergers and rotation
• The secular evolution of the Kuiper belt after a close stellar encounter
• Globular clusters as tracers of the host galaxy mass distribution: the Fornax dSph test case
• Intermediate mass black holes in globular clusters: effects on jerks and jounces of millisecond pulsars
• A connected component-based method for efficiently integrating multiscale $N$-body systems
• A New Hybrid Technique for Modeling Dense Star Clusters
• Future merger of the Milky Way with the Andromeda galaxy and the fate of their supermassive black holes
• Are we observing a NSC in course of formation in the NGC 4654 galaxy?
• Stellar hardening of massive black hole binaries: the impact of the host rotation
• Rapid Mass segregation in small stellar clusters
• Astrophysics with the Laser Interferometer Space Antenna
• GraviDy, a GPU modular, parallel direct-summation $N-$body integrator: Dynamics with softening
• Globular Clusters in the Galactic Center Region: expected behavior in the infalling and merger scenario
• Dynamical Processes in Globular Clusters
• The collision between the Milky Way and Andromeda and the fate of their Supermassive Black Holes
• Direct $N$-body code on low-power embedded ARM GPUs
• Population Synthesis of Gravitational Wave Sources
• The interaction between supermassive black holes and globular clusters
• Direct N-body application on low-power and energy-efficient parallel architectures
• KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics
• Sapporo2: A versatile direct $N$-body library
• Selfgravitating disks in binary systems: an SPH approach -- I. Implementation of the code and reliability tests
• Relationship between low-discrepancy sequence and static solution to multi-bodies problem
• A heterogeneous many-core platform for experiments on scalable custom interconnects and management of fault and critical events, applied to many-process applications: Vol. II, 2012 technical report