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
