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.

GRAPE-5: A Special-Purpose Computer for N-body Simulation

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

GRAPESPH: Cosmological SPH simulations with the special purpose hardware GRAPE

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

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

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

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

Performance and accuracy of a GRAPE-3 system for collisionless N-body simulations

The PCI Interface for GRAPE Systems: PCI-HIB

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

The Formation of Galaxies: A Challenge for Supercomputers -- a Simple Task for GRAPE ?

Collisional Stellar Dynamics, Gas Dynamics and Special Purpose Computing

AHA-GRAPE: Adaptive Hydrodynamic Architecture - GRAvity PipE

Heterogeneous Programming and Optimization of Gyrokinetic Toroidal Code and Large-Scale Performance Test on TH-1A.

Hardware acceleration of tensor-structured multilevel ewald summation method on MDGRAPE-4A, a special-purpose computer system for molecular dynamics simulations

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

Astrophysical N-Body Simulations: Algorithms and Challenges

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

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

VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods

Acceleration Method of Neighbor Search with GRAPE and Morton-ordering