Solving the Laplace equation on hp-adaptive meshes with a MatrixFree hybrid
multigrid approach on thousands of processors.
