Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

TL;DR

This paper discusses the evolution of GROMACS, a molecular dynamics simulation package, highlighting its parallelization strategies, acceleration techniques, and challenges faced in achieving exascale performance on supercomputers.

Contribution

It presents the advancements in GROMACS' parallelization, acceleration methods, and algorithmic redesigns to enable exascale molecular dynamics simulations.

Findings

GROMACS 4.6 uses SIMD, GPU offloading, OpenMP, and MPI for high performance.

Revisiting fundamental algorithms was necessary for acceleration.

Future challenges include implementing fine-grained task parallelism.

Abstract

GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.