Numerical Linear Algebra
where the algebra stops being the arithmetic
Concept

Circulant matrix — where it appears

2 essays name this object, across one field. What follows is each of them, and the objects they name alongside it.

Also named here as discrete fourier transform — the same set of essays touches all of them, so they are one junction rather than several.

024681012141610⁻¹1index kλthe transformthe eigensolverC = F* Λ F is a factorisationworst relative disagreement2.6·10⁻¹⁵‖Cx − b‖/‖b‖ from the transform solve4.7·10⁻¹⁶imaginary part of a real spectrum1.2·10⁻¹⁶n = 16, and the whole matrix is 16 numberseigenvectors known in advance

The matrix that is one row

A circulant of size 16 is sixteen numbers, has no zero entry anywhere, and hands over its entire spectrum in closed form — the discrete Fourier transform of its first column, exactly. An eigensolver spends a sweep of Jacobi rotations over 256 entries arriving at the same answer, and agrees to 1.2·10⁻¹⁵.

structure · circulant 10¹10²10²size ncondition numberlimit 81measureda limit, as a fraction of itselfreached at n = 1280.99still to go0.011κ at n = 8, as a fraction0.52every point is below the line and none of them is on itthe limit is not a value

A limit the matrix never reaches

Szegő's theorem gives a Toeplitz family's condition number in closed form — ((1+ρ)/(1−ρ))², which is 81 at ρ = 0.8. The 8×8 section reaches 52% of it, the 128×128 reaches 98.9%, and none of them ever arrives. A statement about a family is not a statement about the matrix in front of you.

structure · toeplitz

Named alongside it

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

Condition numberDiscrete fourier transformAsymptotic analysisEigenvaluesExact inverseFactorisationFast fourier transformMatrix structureResidualSymbolSzego theoremToeplitz matrix

All concepts