Work

p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures

Year: 2025

Type: article

Abstract: Heterogeneous computing and exploiting integrated CPU–GPU architectures has become a clear current trend since the flattening of Moore’s Law. In this work, we propose a numerical and algorithmic re-design of a p-adaptive quadrature-free discontinuous Galerkin (DG) method for the shallow water equations. Our new approach separates the computations of the non-adaptive (lower-order) and adaptive (higher-order) parts of the discretization from each other. Thereby, we can overlap computations of the lower-order and the higher-order DG solution components. Furthermore, we investigate execution times of main computational kernels and use automatic code generation to optimize their distribution between the CPU and GPU. Several setups, including a prototype of a tsunami simulation in a tide-driven flow scenario, are investigated, and the results show that significant performance improvements can be achieved in suitable setups.

Source: GEM - International Journal on Geomathematics

DOI: 10.1007/s13137-025-00267-2

Authors: Sara Faghih-Naini, Vadym Aizinger, Sebastian Kuckuk, Richard Angersbach, Harald Köstler

Institutions: Friedrich-Alexander-Universität Erlangen-Nürnberg, University of Bayreuth

Language: English

Topic: Advanced Numerical Methods in Computational Mathematics

Subfield: Computational Mechanics

Field: Engineering

Domain: Physical Sciences

Open Access status: hybrid

Estimated APC paid: $2,490

Funders: Friedrich-Alexander-Universität Erlangen-Nürnberg, Deutsche Forschungsgemeinschaft, Mercator Research Center Ruhr