Pressure Mass Matrix

I updated the code to produce a (stokes) pressure mass matrix which is
used as a preconditioner for the Schur complement.