The determinant can be computed from the LU factorization of a matrix
(http://en.wikipedia.org/wiki/Determinant). We don't have code that can
automatically pull out the appropriate values but it shouldn't be too hard to
write. If you are running sequentially then you can just use PETS's LU
factorization and then access the underlying data structure to pull out the
diagonal entries. If you want in parallel you could use MUMPS or SuperLU_Dist()
and then pull out the diagonals from their data structure (a bit of work to
figure out how to do that).
Is this what you are looking for?
Barry
On Jul 28, 2010, at 10:52 AM, ?????? ???? (Nemanja Ilic) wrote:
> Hello,
>
> Is there an easy way to calculate determinant in PETSc? I thought it could be
> done with MatNorm(), but I get bad results. Is there another way?
>
> Thank you in advance,
> Best regards,
> Nemanja Ilic
