[petsc-dev] Adding support for new sparse BLAS inspector-executor routines in MKL

[Barry Smith]

> On Sep 30, 2015, at 9:58 AM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Wed, Sep 30, 2015 at 9:48 AM, Jed Brown <jed at jedbrown.org> wrote:
> Matthew Knepley <knepley at gmail.com> writes:
> > Richard, this is somewhat subtle since PETSc's inheritance model is sucky
> > and somewhat broken. Look at
> > the factorization classes, AIJMUMPS etc., before coding anything.
> 
> Richard wants to use the MKL implementations of stuff like MatMult and
> MatAXPY, not just factorization.  So it doesn't look like AIJMUMPS.
> 
> It will look like that in the sense that he wants to keep some data structures and
> not others, and override some functions and not others. AIJMUMPS shows you
> how to keep things in sync since there are too many places where they can go
> wrong. Factorization is beside the point.

   No it isn't. We have a special infrastructure for factorizations that is distinct from making subclasses of mat subclasses.

  Barry

> 
>    Matt
> 
> -- 
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
> -- Norbert Wiener