Performance analysis of parallel gravitational $N$-body codes on large GPU cluster

Abstract

We compare the performance of two very different parallel gravitational $N$-body codes for astrophysical simulations on large GPU clusters, both pioneer in their own fields as well as in certain mutual scales - NBODY6++ and Bonsai. We carry out the benchmark of the two codes by analyzing their performance, accuracy and efficiency through the modeling of structure decomposition and timing measurements. We find that both codes are heavily optimized to leverage the computational potential of GPUs as their performance has approached half of the maximum single precision performance of the underlying GPU cards. With such performance we predict that a speed-up of $200-300$ can be achieved when up to 1k processors and GPUs are employed simultaneously. We discuss the quantitative information about comparisons of two codes, finding that in the same cases Bonsai adopts larger time steps as well as relative energy errors than NBODY6++, typically ranging from $10-50$ times larger, depending on the chosen parameters of the codes. While the two codes are built for different astrophysical applications, in specified conditions they may overlap in performance at certain physical scale, and thus allowing the user to choose from either one with finetuned parameters accordingly.

15 pages, 7 figures, 3 tables, accepted for publication in Research in Astronomy and Astrophysics (RAA)

References in corpus (6)

• Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe
• Introducing the Illustris Project: the evolution of galaxy populations across cosmic time
• NBODY6++GPU: Ready for the gravitational million-body problem
• Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
• 6th and 8th Order Hermite Integrator for N-body Simulations
• 24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs

Cited by in corpus (5)

• The DRAGON simulations: globular cluster evolution with a million stars
• Extreme-mass-ratio burst detection with TianQin
• Revisit the rate of tidal disruption events: the role of the partial tidal disruption event
• New features of parallel implementation of N-body problems on GPU
• Is the overconcentration of pristine populations in Galactic globular clusters real? An N-body approach to the problem