Eqs. (3.21)–(3.22) — recursion for the expected suboptimality
ProvedStochQuasiNewton.SQN.suboptimality_recursionconvergence-ratep2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1stochastic-optimization
In the setting of Eq. (3.18) (a objective with , a Newton-like iteration with step lengths and constants , ), let minimize and put
Then for every , and are integrable and
This is the recursion from which the rate of Theorem 3.2 is derived.
Formalization Note The step lengths are general (); the recursion does not use . The page's intermediate inequality (3.20) contains the typo and is not restated.
Preamble
import Mathlib import Definitions.Def_StochQuasiNewton_SQN_FiniteSum import Definitions.Def_StochQuasiNewton_SQN_NewtonLike open scoped RealInnerProductSpace open MeasureTheory
Formal statement
namespace StochQuasiNewton.SQN
/-- Eqs. (3.21)–(3.22): in the setting of (3.18), with a minimizer `w*` of `F` and
`φ_k = E[F(w^k) − F(w*)]`, for every `k ≥ 1` the functions `F(w^k)`, `F(w^{k+1})` are integrable
and `φ_{k+1} ≤ (1 − 2 α^k μ₁ λ) φ_k + (Λ/2)(α^k μ₂)² γ²`. -/
theorem suboptimality_recursion {n : ℕ} {Ω Ξ : Type*} [mΩ : MeasurableSpace Ω]
[MeasurableSpace Ξ] (P : Measure Ω) [IsProbabilityMeasure P] (ℱ : Filtration ℕ mΩ)
(F : EuclideanSpace ℝ (Fin n) → ℝ) (hF : ContDiff ℝ 2 F) (lam Lam : ℝ) (hlam : 0 < lam)
(hLam : 0 < Lam)
(hHess : ∀ w : EuclideanSpace ℝ (Fin n),
StrictLoewnerBounds lam Lam (fderiv ℝ (gradient F) w))
(wstar : EuclideanSpace ℝ (Fin n)) (hwstar : IsMinOn F Set.univ wstar)
(G : EuclideanSpace ℝ (Fin n) → Ξ → EuclideanSpace ℝ (Fin n)) (γ : ℝ) (α : ℕ → ℝ)
(hα : ∀ k, 1 ≤ k → 0 < α k) (μ₁ μ₂ : ℝ) (hμ₁ : 0 < μ₁) (hμ : μ₁ ≤ μ₂)
(w1 : EuclideanSpace ℝ (Fin n)) (w : ℕ → Ω → EuclideanSpace ℝ (Fin n)) (ξ : ℕ → Ω → Ξ)
(H : ℕ → Ω → Matrix (Fin n) (Fin n) ℝ)
(hrun : NewtonLikeIteration P ℱ F G γ α μ₁ μ₂ w1 w ξ H) (k : ℕ) (hk : 1 ≤ k) :
Integrable (fun ω => F (w k ω)) P ∧ Integrable (fun ω => F (w (k + 1) ω)) P ∧
∫ ω, (F (w (k + 1) ω) - F wstar) ∂P ≤
(1 - 2 * α k * μ₁ * lam) * ∫ ω, (F (w k ω) - F wstar) ∂P +
Lam / 2 * (α k * μ₂) ^ 2 * γ ^ 2 := by sorry
end StochQuasiNewton.SQN
Source
Byrd, Hansen, Nocedal, Singer, A Stochastic Quasi-Newton Method for Large-Scale Optimization, SIAM J. Optim. 26(2) (2016), p. 1017, Eqs. (3.21)–(3.22)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.