how to inverse a sparse matrix in Petsc?

Yujie recrusader at
Mon Feb 4 12:20:28 CST 2008

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.


On 2/4/08, Barry Smith <bsmith at> 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: <>

More information about the petsc-users mailing list