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

[Barry Smith]

> On Sep 30, 2015, at 1:25 PM, Richard Mills <richardtmills at gmail.com> wrote:
> 
> Ha!  Thanks, Barry -- your way with words is always appreciated. =)
> 
> Matt, the inheritance model *is* kinda sucky and broken but at least you guys didn't use C++, so I'll take it.  I will take a look at the different subclasses (including the one that I wrote =) ) to refresh my memory about how to do this correctly.
> 
> I need to get an initial implementation in place first, but I'll go ahead and ask about how I'll eventually need to set things up in BuildSystem for this.  Perhaps I should follow a similar model to the MKL Pardiso support, i.e., make the user explicitly ask for this?  Something like '--with-mkl_spmv2-dir=LOCATION_OF_INTEL_MKL'?
> 
> What is the workflow for making changes to BuildSystem these days?

  BuildSystem is just in the PETSc repository so you just make changes there with a branch, test etc in the same way. 


>  Do we have this documented somewhere?  I ask because I'll need to add something for configuring the MKL inspector-executor stuff, and also because, in a separate issue, I need to add support for the memkind library (https://github.com/memkind/memkind).  (Yes, Jed, I know you don't like it but it's what we've got for now for using KNL MCDRAM.  I hope someone comes up with an intelligent middleware approach to dealing with this at some point.)
> 
> --Richard
> 
> On Wed, Sep 30, 2015 at 8:06 AM, Jed Brown <jed at jedbrown.org> wrote:
> Barry Smith <bsmith at mcs.anl.gov> writes:
> 
> >> On 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.
> >
> >   Jed and Matt,
> >
> >    Remember Richard actually wrote the subclasses of AIJ for the Cray vector machine so you guys are teaching your grandmother to suck eggs.
> 
> Where is @shitbsmithsays when you need it?
> 
> >   Barry