§5.3.2, proof of Theorem 5.3, p. 321 — ∫₀¹ ‖∇²f(x_k) − ∇²f(x* + s(x_k − x*))‖ ds ≤ (M/2)‖x_k − x*‖
OpenConvexOptAlg.Newton.lipschitz_integral_boundconvex-optimizationhessiannewton-methodp2o-batch-pfp2ap2o-gran-per-chapterp2o-plan-bookp2o-v1
Let be -Lipschitz in operator norm, i.e. for all (for instance, the Hessian of a function with -Lipschitz Hessian). Then for all ,
In the analysis of Newton's method, with and , this bounds the integral in the error representation of a Newton step by , which produces the quadratic rate.
Formalization Note The book states the bound at an iterate for the Hessian of ; it uses only the Lipschitz property, so it is stated for any -Lipschitz map and any point . No sign condition on is needed (a negative makes the hypothesis unsatisfiable unless ). is EuclideanSpace ℝ (Fin n) and on linear maps is the operator norm.
Preamble
import Mathlib import Definitions.Def_ConvexOptAlg_Newton_Defs
Formal statement
namespace ConvexOptAlg.Newton
/-- The Lipschitz integral bound in the proof of Theorem 5.3 (Bubeck, arXiv:1405.4980v2, §5.3.2,
p. 321, fifth display of the proof): if the Hessian map `H` is `M`-Lipschitz in operator norm, then
for all points `x∗, y ∈ ℝⁿ`,
`∫₀¹ ‖∇²f(y) − ∇²f(x∗ + s(y − x∗))‖ ds ≤ (M/2)‖y − x∗‖`.
The book states it at `y = x_k`; it uses only the Lipschitz property, so it is stated for any
`M`-Lipschitz map `H : ℝⁿ → L(ℝⁿ, ℝⁿ)` and any `y`. -/
theorem lipschitz_integral_bound {n : ℕ}
(H : EuclideanSpace ℝ (Fin n) → (EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n)))
(M : ℝ) (hHL : IsLipschitzHessian H M) (xstar y : EuclideanSpace ℝ (Fin n)) :
∫ s in (0 : ℝ)..1, ‖H y - H (xstar + s • (y - xstar))‖ ≤ M / 2 * ‖y - xstar‖ := by sorry
end ConvexOptAlg.Newton
Source
Bubeck, arXiv:1405.4980v2, §5.3.2, proof of Theorem 5.3, p. 321, fifth display