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? I added implementation of these routines. The tests run well. The update code and ex104.c have been pushed to petsc-dev. Thanks for the contribution, Hong On Mon, 29 Aug 2005, Simon Burton wrote:

On Sun, 28 Aug 2005 09:29:26 -0500 (CDT) Hong Zhang <[email protected]> wrote:

...

...

I sent through a patch to implement MatMatMult* for seqdense matrices.

Where is this patch? through bk or an attached file? I don't see MatMatMult_seqdense from petsc-dev.

Also, do you have testing code for it? I can test it once I see the MatMatMult_seqdense.

We tested it using python. I see also there is a test in src/mat/examples/tests/ex94.c .

Simon.

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