<div dir="ltr"><div><div>I found one issue on my side. The preallocation is not right for the BAIJ matrix.  Will this slow down MatLUFactor and MatSolve?<br><br></div>How to converge AIJ to BAIJ using a command-line option?<br><br></div>Fande,<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 7, 2017 at 3:26 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">"Kong, Fande" <<a href="mailto:fande.kong@inl.gov">fande.kong@inl.gov</a>> writes:<br>
<br>
> On Tue, Mar 7, 2017 at 3:16 PM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
><br>
>> Hong <<a href="mailto:hzhang@mcs.anl.gov">hzhang@mcs.anl.gov</a>> writes:<br>
>><br>
>> > Fande,<br>
>> > Got it. Below are what I get:<br>
>><br>
>> Is Fande using ILU(0) or ILU(k)?  (And I think it should be possible to<br>
>> get a somewhat larger benefit.)<br>
>><br>
><br>
><br>
> I am using ILU(0). Will it be much better to use ILU(k>0)?<br>
<br>
</span>It'll be slower, but might converge faster.  You asked about ILU(k) so I<br>
assumed you were interested in k>0.<br>
</blockquote></div><br></div>