Hi, How does PETSc handle the solution of a singular system Ax=b which also happens to be inconsistent? More specifically: is b orthogonalized against Null(A) internally by PETSc? Thanks Håkan Jakobsson