Theorem 6.1, p. 24 — Nesterov's method (Algorithm 1) has 0 < α_k ≤ 1, ℒ_{k+1} ≤ (1 − α_k)ℒ_k, (99) and (100)
ProvedNAGFlow.Nesterov.theorem_6_1Let be a real Hilbert space and with , that is, is continuously differentiable, -convex, and has -Lipschitz gradient. Let be a global minimiser of . Let be a run of Nesterov's accelerated gradient method (Algorithm 1): , and for every , solves , (so ), , , and . Let
and (89). Then for every ,
for all
and for all
The theorem gives a Lyapunov-function proof of the optimal accelerated rates of Nesterov's method, simultaneously for the convex case (, rate ) and the strongly convex case (, linear rate once ), interpreting the method as a corrected semi-implicit discretization of the NAG flow.
Formalization Note. and are the inner product and norm of (Riesz); is an explicit gradient map. Algorithm 1 writes ; the run records only (the positive root of step 2), and is concluded. Step 5 is an inequality, so the theorem covers every sufficient-decrease choice of . The theorem's hypothesis "" follows from steps 2–3 and is not restated. is the value produced by the run. is written out inline.
import Mathlib import Definitions.Def_NAGFlow_Nesterov_Setting
namespace NAGFlow.Nesterov
/-- Theorem 6.1 (Luo & Chen, arXiv:1909.03145v4, p. 24). 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
Algorithm 1 (equivalently, of (96)–(97) with `Lα_k² = γ_{k+1}`). Then, with `ℒ_k` the Lyapunov
function (74):
* `0 < α_k ≤ 1` for every `k`;
* (98) `ℒ_{k+1} ≤ (1 − α_k)ℒ_k` for every `k ∈ ℕ`;
* (99) `ℒ_k ≤ ℒ₀ · min{4L/(√γ₀ k + 2√L)², (1 − √(min{γ₁, μ}/L))^k}` for every `k ≥ 0`;
* (100) `ℒ_k ≤ C_{γ₀,L} · min{4/k², (1 − √(min{γ₁, μ}/L))^{k−1}}` for every `k ≥ 1`, where
(89) `C_{γ₀,L} = (L/γ₀)(f(x₀) − f(x*)) + (L/2)‖v₀ − x*‖²`. -/
theorem theorem_6_1 {V : Type*} [NormedAddCommGroup V] [InnerProductSpace ℝ V] [CompleteSpace V]
(f : V → ℝ) (gradf : V → V) (μ L : ℝ) (hf : NAGFlow.PredCorr.IsS11 f gradf μ L)
(xstar : V) (hxstar : ∀ z, f xstar ≤ f z)
(α γ : ℕ → ℝ) (x y v : ℕ → V) (hrun : IsNAGRun f gradf μ L α γ x y v) :
(∀ k : ℕ, 0 < α k ∧ α k ≤ 1) ∧
(∀ k : ℕ, lyap f xstar x v γ (k + 1) ≤ (1 - α k) * lyap f xstar x v γ k) ∧
(∀ k : ℕ, lyap f xstar x v γ k ≤
lyap f xstar x v γ 0 *
min (4 * L / (Real.sqrt (γ 0) * k + 2 * Real.sqrt L) ^ 2)
((1 - Real.sqrt (min (γ 1) μ / L)) ^ k)) ∧
(∀ k : ℕ, 1 ≤ k → 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 (γ 1) μ / L)) ^ (k - 1))) := by sorry
end NAGFlow.Nesterov
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.