Concept

Contour integral — where it appears

An integral taken around a closed curve in the complex plane and evaluated by quadrature at points on it. For an analytic integrand the trapezoidal rule converges geometrically rather than algebraically, which is what makes contour methods for eigenvalues practical.

Named by 8 essays across 3 fields — each of them below, with the objects they name alongside it.

-2-101234-5-3-1135real partimaginary part4 insidea countable spectruminside the contour4drawn12existing∞worst branch residual1.6·10⁻¹⁵there is no last eigenvalueso the question has to change

A problem with infinitely many eigenvalues

Let the matrix depend on λ through something that is not a polynomial and three things stop being true at once. There is no linearisation, there is no characteristic polynomial, and "compute the spectrum" is not a request that can be granted — the only finite question is how many eigenvalues are inside this circle.

polynomial · Nonlinear eigenvalue
11.31.61.92.210⁻¹⁶10⁻¹³10⁻¹⁰10⁻⁷10⁻⁴10⁻¹10²log₁₀ quadrature pointsdistance from the true counthalf an eigenvaluean integer, eventuallytrue count2at 4 points2at 128 points2finest error1.4·10⁻¹³the integral is an integerand a rounding hides how far it was

The last digit is the cheapest

Every cost curve measured here has the same shape: the first digits are cheap and the last ones are not. One method inverts it. Doubling the work buys twice as many digits as the previous doubling did, so the price of a digit halves every time it is paid.

cost · Nonlinear cost
11.31.61.92.210⁻¹⁶10⁻¹³10⁻¹⁰10⁻⁷10⁻⁴10⁻¹10²log₁₀ quadrature pointsdistance from the true counthalf an eigenvaluean integer, eventuallytrue count2at 4 points2at 128 points2finest error1.4·10⁻¹³the integral is an integerand a rounding hides how far it was

Counting what is inside a circle

A trace of a matrix nobody wants to form, integrated around a contour, gives an integer — how many eigenvalues are inside. It converges exponentially, it is estimated with random probes, and the probe block is a ceiling that the answer does not mention.

randomised · Trace estimation
123456024681012moments Krank of the block Hankel12 eigenvalues insidethe ceiling K·ℓthe ranka ceiling with a knob on itprobes2eigenvalues inside12rank at K = 12rank at K = 612solves, at every K512one ceiling per probeand K of them per moment

A ceiling with a knob on it

A contour method returns at most as many eigenvalues as its probe block has columns, and the object that comes back does not distinguish that from having found everything. One line of the derivation multiplies the ceiling by a number the caller chooses, and it costs no extra solves at all.

polynomial · Nonlinear eigenvalue
123456110²10⁴10⁶moments Kκ of the block Hankelintegrated in zin (z − c)/ρone division per pointradius6unscaled, K = 62.3·10⁶scaled, K = 61750the factor between1291the ceiling rises by Kand so did the conditioning

The conditioning that rises with the ceiling

Higher moments multiply a contour method's ceiling by K and grade its block Hankel over ρ to the 2K, so the two knobs are the same knob. One division per quadrature point separates them, and the measurement of what it is worth grows from twenty to twenty thousand.

polynomial · Nonlinear eigenvalue
-6-4-20246-6-4-20246real partimaginary partan infinite spectrum, finitely askedinside the contour12probes4moments used3values returned24worst against the closed form6.2·10⁻¹⁵infinitely many eigenvaluesand a question with an answer

Where a contour's budget should go

A contour method's ceiling is the number of probes times the number of moments, and the moments are free in solves while the probes are not. Four ways of reaching one ceiling come out four orders apart, the ordering is not monotone, and what separates the best two is not the usual draw but the unlucky one.

polynomial · Nonlinear eigenvalue
123456700.10.20.30.4circle radius Rdigits per quadrature point27 points41 points78 pointseight eigenvalues, 4.19 to 4.56line: the annulus · dots: measured · large: ten digits at the geometric meanon a shared ray the balance is a cancellation

The circle between two eigenvalues

A contour count's error is not approximately governed by the nearest eigenvalue; it is exactly one closed-form term per eigenvalue, and summing those terms reproduces the quadrature to a millionth at 1,720 radius and point pairs. Three things follow. The rate is the ratio of the two moduli the circle sits between, not a distance, so two circles 0.2 from their nearest eigenvalue converge three times apart. Ten digits cost about 21 points divided by log₁₀ of that ratio, 71 points with twelve eigenvalues inside and 3,476 with six. And the best circle is not halfway: at the geometric mean of two eigenvalues on one ray their two terms are equal and opposite, and ten digits cost 27 points where the midpoint needs 84.

cost · Nonlinear cost
1020501002005000.822 to 1.58741220.822 to 2.68678210.822 to 3.644341321.587 to 2.686381.587 to 3.644202.686 to 3.64426points for ten digitsbest circle about the originbest circle about any real centrethe four together: 341 about the origin, 32 about 2.3the centre is worth most where the rays differ

A centre that chooses which pair cancels

A contour count's error is one closed-form term per eigenvalue, and two real eigenvalues on one ray cancel at their geometric mean — about the origin. Moved along the real axis, the centre was predicted to buy a further factor of two by putting real eigenvalues on two rays, the far one cancelling with a sign flip. The flip is parity: an eigenvalue behind the centre changes sign only at odd point counts. The factor is not two. Counting the delay problem's first two and three real eigenvalues, the best centre costs 22 and 21 points against the origin's 41 and 78; counting all four, 32 against 341, because about the origin the circle has to pass between a real eigenvalue and a complex pair and about 2.3 it does not. The centre does not add a cancellation. It chooses which pair cancels and pushes everything else away.

cost · Nonlinear cost

Named alongside it

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

Nonlinear eigenvalue problemQuadratureContour eigensolverExact ground truthRankArithmetic costComplex arithmeticConvergence rateMomentNumerical rankArgument principleCondition number

All concepts