In theory it is all there. Just do MatMatMult() with your dense matrices Barry On Dec 29, 2008, at 1:55 PM, Yujie wrote: > Hi, > > I am wondering whether it have been realized in PETSc to do > MPIDense*MPIDense? thanks a lot. Happy holiday:). > > Regards, > > Yujie >