Fast and scalable solvers for high-order hybridized discontinuous Galerkin methods with applications to fluid dynamics and magnetohydrodynamics