All Topics
All Topics
Technology
Technology
Design
Design
Programming
Programming
Science
Science
News
News
Gaming
Gaming
Entertainment
Entertainment
Business
Business
Finance
Finance
Sports
Sports
Health
Health
Food
Food
Travel
Travel
Art
Art
Music
Music
Books
Books
Education
Education
Politics
Politics
Personal
Personal
No algorithm. No AI slop. No ads. Just RSS. Pro-human. Indie writers. Real journalism. Open web. Chronological. Hand toasted.

Deriving the Sparse Cholesky Elimination Tree for Matrix Factorization

By

selimthegrim

21d ago· 6 min readenInsight

Summary

This article provides a technical derivation of the elimination tree for the right-looking sparse Cholesky algorithm (A = LL^T) for sparse matrices. The author explains how this tree serves as the foundation for most sparse factorization software, revealing two key insights: where fill-in nonzeros appear in matrix L that weren't in the original A, and the task dependency graph for the factorization process. Unlike traditional presentations that build from sparse triangular solves, the author derives the elimination tree directly from the Cholesky factorization itself.

Key quotes

· 3 pulled
This tree forms the foundation for most sparse factorization software, even when the underlying assumptions of Cholesky (symmetric and definite) do not apply.
Ultimately this tree tells us two things: 1. where nonzeros appear in the matrix L even if not present in the original A (i.e. 'fill-in') and 2. the task dependency graph of our resulting factorization.
I wanted to instead work directly from a Cholesky factorization, which is what I do below.
Snippet from the RSS feed
Here I derive the elimination tree for the (right-looking) sparse Cholesky algorithm for computing A = LL^T for lower triangular L and sparse matrices A. This tree forms the foundation for most sparse factorization software, even when the underlying assum

You might also wanna read