patch for seqdense

[Simon Burton]

On Fri, 2 Sep 2005 15:41:50 -0500 (CDT)
Hong Zhang <hzhang at mcs.anl.gov> 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.

bye,

Simon.

-- 
Simon Burton, B.Sc.
Licensed PO Box 8066
ANU Canberra 2601
Australia
Ph. 61 02 6249 6940
http://arrowtheory.com