Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Remark 4(1) — Theorem 3 when LβL_\betaLβ​ is a KL function

Proved
NonconvexSplitting.ADMMKL.remark4_kl_convergence

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).

Suppose that LβL_\betaLβ​, as a function of (x,y,z)(x,y,z)(x,y,z) on Rn×Rm×Rm\mathbb R^n\times\mathbb R^m\times\mathbb R^mRn×Rm×Rm, is a KL function, and that (xt,yt,zt)(x^t,y^t,z^t)(xt,yt,zt) has a cluster point (x∗,y∗,z∗)(x^*,y^*,z^*)(x∗,y∗,z∗). Then (xt,yt,zt)→(x∗,y∗,z∗)(x^t,y^t,z^t)\to(x^*,y^*,z^*)(xt,yt,zt)→(x∗,y∗,z∗), x∗x^*x∗ is a stationary point of min⁡xh(x)+P(Mx)\min_x h(x)+P(\mathcal Mx)minx​h(x)+P(Mx), and

∑t=0∞∥xt+1−xt∥<∞.\sum_{t=0}^{\infty}\|x^{t+1}-x^t\|<\infty .t=0∑∞​∥xt+1−xt∥<∞.

This is the analytic core of Theorem 3: the semi-algebraic hypothesis is used only to guarantee the KL property 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. Finite length is Summable (fun t => ‖x (t+1) - x t‖); summability does not depend on the first term, so this is the paper's (34) summed from t=1t=1t=1.

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_Shared_IsStationary
import Definitions.Def_NonconvexSplitting_ADMMKL_AugLagProd
import Definitions.Def_NonconvexSplitting_ADMMKL_KLProperty
open NonconvexSplitting.Shared

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

/-- Remark 4(1) of Li–Pong (p. 16): the conclusion of Theorem 3 holds when the semi-algebraicity
of `h` and `P` is replaced by the assumption that the augmented Lagrangian `L_β` (a function of
the triple `(x, y, z)`) is a KL function: the ADMM sequence converges to its cluster point
`(x*, y*, z*)`, `x*` is stationary, and `∑ ‖x^{t+1} - x^t‖ < ∞`. -/
theorem remark4_kl_convergence {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)
    (hKL : IsKLFunction (augLagX h P M β))
    (xs : EuclideanSpace ℝ (Fin n)) (ys zs : EuclideanSpace ℝ (Fin m))
    (hclu : ∃ s : ℕ → ℕ, StrictMono s ∧
      Tendsto (fun i => (x (s i), y (s i), z (s i))) atTop (𝓝 (xs, ys, zs))) :
    Tendsto x atTop (𝓝 xs) ∧ Tendsto y atTop (𝓝 ys) ∧ Tendsto z atTop (𝓝 zs) ∧
      IsStationary h P M xs ∧ Summable (fun t => ‖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. 16, Remark 4(1)
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