On Sep 30, 2015, at 9:48 AM, Jed Brown <[email protected]> wrote:

Matthew Knepley <[email protected]> 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. Barry