<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">All I know is that we use SPARSKIT2's iluk.<div><br></div><div>I am aware that the code snippet I gave only works in serial.<br><div>
<div style="color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>--</div><div> Hugo Gagnon</div></div>
</div>
<br><div><div>On 2013-04-15, at 5:15 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">Hugo Gagnon <<a href="mailto:opensource.petsc@user.fastmail.fm">opensource.petsc@user.fastmail.fm</a>> writes:<br><br><blockquote type="cite">For the problem I'm describing my serial in-house solver does not work<br>with ILU(0) but works with ILU(3). I have no option to run Jacobi.<br>When I apply the same problem to PETSc's PC solver with ILU(3) in<br>serial I get KSP_DIVERGED_INDEFINITE_PC<br></blockquote><br>Does your in-house ILU(3) use a different ordering? What shift scheme<br>does it use?<br><br><blockquote type="cite">on the first iteration (in MPI the solution somewhat converges but<br>very slowly).<br><br>call KSPGetPC(Pksp,Ppc,Pierr)<br>call PCSetType(Ppc,PCILU,Pierr)<br>call PCFactorSetLevels(Ppc,3,Pierr)<br><br>This effectively changes the fill level from 0 to 3, right?<br></blockquote><br>This only works in serial. Check the -ksp_view output to see what is<br>done. You should just call KSPSetFromOptions() and use run-time options<br>to configure the solver. You can do it from code later, but writing<br>code is slow to figure out what works.<br></blockquote></div><br></div></body></html>