[petsc-users] Incomplete Cholesky Factorization in petsc

[Wenlong Gong]

Hello,

I'm trying to use the Incomplete Cholesky Factorization for a sparse matrix
in petsc. I started with a 10*10 matrix and used ksp and pc in order to get
the ICC(0) factor, with no fill-in, natural ordering. However, the returned
factor matrix does not match with the answer I got from matlab ichol()
function.

The code with the hard-coded data is attached here. I would appreciate if
anyone could help check if I did anything wrong.Please let me know if there
is easier way to get this incomplete cholesky factor. Thanks!

Best regards,
Wendy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171204/47bce13e/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: repex.c
Type: text/x-csrc
Size: 2982 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20171204/47bce13e/attachment.bin>