patch for seqdense

Simon Burton simon at
Tue Sep 6 19:45:17 CDT 2005

On Fri, 2 Sep 2005 15:41:50 -0500 (CDT)
Hong Zhang <hzhang at> wrote:

> Simon,
> I tested MatMatMult* using a new file
> petsc-dev/src/mat/examples/tests/ex104.c.
> Why you did not implement
> MatMatMultSymbolic_SeqDense_SeqDense()
> and
> MatMatMultTransposeSymbolic_SeqDense_SeqDense()?
> Did you use inplace mat product in your python test?

We ended up wrapping BLASgemm_ and using it directly.
This call provides a kind of MatMatMultAdd (except the result
gets stored in the add matrix) which we need. But obviously
this will not work for the other matrix types.



Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Ph. 61 02 6249 6940 

More information about the petsc-dev mailing list