<div dir="ltr">Thanks, Barry.<div><br></div><div>More questions. Is bs usually the number of unknowns associated with a grid point. Do I need to change MatSetValues if use a <span style="font-size:13px">MATSEQBAIJ </span><span style="font-size:13px">matrix?</span></div><div><span style="font-size:13px"><br></span></div><div>Fande,</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 1, 2015 at 2:36 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:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Fande,<br>
<br>
    The MATSEQBAIJ matrix automatically uses a point block ILU; with the block size of the matrix. So if your matrix naturally has a block size bs if you use the MATSEQBAIJ matrix type it will be point block ILU of block size bs.<br>
<br>
    We don't have a variable point block ILU in PETSc for SeqAIJ. It is possible to write but would be rather involved.<br>
<span class="HOEnZb"><font color="#888888"><br>
   Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
> On Mar 1, 2015, at 1:47 PM, Fande Kong <<a href="mailto:fd.kong@siat.ac.cn">fd.kong@siat.ac.cn</a>> wrote:<br>
><br>
> Hi all,<br>
><br>
> I have a code with ILU as subdomain solver, and it works well. Now I want to try a block version of ILU, there are some easy ways to implement this?<br>
><br>
> Fande,<br>
<br>
<br>
</div></div></blockquote></div><br></div>