> We are stoing the diagonal block and offdiagonal > block > separately. However both blocks are on the same > processor. i.e each > processor stores m*N values - in 2 submatrices m*n, > m*(N-n). To > understand this better - check manpage for > MatCreateMPIAIJ(). Thanks. But this is completely different from what I read from the PETSC mannual. http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/Mat/MatCreateMPIAIJ.html Here is says: "The DIAGONAL portion of the local submatrix of a processor can be defined as the submatrix which is obtained by extraction the part corresponding to the rows r1-r2 and columns r1-r2 of the global matrix, where r1 is the first row that belongs to the processor, and r2 is the last row belonging to the this processor. This is a square mxm matrix. The remaining portion of the local submatrix (mxN) constitute the OFF-DIAGONAL portion."
So the two matrices are mxm and mx(N-m) instead of what you said: mxn and mx(N-n) However, the code seems to act like what you described. They are equivalent for square matrices but not the same for non-square matrices like the one I showed. Could you please clarify whether the manual is accurate or not? Thank you very much. Shi ____________________________________________________________________________________ Sucker-punch spam with award-winning protection. Try the free Yahoo! Mail Beta. http://advision.webevents.yahoo.com/mailbeta/features_spam.html
