Firedrake: automating the finite element method by composing abstractions

F Rathgeber, DA Ham, L Mitchell, M Lange… - ACM Transactions on …, 2016 - dl.acm.org
Firedrake is a new tool for automating the numerical solution of partial differential equations.
Firedrake adopts the domain-specific language for the finite element method of the FEniCS …

Modular and flexible spectral-element waveform modelling in two and three dimensions

M Afanasiev, C Boehm, M van Driel… - Geophysical Journal …, 2019 - academic.oup.com
In this paper, we present a series of mathematical abstractions for seismologically relevant
wave equations discretized using finite-element methods, and demonstrate how these …

The PetscSF scalable communication layer

J Zhang, J Brown, S Balay… - … on Parallel and …, 2021 - ieeexplore.ieee.org
PetscSF, the communication component of the Portable, Extensible Toolkit for Scientific
Computation (PETSc), is designed to provide PETSc's communication infrastructure suitable …

Monolithic multigrid for implicit Runge–Kutta discretizations of incompressible fluid flow

R Abu-Labdeh, S MacLachlan, PE Farrell - Journal of Computational …, 2023 - Elsevier
Most research on preconditioners for time-dependent PDEs has focused on implicit multi-
step or diagonally-implicit multi-stage temporal discretizations. In this paper, we consider …

A dynamically load-balanced parallel p-adaptive implicit high-order flux reconstruction method for under-resolved turbulence simulation

L Wang, MK Gobbert, M Yu - Journal of Computational Physics, 2020 - Elsevier
We present a dynamically load-balanced parallel p-adaptive implicit high-order flux
reconstruction method for under-resolved turbulence simulation. The high-order explicit first …

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

GT Bercea, ATT McRae, DA Ham… - Geoscientific Model …, 2016 - gmd.copernicus.org
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 …

Monolithic multigrid methods for magnetohydrodynamics

JH Adler, TR Benson, EC Cyr, PE Farrell… - SIAM Journal on …, 2021 - SIAM
The magnetohydrodynamics equations model a wide range of plasma physics applications
and are characterized by a nonlinear system of partial differential equations that strongly …

Fully parallel mesh I/O using PETSc DMPlex with an application to waveform modeling

V Hapla, MG Knepley, M Afanasiev, C Boehm… - SIAM Journal on …, 2021 - SIAM
Large-scale PDE simulations using high-order finite-element methods on unstructured
meshes are an indispensable tool in science and engineering. The widely used open …

Conservative projection between finite element and particle bases

JV Pusztay, MG Knepley, MF Adams - SIAM Journal on Scientific Computing, 2022 - SIAM
Particle-in-cell methods employ particle representations of unknown fields but also employ
continuum fields for other parts of the problem. Thus projection between particle and …

Irksome: Automating Runge–Kutta time-step** for finite element methods

PE Farrell, RC Kirby… - ACM Transactions on …, 2021 - dl.acm.org
While implicit Runge–Kutta (RK) methods possess high order accuracy and important
stability properties, implementation difficulties and the high expense of solving the coupled …