Spectral-element simulations on quadrilaterals and hexahedra rely on the Gauss–Lobatto (GL) integration rule to enable explicit simulations with optimal spatial convergence rates. In this work, it is ...