New features of parallel implementation of N-body problems on GPU

DOI10.14529/MMP180111MaRDI QIDQ4692150zbMATH OpenOpenAlexFDO

Authors S. S. Khrapov, S. A. Khoperskov, A. V. Khoperskov

Publication date 26 October 2018

Published in Bulletin of the South Ural State University. Series "Mathematical Modelling, Programming and Computer Software" (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1803.01190, http://mathnet.ru/eng/vyuru423

zbMATH Keywords

\(N\)-body problemcollisionless systemgravitational instabilitymulti-GPUsingle and double precision numerical simulation

Mathematics Subject Classification ID

Simulation of dynamical systems (37M05) Modes of computation (nondeterministic, parallel, interactive, probabilistic, etc.) (68Q10) (n)-body problems (70F10)

Abstract: This paper focuses on the parallel implementation of a direct $N$-body method~(particle-particle algorithm) and the application of multiple GPUs for galactic dynamics simulations. Application of a hybrid OpenMP-CUDA technology is considered for models with a number of particles $Nsim10^{5}div10^7$. By means of $N$-body simulations of gravitationally unstable stellar galactic we have investigated the algorithms parallelization efficiency for various Nvidia Tesla graphics processors~(K20, K40, K80). Particular attention was paid to the parallel performance of simulations and accuracy of the numerical solution by comparing single and double floating-point precisions~(SP and DP). We showed that the double-precision simulations are slower by a factor of~$1.7$ than the single-precision runs performed on Nvidia Tesla K-Series processors. We also claim that application of the single-precision operations leads to incorrect result in the evolution of the non-axisymmetric gravitating $N$-body systems. In particular, it leads to significant quantitative and even qualitative distortions in the galactic disk evolution. For instance, after $10^4$ integration time steps for the single-precision numbers the total energy, momentum, and angular momentum of a system with $N=2^{20}$ conserve with accuracy of $10^{-3}$, $10^{-2}$ and $10^{-3}$ respectively, in comparison to the double-precision simulations these values are $10^{-5}$, $10^{-15}$ and $10^{-13}$, respectively. Our estimations evidence in favour of usage of the second-order accuracy schemes with double-precision numbers since it is more efficient than in the fourth-order schemes with single-precision numbers.

Recommendations

• GANBISS: a new GPU accelerated \(N\)-body code for binary star systems
• A sparse octree gravitational N-body code that runs entirely on the GPU processor
• scientific article; zbMATH DE number 1728464
• Direct \(N\)-body simulations
• \(N\)-body simulations of massively parallel architectures

Cited in

(5)

• Numerical modeling of the collisions of spheroidal galaxies: mass loss efficiency by baryon components
• A Stochastic-Based Optimized Schwarz Method for the Gravimetry Equations on GPU Clusters
• Analysis and parallel implementation of a forced N-body problem
• High-precision numerical simulations on a CUDA GPU: Kerr black hole tails
• GANBISS: a new GPU accelerated \(N\)-body code for binary star systems

Describes a project that uses

Uses Software

• CUDA

This page was built for publication: New features of parallel implementation of \(N\)-body problems on GPU

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4692150)