how to inverse a sparse matrix in Petsc?

[Waad Subber]

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>