Towards optimal parallel PM N-body codes: PMFAST
Hugh Merz, Ue-Li Pen, Hy Trac
Abstract
We present a new parallel PM N-body code named PMFAST that is freely available to the public. PMFAST is based on a two-level mesh gravity solver where the gravitational forces are separated into long and short range components. The decomposition scheme minimizes communication costs and allows tolerance for slow networks. The code approaches optimality in several dimensions. The force computations are local and exploit highly optimized vendor FFT libraries. It features minimal memory overhead, with the particle positions and velocities being the main cost. The code features support for distributed and shared memory parallelization through the use of MPI and OpenMP, respectively. The current release version uses two grid levels on a slab decomposition, with periodic boundary conditions for cosmological applications. Open boundary conditions could be added with little computational overhead. We present timing information and results from a recent cosmological production run of the code using a 37123 mesh with 6.4 x 109 particles. PMFAST is cost-effective, memory-efficient, and is publicly available.
Create a lesson
Related papers
On binary pulsars and the force of gravity
Davor Palle
Tidal torques. A critical review of some techniques
Michael Efroimsky, James G. Williams
Dynamics of a Spherical Accretion Shock with Neutrino Heating and Alpha-Particle Recombination
Rodrigo Fernández, Christopher Thompson
Asymptotically FRW black holes
J. T. Firouzjaee, Reza Mansouri
Reaction of Accretion Disks to Abrupt Mass Loss During Binary Black Hole Merger
Sean M. O'Neill, M. Coleman Miller, Tamara Bogdanovic et al.
A Gamma-Ray Burst/Pulsar for Cosmic-Ray Positrons with a Dark Matter-like Spectrum
Kunihito Ioka