<div dir="ltr"><div>Yes, I was suprised that I can still use MatSetValuesBlocked after I did MatConvert MATBAIJ to MATAIJ.</div><div> </div><div>Then I just found the mat type is still seqbaij which means the MatConvert didn't work.</div><div> </div><div>This is how I did MatConvert</div><div> </div><div>MatConvert(J, MATAIJ, MAT_REUSE_MATRIX, & J);</div><div> </div><div>where J is my Jacobian matrix.</div><div> </div><div>Thanks,</div><div> </div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Mar 18, 2015 at 2:24 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote style="margin:0px 0px 0px 0.8ex;padding-left:1ex;border-left-color:rgb(204,204,204);border-left-width:1px;border-left-style:solid" class="gmail_quote"><br>
  Are you always using MatSetValuesBlocked()?<br>
<br>
  Did you look at the three matrices (before passing to the solve), do they all have the exact same nonzero structure and entries?<br>
<span><font color="#888888"><br>
  Barry<br>
</font></span><div><div><br>
> On Mar 18, 2015, at 1:45 PM, Chung-Kan Huang <<a href="mailto:ckhuangf@gmail.com" target="_blank">ckhuangf@gmail.com</a>> wrote:<br>
><br>
> So go back to my original question,<br>
><br>
> I found the resuls look like<br>
> MATBAIJ == MATBAIJ2MATAIJ != MATAIJ<br>
> as I was expected to see<br>
> MATBAIJ != MATBAIJ2MATAIJ == MATAIJ<br>
><br>
> because I thought MATBAIJ2MATAIJ & MATAIJ should use the same ILU factorization.<br>
><br>
> On Wed, Mar 18, 2015 at 1:33 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
>    As I think I said earlier: ILU will behave differently with BAIJ and AIJ format. With BAIJ it uses a point block factorization (doing an ILU treating each block in the matrix as a single point (and doing a dense factorization for the inverse) while AIJ ignores the block size and just uses a point ILU. So in general the convergence with BAIJ vs AIJ and ILU will be different; usually BAIJ will converge better but it need not.<br>
><br>
>    Barry<br>
><br>
> > On Mar 18, 2015, at 1:23 PM, Chung-Kan Huang <<a href="mailto:ckhuangf@gmail.com" target="_blank">ckhuangf@gmail.com</a>> wrote:<br>
> ><br>
> > Sorry I didn't state my problem more accurate.  The entries in the input matrixes and vectors​ are correct -- they have right values and correct nonzero structure.<br>
> > The KSP was setup so I can solve the lienar systems for an unsteady state convective problem.<br>
> > I use KSPBCGS + PCILU and so far I only try sequential on a very small problem.<br>
> > I wonder if blocked ilu is actually used in  MATBAIJ2MATAIJ case.<br>
> ><br>
> > Thanks,<br>
> ><br>
> ><br>
><br>
><br>
><br>
><br>
> --<br>
> Cheers<br>
><br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div><p><strong>Cheers</strong></p></div>
</div></div>