[petsc-users] ILUTP

Torquil Macdonald Sørensen torquil at gmail.com
Sat Oct 26 15:05:04 CDT 2013


I have a linear problem Ax=b, with complex values, that is solved very well
using the ILUTP and GMRES implementations in GMM++. None of the other
preconditioners available in GMM++ would work (e.g. ILU and ILUT) for this

Then I tried the same problem using ILU and GMRES in PETSc, with no
success, despite a lot of attempts at adjusting ILU settings. I always end
up with gigantic residual norm values. The same PETSc program works well
when I apply it to a different matrix A.

I'm now suspecting that the ILU options cannot be set so as to obtain ILUTP

What would be the recommended method to access an ILUTP preconditioner from

According to the PETSc website, a preconditioner named ILUDT is available
by using the external package Hypre, but I had to deselect Hypre during the
PETSc build due to my use of complex numbers... So do you guys think that I
should transform everything to a real representation and try Hypre/ILUDT?

Best regards
Torquil Sørensen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20131026/b1bdd4fc/attachment.html>

More information about the petsc-users mailing list