<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 7, 2017 at 7:37 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class="gmail-"><br>
> On Mar 7, 2017, at 4:35 PM, Kong, Fande <<a href="mailto:fande.kong@inl.gov">fande.kong@inl.gov</a>> wrote:<br>
><br>
> 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>
</span>  No, but you should still fix it.<br>
<span class="gmail-"><br>
><br>
> How to converge AIJ to BAIJ using a command-line option?<br>
<br>
</span>   Instead of using MatCreateSeq/MPIAIJ() at the command line you would use<br>
<br>
   MatCreate()<br>
   MatSetSizes()<br>
   MatSetBlockSize()<br>
   MatSetFromOptions()<br></blockquote><div><br></div><div> MatSetFromOptions() has to be called before "MatXXXXSetPreallocation"?  What happens if I call MatSetFromOptions() right after "MatXXXXSetPreallocation"?</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
   MatMPIAIJSetPreallocation()<br>
   MatMPIBAIJSetPreallocation() and any other preallocations you want<br>
   MatSetValues.....<wbr>MatAssemblyBegin/End()<br>
<br>
   Then you can use -mat_type baij or aij to set the type.<br>
<span class="gmail-HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="gmail-HOEnZb"><div class="gmail-h5"><br>
><br>
> Fande,<br>
><br>
> On Tue, Mar 7, 2017 at 3:26 PM, Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br>
> "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>
> It'll be slower, but might converge faster.  You asked about ILU(k) so I<br>
> assumed you were interested in k>0.<br>
><br>
<br>
</div></div></blockquote></div><br></div></div>