Numerical Linear Algebra
where the algebra stops being the arithmetic
Concept

Rotation — where it appears

An orthogonal matrix with determinant +1, which turns space without reflecting it. The rotations are one of the two components of the orthogonal matrices, so a nearest-orthogonal-matrix computation can return a reflection instead unless the determinant is constrained.

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

how far is it from A to an orthogonal matrix?smaller is nearer · the polar factor minimises this in every unitarily invariant normpolar factor U1.8554QR, signs fixed2.1265QR as returned3.8226200 drawn at randomκ = 10polar factor1.9QR, signs fixed2.1QR as returned3.8best of 200 random2.7‖A − QR‖ is the same either wayand ‖A − Q‖ is not

The nearest orthogonal matrix

Every field that has to clean up a drifted rotation reaches for QR, and QR does not answer the question. The nearest orthogonal matrix is the orthogonal factor of the polar decomposition — nearer by about a tenth, and, more to the point, the same matrix whatever order the columns were written in. QR's answer changes completely.

orthogonality · Polar decomposition 110¹10²01020304050noise ÷ thicknesstrials mirrored, %a coin: 50%t = 10⁻²t = 10⁻³t = 10⁻⁴per cent mirroredσ/t = 3, mean of three7.7σ/t = 10, mean of three34σ/t = 100, mean of three4720 points, 400 trials a stop, one seed per thicknessthe ratio decides, not the thinness

A rotation that comes back mirrored

Align twenty noisy points and the nearest orthogonal matrix to the answer is a reflection in 7.7 per cent of trials at noise three times the set's thickness and a third of them at ten — at thicknesses of 10⁻², 10⁻³ and 10⁻⁴ alike. The determinant fix is never a small correction. It moves the answer by exactly 2, it costs exactly 4σ₃ of residual, and it leaves the rotation's error at half the noise however thin the set becomes.

orthogonality · Polar decomposition

Named alongside it

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

OrthogonalityPolar decompositionCondition numberDeterminantFrobenius normGram–SchmidtHouseholderInvarianceQR factorisationReflectionSensitivitySingular value decomposition

All concepts