Abstract
We present abacus, a fast and accurate cosmological N-body code based on a new method for calculating the gravitational potential from a static multipole mesh. The method analytically separates the near-and far-field forces, reducing the former to direct 1/r2 summation and the latter to a discrete convolution over multipoles. The method achieves 70 million particle updates per second per node of the Summit supercomputer, while maintaining a median fractional force error of 10-5. We express the simulation time-step as an event-driven 'pipeline', incorporating asynchronous events such as completion of co-processor work, input/output, and network communication. abacus has been used to produce the largest suite of N-body simulations to date, the abacussummit suite of 60 trillion particles, incorporating on-The-fly halo finding. abacus enables the production of mock catalogues of the volume and resolution required by the coming generation of cosmological surveys.
Original language | English (US) |
---|---|
Pages (from-to) | 575-596 |
Number of pages | 22 |
Journal | Monthly Notices of the Royal Astronomical Society |
Volume | 508 |
Issue number | 1 |
DOIs | |
State | Published - Nov 1 2021 |
Keywords
- cosmology: Theory
- methods: numerical
ASJC Scopus subject areas
- Astronomy and Astrophysics
- Space and Planetary Science