<div dir="ltr">I guess I missed out an important piece of information, due to legacy code issues, I am working with Petsc version 2.3.2. Unfortunately, it doesn't have <span style="font-size:12.8px">PCFactorSetMatOrderingType(). Thanks for the advice though. I still don't understand, why would it ask for so much memory, it is almost 16 Petabytes! Any other advice would be greatly appreciated.</span><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Regards,</span></div><div><span style="font-size:12.8px">Rachit Prasad</span></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Sep 28, 2017 at 5:13 PM, Kong, Fande <span dir="ltr"><<a href="mailto:fande.kong@inl.gov" target="_blank">fande.kong@inl.gov</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">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.<span class="HOEnZb"><font color="#888888"><br><div><br></div><div>Fande,</div></font></span></div><div class="HOEnZb"><div class="h5"><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 -<a href="tel:(214)%20748-3648" value="+12147483648" target="_blank">2147483648</a> Memory used by process -<a href="tel:(214)%20748-3648" value="+12147483648" target="_blank">2147483648</a></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>
</div></div></blockquote></div><br></div>