<div dir="ltr"><div><div><div>Hi!<br><br></div>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 problem.<br>
<br></div>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.<br>
<br></div><div>I'm now suspecting that the ILU options cannot be set so as to obtain ILUTP behaviour.<br></div><div><br></div><div>What would be the recommended method to access an ILUTP preconditioner from PETSc?<br>
<br>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?<br>
<br></div><div>Best regards<br>Torquil Sørensen<br><br></div></div>