Devito: a scalable and portable stencil DSL and compiler

Devito is a Python package to implement optimized stencil computation (e.g., finite differences, image processing, machine learning) from high-level symbolic problem definitions. Devito builds on SymPy and employs automated code generation and just-in-time compilation to execute optimized computational kernels on several computer platforms, including CPUs, GPUs, and clusters thereof.
Symbolic computation is a powerful tool that allows users to:
Devito is primarily designed to create wave propagation kernels for use in seismic inversion problems. A tutorial for the generation of a modelling operator using an acoustic wave equation can be found here. A paper outlining the symbolic API and the verification procedures of the acoustic operator can be found here (Louboutin et al. 2019) and the internals of the compiler are described in this paper (Luporini et al. 2020).
You can find below an introductory tutorial to Devito.