Remark 4(1) — Theorem 3 when is a KL function
ProvedNonconvexSplitting.ADMMKL.remark4_kl_convergenceLet be twice continuously differentiable with bounded Hessian, proper and closed, linear and , with augmented Lagrangian . Suppose Assumption 1 holds with (hence ), with witnesses , and let be generated by the ADMM
Suppose that , as a function of on , is a KL function, and that has a cluster point . Then , is a stationary point of , and
This is the analytic core of Theorem 3: the semi-algebraic hypothesis is used only to guarantee the KL property of .
Formalization Note The ADMM is the relation IsProxADMMSeq with (the Bregman term vanishes); Assumption 1 is instantiated with and , and is left as a parameter (Assumption 1 forces ). Values of and live in EReal; iterates are indexed from , and 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 .
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
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
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.