Add function asInverseMatrix
There is a function asTransformedMatrix in the cpppetsc library. This function allows to lazily apply a transform to the matrix when it is multiplied by a vector. A similar function asInverseMatrix should allow to compute A^-1 x using PETSc's KSP based on matrix-by-vector multiplications.