Towards generating efficient flow solvers with the ExaStencils approach

Abstract

Summary ExaStencils aims at providing intuitive interfaces for the specification of numerical problems and resulting solvers, particularly those from the class of (geometric) multigrid methods. It envisions a multi‐layered domain‐specific language and a sophisticated code generation framework ultimately emitting source code in a chosen target language. We present our recent advances in fully generating solvers applied to 3D fluid mechanics for nonisothermal/non‐Newtonian flows. In detail, a system of time‐dependent, nonlinear partial differential equations is discretized on a cubic, nonuniform, and staggered grid using finite volumes. We examine the contained problem of coupled Navier‐Stokes and temperature equations, which are linearized and solved using the SIMPLE algorithm and geometric multigrid solvers, as well as the incorporation of non‐Newtonian properties. Furthermore, we provide details on necessary extensions to our domain‐specific language and code generation framework, in particular, those concerning the handling of boundary conditions, support for nonequidistant staggered grids, and supplying specialized functions to express operations reoccurring in the scope of finite volume discretizations. Many of these enhancements are generalizable and thus suitable for utilization in similar projects. Lastly, we demonstrate the applicability of our code generation approach by providing convincing performance results for fully generated and automatically parallelized solvers.

Developing performance portable plasma edge simulations: A survey

“…OPS [89] is a multi-block structured-mesh DSL embedded in both Fortran and C/C++, targeting CPUs, GPUs and clusters of CPUs/GPUs -it uses a source-to-source translation strategy to generate code for a variety of parallelisations. ExaSlang [90] is part of a larger European project, ExaStencils [91], which allows the description of PDE computations at many levelsincluding at the level of structured-mesh stencil algorithms. It is embedded in Scala and targets MPI and CPUs, with limited GPU support.…”

Productivity, performance, and portability for computational fluid dynamics applications

“…OPS [58] (active) is a multi-block structured mesh DSL embedded in both Fortran and C/C++, targeting CPUs, GPUs and MPI -it uses a source-to-source translation strategy to generate code for a variety of parallelisations. ExaSlang [59] (active) is part of a larger European project, Exastencils, which allows the description of PDE computations at many levels -including at the level of structured-mesh stencil algorithms. It is embedded in Scala, and targets MPI and CPUs, with limited GPU support.…”

Code generation approaches for parallel geometric multigrid solvers

“…Examples for the latter are address pre-calculation, loop transformation, a sophisticated common sub-expression elimination (CSE) [34] and vectorization. The emitted code can be executed on various platforms ranging from traditional CPU [36,35] and hybrid CPU-GPU [37] clusters, ARM-based architectures [39] to reconfigurable hardware [51,50].…”