what is the difference between sequantial and parallel AIJ matrix? Assuming there is a matrix A, if I partitaion this matrix into A1, A2, Ai... An. A is a parallel AIJ matrix at the whole view, Ai is a sequential AIJ matrix? I want to operate Ai at each node. In addition, whether is it possible to get general inverse using MatMatSolve() if the matrix is not square? Thanks a lot.
Regards, Yujie On 2/4/08, Barry Smith <bsmith at mcs.anl.gov> wrote: > > > For sequential AIJ matrices you can fill the B matrix with the > identity and then use > MatMatSolve(). > > Note since the inverse of a sparse matrix is dense the B matrix is > a SeqDense matrix. > > Barry > > On Feb 4, 2008, at 12:37 AM, Yujie wrote: > > > Hi, > > Now, I want to inverse a sparse matrix. I have browsed the manual, > > however, I can't find some information. could you give me some advice? > > > > thanks a lot. > > > > Regards, > > Yujie > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080204/4dd1ff62/attachment.htm>
