Theorem 3.7 — proximal quasi-Newton methods converge q-superlinearly under the Dennis–Moré criterion
ProvedProxNewton.Exact.prox_quasi_newton_superlinearLet be twice continuously differentiable and strongly convex with constant , with and Lipschitz continuous with constants and . Let be a proper closed convex function (possibly taking the value ) with domain , and let be an optimal solution of
Fix , a backtracking factor , and . Let be a run of Algorithm 1 (a proximal quasi-Newton method) from any such that for every and satisfies the Dennis–Moré criterion
Then converges to q-superlinearly:
This is the composite counterpart of the Dennis–Moré characterization of superlinear convergence of quasi-Newton methods, and it covers proximal BFGS and L-BFGS-type methods whenever their Hessian approximations satisfy the criterion.
Formalization Note The little- is stated without quotients: for every , eventually . The line search is backtracking from the unit step with factor (the least with satisfying (2.19)); a line search allowed to return any admissible step would make the theorem false. The same constant bounds 's strong convexity and from below, as printed.
import Mathlib import Definitions.Def_ProxNewton_Exact_Basic open scoped RealInnerProductSpace Topology open Filter
namespace ProxNewton.Exact
/-- Theorem 3.7, arXiv:1206.1623v13, p. 14. Let `g` be twice continuously differentiable and
strongly convex with constant `m > 0`, with `∇g` and `∇²g` Lipschitz with constants `L1`, `L2`;
let `h` be proper closed convex with domain `D`, and `xstar` an optimal solution of (1.1). Let
`(x, H, Δ, t)` be a run of Algorithm 1 (a proximal quasi-Newton method) with `α ∈ (0, 1/2)` and
backtracking factor `β ∈ (0, 1)`, from any `x 0 ∈ D`. If `mI ⪯ H k ⪯ MI` for all `k`, with
`0 < m ≤ M`, and `{H k}` satisfies the Dennis–Moré criterion (3.2), then `x k → xstar`
q-superlinearly: `x k → xstar` and, for every `ε > 0`, eventually
`‖x (k+1) − xstar‖ ≤ ε ‖x k − xstar‖`. -/
theorem prox_quasi_newton_superlinear {n : ℕ} (g : EuclideanSpace ℝ (Fin n) → ℝ)
(D : Set (EuclideanSpace ℝ (Fin n))) (h : EuclideanSpace ℝ (Fin n) → ℝ)
(m M L1 L2 α β : ℝ)
(hg : ContDiff ℝ 2 g) (hm : 0 < m) (hsc : StronglyConvexWith g m) (hL1 : 0 ≤ L1)
(hgL : ∀ x y : EuclideanSpace ℝ (Fin n), ‖gradient g x - gradient g y‖ ≤ L1 * ‖x - y‖)
(hL2 : 0 ≤ L2)
(hHL : ∀ x y : EuclideanSpace ℝ (Fin n), ‖hessian g x - hessian g y‖ ≤ L2 * ‖x - y‖)
(hD : IsProperClosedConvex D h)
(xstar : EuclideanSpace ℝ (Fin n)) (hstar : IsMinimizer g D h xstar)
(hα : 0 < α) (hα2 : α < 1 / 2) (hβ : 0 < β) (hβ1 : β < 1) (hmM : m ≤ M)
(x : ℕ → EuclideanSpace ℝ (Fin n))
(H : ℕ → EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n))
(Δ : ℕ → EuclideanSpace ℝ (Fin n)) (t : ℕ → ℝ)
(hrun : IsProxNewtonTypeRun g D h α β x H Δ t)
(hHb : ∀ k : ℕ, IsBoundedBetween (H k) m M)
(hDM : DennisMore g xstar x H) :
Tendsto x atTop (𝓝 xstar) ∧
∀ ε : ℝ, 0 < ε → ∀ᶠ k in atTop, ‖x (k + 1) - xstar‖ ≤ ε * ‖x k - xstar‖ := by sorry
end ProxNewton.Exact
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.