https://hgpu.org/?p=11080
Multilayered Abstractions for Partial Differential Equations