4
$\begingroup$

Inverse of a sparse matrix could be dense, but what if I'm only interested in the main diagonal of the result? Is there a method that is more efficient than computing the full inverse?

3 Answers 3

2

You may be interested in work of Jok M. Tang and Yousef Saad. Their recent paper "Domain-decomposition-type methods for computing the diagonal of a matrix inverse" has the following abstract:

alt text

  • 0
    I was hoping for something that I could quickly implement using standard linear algebra routines, I guess it's not possible2011-01-23
1

This(explanation of how to find diagonal of an inverse of a sparse matrix) should help.

1

We recently implemented in our solver PARDISO such as option where you can these results. The software is available at www.pardiso-project.org.

It is not yet documented in the manual b

Regards, Olaf