papers

Publications (22)

cs.GL2020

An Environment for Sustainable Research Software in Germany and Beyond: Current State, Open Challenges, and Call for Action

Hartwig Anzt, Felix Bach, Stephan Druskat +40

Research software has become a central asset in academic research. It optimizes existing and enables new research methods, implements and embeds research knowledge, and constitutes…

cs.MS2020

Vectorization and Minimization of Memory Footprint for Linear High-Order Discontinuous Galerkin Schemes

Jean-Matthieu Gallard, Leonhard Rannabauer, Anne Reinarz +1

We present a sequence of optimizations to the performance-critical compute kernels of the high-order discontinuous Galerkin solver of the hyperbolic PDE engine ExaHyPE -- successiv…

math.NA2018

Efficient implementation of ADER discontinuous Galerkin schemes for a scalable hyperbolic PDE engine

Michael Dumbser, Francesco Fambri, Maurizio Tavelli +2

In this paper we discuss a new and very efficient implementation of high order accurate ADER discontinuous Galerkin (ADER-DG) finite element schemes on modern massively parallel su…

math.NA2025

Mixed-Precision in High-Order Methods: the Impact of Floating-Point Precision on the ADER-DG Algorithm

Marc Marot-Lassauzaie, Michael Bader

We present a mixed-precision implementation of the high-order discontinuous Galerkin method with ADER time stepping (ADER-DG) for solving hyperbolic systems of partial differential…

physics.comp-ph2021

3D Acoustic-Elastic Coupling with Gravity: The Dynamics of the 2018 Palu, Sulawesi Earthquake and Tsunami

Lukas Krenz, Carsten Uphoff, Thomas Ulrich +4

We present a highly scalable 3D fully-coupled Earth & ocean model of earthquake rupture and tsunami generation. We model seismic, acoustic and surface gravity wave propagation in e…

cs.MS2021

High Performance Uncertainty Quantification with Parallelized Multilevel Markov Chain Monte Carlo

Linus Seelinger, Anne Reinarz, Leonhard Rannabauer +3

Numerical models of complex real-world phenomena often necessitate High Performance Computing (HPC). Uncertainties increase problem dimensionality further and pose even greater cha…

cs.MS2019

Yet Another Tensor Toolbox for discontinuous Galerkin methods and other applications

Carsten Uphoff, Michael Bader

The numerical solution of partial differential equations is at the heart of many grand challenges in supercomputing. Solvers based on high-order discontinuous Galerkin (DG) discret…

cs.DC2020

TeaMPI -- Replication-based Resilience without the (Performance) Pain

Philipp Samfass, Tobias Weinzierl, Benjamin Hazelwood +1

In an era where we can not afford to checkpoint frequently, replication is a generic way forward to construct numerical simulations that can continue to run even if hardware parts…

physics.comp-ph2019

Influence of A-Posteriori Subcell Limiting on Fault Frequency in Higher-Order DG Schemes

Anne Reinarz, Jean-Matthieu Gallard, Michael Bader

Soft error rates are increasing as modern architectures require increasingly small features at low voltages. Due to the large number of components used in HPC architectures, these…

cs.DC2020

Lightweight Task Offloading Exploiting MPI Wait Times for Parallel Adaptive Mesh Refinement

Philipp Samfass, Tobias Weinzierl, Dominic E. Charrier +1

Balancing the workload of sophisticated simulations is inherently difficult, since we have to balance both computational workload and memory footprint over meshes that can change a…

math.NA2021

A stable discontinuous Galerkin method for linear elastodynamics in 3D geometrically complex media using physics based numerical fluxes

Kenneth Duru, Leonhard Rannabauer, Alice-Agnes Gabriel +3

High order accurate and explicit time-stable solvers are well suited for hyperbolic wave propagation problems. As a result of the complexities of real geometries, internal interfac…

math.NA2020

A stable discontinuous Galerkin method for the perfectly matched layer for elastodynamics in first order form

Kenneth Duru, Leonhard Rannabauer, Alice-Agnes Gabriel +2

We present a stable discontinuous Galerkin (DG) method with a perfectly matched layer (PML) for three and two space dimensional linear elastodynamics, in velocity-stress formulatio…

cs.DC2022

An Efficient ADER-DG Local Time Stepping Scheme for 3D HPC Simulation of Seismic Waves in Poroelastic Media

Sebastian Wolf, Martin Galis, Carsten Uphoff +4

Many applications from geosciences require simulations of seismic waves in porous media. Biot's theory of poroelasticity describes the coupling between solid and fluid phases and i…

physics.geo-ph2023

Numerical simulations of seismo-acoustic nuisance patterns from an induced M1.8 earthquake in the Helsinki, southern Finland, metropolitan area

Lukas Krenz, Sebastian Wolf, Gregor Hillers +2

Seismic waves can couple with the atmosphere and generate sound waves. The influence of faulting mechanisms on earthquake sound patterns provides opportunities for earthquake sourc…

math.NA2018

A simple diffuse interface approach on adaptive Cartesian grids for the linear elastic wave equations with complex topography

Maurizio Tavelli, Michael Dumbser, Dominic Etienne Charrier +3

In most classical approaches of computational geophysics for seismic wave propagation problems, complex surface topography is either accounted for by boundary-fitted unstructured m…

physics.comp-ph2019

A High-Order Discontinuous Galerkin Solver with Dynamic Adaptive Mesh Refinement to Simulate Cloud Formation Processes

Lukas Krenz, Leonhard Rannabauer, Michael Bader

We present a high-order discontinuous Galerkin (DG) solver of the compressible Navier-Stokes equations for cloud formation processes. The scheme exploits an underlying parallelized…

cs.MS2020

Role-Oriented Code Generation in an Engine for Solving Hyperbolic PDE Systems

Jean-Matthieu Gallard, Lukas Krenz, Leonhard Rannabauer +2

The development of a high performance PDE solver requires the combined expertise of interdisciplinary teams with respect to application domain, numerical scheme and low-level optim…

cs.SE2021

Doubt and Redundancy Kill Soft Errors -- Towards Detection and Correction of Silent Data Corruption in Task-based Numerical Software

Philipp Samfass, Tobias Weinzierl, Anne Reinarz +1

Resilient algorithms in high-performance computing are subject to rigorous non-functional constraints. Resiliency must not increase the runtime, memory footprint or I/O demands too…

cs.MS2020

ExaHyPE: An Engine for Parallel Dynamically Adaptive Simulations of Wave Problems

Anne Reinarz, Dominic E. Charrier, Michael Bader +13

ExaHyPE ("An Exascale Hyperbolic PDE Engine") is a software engine for solving systems of first-order hyperbolic partial differential equations (PDEs). Hyperbolic PDEs are typicall…

cs.MS2019

Studies on the energy and deep memory behaviour of a cache-oblivious, task-based hyperbolic PDE solver

Dominic E. Charrier, Benjamin Hazelwood, Ekaterina Tutlyaeva +5

We study the performance behaviour of a seismic simulation using the ExaHyPE engine with a specific focus on memory characteristics and energy needs. ExaHyPE combines dynamically a…

cs.DC2010

Fast GPGPU Data Rearrangement Kernels using CUDA

Michael Bader, Hans-Joachim Bungartz, Dheevatsa Mudigere +2

Many high performance-computing algorithms are bandwidth limited, hence the need for optimal data rearrangement kernels as well as their easy integration into the rest of the appli…

astro-ph.IM2018

Exploiting the Space Filling Curve Ordering of Particles in the Neighbour Search of Gadget3

Antonio Ragagnin, Nikola Tchipev, Michael Bader +2

Gadget3 is nowadays one of the most frequently used high performing parallel codes for cosmological hydrodynamical simulations. Recent analyses have shown t\ hat the Neighbour Sear…