Performance analysis of direct N-body algorithms for astrophysical simulations on distributed systems

Authors:
Alessia Gualandris;Simon Portegies Zwart;Alfredo Tirado-Ramos
Affiliations:
Astronomical Institute and Section Computational Science, University of Amsterdam, The Netherlands;Astronomical Institute and Section Computational Science, University of Amsterdam, The Netherlands;Section Computational Science, University of Amsterdam, The Netherlands
Venue:
Parallel Computing
Year:
2007

Citing 4
Cited 4

Systolic and hyper-systolic algorithms for the gravitational N-body problem, with an application to Brownian motion

Journal of Computational Physics
Time Warp Cancellation Optimisations on High Latency Networks

DS-RT '03 Proceedings of the Seventh IEEE International Symposium on Distributed Simulation and Real-Time Applications
Unlocking the Performance of the BlueGene/L Supercomputer

Proceedings of the 2004 ACM/IEEE conference on Supercomputing
The Anatomy of the Grid: Enabling Scalable Virtual Organizations

International Journal of High Performance Computing Applications

Regular Paper: Interactive N-Body Simulations On the Grid: HLA Versus MPI

International Journal of High Performance Computing Applications
Simulating N-Body Systems on the Grid Using Dedicated Hardware

ICCS '08 Proceedings of the 8th international conference on Computational Science, Part I
Parallel Discrete Event N-Body Dynamics

PADS '11 Proceedings of the 2011 IEEE Workshop on Principles of Advanced and Distributed Simulation
A grid service for management of multiple HLA federate processes

PPAM'05 Proceedings of the 6th international conference on Parallel Processing and Applied Mathematics

Quantified Score

Hi-index	0.00

Visualization

Abstract

We discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an O(N^2) scaling with the number of particles. They can be applied to a variety of astrophysical problems, like the evolution of star clusters, the dynamics of black holes, the formation of planetary systems, and cosmological simulations. The simulation of realistic star clusters with sufficiently high accuracy cannot be performed on a single workstation but may be possible on parallel computers or grids. We have implemented two parallel schemes for a direct N-body code and we study their performance on general purpose parallel computers and large computational grids. We present the results of timing analyzes conducted on the different architectures and compare them with the predictions from theoretical models. We conclude that the simulation of star clusters with up to a million particles will be possible on large distributed computers in the next decade. Simulating entire galaxies however will in addition require new hybrid methods to speedup the calculation.