<div dir="ltr">I do not mind accepting Frank's MR. However, if this is a 2x improvement for something that is often done, should we think about having<div>a nicer way to add block sizes? Jed has code for different blocks using the preprocessor (I think) in sfbasic, and we have top-level switches</div><div>in PBJACOBI and (I think) MatInvertDiagonal. Should we agree on a good strategy and use it everywhere?</div><div><br></div><div>  Matt</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 10:32 AM Mark Adams <<a href="mailto:mfadams@lbl.gov">mfadams@lbl.gov</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"><div dir="ltr">Yes, this is a useful thing to have. AIJ assembly is always a problem and this seems like a good optimization to have in there.<div>Thanks,</div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 9:47 AM Frank Bramkamp <<a href="mailto:bramkamp@nsc.liu.se" target="_blank">bramkamp@nsc.liu.se</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"><br>
Dear Mark,<br>
<br>
then I can do this.<br>
I will add it to the sequential and parallel part.<br>
Probably useful for everyone.<br>
<br>
Greetings, Frank<br>
<br>
<br>
<br>
</blockquote></div>
</blockquote></div><br clear="all"><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>