A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake

Gheorghe-Teodor Bercea, Andrew T T McRae, David A Ham, Lawrence Mitchell, Florian Rathgeber, Luigi Nardi, Fabio Luporini, Paul H J Kelly

Research output: Contribution to journalArticlepeer-review

21 Citations (SciVal)
177 Downloads (Pure)

Abstract

We present a generic algorithm for numbering and then efficiently iterating over the data values attached to an extruded mesh. An extruded mesh is formed by replicating an existing mesh, assumed to be unstructured, to form layers of prismatic cells. Applications of extruded meshes include, but are not limited to, the representation of 3D high aspect ratio domains employed by geophysical finite element simulations. These meshes are structured in the extruded direction. The algorithm presented here exploits this structure to avoid the performance penalty traditionally associated with unstructured meshes. We evaluate the implementation of this algorithm in the Firedrake finite element system on a range of low compute intensity operations which constitute worst cases for data layout performance exploration. The experiments show that having structure along the extruded direction enables the cost of the indirect data accesses to be amortized after 10-20 layers as long as the underlying mesh is well-ordered. We characterise the resulting spatial and temporal reuse in a representative set of both continuous-Galerkin and discontinuous-Galerkin discretisations. On meshes with realistic numbers of layers the performance achieved is between 70% and 90% of a theoretical hardware-specific limit.
Original languageEnglish
Pages (from-to)3803-3815
JournalGeoscientific Model Development
Volume9
Issue number10
DOIs
Publication statusPublished - 27 Oct 2016

Fingerprint

Dive into the research topics of 'A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake'. Together they form a unique fingerprint.

Cite this