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>
