mgw-convergence
At its defaults it draws minres on a saddle-point system under three schur-complement approximations, and under none. The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 4.109·10⁻¹⁵ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 11 steps and replacing it by a scaled AAᵀ costs 11; the unpreconditioned system takes 22. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
mgw-convergence is one function in lib/figures/blockprec.js —
preconditioning a saddle point — three eigenvalues, and a spectrum with no constraint in it. Everything below came out of it during this build, at
arguments taken from the essays rather than invented for this page. A figure here is the
figure a reader meets in an essay, and if the generator changes, this page changes with it.
At its defaults
Drawn even though every essay passes arguments — which on this site is every essay, at 100% of placements since the standard pass. A default nothing exercises is a trap for the next essay to call this with none, and this is the page where a default that has drifted from the figures around it becomes visible.
The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 4.109·10⁻¹⁵ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 11 steps and replacing it by a scaled AAᵀ costs 11; the unpreconditioned system takes 22. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
kappaH: 10000
The arguments are the ones A preconditioner that need not know the constraint passes. A value drawn at the generator's defaults instead would be a picture no essay asked for and no assertion has been run against.
The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 2.947·10⁻¹⁴ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 13 steps and replacing it by a scaled AAᵀ costs 13; the unpreconditioned system takes 25. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
kappaH: 1
The arguments are the ones A preconditioner that need not know the constraint passes. A value drawn at the generator's defaults instead would be a picture no essay asked for and no assertion has been run against.
The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 6.534·10⁻¹⁶ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 3 steps and replacing it by a scaled AAᵀ costs 3; the unpreconditioned system takes 13. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
m: 3
The arguments are the ones The zero that is not a missing entry passes. A value drawn at the generator's defaults instead would be a picture no essay asked for and no assertion has been run against.
The same system at 12 unknowns and 3 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 4.31·10⁻¹⁵ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 7 steps and replacing it by a scaled AAᵀ costs 7; the unpreconditioned system takes 18. The exact preconditioner is unaffordable — forming S costs 3 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
kappaH: 100
The arguments are the ones Three eigenvalues, and two are the golden ratio passes. A value drawn at the generator's defaults instead would be a picture no essay asked for and no assertion has been run against.
The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 4.109·10⁻¹⁵ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 11 steps and replacing it by a scaled AAᵀ costs 11; the unpreconditioned system takes 22. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
kappaH: 100000000
The arguments are the ones Three eigenvalues, and two are the golden ratio passes. A value drawn at the generator's defaults instead would be a picture no essay asked for and no assertion has been run against.
The same system at 12 unknowns and 5 constraints, solved four ways. With the exact Schur complement the preconditioned matrix has three distinct eigenvalues and the residual falls to 5.317·10⁻¹⁰ in three steps, after which nothing is left to remove. Replacing S by A diag(H)⁻¹Aᵀ costs 15 steps and replacing it by a scaled AAᵀ costs 17; the unpreconditioned system takes 42. The exact preconditioner is unaffordable — forming S costs 5 solves with H and a decomposition — so its value is as the statement the cheap ones are measured against, and the measurement needs no reference solution: the distance from {1 − φ, 1, φ} is a property of the approximation alone.
What it checked while drawing
Every figure above asserted its own claims on the way to being drawn, and a claim that failed
would have failed the build rather than drawn a wrong picture. Those assertions used to leave
no trace at all: a passing one returned true and the only evidence the figure had
checked anything was that nothing crashed. The list below is what they actually said, collected
by running this generator with an observer installed — not a description of
what it is believed to check.
10 distinct claims across 6 sets of arguments, grouped below by shape — because most of them are one sentence with a different number in it, and how many separate times that sentence was put to the test is the informative part.
a constraint no larger than the problem
a positive definite preconditioner, since MINRES needs one
a problem with something left to minimise
a Schur approximation this file implements
and the same system unpreconditioned needs more
fewer constraints than unknowns, so something is left to minimise
Jacobi needs a symmetric matrix
LU is for square matrices
matmul shapes agree
three steps with the exact Schur complement take the residual below 10⁻⁸
Against the rule
It draws a decomposition and prints its residual. It calls
minres, blockDiagonalRun,
and every figure above carries the badge — which residualcheck verifies by looking
for it in the emitted SVG rather than by finding the call that builds one. A badge that is
constructed and then left out of the body is the failure that check exists for.
Across the library: the rule bites on 165
of 306 generators —
150 print a residual and
15 are exempt with a published reason;
141 factorise nothing.
Read from lib/residual-rule.js, which is the same body the gate enforces from,
and the gate's last check fails the build if this page and it disagree about any generator.
Where it is called
Changing this generator changes every figure on this list. That is what makes the list worth publishing rather than keeping in a check script.
A preconditioner that need not know the constraint
Keep the constraint block exactly and replace the objective block by anything positive definite on the null space. The preconditioned matrix then has 2m eigenvalues at exactly one, and its remaining n − m are the generalised eigenvalues of a pencil in which the constraint does not appear. Sweep its condition number over six decades and they do not move in six digits.
The matrix a constraint makesThe zero that is not a missing entry
A constrained minimisation produces a matrix with a zero block, and the zero is a theorem rather than a sparsity pattern. No pivot order makes it positive definite, no precision changes that, and Cholesky does not fail somewhere on it — it fails at the first constraint row, on a number the problem already contained.
The matrix a constraint makesThree eigenvalues, and two are the golden ratio
Precondition a saddle-point system by the block diagonal of its own two definite pieces and the preconditioned matrix has exactly three distinct eigenvalues — 1, and the two roots of λ² − λ − 1. A minimal polynomial of degree three means three steps, at every conditioning, and the preconditioner nobody can afford turns out to be the statement the affordable ones are measured against.