Theorem 6.1, proof — when
ProvedTraceEstimation.Rayleigh.quadratic_form_rangeeigenvaluesp2o-batch-p200ap2o-gran-per-chapterp2o-plan-paperp2o-v1trace-estimation
Let be a nonzero symmetric positive semi-definite matrix with largest eigenvalue , and let be the ratio between its largest and smallest nonzero eigenvalue. For every vector with ,
Applied to each sample vector of a normalized Rayleigh-quotient estimator, this says every summand lies in the fixed interval , which is what makes Hoeffding's inequality applicable.
Formalization Note The statement is deterministic (a single vector ) and records all four links of the displayed chain as a conjunction. guarantees , so the division by is genuine.
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), second display: for a symmetric positive
semi-definite `A ≠ 0` with largest eigenvalue `λ_n`, and every `z ∈ ℝⁿ` with `zᵀz = n`,
`0 ≤ zᵀAz ≤ λ_n zᵀz = n λ_n ≤ (n / rank(A)) · trace(A) · κ_f(A)`.
This is the deterministic range of one sample of a normalized Rayleigh-quotient estimator. -/
theorem quadratic_form_range {n : ℕ} (A : Matrix (Fin n) (Fin n) ℝ) (hA : A.PosSemidef)
(hA0 : A ≠ 0) (z : Fin n → ℝ) (hz : z ⬝ᵥ z = (n : ℝ)) :
0 ≤ z ⬝ᵥ (A *ᵥ z) ∧
z ⬝ᵥ (A *ᵥ z) ≤ lambdaMax hA.isHermitian * (z ⬝ᵥ z) ∧
lambdaMax hA.isHermitian * (z ⬝ᵥ z) = (n : ℝ) * lambdaMax hA.isHermitian ∧
(n : ℝ) * lambdaMax hA.isHermitian ≤
(n : ℝ) / (A.rank : ℝ) * 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 2
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.