A few random fixes (!218)

The only notable change is the use use of the Lagrange derivative matrix for all DG kinds.
That's done because not all DG kinds have the derivative defined on their nodal points (`:modal_bspline1`).