Re: [petsc-dev] How to check that MatMatMult is available

2019-09-19 Thread Pierre Jolivet via petsc-dev
> On 20 Sep 2019, at 7:36 AM, Jed Brown > wrote: > > Pierre Jolivet via petsc-dev > writes: > >> Hello, >> Given a Mat A, I’d like to know if there is an implementation available for >> doing C=A*B >> I was previously using MatHasOperati

Re: [petsc-dev] How to check that MatMatMult is available

2019-09-19 Thread Jed Brown via petsc-dev
Pierre Jolivet via petsc-dev writes: > Hello, > Given a Mat A, I’d like to know if there is an implementation available for > doing C=A*B > I was previously using MatHasOperation(A, MATOP_MATMAT_MULT, &hasMatMatMult) > but the result is not correct in at least two cases: Do you want MATOP_MAT_

[petsc-dev] How to check that MatMatMult is available

2019-09-19 Thread Pierre Jolivet via petsc-dev
Hello, Given a Mat A, I’d like to know if there is an implementation available for doing C=A*B I was previously using MatHasOperation(A, MATOP_MATMAT_MULT, &hasMatMatMult) but the result is not correct in at least two cases: 1) A is a MATTRANSPOSE and the underlying Mat B=A^T has a MatTransposeMa