Numerical Linear Algebra
where the algebra stops being the arithmetic
Concept

Range finder — where it appears

The randomised procedure that produces an orthonormal basis for the column space of a matrix from a few products with a random one. It is what turns a matrix that can only be applied into a matrix that can be factorised.

Named by 2 essays across one field — each of them below, with the objects they name alongside it.

012345610⁻¹⁶10⁻¹³10⁻¹⁰10⁻⁷10⁻⁴10⁻¹round‖AΩ‖ ⁄ ‖A‖ at that roundredrawn each roundhalf redrawnkept: the sketch is now the zero matrixthe randomness is spentkept, round 12.6kept, round 29·10⁻¹⁵redrawn, round 23.1kept, final error0.28redrawn, final error0.11a random matrix used twiceis not random the second time

The sketch that is spent

Every other object a sequence carries has a shelf life. A random sketch has one use. Deflate what its first round found and apply it again, and it returns the zero matrix — 9.0·10⁻¹⁵ where the first round saw 2.62 — because the input has been made orthogonal to the very draw the guarantee is over.

randomised · sketching 567891010¹10²10³10⁴10⁵10⁶log₂ nnumbers the construction touchedentries, n²products with the operatoran operator, applied a few hundred timesproducts at n = 512256entries at n = 5122.6·10⁵per doubling48‖A − A_H‖ ⁄ ‖A‖4·10⁻⁷excess over the compression7.3no entry of the matrixwas ever read

Built from products alone

A 512-square hierarchical representation, at a relative error of 4·10⁻⁷, from 256 applications of an operator that is never assembled. The compression route reads 262,144 entries; this one reads none, and pays for it with a factor of seven against the representation the entries would have given.

randomised · sketching

Named alongside it

The objects these essays reach for when they reach for this one.

Randomised SVDAdaptivityAdmissibilityBlack box constructionDeflationHierarchical matrixLow-rank approximationMatrix-freeOff diagonal rankOrthogonal projectionOversamplingProbability bound

All concepts