Theorem 6.1, proof —
ProvedTraceEstimation.Rayleigh.rank_mul_lambdaMax_lecondition-numbereigenvaluesp2o-batch-p200ap2o-gran-per-chapterp2o-plan-paperp2o-v1trace-estimation
Let be a nonzero symmetric positive semi-definite matrix, let be its largest eigenvalue, and let be the ratio between its largest and smallest nonzero eigenvalue. Then
This linear-algebra inequality converts the largest eigenvalue, which bounds a single sample , into a multiple of the trace, so that the range of each sample is measured relative to the quantity being estimated.
Formalization Note The page derives the inequality from an indexing with and writes ; with that , , whereas needs to be the smallest nonzero eigenvalue (an off-by-one slip). Only the inequality is stated, with as defined in kappaF. The hypothesis excludes the placeholder value of .
Preamble
import Mathlib import Definitions.Def_TraceEstimation_Rayleigh_kappaF
Formal statement
namespace TraceEstimation.Rayleigh
open Matrix
/-- Avron–Toledo, proof of Theorem 6.1 (p. 8:10), first display: for a symmetric positive
semi-definite `A ≠ 0` with largest eigenvalue `λ_n` and `κ_f(A)` the ratio of its largest to its
smallest nonzero eigenvalue,
`trace(A) · κ_f(A) ≥ rank(A) · λ_n`.
(The page indexes the eigenvalues as `0 = λ_1 = ⋯ = λ_k ≤ ⋯ ≤ λ_n` with `k = n − rank(A) + 1`,
which makes `λ_k = 0`; the intended `λ_k` is the smallest nonzero eigenvalue. Only the
inequality is stated.) -/
theorem rank_mul_lambdaMax_le {n : ℕ} (A : Matrix (Fin n) (Fin n) ℝ) (hA : A.PosSemidef)
(hA0 : A ≠ 0) :
(A.rank : ℝ) * lambdaMax hA.isHermitian ≤ A.trace * kappaF hA.isHermitian := by sorry
end TraceEstimation.Rayleigh
Source
Avron and Toledo, Randomized algorithms for estimating the trace of an implicit symmetric positive semi-definite matrix, J. ACM 58(2), Article 8 (2011), p. 8:10, Section 6, proof of Theorem 6.1, display 1
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.