How well do STARLAB and NBODY compare? II: Hardware and accuracy
arXiv:1201.5692 · doi:10.1111/j.1365-2966.2012.20581.x
Abstract
Most recent progress in understanding the dynamical evolution of star clusters relies on direct N-body simulations. Owing to the computational demands, and the desire to model more complex and more massive star clusters, hardware calculational accelerators, such as GRAPE special-purpose hardware or, more recently, GPUs (i.e. graphics cards), are generally utilised. In addition, simulations can be accelerated by adjusting parameters determining the calculation accuracy (i.e. changing the internal simulation time step used for each star). We extend our previous thorough comparison (Anders et al. 2009) of basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY6. Here we focus on differences arising from using different hardware accelerations (including the increasingly popular graphic card accelerations/GPUs) and different calculation accuracy settings. We use the large number of star cluster models (for a fixed stellar mass function, without stellar/binary evolution, primordial binaries, external tidal fields etc) already used in the previous paper, evolve them with STARLAB/KIRA (and NBODY6, where required), analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we apply the bootstrap algorithm for functional dependencies developed in our previous study. In general we find very high comparability of the simulation results, independent of the used computer hardware (including the hardware accelerators) and the used N-body code. For the tested accuracy settings we find that for reduced accuracy (i.e. time step at least a factor 2.5 larger than the standard setting) most simulation results deviate significantly from the results using standard settings. The remaining deviations are comprehensible and explicable.
14 pages incl. 3 pages with figures and 4 pages of tables (analysis results), MNRAS in press
References in corpus (13)
- GALEV evolutionary synthesis models - I. Code, input physics and web-interface
- A multiphysics and multiscale software environment for modeling astrophysical systems
- The Origin of OB Runaway Stars
- Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg
- Teraflop per second gravitational lensing ray-shooting using graphics processing units
- Hyperfast pulsars as the remnants of massive stars ejected from young star clusters
- A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
- Parallel Algorithm for Solving Kepler's Equation on Graphics Processing Units: Application to Analysis of Doppler Exoplanet Searches
- A parallel gravitational N-body kernel
- MYRIAD: A new N-body code for simulations of Star Clusters
- Clearing the Dust from Globular Clusters
- High Performance Gravitational N-body Simulations on a Planet-wide Distributed Supercomputer