N-body simulation for self-gravitating collisional systems with a new SIMD instruction set extension to the x86 architecture, Advanced Vector eXtensions

Abstract

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of 20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved 90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N 105 on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.

14 pages, 9 figures, 3 tables, accepted for publication in New Astronomy. The code is publicly available at http://code.google.com/p/phantom-grape/

References in corpus (1)

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

Cited by in corpus (29)

• Simulating cosmic structure formation with the GADGET-4 code
• The Uchuu Simulations: Data Release 1 and Dark Matter Halo Concentrations
• Accelerating NBODY6 with Graphics Processing Units
• Implementation and performance of FDPS: A Framework Developing Parallel Particle Simulation Codes
• Hierarchical Formation of Dark Matter Halos and the Free Streaming Scale
• Enrichment of r-process elements in dwarf spheroidal galaxies in chemo-dynamical evolution model
• Dynamical evolution of stellar-mass black holes in dense stellar clusters: estimate for merger rate of binary black holes originating from globular 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
• Double-Detonation Models for Type Ia Supernovae: Trigger of Detonation in Companion White Dwarfs and Signatures of Companions' Stripped-off Materials
• Phantom-GRAPE: numerical software library to accelerate collisionless $N$-body simulation with SIMD instruction set on x86 architecture
• The Abundance and Structure of Subhaloes near the Free Streaming Scale and Their Impact on Indirect Dark Matter Searches
• 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
• On the evolution of gas clouds exposed to AGN radiation. I. Three-dimensional radiation hydrodynamic simulations
• Accelerated orbital decay of supermassive black hole binaries in merging nuclear star clusters
• 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 ?
• A review of High Performance Computing foundations for scientists
• Adaptive fast multipole methods on the GPU
• 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
• 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
• A pilgrimage to gravity on GPUs
• First Star Survivors as Metal-Rich Halo Stars that Experienced Supernova Explosions in Binary Systems
• Acceleration of the tree method with SIMD instruction set
• Sapporo2: A versatile direct $N$-body library