## Inverse of Sparse Matrix

mlimberger
Registered Member
Posts
9
Karma
0
OS

### Inverse of Sparse Matrix

Tue Nov 13, 2012 12:03 pm
Hello.

I know, there are already two posts concerning this topic, but unfortunately they doesn't help on.
ALthough it is not recommended, I want to inverse a sparse matrix. In particular, I want to calculate
the covariance matrix from a given sparse normal equations matrix. For solving my least squares system
I use a cholesky decomposition and solve it for the unknowns. This works quite well.

Eigen::SparseMatrix A(n,n); //n=17495
Eigen::VectorXd b(n);
Ax = b;
Eigen::SimplicialCholesky<Eigen::SparseMatrix<double> > chol(A);
x = chol.solve(b);

But how to compute A^(-1)

Marco
ggael
Moderator
Posts
2147
Karma
15
OS

### Re: Inverse of Sparse Matrix

Tue Nov 13, 2012 8:17 pm
call solve with a sparse matrix representing the identity matrix.
mlimberger
Registered Member
Posts
9
Karma
0
OS

### Re: Inverse of Sparse Matrix

Thu Nov 15, 2012 10:02 am
Thank you! My approach was to use the Lower triangular matrix from the Cholesky LLT. But you are completely right, using the Identity matrix works fine!

Again, thanks!

## Who is online

Registered users: alake, ant, Baidu [Spider], bilbo, Bing [Bot], Charblaze, dequire, DestinTheGreat, einar, Exabot [Bot], fences, ggael, Google [Bot], Hei Ku, kainz.a, kdeuserk, La Ninje, lazyit, Leroux, Majestic-12 [Bot], prowlerxpla, rzeta, scummos, SecretCode, sfepa, sir_herrbatka, Sogatori, Sogou [Bot], trebor, Uri_Herrera, Yahoo [Bot]