On Wed, Aug 8, 2018 at 3:51 AM, Barna Becsek <barnabec...@gmail.com> wrote:

> Hello,
>
> is there a way to compute a SparseMatrix’s determinant? I would like to
> solve a linear system but my final residual always equals exactly 0 after
> just one iteration (where the solve terminates). Hence, I would like to
> check the matrix for singularity.
>

If you are using PETSc, for a small enough problem you could run with the
"-pc_type svd -pc_svd_monitor" command line options and look at the
singular values. If any of them are zero or nearly zero, the system is
singular.

-- 
John
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to