Mohammad Mirzadeh <[email protected]> writes: > Yes. To be precise this is the set of functions I call: > > ierr = MatNullSpaceCreate(mpicomm, PETSC_FALSE, 1, &null_space, > &A_null_space); CHKERRXX(ierr); > > ierr = MatSetNullSpace(A, A_null_space); CHKERRXX(ierr); > > ierr = KSPSetNullSpace(ksp, A_null_space); CHKERRXX(ierr); > > ierr = MatNullSpaceRemove(A_null_space, rhs_, NULL); CHKERRXX(ierr);
Is the matrix symmetric? If not, the right and left null spaces could be different, in which case this system might be inconsistent.
pgpw51YexNm9E.pgp
Description: PGP signature
