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
2195
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: Baidu [Spider], Bing [Bot], Dimitrios, drnn1076, Exabot [Bot], Google [Bot], google01103, Hans, Horus, koriun, Majestic-12 [Bot], pbCyanide, pedrorodriguez, random_fan, scummos, SeaJey, Sogatori, Sogou [Bot], TheraHedwig, vascobasque, whatthefunk, Yahoo [Bot], yurchor, z-uo