English

Scalable communication for high-order stencil computations using CUDA-aware MPI

Distributed, Parallel, and Cluster Computing 2022-05-11 v2 Computational Physics Fluid Dynamics

Abstract

Modern compute nodes in high-performance computing provide a tremendous level of parallelism and processing power. However, as arithmetic performance has been observed to increase at a faster rate relative to memory and network bandwidths, optimizing data movement has become critical for achieving strong scaling in many communication-heavy applications. This performance gap has been further accentuated with the introduction of graphics processing units, which can provide by multiple factors higher throughput in data-parallel tasks than central processing units. In this work, we explore the computational aspects of iterative stencil loops and implement a generic communication scheme using CUDA-aware MPI, which we use to accelerate magnetohydrodynamics simulations based on high-order finite differences and third-order Runge-Kutta integration. We put particular focus on improving intra-node locality of workloads. Our GPU implementation scales strongly from one to 6464 devices at 50%50\%--87%87\% of the expected efficiency based on a theoretical performance model. Compared with a multi-core CPU solver, our implementation exhibits 2020--60×60\times speedup and 99--12×12\times improved energy efficiency in compute-bound benchmarks on 1616 nodes.

Keywords

Cite

@article{arxiv.2103.01597,
  title  = {Scalable communication for high-order stencil computations using CUDA-aware MPI},
  author = {Johannes Pekkilä and Miikka S. Väisälä and Maarit J. Käpylä and Matthias Rheinhardt and Oskar Lappi},
  journal= {arXiv preprint arXiv:2103.01597},
  year   = {2022}
}

Comments

15 pages, 15 figures. Updated with the accepted manuscript. More extensive tests added and wording clarified in several places. Please refer to the published article for the most polished version

R2 v1 2026-06-23T23:39:13.283Z