Everything uses a mask now in order to weed out Dirichlet degrees of freedom. However, matvec can still be used with optional arguments specifying the range of rows and columns to be multiplied. This will be needed later for alternating Schwarz preconditioners.