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

[Barry Smith]

> On Sep 30, 2015, at 10:39 AM, Matthew Knepley <knepley at gmail.com> wrote:
> 
> On Wed, Sep 30, 2015 at 10:29 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
> > 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.
> 
> What matrix subclasses do we have that are not factorizations?

MatCreate_MPIAIJCR
MatCreate_MPIAIJPERM
MatCreate_MPIAIJCUSP
MatCreate_MPIAIJCUSPARSE
MatCreate_MPIAIJViennaCL



> 
>   Matt
>  
>   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
> 
> 
> 
> 
> -- 
> 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