<div dir="ltr"><div dir="ltr">On Wed, May 27, 2020 at 7:09 PM Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><br><div><br></div><div>Matt: Could I use BAIJ with Plex?</div></div></blockquote><div><br></div><div>Plex does this automatically if you have blocks.</div></div></div></blockquote><div><br></div><div>I use DMCreateMatrix with forest or plex and I seem to get AIJ matrices. Where does Plex get the block size? </div><div><br></div><div>I have not yet verified that bs is set in this matrix.</div></div></div>
</blockquote></div><br clear="all"><div>I think I may know what your problem is. Plex evaluates the blocksize by looking for an equal number of dofs</div><div>on each point. This is sufficient, but not necessary. If you are using higher order methods, there is block structure</div><div>there that I will not see.</div><div><br></div><div>Jed, is there an obvious way to see that structure that I am missing?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><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><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>