<div dir="ltr">Thanks, Matt.<div><br></div><div>Could bs be different from rows. I am solving a fluid-structure problem, which has different bs in fluid and solid parts.</div><div><br></div><div>Fande,</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Mar 1, 2015 at 3:57 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</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"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Sun, Mar 1, 2015 at 4:55 PM, Fande Kong <span dir="ltr"><<a href="mailto:fd.kong@siat.ac.cn" target="_blank">fd.kong@siat.ac.cn</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">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></blockquote><div><br></div></span><div>No, but it does enable you to use MatSetValuesBlock() if you want.</div><div><br></div><div>   Matt</div><span class=""><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><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><font color="#888888"><br>
   Barry<br>
</font></span><div><div><br>
> On Mar 1, 2015, at 1:47 PM, Fande Kong <<a href="mailto:fd.kong@siat.ac.cn" target="_blank">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>
</blockquote></span></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</font></span></div></div>
</blockquote></div><br></div>