Applying Fast Fourier Transforms to Accelerate Spatially and Temporally Inhomogeneous Stencil Computations

Stencil computations are essential for simulating the evolution of physical systems across multi-dimensional grids over multiple timesteps. State-of-the-art techniques in this field fall into three major groups: tiled looping algorithms, divide-and-conquer trapezoidal algorithms, and Krylov subspace methods.