Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Eq. (35) — a subgradient of LβL_\betaLβ​ bounded by C∥xt+1−xt∥C\|x^{t+1}-x^t\|C∥xt+1−xt∥

Proved
NonconvexSplitting.ADMMKL.eq35_subgradient_bound

by mikedeng1 · Sep 29, 2026 · Mathlib 0df444a (Lean v4.33.1)

admmkurdyka-lojasiewicznonconvex-optimizationp2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1

Let h:Rn→Rh:\mathbb R^n\to\mathbb Rh:Rn→R be twice continuously differentiable with bounded Hessian, P:Rm→(−∞,+∞]P:\mathbb R^m\to(-\infty,+\infty]P:Rm→(−∞,+∞] proper and closed, M:Rn→Rm\mathcal M:\mathbb R^n\to\mathbb R^mM:Rn→Rm linear and β>0\beta>0β>0, with augmented Lagrangian Lβ(x,y,z)=h(x)+P(y)−⟨z,Mx−y⟩+β2∥Mx−y∥2L_\beta(x,y,z)=h(x)+P(y)-\langle z,\mathcal Mx-y\rangle+\frac\beta2\|\mathcal Mx-y\|^2Lβ​(x,y,z)=h(x)+P(y)−⟨z,Mx−y⟩+2β​∥Mx−y∥2. Suppose Assumption 1 holds with T1=0\mathcal T_1=0T1​=0 (hence ϕ=0\phi=0ϕ=0), with witnesses σ,Q1,Q2,δ,Q3,γ\sigma,\mathcal Q_1,\mathcal Q_2,\delta,\mathcal Q_3,\gammaσ,Q1​,Q2​,δ,Q3​,γ, and let (xt,yt,zt)t≥0(x^t,y^t,z^t)_{t\ge0}(xt,yt,zt)t≥0​ be generated by the ADMM

yt+1∈Arg min⁡yLβ(xt,y,zt),xt+1∈Arg min⁡xLβ(x,yt+1,zt),zt+1=zt−β(Mxt+1−yt+1).y^{t+1}\in\operatorname*{Arg\,min}_y L_\beta(x^t,y,z^t),\qquad x^{t+1}\in\operatorname*{Arg\,min}_x L_\beta(x,y^{t+1},z^t),\qquad z^{t+1}=z^t-\beta(\mathcal Mx^{t+1}-y^{t+1}).yt+1∈yArgmin​Lβ​(xt,y,zt),xt+1∈xArgmin​Lβ​(x,yt+1,zt),zt+1=zt−β(Mxt+1−yt+1).

Then there is a constant C>0C>0C>0 such that for every t≥1t\ge1t≥1 there is w∈∂Lβ(xt+1,yt+1,zt+1)w\in\partial L_\beta(x^{t+1},y^{t+1},z^{t+1})w∈∂Lβ​(xt+1,yt+1,zt+1), the limiting subdifferential of LβL_\betaLβ​ as a function of (x,y,z)(x,y,z)(x,y,z), with

∥w∥≤C ∥xt+1−xt∥.\|w\|\le C\,\|x^{t+1}-x^t\| .∥w∥≤C∥xt+1−xt∥.

In particular dist⁡(0,∂Lβ(xt+1,yt+1,zt+1))≤C∥xt+1−xt∥\operatorname{dist}(0,\partial L_\beta(x^{t+1},y^{t+1},z^{t+1}))\le C\|x^{t+1}-x^t\|dist(0,∂Lβ​(xt+1,yt+1,zt+1))≤C∥xt+1−xt∥, which is (35).

This is the "relative error" half of the KL argument: small primal steps force nearly stationary points of LβL_\betaLβ​.

Formalization Note The ADMM is the relation IsProxADMMSeq with ϕ=0\phi = 0ϕ=0 (the Bregman term vanishes); Assumption 1 is instantiated with ϕ=0\phi = 0ϕ=0 and T1=0\mathcal T_1 = 0T1​=0, and T2\mathcal T_2T2​ is left as a parameter (Assumption 1 forces 0⪯T2⪯T1=00 \preceq \mathcal T_2 \preceq \mathcal T_1 = 00⪯T2​⪯T1​=0). Values of LβL_\betaLβ​ and PPP live in EReal; iterates are indexed from 000, and y0y^0y0 is not constrained. The existential form (an explicit subgradient) implies the distance form. The range t≥1t\ge1t≥1 is where the paper's (14) applies; the paper leaves it implicit.

Preamble
import Mathlib
import Definitions.Def_NonconvexSplitting_Shared_StandingAssumptions
import Definitions.Def_NonconvexSplitting_Shared_Assumption1
import Definitions.Def_NonconvexSplitting_Shared_IsProxADMMSeq
import Definitions.Def_NonconvexSplitting_Shared_LimitingSubdiff
import Definitions.Def_NonconvexSplitting_ADMMKL_AugLagProd
open NonconvexSplitting.Shared

open Filter Topology
open scoped InnerProductSpace
Formal statement
namespace NonconvexSplitting.ADMMKL

/-- Eq. (35) of Li–Pong (p. 14): along the ADMM (proximal ADMM with `φ = 0`) under Assumption 1
with `T₁ = 0`, the limiting subdifferential of `L_β` (on the triple space) at the new iterate
contains a vector of norm at most `C ‖x^{t+1} - x^t‖`, for a constant `C > 0` and all `t ≥ 1`. -/
theorem eq35_subgradient_bound {n m : ℕ} (h : EuclideanSpace ℝ (Fin n) → ℝ) (P : EuclideanSpace ℝ (Fin m) → EReal)
    (M : EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin m)) (β : ℝ)
    (hstd : StandingAssumptions h P (fun _ => 0) β)
    (σ δ γ : ℝ) (Q1 Q2 T2 Q3 : EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n))
    (hA : Assumption1 h (fun _ => 0) M β σ Q1 Q2 0 T2 δ Q3 γ)
    (x : ℕ → EuclideanSpace ℝ (Fin n)) (y z : ℕ → EuclideanSpace ℝ (Fin m))
    (hseq : IsProxADMMSeq h P (fun _ => 0) M β x y z) :
    ∃ C : ℝ, 0 < C ∧ ∀ t : ℕ, 1 ≤ t →
      ∃ w ∈ LimitingSubdiff (augLagX h P M β) (pack (x (t + 1)) (y (t + 1)) (z (t + 1))),
        ‖w‖ ≤ C * ‖x (t + 1) - x t‖ := by sorry

end NonconvexSplitting.ADMMKL
Source
Li & Pong, Global Convergence of Splitting Methods for Nonconvex Composite Optimization, arXiv:1407.0753v6, p. 14, Eq. (35)
Human review
  • Endorsed by Shuze Chen · Oct 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by mikedeng1 · Oct 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me