English

Domain-specific implementation of high order Discontinuous Galerkin methods in spherical geometry

Numerical Analysis 2023-05-23 v3 Distributed, Parallel, and Cluster Computing Numerical Analysis Programming Languages

Abstract

We assess two domain-specific languages included in the GridTools ecosystem as tools for implementing a high-order Discontinuous Galerkin discretization of the shallow water equations. Equations in spherical geometry are considered, thus providing a blueprint for the application of domain-specific languages to the development of global atmospheric models. The results demonstrate that domain-specific languages designed for finite difference/volume methods can be successfully extended to implement a Discontinuous Galerkin solver.

Keywords

Cite

@article{arxiv.2303.11767,
  title  = {Domain-specific implementation of high order Discontinuous Galerkin methods in spherical geometry},
  author = {Kalman Szenes and Niccolò Discacciati and Luca Bonaventura and William Sawyer},
  journal= {arXiv preprint arXiv:2303.11767},
  year   = {2023}
}

Comments

24 pages; will be subsequently submitted to Computer Physics Communications