Theorem 5.2, p. 22 — with Lα_k² = γ_k(1 + α_k), the gradient-corrected scheme (91) has ℒ_{k+1} ≤ ℒ_k/(1 + α_k), (87) and (88)
ProvedNAGFlow.GradCorr.theorem_5_2Let be a real Hilbert space and with , and let be a global minimiser of . Let satisfy the parameter equation (73),
with the step sizes chosen by , and let be a run of the corrected scheme (91),
Let (74). Then
- for every ,
- for every ,
- for every ,
where
The corrected scheme evaluates the gradient once per step, at , and attains the same accelerated rates as the predictor–corrector scheme of Theorem 5.1: for convex and the linear rate for strongly convex , in one statement covering and .
Formalization Note. is a real Hilbert space, and the dual norm is its norm (Riesz). The gradient is an explicit map with HasGradientAt at every point. The scheme and (73) are hypotheses on given sequences; is any positive sequence satisfying , not a formula. The positivity and are those of (73) and of the scheme; is part of the class. The paper states (87) and (88) by reference ("both the two estimates (87) and (88) hold true here"); they are written out in full, with the constants of p. 21. The negative and powers are integer powers of a number . When the second entry of each minimum equals .
import Mathlib import Definitions.Def_NAGFlow_GradCorr_Setting
namespace NAGFlow.GradCorr
/-- Theorem 5.2 (Luo & Chen, arXiv:1909.03145v4, p. 22). Let `f ∈ S^{1,1}_{μ,L}` with
`0 ≤ μ ≤ L < ∞`, let `x*` be a global minimiser of `f`, and let `(x, y, v)` be a run of the corrected
scheme (91) together with (73) (`γ₀ > 0`, `α_k > 0`), with step sizes chosen by
`Lα_k² = γ_k(1 + α_k)`. Then, with `ℒ_k` the Lyapunov function (74):
(92) `ℒ_{k+1} ≤ ℒ_k/(1 + α_k)` for every `k ∈ ℕ`;
(87) `ℒ_k ≤ ℒ₀ · min{4L/(√γ₀ k + 2√L)², (1 + √(min{γ₀, μ}/L))^{−k}}` for every `k ≥ 0`;
(88) `ℒ_k ≤ C_{γ₀,L} · min{4/k², (1 + √(min{γ₀, μ}/L))^{1−k}}` for every `k ≥ 1`, where
(89) `C_{γ₀,L} = (L/γ₀)(f(x₀) − f(x*)) + (L/2)‖v₀ − x*‖²`. -/
theorem theorem_5_2 {V : Type*} [NormedAddCommGroup V] [InnerProductSpace ℝ V] [CompleteSpace V]
(f : V → ℝ) (gradf : V → V) (μ L : ℝ) (hf : NAGFlow.PredCorr.IsS11 f gradf μ L)
(xstar : V) (hxstar : ∀ y, f xstar ≤ f y)
(α γ : ℕ → ℝ) (hγ : NAGFlow.PredCorr.IsGammaRun μ α γ)
(hstep : ∀ k, L * α k ^ 2 = γ k * (1 + α k))
(x y v : ℕ → V) (hrun : IsGCRun gradf μ L α γ x y v) :
(∀ k : ℕ, NAGFlow.PredCorr.lyap f xstar x v γ (k + 1) ≤ NAGFlow.PredCorr.lyap f xstar x v γ k / (1 + α k)) ∧
(∀ k : ℕ, NAGFlow.PredCorr.lyap f xstar x v γ k ≤
NAGFlow.PredCorr.lyap f xstar x v γ 0 *
min (4 * L / (Real.sqrt (γ 0) * k + 2 * Real.sqrt L) ^ 2)
((1 + Real.sqrt (min (γ 0) μ / L)) ^ (-(k : ℤ)))) ∧
(∀ k : ℕ, 1 ≤ k → NAGFlow.PredCorr.lyap f xstar x v γ k ≤
(L / γ 0 * (f (x 0) - f xstar) + L / 2 * ‖v 0 - xstar‖ ^ 2) *
min (4 / (k : ℝ) ^ 2) ((1 + Real.sqrt (min (γ 0) μ / L)) ^ (1 - (k : ℤ)))) := by sorry
end NAGFlow.GradCorr
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.