Theorem 15, p. 29 — (FDM) converges linearly in function values on F_{L_f,κ_f} with rate 1/(1 + Lκ_f/(4(L_f + L̄_f + βL̄_f)²))
ProvedNecoaraNG.FDM.theorem_15Consider problem : , where is closed and convex, is convex on and differentiable at every point of with gradient Lipschitz on with constant ,
and the optimal set is nonempty, with and . Assume has quadratic functional growth with constant : for every and every nearest point , . Together these say belongs to the class .
Let and let , , be a run of the feasible descent method (FDM): and, for every ,
Then for every
The theorem shows that the whole family of feasible descent methods, which includes proximal point, coordinate descent, extragradient and matrix splitting methods, converges linearly in function values under quadratic functional growth, without strong convexity and without assuming an error bound.
Formalization Note is EuclideanSpace ℝ (Fin n); is defined on all of and only its values and gradients at points of matter. The projection is the predicate IsNearest, so a run is any pair of sequences satisfying the conditions, and the errors are quantified together with the iterates. "Simple" (cheap projection) and "closed" for are dropped: the first is not a mathematical condition and the second follows from continuity on . No upper bound on is assumed, as on the page.
import Mathlib import Definitions.Def_NecoaraNG_FDM_Setting
namespace NecoaraNG.FDM
open scoped InnerProductSpace
/-- Theorem 15, p. 29, (64): for `f` in the class `F_{L_f, κ}(X)` (convex, (1), (22)), every run of
the feasible descent method (FDM) satisfies
`f(x^k) - f* ≤ (1 / (1 + L κ / (4 (L_f + L̄_f + β L̄_f)²)))^k (f(x⁰) - f*)`, with `f* = f x*`. -/
theorem theorem_15 {n : ℕ} (X : Set (NecoaraNG.Chain.E n)) (hXc : IsClosed X) (hXconv : Convex ℝ X)
(f : NecoaraNG.Chain.E n → ℝ) (hf : ConvexOn ℝ X f) (hdiff : ∀ x ∈ X, DifferentiableAt ℝ f x)
(Lf : ℝ) (hLf : 0 < Lf)
(hL : ∀ x ∈ X, ∀ y ∈ X, ‖gradient f x - gradient f y‖ ≤ Lf * ‖x - y‖)
(xstar : NecoaraNG.Chain.E n) (hxstar : xstar ∈ NecoaraNG.Chain.optSet X f)
(κ : ℝ) (hκ : 0 < κ) (hF : NecoaraNG.Chain.QuadFunGrowth X f κ)
(β L Lbar : ℝ) (hβ : 0 < β) (hLpos : 0 < L) (hLbar : 0 < Lbar)
(α : ℕ → ℝ) (e x : ℕ → NecoaraNG.Chain.E n) (hrun : IsFDMRun X f β L Lbar α e x) :
∀ k : ℕ, f (x k) - f xstar ≤
(1 / (1 + L * κ / (4 * (Lf + Lbar + β * Lbar) ^ 2))) ^ k * (f (x 0) - f xstar) := by sorry
end NecoaraNG.FDM
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.