<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sun, Mar 4, 2018 at 6:57 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">"Zhang, Hong" <<a href="mailto:hongzhang@anl.gov">hongzhang@anl.gov</a>> writes:<br>
<br>
> We have thought about subclassing long before. As Barry suggested, an<br>
> ideal solution is to make a base class (XAIJ ?) for AIJ, BAIJ, SBAIJ,<br>
> and SELL. The base class will have diagonal and off diagonal parts<br>
> without binding their storage format to CSR-type. This requires a lot<br>
> of efforts to refactor the current AIJ. But I think we should<br>
> eventually do it.<br>
<br>
It sounds like we're mixing the concerns of parallel decomposition (for<br>
which "Seq*" isn't relevant) with the sequential format as some CSR-like<br>
variant.  And apart from a couple abstraction-breaking performance<br>
optimizations (that may not be necessary and could be done differently),<br>
they are nearly independent.<br>
<br>
Now there are some good reasons why we might want SeqAIJ to run inside a<br>
different parallel decomposition.  MATIS is one such format.  If we're<br>
serious about working with irregular problems, we'll eventually need a<br>
2D sparse decomposition.  We could get it from an existing package such<br>
as CombBLAS, but if we do preconditioners that depend on matrix format,<br>
we'll likely want it to work on some native format.<br>
<br>
Also, we need a better way to do small-size specialization.  I<br>
personally like doing it via inline functions instead of macros and<br>
copy-pasta, though the XL compiler can be offensively negligent.<br>
</blockquote></div><br>I have missed the larger point. What do you think should be done?</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><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.caam.rice.edu/~mk51/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div>
</div></div>