<div>Dear PETSc users,<br></div><div><br></div><div>I would like to know if it is possible to provide PETSc with a ILU-preconditioner matrix computed externally beforehand.<br></div><div><div>I tried and built a PCSHELL, to which I pass the externally-computed preconditioner matrix as "Pmat" using the KSPSetOperators function. Then I wanted to use that Pmat in PCApply by calling MatSolve as it seems to be done in the ILU case. Though, this fails since the mat->ops->solve (with mat being my PC Matrix) is a null pointer.<br></div><div><br></div><div>I guess the way I set the matrix (as a MATSEQAIJ) is not sufficient for PETSc to know what function to use as MatSolve.<br></div><div><br></div><div>How could I achieve providing my own ILU-decomposed matrix and feed PETSc's PCG with it? Is it actually possible?<br></div><div><br></div><div>Thanks for your help,<br><br></div></div><div class="protonmail_signature_block"><div class="protonmail_signature_block-user"><div><span class="colour" style="color:rgb(107,107,107);">Milan Pelletier</span><br></div></div><div class="protonmail_signature_block-proton protonmail_signature_block-empty"></div></div><div><br></div>