how to inverse a sparse matrix in Petsc?

Dave May dave.mayhem23 at gmail.com
Mon Feb 4 08:04:17 CST 2008


Hi,
    Does anyone know how much faster (approximately) using MatMatSolve is
compared
to using PCComputeExplicitOperator(), when the PC in the latter function is
defined to be LU?

Cheers,
    Dave.


On Feb 5, 2008 12:10 AM, 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
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080205/44383e32/attachment.htm>


More information about the petsc-users mailing list