High-Order Discontinuous Galerkin Methods by GPU Metaprogramming