<div dir="ltr">Calling PCFactorSetMatOrderingType() (or command line option: -pc_factor_mat_ordering_type) to use RCM or 1WD usually helps me a lot. My application is based on highly-unstructured meshes.<br><div><br></div><div>Fande,</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 28, 2017 at 2:53 PM, Rachit Prasad <span dir="ltr"><<a href="mailto:rachitp@vt.edu" target="_blank">rachitp@vt.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi,<div><br></div><div>I am trying to solve a highly ill-conditioned matrix and am unable to get convergence when using ILU(0) or Jacobi and SOR preconditioners. I tried to implement a higher order ILU by increasing the level of fill-in to 1 by calling</div><div><br></div><div><div>call PCFactorSetLevels(pc,1,ierr)</div></div><div><br></div><div>However, when I do that I get the following memory error:</div><div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------<wbr>------</div><div>[0]PETSC ERROR: Out of memory. This could be due to allocating</div><div>[0]PETSC ERROR: too large an object or bleeding by not properly</div><div>[0]PETSC ERROR: destroying unneeded objects.</div><div>[0]PETSC ERROR: Memory allocated -2147483648 Memory used by process -2147483648</div><div>[0]PETSC ERROR: Try running with -malloc_dump or -malloc_log for info.</div><div>[0]PETSC ERROR: Memory requested 18446744071863943168!</div></div><div><br></div><div>While I do understand that increasing the level of fill-in will require higher memory, the memory requested seems to be way too high. Am I implementing the higher order ILU correctly? Is there any other subroutine which I need to call?</div><div><br></div><div>Regards,</div><div>Rachit</div><div><br></div></div>
</blockquote></div><br></div>