Phantom-GRAPE: numerical software library to accelerate collisionless $N$-body simulation with SIMD instruction set on x86 architecture

Abstract

(Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed. Using an Intel Core i7--2600 processor, we measure the performance of our library for both the forces. In the case of Newton's forces, we achieve 2 x 10^9 interactions per second with 1 processor core, which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With 4 processor cores, we obtain the performance of 8 x 10^9 interactions per second. In the case of the arbitrarily shaped forces, we can calculate 1 x 10^9 and 4 x 10^9 interactions per second with 1 and 4 processor cores, respectively. The performance with 1 processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend weakly on the number of particles. It is good contrast with the fact that the performance of force calculations accelerated by GPUs depends strongly on the number of particles. Substantially weak dependence of the performance on the number of particles is suitable to collisionless N-body simulations, since these simulations are usually performed with sophisticated N-body solvers such as Tree- and TreePM-methods combined with an individual timestep scheme. Collisionless N-body simulations accelerated with our library have significant advantage over those accelerated by GPUs, especially on massively parallel environments.

19 pages, 11 figures, 4tables, accepted for publication in New Astronomy

References in corpus (3)

• 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
• A parallel gravitational N-body kernel

Cited by in corpus (38)

• Simulating cosmic structure formation with the GADGET-4 code
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• Dynamics of Non-Steady Spiral Arms in Disk Galaxies
• Implementation and performance of FDPS: A Framework Developing Parallel Particle Simulation Codes
• Hierarchical Formation of Dark Matter Halos and the Free Streaming Scale
• Formation of Intermediate-Mass Black Holes through Runaway Collisions in the First Star Clusters
• Where are the Low-Mass Population III Stars ?
• SILVERRUSH. VI. A simulation of Ly$\alpha$ emitters in the reionization epoch and a comparison with Subaru Hyper Suprime-Cam survey early data
• Short-term dynamical evolution of grand-design spirals in barred galaxies
• Double-Detonation Models for Type Ia Supernovae: Trigger of Detonation in Companion White Dwarfs and Signatures of Companions' Stripped-off Materials
• Age Dating the Galactic Bar with the Nuclear Stellar Disc
• The Abundance and Structure of Subhaloes near the Free Streaming Scale and Their Impact on Indirect Dark Matter Searches
• Growth of intermediate mass black holes by tidal disruption events in the first star clusters
• Three-Dimensional Simulation of Double-Detonations in the Double-Degenerate Model for Type Ia Supernovae and Interaction of Ejecta with a Surviving White Dwarf Companion
• Eventful Evolution of Giant Molecular Clouds in Dynamically Evolving Spiral Arms
• On the evolution of gas clouds exposed to AGN radiation. I. Three-dimensional radiation hydrodynamic simulations
• Gas velocity patterns in simulated galaxies: Observational diagnostics of spiral structure theories
• Imprint of anisotropic primordial non-Gaussianity on halo intrinsic alignments in simulations
• Radial distributions of arm-gas offsets as an observational test of spiral theories
• A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics
• Tidal Disruption of a White Dwarf by a Black Hole: The Diversity of Nucleosynthesis, Explosion Energy, and the Fate of Debris Streams
• Does Explosive Nuclear Burning occur in Tidal Disruption Events of White Dwarfs by Intermediate Mass Black Holes ?
• Properties of thick discs formed in clumpy galaxies
• GOTHIC: Gravitational oct-tree code accelerated by hierarchical time step controlling
• SIRIUS project II: a new tree-direct hybrid code for smoothed particle hydrodynamics/N-body simulations of star clusters
• Simulating the Smallest Ring World of Chariklo
• Constraining Self-Interacting Dark Matter with Dwarf Spheroidal Galaxies and High-resolution Cosmological $N$-body Simulations
• Santa Barbara Cluster Comparison Test with DISPH
• High-resolution hydrodynamic simulation of tidal detonation of a helium white dwarf by an intermediate mass black hole
• Global N-Body Simulation of Galactic Spiral Arms
• Imprints of Zero-Age Velocity Dispersions and Dynamical Heating on the Age-Velocity dispersion Relation
• A pilgrimage to gravity on GPUs
• Emergence of a stellar cusp by a dark matter cusp in a low-mass compact ultra-faint dwarf galaxy
• First Star Survivors as Metal-Rich Halo Stars that Experienced Supernova Explosions in Binary Systems
• Acceleration of the tree method with SIMD instruction set
• Age distribution of stars in boxy/peanut/X-shaped bulges formed without bar buckling
• A 400 Trillion-Grid Vlasov Simulation on Fugaku Supercomputer: Large-Scale Distribution of Cosmic Relic Neutrinos in a Six-dimensional Phase Space
• Cosmological Calculations on the GPU