<div dir="ltr"> I have DM_BOUNDARY_NONE.<br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 27, 2020 at 7:22 PM Jed Brown <<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> writes:<br>
<br>
> On Wed, May 27, 2020 at 7:09 PM Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br>
><br>
>><br>
>><br>
>>>><br>
>>>> Matt: Could I use BAIJ with Plex?<br>
>>>><br>
>>><br>
>>> Plex does this automatically if you have blocks.<br>
>>><br>
>><br>
>> I use DMCreateMatrix with forest or plex and I seem to get AIJ matrices.<br>
>> Where does Plex get the block size?<br>
>><br>
>> I have not yet verified that bs is set in this matrix.<br>
>><br>
><br>
> I think I may know what your problem is. Plex evaluates the blocksize by<br>
> looking for an equal number of dofs<br>
> on each point. This is sufficient, but not necessary. If you are using<br>
> higher order methods, there is block structure<br>
> there that I will not see.<br>
><br>
> Jed, is there an obvious way to see that structure that I am missing?<br>
<br>
gcd of all the point sizes.<br>
<br>
Note that you don't have constant block sizes if you eliminate some<br>
fields/components in boundary conditions.<br>
</blockquote></div>