Newton-decrement contraction of the pure step
ProvedConvexOptimization.sc_newton_decrement_contractionQuadratic contraction of the Newton decrement — inequality (9.55) of Boyd & Vandenberghe.
Let be open and convex, let be self-concordant on with positive definite Hessian and with closed sublevel sets . Let , let solve , and let be the Newton decrement at , . If , then the full Newton step remains in the domain, , and the decrement at the new point satisfies
Two things are asserted at once, and both matter. The domain statement is not automatic — for a barrier objective, leaving means violating a constraint — and self-concordance is precisely what guarantees the unit step is safe once . The contraction then gives quadratic convergence beyond an absolute threshold: for , say, the decrement is at most squared at each step, so the number of iterations to reach any accuracy is bounded without reference to problem constants.
Formalization Note The conclusion quantifies over any Newton direction and any nonnegative with , rather than asserting uniqueness of the decrement. The closed-sublevel-set hypothesis is the standing assumption of B&V §9.6.4. Source: B&V §9.6.4, p. 505, eq. (9.55) / exercise 9.18; proof in Nesterov, Lectures on Convex Optimization, Theorem 5.2.2.
import Mathlib import Definitions.Def_ConvexOptimization_selfConcordance open scoped RealInnerProductSpace ENNReal open MeasureTheory
theorem ConvexOptimization.sc_newton_decrement_contraction {n : ℕ}
(Ω : Set (EuclideanSpace ℝ (Fin n))) (hΩo : IsOpen Ω) (hΩc : Convex ℝ Ω)
(f : EuclideanSpace ℝ (Fin n) → ℝ) (hsc : IsSelfConcordantOn Ω f)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(hg : ∀ x ∈ Ω, HasGradientAt f (g x) x)
(H : EuclideanSpace ℝ (Fin n) →
EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n))
(hH : ∀ x ∈ Ω, HasFDerivAt g (H x) x)
(hHpd : ∀ x ∈ Ω, ∀ v, v ≠ 0 → 0 < ⟪H x v, v⟫)
(hclosed : ∀ c : ℝ, IsClosed {y | y ∈ Ω ∧ f y ≤ c})
(x : EuclideanSpace ℝ (Fin n)) (hx : x ∈ Ω)
(Δ : EuclideanSpace ℝ (Fin n)) (hΔ : H x Δ = -g x)
(lam : ℝ) (hlam0 : 0 ≤ lam) (hlam : lam ^ 2 = ⟪g x, -Δ⟫) (hlt : lam < 1) :
x + Δ ∈ Ω ∧
∀ (Δ' : EuclideanSpace ℝ (Fin n)) (lam' : ℝ),
H (x + Δ) Δ' = -g (x + Δ) → 0 ≤ lam' → lam' ^ 2 = ⟪g (x + Δ), -Δ'⟫ →
lam' ≤ (lam / (1 - lam)) ^ 2 := by
sorry
Read-back
What the Lean code literally says, in plain math · claude-fable-5
Fix a natural number , a set that is open and convex (two hypotheses), and a total function that is self-concordant on , meaning: is convex and for , , ; is three times continuously differentiable on (within-set sense); and for every , , the line restriction satisfies at (total, possibly junk, derivatives; real power with for ). Further hypotheses: with differentiable at each with gradient ; assigning to each point a continuous linear operator, with Fréchet-differentiable at each with derivative ; for every and , (symmetry not assumed); and for EVERY real , the sublevel set is closed. (No minimizer is hypothesized in this theorem.) Given: ; with ; and a real with , (forcing ; vacuous otherwise), and (strict). Conclusion, a conjunction of two claims: (1) — the full, undamped Newton step from lands in ; and (2) for every vector and real such that , , and , one has
Note carefully: the bound is on the first power of , compared against the square of the ratio — it is , not a statement about . If no pair satisfies the three conditions at , claim (2) is vacuously true.
Confirmed by the mission captain (proposal self-audit).