Deriving the Sparse Cholesky Elimination Tree for Matrix Factorization
By
selimthegrim
Hand-rolled, kettle-boiled, baked to perfection. Worth every minute at the bakery.
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 pulledThis 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.
You might also wanna read
Research Shows All 23-Bit Still Life Patterns in Conway's Game of Life Are Glider Constructible
The article discusses research in Conway's Game of Life about which still life patterns can be constructed by colliding gliders. It explains
GitHub Repository Claims Solution to Navier-Stokes Millennium Prize Problem Using Limit Functional Theorem and LLM Collaboration
A GitHub repository claims to present a solution to the Navier-Stokes Clay Math Millennium Prize Problem using a theorem about extending lim
Exploring Strange Attractors and Chaos Theory with Three.js Visualizations
The article describes the author's personal experience discovering strange attractors while working with Three.js, leading to a deep fascina

The Busy Beaver Problem: Recent BB(5) Breakthrough and the Antihydra Challenge
The article explores the busy beaver problem in theoretical computer science, focusing on the recent breakthrough in determining BB(5) as 47
benbrubaker.com·7mo agoUnderstanding Wilson's Algorithm for Uniform Spanning Trees
The article explains Wilson's Algorithm, a method for generating uniform spanning trees using loop-erased random walks. It describes how the
An Introduction to the Lambda Calculus: Functions, Abstraction, and Computation
The λ-calculus is a simple yet powerful notation for functions and application, centered on applying functions to arguments and forming func
