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

Matthew Knepley knepley at gmail.com
Wed Sep 30 10:39:37 CDT 2015


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?

  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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20150930/7812a6b5/attachment.html>


More information about the petsc-dev mailing list