> You can use -pc_type lu to test for a rank deficient preconditioning matrix. > > Matt I get: [0]PETSC ERROR: No support for this operation for this object type! [0]PETSC ERROR: Matrix format mpiaij does not have a built-in PETSc LU! Dominik