Fast and Accurate Finite-Element Multigrid Solvers for PDE Simulations on GPU Clusters