A Particle-Mesh Integrator for Galactic Dynamics Powered by GPGPUs

Authors:
Dominique Aubert;Mehdi Amini;Romaric David
Affiliations:
Observatoire Astronomique, Universite de Strasbourg, France;Direction Informatique, Universite de Strasbourg, France;Direction Informatique, Universite de Strasbourg, France
Venue:
ICCS '09 Proceedings of the 9th International Conference on Computational Science: Part I
Year:
2009

Citing 4
Cited 0

Computer simulation using particles

Computer simulation using particles
Numerical Recipes in C++: the art of scientific computing

Numerical Recipes in C++: the art of scientific computing
Performance and accuracy of hardware-oriented native-, emulated-and mixed-precision solvers in FEM simulations

International Journal of Parallel, Emergent and Distributed Systems
Fast parallel Particle-To-Grid interpolation for plasma PIC simulations on the GPU

Journal of Parallel and Distributed Computing

Quantified Score

Hi-index	0.00

Visualization

Abstract

We present a particle-mesh N-body integrator running on GPU using CUDA. Relying on a grid-based description of the gravitational potential, it can simulate the evolution of self-interacting `stars' in order to model e.g. galaxies. All the steps of the application have been ported on the GPU , namely 1/ an histogramming algorithm with CUDPP, 2/ of the resolution of the Poisson equation by means of FFT with CUFFT and multi-grid relaxation, 3/ of an optimized finite difference scheme to compute the accelerations of stars and 4/ of an update procedure for positions and velocities. We present several tests at different resolution, and reach a speedup from 2 to 50 depending on the resolution and on the test case.