https://hgpu.org/?p=3638
A pseudospectral matrix method for time-dependent tensor fields on a spherical shell