Interface to Petsc data structures

The coefficient matrix used in the preconditioner setup (section 2.2), and the vectors involved in the preconditioner solve need to be in a suitable form. The vector can most likely be transformed in a single (essentially costless) instruction; for the matrix it takes some work, and in fact an extra copy of the matrix will be built.