I have further checked this function. In MatMatMult(Mat A,Mat B,MatReuse scall,PetscReal fill,Mat *C) I am wondering why the type of C is MATAIJ when the types of A and B are MATAIJ. although A and B is MATAIJ, C should be dense. If C uses MATAIJ type, it should take more memory, is it right?
thanks a lot. Regards, Yujie -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080501/b327f602/attachment.htm>
