Efficient mesh management in Firedrake using PETSc-DMPlex

Abstract

The use of composable abstractions allows the application of new and established algorithms to a wide range of problems while automatically inheriting the benefits of well-known performance optimisations. This work highlights the composition of the PETSc DMPlex domain topology abstraction with the Firedrake automated finite element system to create a PDE solving environment that combines expressiveness, flexibility and high performance. We describe how Firedrake utilises DMPlex to provide the indirection maps required for finite element assembly, while supporting various mesh input formats and runtime domain decomposition. In particular, we describe how DMPlex and its accompanying data structures allow the generic creation of user-defined discretisations, while utilising data layout optimisations that improve cache coherency and ensure overlapped communication during assembly computation.

12 pages, 6 figures, submitted to SISC CSE Special Issue

Cited by in corpus (10)

• Firedrake: automating the finite element method by composing abstractions
• A dynamically load-balanced parallel $p$-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation
• A structure-exploiting numbering algorithm for finite elements on extruded meshes, and its performance evaluation in Firedrake
• Fully Parallel Mesh I/O using PETSc DMPlex with an Application to Waveform Modeling
• Monolithic multigrid for implicit Runge-Kutta discretizations of incompressible fluid flow
• Composable block solvers for the four-field double porosity/permeability model
• Droplet formation simulation using mixed finite elements
• Conservative Projection Between Finite Element and Particle Bases
• Efficient N-to-M Checkpointing Algorithm for Finite Element Simulations
• Structure Matters: A Scale-Resolved Numerical Operando Approach for Lithium-Sulfur Batteries