Suboptimality bound from strong convexity
ProvedConvexOptimization.strong_convexity_quadratic_lower_boundSuboptimality is controlled by the gradient norm — inequality (9.9) of Boyd & Vandenberghe, the standard stopping criterion for unconstrained minimization.
Let have gradient field , and let be such that satisfies the strong-convexity lower bound
Let be a global minimizer of and write for the optimal value. Then for every
The bound converts a computable quantity, the gradient norm at the current iterate, into a certificate of suboptimality: already guarantees . It is what turns the gradient contraction of Newton's quadratically convergent phase into a bound on the objective error, and hence is used directly in the mission's goal theorem.
Formalization Note The gradient is an explicit field g with ∀ x, HasGradientAt f (g x) x; the minimizer is stated as IsMinOn f Set.univ xstar, and appears as f xstar. Strong convexity enters as the displayed inequality for all rather than through a Hessian hypothesis, which keeps the statement usable for functions that are not twice differentiable. Source: B&V §9.1.2 p. 460, eq. (9.9).
import Mathlib open scoped RealInnerProductSpace ENNReal open MeasureTheory
theorem ConvexOptimization.strong_convexity_quadratic_lower_bound {n : ℕ} (m : ℝ) (hm : 0 < m)
(f : EuclideanSpace ℝ (Fin n) → ℝ)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(hg : ∀ x, HasGradientAt f (g x) x)
(hsc : ∀ x y : EuclideanSpace ℝ (Fin n),
f x + ⟪g x, y - x⟫ + m / 2 * ‖y - x‖ ^ 2 ≤ f y)
(xstar : EuclideanSpace ℝ (Fin n)) (hstar : IsMinOn f Set.univ xstar)
(x : EuclideanSpace ℝ (Fin n)) :
f x - f xstar ≤ ‖g x‖ ^ 2 / (2 * m) := by
sorry