Random C++ Part 2: Sparse partial inverses in Eigen
Read OriginalThis technical article details the implementation of an algorithm to compute the non-zero entries of the inverse of a sparse symmetric positive definite matrix, specifically the elements corresponding to the original matrix's non-zero pattern. It focuses on using the Takahashi recursions with a Cholesky factorization in C++ via the Eigen library, discussing motivation from statistical applications and the challenges of working with sparse matrix internals.
Comments
No comments yet
Be the first to share your thoughts!
Browser Extension
Get instant access to AllDevBlogs from your browser
Top of the Week
1
Quoting Thariq Shihipar
Simon Willison
•
2 votes
2
Using Browser Apis In React Practical Guide
Jivbcoop
•
2 votes
3
Top picks — 2026 January
Paweł Grzybek
•
1 votes
4
In Praise of –dry-run
Henrik Warne
•
1 votes
5
Deep Learning is Powerful Because It Makes Hard Things Easy - Reflections 10 Years On
Ferenc Huszár
•
1 votes
6
Vibe coding your first iOS app
William Denniss
•
1 votes
7
AGI, ASI, A*I – Do we have all we need to get there?
John D. Cook
•
1 votes
8
Dew Drop – January 15, 2026 (#4583)
Alvin Ashcraft
•
1 votes