Before revision:  We observe speedups of 2x on the Cloverleaf 2D/3D proxy application, which contain 83/141 loops respectively  .
After revision:  We observe speedups of 2x on the Cloverleaf 2D/3D proxy application, which contain 83/141 loops respectively , 3.5\times on the linear solver TeaLeaf, and 1.7\times on the compressible Navier-Stokes solver OpenSBLI. We demonstrate strong and weak scalability up to 4608 cores of CINECA's Marconi supercomputer. We also evaluate our algorithms on Intel's Knights Landing, demonstrating maintained throughput as the problem size grows beyond 16GB, and we do scaling studies up to 8704 cores .