<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Sep 30, 2015 at 9:48 AM, 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">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
> Richard, this is somewhat subtle since PETSc's inheritance model is sucky<br>
> and somewhat broken. Look at<br>
> the factorization classes, AIJMUMPS etc., before coding anything.<br>
<br>
Richard wants to use the MKL implementations of stuff like MatMult and<br>
MatAXPY, not just factorization.  So it doesn't look like AIJMUMPS.<br>
</blockquote></div><br>It will look like that in the sense that he wants to keep some data structures and</div><div class="gmail_extra">not others, and override some functions and not others. AIJMUMPS shows you</div><div class="gmail_extra">how to keep things in sync since there are too many places where they can go</div><div class="gmail_extra">wrong. Factorization is beside the point.</div><div class="gmail_extra"><br></div><div class="gmail_extra">   Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature">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></div>