Theorem 9 — rate of the accelerated random method (Eq. (62)) — goal theorem
ProvedRandomGradFree.Accelerated.accelerated_random_method_rateLet be a real inner product space of dimension . Let be differentiable with -Lipschitz gradient, , and strongly convex with parameter :
( is allowed and means is convex). Let be a minimizer of , and write and . Let , and set
Consider a run of the accelerated random method (Eq. (60)) with these , starting point , sequences with , , and i.i.d. standard Gaussian directions. Let , and let and (, ) be as in the proof. Then for every ,
where
- ;
- ;
- ;
- if , then .
With small this gives the rate of an accelerated method that uses only function values, and linear convergence with ratio in the strongly convex case.
Formalization Note The paper prints . This statement uses : the proof (pp. 549–550) needs , and , all of which hold only with ; with the printed value the first step of the proof fails. The paper's bounds are stated as separate conjuncts; the bound is guarded by because at the paper's value is . Convexity, solvability and are the standing assumptions of problem (53) in Section 5. The oracle at is the limiting oracle . is the Bochner integral of ; under the hypotheses it is integrable.
import Mathlib import Definitions.Def_RandomGradFree_Accelerated_psi import Definitions.Def_RandomGradFree_Accelerated_C import Definitions.Def_RandomGradFree_Accelerated_IsAcceleratedRandomRun open MeasureTheory ProbabilityTheory
namespace RandomGradFree.Accelerated
theorem accelerated_random_method_rate {E : Type*} [NormedAddCommGroup E]
[InnerProductSpace ℝ E] [FiniteDimensional ℝ E] [MeasurableSpace E] [BorelSpace E]
(hdim : 2 ≤ Module.finrank ℝ E)
(f : E → ℝ) (L₁ : ℝ) (hL₁ : 0 < L₁) (hdiff : Differentiable ℝ f)
(hgrad : ∀ x y, ‖gradient f x - gradient f y‖ ≤ L₁ * ‖x - y‖)
(τ : ℝ) (hτ : 0 ≤ τ)
(hsc : ∀ x y, f y ≥ f x + inner ℝ (gradient f x) (y - x) + τ / 2 * ‖y - x‖ ^ 2)
(xstar : E) (hopt : ∀ y, f xstar ≤ f y)
(μ : ℝ) (hμ : 0 ≤ μ)
(θ : ℝ) (hθ : θ = 1 / (16 * ((Module.finrank ℝ E : ℝ) + 4) ^ 2 * L₁))
(h : ℝ) (hh : h = 1 / (4 * ((Module.finrank ℝ E : ℝ) + 4) * L₁))
(x₀ : E) (γ α : ℕ → ℝ)
{Ω : Type*} [MeasurableSpace Ω] (P : Measure Ω) [IsProbabilityMeasure P]
(u x v : ℕ → Ω → E) (hrun : IsAcceleratedRandomRun P f μ τ θ h x₀ γ α u x v) (k : ℕ) :
(∫ ω, f (x k ω) ∂P) - f xstar
≤ psi α k * (f x₀ - f xstar + γ 0 / 2 * ‖x₀ - xstar‖ ^ 2)
+ μ ^ 2 * L₁ * ((Module.finrank ℝ E : ℝ)
+ 3 * ((Module.finrank ℝ E : ℝ) + 8) / 16 * C α k) ∧
psi α k ≤ (1 - Real.sqrt (τ / L₁) / (4 * ((Module.finrank ℝ E : ℝ) + 4))) ^ k ∧
psi α k ≤ 1 / (1 + k / (8 * ((Module.finrank ℝ E : ℝ) + 4)) * Real.sqrt (γ 0 / L₁)) ^ 2 ∧
C α k ≤ k ∧
(0 < τ → C α k ≤ 4 * ((Module.finrank ℝ E : ℝ) + 4) / Real.sqrt (τ / L₁)) := by sorry
end RandomGradFree.Accelerated
Read-back
What the Lean code literally says, in plain math · claude-opus-5-5
Let be a finite-dimensional real inner-product space of dimension . carries its Borel -algebra and is assumed. The following data and hypotheses are given:
- Objective. is differentiable everywhere. Its gradient is -Lipschitz for a constant :
- Strong convexity. is a real number, and for all :
When this is just the first-order convexity inequality.
- Minimizer. is a global minimizer: for every .
- Parameters. is a real number. The reals and are fixed to
Both are strictly positive because .
- Starting point and sequences. is a point. and are arbitrary real sequences. No sign or size condition is imposed on them in this statement itself.
- Probability space. is a measurable space with a probability measure .
- Random sequences. for are three sequences of functions. No measurability or integrability of them is assumed in this statement itself.
- Run hypothesis. The hypothesis
IsAcceleratedRandomRunholds for the tuple . This predicate is defined in an imported file whose code I was not given. Whatever it requires of , , , and (for example a recursion, the distribution of random directions, measurability, or conditions on and ), this statement assumes exactly that and nothing more. I cannot say what it contains or whether it can be satisfied. - Iteration index. is arbitrary.
The statement also uses two imported functions whose definitions I was not given: (psi α k) and (C α k). Both are real numbers that depend only on the sequence and the index . Beyond the inequalities below, the statement asserts nothing about them.
Conclusion. All five of the following hold simultaneously:
- The expected value of at the -th iterate satisfies
- A geometric bound:
- A polynomial bound:
- .
- If , then
Degenerate cases.
- The integral in (1). It is the Lebesgue integral of . If that function is not integrable (or not measurable), the integral counts as , and (1) then reads the right-hand side. Whether
IsAcceleratedRandomRunrules this out cannot be seen here. - . Conjunct (2) becomes , and (3) also becomes . Conjunct (4) becomes . Conjunct (5) gives the upper bound when .
- . Conjunct (2) reduces to , and (5) is vacuous.
- Size of . The Lipschitz-gradient hypothesis together with the strong-convexity inequality forces on a nonzero space, and is nonzero because . Hence , and the base in (2) lies between and .
- Negative . No sign is imposed on here. If , the real square root returns , so (3) becomes . In that case the term in (1) is non-positive.
- . The additive term in (1) vanishes.
- Satisfiability. All hypotheses other than
IsAcceleratedRandomRuncan be satisfied jointly, for example by a quadratic . Whether the run hypothesis is satisfiable, and hence whether the theorem is vacuous, depends on the unseen definition.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.