A GPU accelerated Barnes-Hut Tree Code for FLASH4

Abstract

We present a GPU accelerated CUDA-C implementation of the Barnes Hut (BH) tree code for calculating the gravitational potential on octree adaptive meshes. The tree code algorithm is implemented within the FLASH4 adaptive mesh refinement (AMR) code framework and therefore fully MPI parallel. We describe the algorithm and present test results that demonstrate its accuracy and performance in comparison to the algorithms available in the current FLASH4 version. We use a MacLaurin spheroid to test the accuracy of our new implementation and use spherical, collapsing cloud cores with effective AMR to carry out performance tests also in comparison with previous gravity solvers. Depending on the setup and the GPU/CPU ratio, we find a speedup for the gravity unit of at least a factor of 3 and up to 60 in comparison to the gravity solvers implemented in the FLASH4 code. We find an overall speedup factor for full simulations of at least factor 1.6 up to a factor of 10

For further information see: http://www.hs.uni-hamburg.de/gpubh

References in corpus (4)

• Modeling Collapse and Accretion in Turbulent Gas Clouds: Implementation and Comparison of Sink Particles in AMR and SPH
• A Direct Multigrid Poisson Solver for Oct-Tree Adaptive Meshes
• An Improved Multipole Approximation for Self-Gravity and Its Importance for Core-Collapse Supernova Simulations
• 24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs

Cited by in corpus (13)

• FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation
• GAMER-2: a GPU-accelerated adaptive mesh refinement code -- accuracy, performance, and scalability
• Tree-based solvers for adaptive mesh refinement code FLASH - I: gravity and optical depths
• On the shape and completeness of the column density probability distribution function of molecular clouds
• Global dynamics of the interstellar medium in magnetised disc galaxies
• The driving of turbulence in simulations of molecular cloud formation and evolution
• Deuterium fractionation and H2D+ evolution in turbulent and magnetized cloud cores
• The Origin of Filamentary Star Forming Clouds in Magnetised Galaxies
• A GPU-Accelerated Fast Summation Method Based on Barycentric Lagrange Interpolation and Dual Tree Traversal
• Fast deuterium fractionation in magnetized and turbulent filaments
• Magnetohydrodynamics with GAMER
• The Spatial Power Spectrum and Derived Turbulent Properties of Isolated Galaxies
• The turbulence driving parameter of molecular clouds in disc galaxies