p. 156 —
ProvedNumStochOpt.Nonstationary.lemma_p156_step_boundp2o-batch-b23ap2o-gran-per-chapterp2o-plan-bookp2o-v1projectionstochastic-optimizationsubgradient-method
Let be a convex compact set, , and let satisfy for all . Let satisfy the projected subgradient recursion (6.41)
If and , then
In the proof of Theorem 6.3 this bounds the distance travelled between the index of a subsequence and the exit time , which converts the decrease of into a decrease proportional to .
Formalization Note The book writes "where is a constant"; the proof yields = the constant of hypothesis (d) of Theorem 6.3, . The hypothesis holds for every (all iterates after the first are projections onto ); in the book is large. is the step-size convention of the chapter, not printed in Theorem 6.3.
Preamble
import Mathlib import Definitions.Def_NumStochOpt_QuasiFejer_ProjectionMethod
Formal statement
namespace NumStochOpt.Nonstationary
/-- Proof of Theorem 6.3, p. 156 (unnumbered display): "in view of the properties of `π_X`",
`‖x^τ - x^a‖ ≤ ∑_{s=a}^{τ-1} ‖x^{s+1} - x^s‖ ≤ C ∑_{s=a}^{τ-1} ρ_s` along the iteration (6.41)
`x^{s+1} = π_X[x^s - ρ_s g_s]`, where `C` is the bound of hypothesis (d), `‖g_s‖ ≤ C`, and
`x^a ∈ X` (which holds for every `a ≥ 1`). -/
theorem lemma_p156_step_bound {n : ℕ}
(X : Set (EuclideanSpace ℝ (Fin n))) (hXconv : Convex ℝ X) (hXcpt : IsCompact X)
(x g : ℕ → EuclideanSpace ℝ (Fin n)) (ρ : ℕ → ℝ) (C : ℝ)
(hrec : ∀ s, x (s + 1) = NumStochOpt.QuasiFejer.projX X (x s - ρ s • g s))
(hρnn : ∀ s, 0 ≤ ρ s) (hbound : ∀ s, ‖g s‖ ≤ C)
(a b : ℕ) (hab : a ≤ b) (hxa : x a ∈ X) :
‖x b - x a‖ ≤ ∑ s ∈ Finset.Ico a b, ‖x (s + 1) - x s‖ ∧
∑ s ∈ Finset.Ico a b, ‖x (s + 1) - x s‖ ≤ C * ∑ s ∈ Finset.Ico a b, ρ s := by sorry
end NumStochOpt.Nonstationary
Source
Yu. Ermoliev, "Stochastic Quasigradient Methods", in Ermoliev & Wets (eds.), Numerical Techniques for Stochastic Optimization, Springer 1988, Ch. 6, p. 156, proof of Theorem 6.3, unnumbered display
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.