[petsc-dev] Questions about MatMatTransposeMult/MatMatMultTranspose

[Dmitry Karpeev]

Is there any reason that there isn't dual dispatch support for these
routines like there is for MatMatMult?
If there is no particular reason, I would like to add it, largely for the
same reasons I'd like symmetric
dual dispatch for MatMatMult: the ability to add new implementations
without screwing around with all of the
constructors.

Similar considerations would apply to the MatMat*Symbolic routines, as well
as MatPtAP.
I still maintain that the right way to handle these is triple dispatch that
depends on the type of the output
Mat as well, but that would only matter in a handful of corner cases.

And a minor question: wouldn't it be better to name these two functions
MatTransposeMatMult and MatMatTransposeMult
respectively, so that the "Transpose" operation sits next to the matrix it
operates on?

Dmitry.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111101/239b1fc3/attachment.html>