Modernizing Checkpointing and Output in the Earthquake Simulation Software SeisSol

Applicant

Prof. Dr. Michael Bader
Scientific Computing in Computer Science – group for Hardware-aware algorithms and software for HPC
Technical University of Munich

Project Overview

SeisSol is a software package for the simulation of earthquake dynamics and seismic wave propagation on domains with complicated geometry (e.g., topography, complex fault systems). It is developed by the groups of Michael Bader (TUM) and Alice-Agnes Gabriel (UC San Diego, LMU). Over the last decade, SeisSol has been re-engineered towards a high-performance, highly scalable solver that exploits a large percentage of the performance offered by current petascale systems. It has been used to realise a range of landmark simulations, especially focusing on earthquakes with highly complex source dynamics, tsunamigenic earthquakes and coupled earthquake-tsunami events. Recent projects have incorporated SeisSol into workflows for physics-based seismic hazard analysis (PSHA) and have developed novel coupled elastic-acoustic models (for simulation of 3D tsunamigenesis as well as of noise patterns of induced earthquakes).
In the proposed project, we will re-redesign, improve and optimise the output and checkpointing functionality of SeisSol. The current implementation of I/O and checkpointing was designed and (besides a few later extensions) implemented in 2018. In the meantime, SeisSol has been substantially extended in functionality, and new requirements from new applications have emerged. Together with new developments in the simulation and visualisation infrastructure, this has led to an extensive list of deficits in the current implementation. In the proposed project, we will tackle these deficits in a combined effort, to re-establish efficient checkpointing and I/O for SeisSol’s entire simulation capabilities, and to modernize and improve output routines to meet the requirements of new use cases and workflows at multi-peta- and exascale.