Theorem 7.2, p. 31 — Semi-APGM (Algorithm 2) has ℒ_{k+1} ≤ ℒ_k/(1 + α_k), (87) and (88)
ProvedNAGFlow.APGM.theorem_7_2Let be a real Hilbert space and consider (110), where with and is proper, closed and convex. Let be a minimiser of . Let be generated by Algorithm 2 (Semi-APGM) with , , and , and let
Then
for all
and for all
The theorem gives an accelerated proximal gradient method with one proximal evaluation per iteration whose rate covers the convex case (, rate ) and the strongly convex case (, linear rate ) in one statement, exactly as Theorem 5.1 does for smooth .
Formalization Note. is encoded by its domain and its real values on ; values off are never used. Hypotheses made explicit: (the page divides by and takes ), and (Algorithm 2's input and step 2), a minimiser of over (argmin , p. 2), and , which is added: if then and (87)–(88) are void on the page. The page's "(110) " is read as , as in (104) and Lemma 7.1, since is nonsmooth. Powers with negative exponent are integer powers of positive reals; is stated for , as on the page.
import Mathlib import Definitions.Def_NAGFlow_APGM_Setting
namespace NAGFlow.APGM
/-- Theorem 7.2 (Luo & Chen, arXiv:1909.03145v4, p. 31). Let `h ∈ S^{1,1}_{μ,L}` with
`0 ≤ μ ≤ L < ∞`, let `g : V → ℝ ∪ {+∞}` be proper, closed and convex (encoded by `(g, D)`,
`D = dom g`), `f = h + g`, and let `x* ∈ dom g` be a minimiser of `f`. For every run
`(α, γ, x, y, w, v)` of Algorithm 2 (Semi-APGM) started at `x₀ ∈ dom g`, with
`ℒ_k = f(x_k) − f(x*) + (γ_k/2)‖v_k − x*‖²`:
(116) `ℒ_{k+1} ≤ ℒ_k/(1 + α_k)` for every `k ∈ ℕ`;
(87) `ℒ_k ≤ ℒ₀ · min{4L/(√γ₀ k + 2√L)², (1 + √(min{γ₀, μ}/L))^{−k}}` for every `k ≥ 0`;
(88) `ℒ_k ≤ C_{γ₀,L} · min{4/k², (1 + √(min{γ₀, μ}/L))^{1−k}}` for every `k ≥ 1`, where
(89) `C_{γ₀,L} = (L/γ₀)(f(x₀) − f(x*)) + (L/2)‖v₀ − x*‖²`. -/
theorem theorem_7_2 {V : Type*} [NormedAddCommGroup V] [InnerProductSpace ℝ V] [CompleteSpace V]
(h : V → ℝ) (gradh : V → V) (μ L : ℝ) (hh : NAGFlow.PredCorr.IsS11 h gradh μ L)
(g : V → ℝ) (D : Set V) (hg : IsProperClosedConvex g D)
(xstar : V) (hxstar : xstar ∈ D) (hmin : ∀ z ∈ D, (h + g) xstar ≤ (h + g) z)
(α γ : ℕ → ℝ) (x y w v : ℕ → V) (hrun : IsAPGMRun gradh g D μ L α γ x y w v)
(hx0 : x 0 ∈ D) :
(∀ k : ℕ, NAGFlow.PredCorr.lyap (h + g) xstar x v γ (k + 1) ≤ NAGFlow.PredCorr.lyap (h + g) xstar x v γ k / (1 + α k)) ∧
(∀ k : ℕ, NAGFlow.PredCorr.lyap (h + g) xstar x v γ k ≤
NAGFlow.PredCorr.lyap (h + g) xstar x v γ 0 *
min (4 * L / (Real.sqrt (γ 0) * k + 2 * Real.sqrt L) ^ 2)
((1 + Real.sqrt (min (γ 0) μ / L)) ^ (-(k : ℤ)))) ∧
(∀ k : ℕ, 1 ≤ k → NAGFlow.PredCorr.lyap (h + g) xstar x v γ k ≤
(L / γ 0 * ((h + g) (x 0) - (h + g) xstar) + L / 2 * ‖v 0 - xstar‖ ^ 2) *
min (4 / (k : ℝ) ^ 2) ((1 + Real.sqrt (min (γ 0) μ / L)) ^ (1 - (k : ℤ)))) := by sorry
end NAGFlow.APGM
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.