Hi There was a discussion between Tim Stitt and petsc developers about matrix inversion, and it was really helpful. That was in last Nov. You can check the emails archive
http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2007/11/threads.html Waad Yujie <recrusader at gmail.com> wrote: 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 > --------------------------------- Looking for last minute shopping deals? Find them fast with Yahoo! Search. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080204/7d617a9b/attachment.htm>
