Multigrid methods for large-scale optimal control problems constrained by partial differential equations