Automatic Generation of Matrix-Free Routines for PDE Solvers with Devito via PETSc

We extend Devito, a domain-specific language (DSL) and compiler for optimised finite differences, by adding matrix-free iterative solver support through an interface with the Portable Extensible Toolkit for Scientific Computing (PETSc). This enhancement expands the capabilities of the symbolic DSL, enabling it to address a wider range of high-performance computing challenges beyond its original scope, including incompressible flow problems in computational fluid dynamics.