An implicit multigrid solver for high-order compressible flow simulations on GPUs