Computational Gravitational Dynamics with Modern Numerical Accelerators
Sterrewacht Leiden, Leiden University, P.O. Box 9513, 2300 RA Leiden, The Netherlands
arXiv:1409.5474 [astro-ph.IM], (18 Sep 2014)
@article{2014arXiv1409.5474P,
author={Portegies Zwart}, S. and {B{‘e}dorf}, J.},
title={"{Computational Gravitational Dynamics with Modern Numerical Accelerators}"},
journal={ArXiv e-prints},
archivePrefix={"arXiv"},
eprint={1409.5474},
primaryClass={"astro-ph.IM"},
keywords={Astrophysics – Instrumentation and Methods for Astrophysics, Astrophysics – Astrophysics of Galaxies, Computer Science – Other Computer Science},
year={2014},
month={sep},
adsurl={http://adsabs.harvard.edu/abs/2014arXiv1409.5474P},
adsnote={Provided by the SAO/NASA Astrophysics Data System}
}
We review the recent optimizations of gravitational N-body kernels for running them on graphics processing units (GPUs), on single hosts and massive parallel platforms. For each of the two main N-body techniques, direct summation and tree-codes, we discuss the optimization strategy, which is different for each algorithm. Because both the accuracy as well as the performance characteristics differ, hybridizing the two algorithms is essential when simulating a large N-body system with high-density structures containing few particles, and with low-density structures containing many particles. We demonstrate how this can be realized by splitting the underlying Hamiltonian, and we subsequently demonstrate the efficiency and accuracy of the hybrid code by simulating a group of 11 merging galaxies with massive black holes in the nuclei.
September 23, 2014 by hgpu
