Theorem 2.13 — Polyak's method contracts by the factor at a sharp minimum
ProvedShorNonsmooth.Fejer.polyak_linear_rate_sharp_minimumLet be convex on with a minimum point , let be such that
and let bound the norm of every subgradient of at every point of the ball . Let and let be generated from by Polyak's method (2.32) with , for any subgradient selection . Then for every
At a sharp minimum Polyak's step converges linearly without any smoothness.
Formalization Note The book assumes " the Lipschitz constant of in the area ", and its proof uses . A Lipschitz constant on the closed ball does not bound subgradients at boundary points of the ball, and the per-step inequality then fails (e.g. on with a steeper slope outside, , ). The statement therefore assumes the subgradient bound the proof uses; it holds whenever is -Lipschitz on an open neighbourhood of the ball. The range is the section's standing range.
import Mathlib import Definitions.Def_ShorNonsmooth_Fejer_PolyakMethod open Filter Topology
namespace ShorNonsmooth.Fejer
/-- Shor (1985), p. 38, Theorem 2.13. Let `f` be convex on `E_n` with a minimum point `x*`, let
`m > 0` be such that `f(x) - f(x*) ≥ m ‖x - x*‖` for all `x`, and let `L` bound the norm of every
subgradient of `f` at every point of the ball `‖x - x*‖ ≤ ‖x₀ - x*‖` (the book: "`L` the Lipschitz
constant of `f` in the area `‖x - x*‖ ≤ ‖x₀ - x*‖`"; see the note below). Then Polyak's method
(2.32) with `c = f(x*)` and `0 < γ < 2`, for any subgradient selection, satisfies
`‖x_{k+1} - x*‖ ≤ q ‖x_k - x*‖` for all `k`, where `q = (1 - γ(2 - γ) m² / L²)^{1/2}`.
Note: the hypothesis is the bound `‖g‖ ≤ L` that the book's proof uses; a Lipschitz constant of
`f` on the closed ball alone does not bound subgradients at boundary points of the ball, and the
per-step inequality can then fail. -/
theorem polyak_linear_rate_sharp_minimum {n : ℕ} (f : EuclideanSpace ℝ (Fin n) → ℝ)
(hf : ConvexOn ℝ Set.univ f) (xstar : EuclideanSpace ℝ (Fin n))
(hmin : ∀ y, f xstar ≤ f y) (m L : ℝ) (hm : 0 < m)
(hsharp : ∀ x, f x - f xstar ≥ m * ‖x - xstar‖)
(γ : ℝ) (hγ0 : 0 < γ) (hγ2 : γ < 2)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(hg : ∀ x, ShorNonsmooth.AlmostDiff.IsSubgradient f x (g x))
(x : ℕ → EuclideanSpace ℝ (Fin n))
(hx : ∀ k, x (k + 1) = polyakStep f g (f xstar) γ (x k))
(hL : ∀ z ∈ Metric.closedBall xstar ‖x 0 - xstar‖, ∀ v, ShorNonsmooth.AlmostDiff.IsSubgradient f z v → ‖v‖ ≤ L) :
∀ k : ℕ, ‖x (k + 1) - xstar‖ ≤
Real.sqrt (1 - γ * (2 - γ) * m ^ 2 / L ^ 2) * ‖x k - xstar‖ := by sorry
end ShorNonsmooth.Fejer
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.