Some changes in petsc_solve_petsc_csr routines:
 petsc_solve_setup_petsc_csr now returns a petsc matrix and petsc_numbering object. Usually this will just be matrix%M and matrix%row_numbering. But this will be different when I have implemented reordering for petsc_csr_matrix type.
Remove some unused variables. In particular this removes an unused argument from setup_ksp_from_options.