Theorem 3.11 — forcing terms (2.25) give q-superlinear convergence
ProvedProxNewton.Inexact.ew_forcing_superlinearAssume the standing assumptions of §3.4 ( twice continuously differentiable and strongly convex with constant , -Lipschitz, -Lipschitz, proper closed convex with domain ), let be the optimal solution of , and let satisfy for all . There is with the following property. Consider any run of the inexact proximal Newton method with unit step lengths under the stopping condition (2.24), whose forcing terms are chosen by
with any first forcing term . If , then q-superlinearly: and, for every , for all sufficiently large .
The forcing terms (2.25) need no tuning: they measure how well the previous model predicted the composite gradient step, and they make the subproblems be solved more accurately exactly when the iterates approach .
Formalization Note (2.25) does not define ; the theorem allows any , the range (2.25) gives the later terms, and does not depend on or on the run. The forcing terms depend on the iterates, unlike in Theorem 3.10. When the quotient is in Lean; this happens only at .
import Mathlib import Definitions.Def_ProxNewton_Inexact_CompositeStep import Definitions.Def_ProxNewton_Inexact_Standing import Definitions.Def_ProxNewton_Inexact_Method
namespace ProxNewton.Inexact
open Filter Topology
/-- Theorem 3.11: with the forcing terms (2.25) (and any first forcing term `η_0 ∈ [0, m/2]`),
the inexact proximal Newton method with unit step lengths converges q-superlinearly to `x⋆` from
every start sufficiently close to `x⋆`. -/
theorem ew_forcing_superlinear {n : ℕ} (g : EuclideanSpace ℝ (Fin n) → ℝ)
(D : Set (EuclideanSpace ℝ (Fin n))) (h : EuclideanSpace ℝ (Fin n) → ℝ) (m L1 L2 M : ℝ)
(xstar : EuclideanSpace ℝ (Fin n))
(hg : SmoothPartAssumptions g m L1 L2) (hh : IsProperClosedConvex D h)
(hxstar : IsMinimizer g D h xstar) (hM : 0 < M) (hHM : HessianLE g M) :
∃ δ : ℝ, 0 < δ ∧ ∀ (η : ℕ → ℝ) (x Δ : ℕ → EuclideanSpace ℝ (Fin n)),
0 ≤ η 0 → η 0 ≤ m / 2 → (∀ k, η (k + 1) = ewForcingTerm g D h M m x k) →
IsInexactProxNewtonRun g D h M η x Δ → ‖x 0 - xstar‖ < δ →
Tendsto x atTop (𝓝 xstar) ∧
∀ ε : ℝ, 0 < ε → ∀ᶠ k in atTop, ‖x (k + 1) - xstar‖ ≤ ε * ‖x k - xstar‖ := by sorry
end ProxNewton.Inexact
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.