High order diamond differencing schemes

We present a generalization if the diamond differencing scheme to high-order spatial orders in one- and two-dimensional Cartesian geometries. Unlike existing variable-order nodal schemes our approach reduces to the conventional low-order diamond differencing scheme in the linear case and feature super convergence characteristics at all orders. This approach is demonstrated on one- and two-dimensional benchmark problems and is compared to spherical harmonics reference solutions.