Eqs. (3.8)–(3.10) — curvature bounds for the correction pairs
ProvedStochQuasiNewton.SQN.curvature_pair_boundsbfgscurvaturep2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1quasi-newton
Let be functions on , let , , and suppose that every subsampled Hessian with samples satisfies for all (Assumption 1(2), Eq. (3.3)). Let be a sample with , let , let , and let as in (3.8). Then
These are the curvature bounds (3.9) and (3.10) for the correction pairs of the SQN method; they bound the scaling and each rank-two term of the L-BFGS update in the proof of Lemma 3.1.
Formalization Note The statement is for an arbitrary nonzero and arbitrary point ; in the paper and .
Preamble
import Mathlib import Definitions.Def_StochQuasiNewton_SQN_FiniteSum open scoped RealInnerProductSpace
Formal statement
namespace StochQuasiNewton.SQN
/-- Eqs. (3.8)–(3.10): if every subsampled Hessian with `b_H` samples satisfies
`λ I ≺ ∇²F_{S_H}(w) ≺ Λ I` (3.3) with `0 < λ` and `0 < Λ`, then a correction pair `y = ∇²F_{S_H}(w̄) s` with `s ≠ 0`
satisfies `λ ‖s‖² ≤ yᵀs ≤ Λ ‖s‖²` (3.9) and `λ ≤ ‖y‖² / (yᵀs) ≤ Λ` (3.10). -/
theorem curvature_pair_bounds {n N : ℕ} (f : Fin N → EuclideanSpace ℝ (Fin n) → ℝ)
(hf : ∀ i, ContDiff ℝ 2 (f i)) (lam Lam : ℝ) (hlam : 0 < lam) (hLam : 0 < Lam) (bH : ℕ)
(hHess : ∀ SH : Finset (Fin N), SH.card = bH → ∀ w : EuclideanSpace ℝ (Fin n),
StrictLoewnerBounds lam Lam (subsampledHessian f SH w))
(SH : Finset (Fin N)) (hSH : SH.card = bH) (wbar s y : EuclideanSpace ℝ (Fin n))
(hs : s ≠ 0) (hy : y = subsampledHessian f SH wbar s) :
(lam * ‖s‖ ^ 2 ≤ ⟪y, s⟫ ∧ ⟪y, s⟫ ≤ Lam * ‖s‖ ^ 2) ∧
(lam ≤ ‖y‖ ^ 2 / ⟪y, s⟫ ∧ ‖y‖ ^ 2 / ⟪y, s⟫ ≤ Lam) := 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. 1015, Eqs. (3.8)–(3.10)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.