A Hybrid MPI+Threads Approach to Particle Group Finding Using Union-Find

03/25/2020
by   James S. Willis, et al.
0

The Friends-of-Friends (FoF) algorithm is a standard technique used in cosmological N-body simulations to identify structures. Its goal is to find clusters of particles (called groups) that are separated by at most a cut-off radius. N-body simulations typically use most of the memory present on a node, leaving very little free for a FoF algorithm to run on-the-fly. We propose a new method that utilises the common Union-Find data structure and a hybrid MPI+threads approach. The algorithm can also be expressed elegantly in a task-based formalism if such a framework is used in the rest of the application. We have implemented our algorithm in the open-source cosmological code, SWIFT. Our implementation displays excellent strong- and weak-scaling behaviour on realistic problems and compares favourably (speed-up of 18x) over other methods commonly used in the N-body community.

READ FULL TEXT
research
07/22/2020

Collectives in hybrid MPI+MPI code: design, practice and performance

The use of hybrid scheme combining the message passing programming model...
research
10/11/2017

Subdomain Deflation and Algebraic Multigrid: Combining Multiscale with Multilevel

The paper proposes a combination of the subdomain deflation method and l...
research
10/11/2017

Subdomain Deflation Combined with Local AMG: a Case Study Using AMGCL Library

The paper proposes a combination of the subdomain deflation method and l...
research
03/09/2020

CUBE – Towards an Optimal Scaling of Cosmological N-body Simulations

N-body simulations are essential tools in physical cosmology to understa...
research
12/13/2020

A GPU-Accelerated Fast Summation Method Based on Barycentric Lagrange Interpolation and Dual Tree Traversal

We present the barycentric Lagrange dual tree traversal (BLDTT) fast sum...
research
08/19/2020

Building Halo Merger Trees from the Q Continuum Simulation

Cosmological N-body simulations rank among the most computationally inte...
research
09/23/2015

A shared memory implementation of pipelined Parareal

The paper introduces an OpenMP implementation of pipelined Parareal and ...

Please sign up or login with your details

Forgot password? Click here to reset