Quadratic (descent-lemma) upper bound for along a trial step inside
ProvedNonmonotoneLS.RLinear.lbar_quadratic_upper_boundline-searchnonmonotonep2o-batch-p100bp2o-gran-per-chapterp2o-plan-paperp2o-v1
Let be , let , be the iterates and directions of a run, and suppose is -Lipschitz on with .
Claim. If and , then
This is the quadratic upper bound (the descent lemma) along the line-search trial step. Every point , , lies in , so Lipschitzness of applies on the whole segment; integrating along the segment gives the claim.
The estimate is what turns the failure of the Armijo test at the next trial step into a lower bound on , and hence into the step used in Eq. (2.11) of the proof of Theorem 2.1 (p. 1048), with the (3.7) argument applied in the setting of Theorem 3.1 (p. 1050).
Preamble
import Mathlib import Definitions.Def_NonmonotoneLS_Shared_Params import Definitions.Def_NonmonotoneLS_Shared_Cost import Definitions.Def_NonmonotoneLS_Shared_Steps import Definitions.Def_NonmonotoneLS_RLinear_Run import Definitions.Def_NonmonotoneLS_RLinear_Regions import Definitions.Def_NonmonotoneLS_RLinear_Constants import Theorems.Thm_NonmonotoneLS_RLinear_step_segment_mem_Lbar open scoped InnerProductSpace NNReal open Filter
Formal statement
namespace NonmonotoneLS.RLinear
/-- Descent/smoothness estimate on the trial segment: if `x k ∈ 𝓛`, `0 ≤ s ≤ μ` and `∇f` is
`L`-Lipschitz on `𝓛̄`, then `f (x k + s • d k) ≤ f (x k) + s ⟨∇f (x k), d k⟩ + (L/2) s² ‖d k‖²`. -/
theorem lbar_quadratic_upper_bound {n : ℕ} (p : Shared.Params)
(f : EuclideanSpace ℝ (Fin n) → ℝ) (hf : ContDiff ℝ 1 f)
(x d : ℕ → EuclideanSpace ℝ (Fin n)) (k : ℕ) (s : ℝ)
(hx : x k ∈ levelSet f (x 0)) (hs0 : 0 ≤ s) (hsμ : s ≤ p.μ)
(L : ℝ≥0) (hLip : LipschitzOnWith L (gradient f) (Lbar p f x d)) :
f (x k + s • d k) ≤ f (x k) + s * ⟪gradient f (x k), d k⟫_ℝ
+ (L : ℝ) / 2 * s ^ 2 * ‖d k‖ ^ 2 := by
sorry
end NonmonotoneLS.RLinearSource
Zhang, Hager, A Nonmonotone Line Search Technique and Its Application to Unconstrained Optimization, SIAM J. Optim. 14 (2004), Lemma 2.1 and Eq. (2.10)-(2.11) (p. 1048), specialised to the mu*d_max region of Theorem 3.1 (p. 1050)