Implemented operator \hat{H} for pressure preconditioner

The operator \hat{H} can be applied to a function and the vertical-diagonal \hat{H}_z = \diag_h{\hat{H}} can be extracted as a banded matrix.