Optimising finite-difference methods for PDEs through parameterised time-tiling in Devito

06/21/2018
by   Nicholas Sim, et al.
0

Finite-difference methods are widely used in solving partial differential equations. In a large problem set, approximations can take days or weeks to evaluate, yet the bulk of computation may occur within a single loop nest. The modelling process for researchers is not straightforward either, requiring models with differential equations to be translated into stencil kernels, then optimised separately. One tool that seeks to speed up and eliminate mistakes from this tedious procedure is Devito, used to efficiently employ finite-difference methods. In this work, we implement time-tiling, a loop nest optimisation, in Devito yielding a decrease in runtime of up to 45 from the acoustic wave equation family, widely used in Devito's target domain of seismic imaging. We present an estimator for arithmetic intensity under time-tiling and a model to predict runtime improvements in stencil computations. We also consider generalisation of time-tiling to imperfect loop nests, a less widely studied problem.

READ FULL TEXT
research
01/20/2022

Finite difference and finite element methods for partial differential equations on fractals

In this paper, we present numerical procedures to compute solutions of p...
research
02/07/2020

FiniteNet: A Fully Convolutional LSTM Network Architecture for Time-Dependent Partial Differential Equations

In this work, we present a machine learning approach for reducing the er...
research
02/14/2017

Symbolic Solutions of Simultaneous First-order PDEs in One Unknown

We propose and implement an algorithm for solving an overdetermined syst...
research
05/06/2022

On the order of accuracy for finite difference approximations of partial differential equations using stencil composition

Stencil composition uses the idea of function composition, wherein two s...
research
04/29/2019

Algorithmic approach to strong consistency analysis of finite difference approximations to PDE systems

For a wide class of polynomially nonlinear systems of partial differenti...
research
10/18/2022

Galois Groups of Linear Difference-Differential Equations

We study the relation between the Galois group G of a linear difference-...
research
10/20/2020

Temporal blocking of finite-difference stencil operators with sparse "off-the-grid" sources

Stencil kernels dominate a range of scientific applications, including s...

Please sign up or login with your details

Forgot password? Click here to reset