how to inverse a sparse matrix in Petsc?

Yujie recrusader at gmail.com
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.

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>


More information about the petsc-users mailing list