Block orthogonalization should use double precision for local sums.
Multigrid preconditioner now uses sloppy precision, this enables the use of mixed-precision.