How statement and proof provenance work
The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.
- Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
- AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
- AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.
These labels describe origin, not correctness: citations and verification chips remain separate evidence.
Direct Matrix Factorisations: LU, Cholesky and QR: Examples and Counterexamples
1 · Prerequisites
- Algebraic Closure, Embeddings, and Separability
- Algebraic Extensions, Extension Degree, and Finite Fields
- Binary Operations, Monoids, Groups and Subgroups
- Compactness in Metric Spaces
- Completeness, Completion, and Uniform Continuity
- Composition Series, the Jordan–Hölder Theorem and Solvable Groups
- Congruences, the Integers Modulo n and the Chinese Remainder Theorem
- Construction of the Natural Numbers
- Construction of the Real Numbers via Cauchy Sequences
- Construction of the Real Numbers via Dedekind Cuts
- Continuity, IVT, EVT, and Uniform Continuity
- Cosets, Index and Lagrange's Theorem
- Countability and Uncountability
- Cyclic Groups and Direct Products
- Determinants of Matrices over a Commutative Ring
- Diagonalisation and the Minimal Polynomial
- Direct Matrix Factorisations: LU, Cholesky and QR
- Divisibility, Euclidean Domains, Principal Ideal Domains and Unique Factorisation
- Divisibility, Greatest Common Divisors and Bézout's Identity
- Dual Spaces, Bilinear and Quadratic Forms, and Sylvester's Law of Inertia
- Eigenvalues, Eigenvectors and the Characteristic Polynomial
- Filters and Ultrafilters
- Finite Counting, Factorials and Binomial Coefficients
- Foundations of the Real Numbers for Analysis
- Group Actions, Orbits, Stabilisers and Cayley's Theorem
- Group Homomorphisms and the Isomorphism Theorems
- Ideals, Quotient Rings and the Isomorphism Theorems for Rings
- Inner Product Spaces, Gram-Schmidt, Projections and Adjoints
- Limits of Real Functions
- Linear Independence, Bases and Dimension
- Linear Transformations, Rank-Nullity and Quotient Spaces
- Matrices, the Matrix of a Linear Map, and Change of Basis
- Matrix Norms, Condition Numbers and Numerical Stability
- Metric Spaces
- Monotone Sequences, Bolzano-Weierstrass, and Cauchy Completeness
- Normal Subgroups and Quotient Groups
- Order, Zorn's Lemma, and the Axiom of Choice
- Polynomial Rings, the Division Algorithm and Roots
- Primes, Euclid's Lemma and the Fundamental Theorem of Arithmetic
- Relations, Functions, and Quotients
- Rings, Subrings, Integral Domains and Fields
- Rⁿ as a Normed Space; Vector-Valued Functions
- Roots, Rational Powers, and Classical Inequalities
- Sequences and Limits
- Simple Field Extensions and the Construction of the Complex Numbers
- Splitting Fields
- Suprema and Infima
- Sylow's Theorems, p-Groups and Nilpotent Groups
- Symmetric Groups, Cycle Decomposition and the Sign Homomorphism
- The Determinant of a Linear Operator, Cofactors and Cramer's Rule
- The Fundamental Theorem of Algebra
- The Fundamental Theorem of Finite Abelian Groups
- The Galois Correspondence
- The Spectral Theorem, Positive Operators and Singular Value Decomposition
- The ZFC Axioms and the Basic Set Constructions
- Topology of ℝ
- Triangularisation, Generalised Eigenspaces and Jordan Canonical Form
- Vector Spaces, Linear Subspaces, Span and Direct Sums
2 · Summary
These examples keep the A-page hypotheses honest. They show an invertible matrix that still fails unpivoted LU, a concrete PLU and LDU computation, a block LU solve through a Schur complement, a worked Cholesky solve, and the exact ways indefinite or merely semidefinite matrices fall outside positive-diagonal Cholesky.
The QR examples compute a short Householder factorization, use Givens rotations to preserve sparsity while zeroing chosen entries, compare reduced QR with the normal equations on a badly scaled least-squares problem, and exhibit fill-in created by Gaussian elimination.
3 · Logical flowchart
4 · Definitions, theorems and proofs
None yet.
5 · Examples, counterexamples and false statements
An invertible matrix can fail unpivoted LU at the first pivot
Statement refuted
Every invertible square matrix has an unpivoted unit-lower LU factorisation.
Facts & Assumptions
Given: The matrix
A square matrix has an unpivoted unit-lower LU factorisation with nonzero pivots exactly when every leading principal minor is nonzero (A square matrix has an unpivoted unit-lower LU factorisation exactly when all leading principal minors are nonzero).
Counterexample
The determinant of is , so is invertible. Its first leading principal minor is the determinant .
By [L1], the vanishing of forbids an unpivoted unit-lower LU factorisation with nonzero pivots. Therefore invertibility alone does not suffice.
Step 2.1 refutes the statement.
A full PLU factorisation can be computed explicitly by hand
Example
For
partial pivoting gives
and these satisfy .
Facts & Assumptions
Given: The displayed matrix and the candidate factors .
Every square matrix admits a PLU factorisation, and partial pivoting records the row swaps in a permutation matrix (Every real or complex square matrix admits a PLU factorisation, and the first failed pivot marks the singular boundary, Permutation matrices, partial pivoting, and the pivot-growth factor).
Verification
The largest entry in modulus in the first column is in row , so the first pivot swap sends that row to the top. The first elimination multipliers are for the new second row and for the new third row, producing the intermediate matrix The second pivot is already the entry in row , so no further swap is needed, and eliminating the entry with multiplier gives the displayed . Recording the two nonzero multipliers in the permuted row order gives the displayed .
Direct multiplication gives Hence the displayed matrices are a correct PLU factorisation.
Steps 1.1-2.1 verify the example.
An LDU factorisation isolates the pivot scalars uniquely
Example
The matrix
has the LDU factorisation
and the diagonal pivots are uniquely determined.
Facts & Assumptions
Given: The displayed matrix and candidate LDU factors.
An LDU factorisation has the displayed unit-lower, diagonal, and unit-upper shape (An LDU factorisation has unit lower-triangular L, diagonal D, and unit upper-triangular U).
LDU factorisations with nonzero diagonal pivots are unique (Normalised LU and LDU factorisations with nonzero pivots are unique).
Verification
Multiplying the right two factors gives and then left multiplication by yields .
The diagonal factor has nonzero entries and , so [L2] applies. Any other LDU factorisation of with nonzero diagonal pivots must therefore have the same diagonal factor and hence the same pivot scalars.
Steps 1.1-2.1 verify the example.
A block LU factorisation turns a linear solve into a Schur-complement solve
Example
For
with the block split, the Schur complement is
and
For the right-hand side , the solve reduces to the Schur-complement solve and yields .
Facts & Assumptions
Given: The displayed matrix , its block split, the right-hand side , and the candidate factorisation.
An invertible leading block yields the block LU factorisation through its Schur complement (An invertible leading block yields block LU through its Schur complement).
Triangular systems are solved by forward and backward substitution (Forward and backward substitution are correct, unique, and quadratic in scalar operations).
Verification
The leading block is , so , which is the displayed Schur complement. The block-LU formula of [L1] gives the displayed factorisation.
Solve : , , . Then solve from the bottom: , , and . This gives , , and .
Steps 1.1-2.1 verify both the factorisation and the Schur-complement solve.
A Cholesky factorisation solves a small positive-definite system efficiently
Example
The symmetric positive-definite matrix has Cholesky factor and for the solve gives .
Facts & Assumptions
Given: The displayed matrix , the candidate factor , and .
Hermitian positive-definite matrices admit a unique Cholesky factorisation with positive diagonal (A matrix admits a Cholesky factorisation with positive diagonal exactly when it is Hermitian positive definite, and that factor is unique).
Cholesky solves by two triangular substitutions (Cholesky solves Hermitian positive-definite systems and has about half the factorisation cost of LU).
Verification
Direct multiplication gives so the displayed matrix is a Cholesky factor.
Solve : and , so and . Then solve : and , hence and .
Steps 1.1-2.1 verify the example.
Indefinite and semidefinite matrices can both fail positive-diagonal Cholesky
Statement refuted
Indefinite matrices, and even positive-semidefinite singular matrices, admit Cholesky factorisations with positive diagonal.
Facts & Assumptions
Given: The matrices
A matrix admits a Cholesky factorisation with positive diagonal exactly when it is Hermitian positive definite (A matrix admits a Cholesky factorisation with positive diagonal exactly when it is Hermitian positive definite, and that factor is unique).
Counterexample
The matrix is Hermitian, but with one has , so is not positive definite. Therefore [L1] forbids a positive-diagonal Cholesky factorisation of .
The matrix is positive semidefinite but singular. If with positive diagonal, then every diagonal entry of is nonzero, so would be invertible and would be invertible as well, a contradiction. Hence also has no such Cholesky factorisation.
Steps 1.1-1.2 refute the statement in both the indefinite and the merely semidefinite cases.
A short dense matrix admits a worked Householder QR factorisation
Example
For
the unit vector gives the Householder reflector
Hence
Facts & Assumptions
Given: The displayed matrix , vector , and reflector .
Successive Householder reflectors produce QR factorisations (Successive Householder or Givens transformations produce full and reduced QR factorisations with the standard dense operation counts).
Verification
The vector has norm , and direct multiplication of gives the displayed matrix . Since , it is an orthogonal reflector.
Multiplying by gives Thus with and upper triangular, exactly as predicted by [L1].
Steps 1.1-2.1 verify the example.
Givens QR can zero selected entries of a sparse matrix one at a time
Example
For first apply a Givens rotation in rows with , , then a second Givens rotation in rows with , . The product zeros one selected entry at a time and yields
Facts & Assumptions
Given: The displayed matrix and the two named Givens rotations.
Givens transformations are unitary and can annihilate a chosen second coordinate (Householder reflectors and Givens transformations are unitary and can annihilate prescribed entries).
Successive Givens transformations produce QR factorizations (Successive Householder or Givens transformations produce full and reduced QR factorisations with the standard dense operation counts).
Verification
The first rotation sends the first column to , so it zeros the entry while preserving the zero in position . Applied to the second column, it gives .
The second rotation acts only on rows , so it keeps the first column fixed and sends to . Hence the displayed upper-triangular matrix.
Since both rotations are unitary, and . This verifies the Givens QR factorisation promised by [L2].
Steps 1.1-3.1 verify the example.
Reduced QR avoids the condition-number squaring seen in the normal equations
Example
Let Then already has reduced QR factorisation with and . The least-squares solution is , while and .
Facts & Assumptions
Given: The displayed matrix , vector , and reduced QR factorisation.
Reduced QR solves full-column-rank least squares through and avoids the condition-number square of the normal equations (Reduced QR over the reals solves full-column-rank least squares without squaring the condition number).
Verification
The columns of are already orthogonal, so and form a reduced QR factorisation. Also .
Solving gives and , hence . The factor has spectral condition number , whereas has spectral condition number . This is exactly the contrast described in [L1].
Steps 1.1-2.1 verify the example.
Sparse Gaussian elimination can create fill-in in the factors
Example
The sparse matrix has a zero in position , but after the first elimination step that entry becomes . Thus Gaussian elimination can create nonzeros that were absent in the original matrix.
Facts & Assumptions
Given: The displayed sparse matrix .
Pivoting language records the elimination process entry by entry (Permutation matrices, partial pivoting, and the pivot-growth factor).
Verification
Use the first pivot . Eliminating the entries below it subtracts row from rows and , giving The entry in position was before elimination and is now .
The new nonzero in step 1.1 is fill-in: it appears after one elimination step even though the corresponding original entry was zero.
Steps 1.1-2.1 verify the example.
Sources
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Example 2.6.1
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Section 2.6
- David Bindel, CS 4220: Numerical Analysis, Blocked LU and Cholesky
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Example 2.9.3
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Section 2.9
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Example 3.4.1
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Givens rotations section
- Tobin A. Driscoll and Richard J. Braun, Fundamentals of Numerical Computation, Section 3.3.3 and Example 3.3.2