[petsc-dev] any way to detect if a matrix is singular?

[Kong, Fande]

Hi Developers,

I am trying to use the native LU in PETSc to compute a matrix (arising from
a fluid-structure interaction problem). I encounter a zero pivot. The
matrix gives the wrong Newton direction when using "-pc_factor_shift_type
 nonzero -pc_factor_shift_amount 1".  I hence guess the matrix is singular.

Any easy way to detect if or not the matrix is singular??


Fande,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20171010/69f96270/attachment.html>