[petsc-users] having issues with nullspace

[Jed Brown]

Mohammad Mirzadeh <mirzadeh at gmail.com> 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.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140306/577e0ef1/attachment.pgp>