[petsc-users] Calculate determinant in PETSc - MatAssemblyEnd breaks

[Немања Илић (Nemanja Ilic)]

Hello,

I did the LU factorization and then pulled out diagonal from the matrix, but the code breaks. It breaks on the MatAssemblyEnd after the followin lines:
"ierr = MatLUFactor(Ap, perm, iperm, &info);CHKERRQ(ierr);
ierr = MatSetUnfactored(Ap);
ierr = MatAssemblyBegin(Ap, MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);
ierr = MatAssemblyEnd(Ap, MAT_FINAL_ASSEMBLY); CHKERRQ(ierr);"

I attached my code.

Thank you in advance for any help.

Best regards,
Nemanja
-------------- next part --------------
A non-text attachment was scrubbed...
Name: math_get_matrix_determinant.c
Type: text/x-csrc
Size: 4435 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100728/01a711e1/attachment.c>