“…This is attributed to the original work of Lax and Wendroff from 1960 where they construct a second-order solver by incorporating the second derivative of the PDE into their method [33]. More recently, higher order (i.e., solvers with order greater than two) versions of these solvers have been investigated for finite volume [26], finite difference [41,32,46,12,14,52], and discontinuous Galerkin discretizations [40,22,36]. A large community centered around Arbitrary DERivative (ADER) discretizations has been very successful with constructing arbitrary order explicit solvers for hyperbolic problems in this category [48,51,15,7], and much of their work relies on symbolic software to generate their code base.…”