Theorem 4: proximal gradient descends and its cluster points are stationary for
ProvedNonconvexSplitting.ProxGrad.theorem4Consider on , where is twice continuously differentiable with a bounded Hessian and is proper and closed with well-defined proximal mappings.
Suppose there exist a twice continuously differentiable convex function and such that for all
Let be generated by the proximal gradient update
with . Then:
- the algorithm is a descent algorithm: for every ;
- every cluster point of , if one exists, is a stationary point: , with the limiting subdifferential.
The constant bounds the curvature of only after its concave part is offset by the convex , so it may be smaller than the Lipschitz modulus of ; the result then permits step sizes beyond the classical range .
Formalization Note Objective values are compared in EReal. A cluster point is a limit of for a strictly increasing . The start is arbitrary, and no boundedness of the sequence is assumed.
import Mathlib import Definitions.Def_NonconvexSplitting_ProxGrad_StandingAssumptions import Definitions.Def_NonconvexSplitting_ProxGrad_HessianSandwich import Definitions.Def_NonconvexSplitting_ProxGrad_IsProxGradSeq import Definitions.Def_NonconvexSplitting_ProxGrad_IsStationary open NonconvexSplitting.Shared open Filter Topology
namespace NonconvexSplitting.ProxGrad
/-- Theorem 4 of Li–Pong (p. 18). Suppose `q` is twice continuously differentiable and convex,
`ℓ > 0`, and `-ℓ I ⪯ ∇²h(x) + ∇²q(x) ⪯ ℓ I` for all `x` (44). Let `{x^t}` be generated by (43)
with `β ∈ (0, 1/ℓ)`. Then the algorithm is a descent algorithm, and every cluster point of
`{x^t}` (if one exists) is a stationary point. -/
theorem theorem4 {n : ℕ} (h q : EuclideanSpace ℝ (Fin n) → ℝ) (P : EuclideanSpace ℝ (Fin n) → EReal) (ℓ β : ℝ)
(x : ℕ → EuclideanSpace ℝ (Fin n))
(hSA : StandingAssumptions h P) (hq : ContDiff ℝ 2 q) (hq_convex : ConvexOn ℝ Set.univ q)
(hℓ : 0 < ℓ) (h44 : HessianSandwich h q ℓ) (hβ : 0 < β)
(hβℓ : β < 1 / ℓ) (hx : IsProxGradSeq h P β x) :
(∀ t, ((h (x (t + 1)) : ℝ) : EReal) + P (x (t + 1)) ≤ ((h (x t) : ℝ) : EReal) + P (x t)) ∧
∀ xs : EuclideanSpace ℝ (Fin n),
(∃ φ : ℕ → ℕ, StrictMono φ ∧ Tendsto (x ∘ φ) atTop (𝓝 xs)) → IsStationary h P xs := by sorry
end NonconvexSplitting.ProxGrad
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.