Eq. (35) — a subgradient of bounded by
ProvedNonconvexSplitting.ADMMKL.eq35_subgradient_boundLet 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
Then there is a constant such that for every there is , the limiting subdifferential of as a function of , with
In particular , which is (35).
This is the "relative error" half of the KL argument: small primal steps force nearly stationary points 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. The existential form (an explicit subgradient) implies the distance form. The range is where the paper's (14) applies; the paper leaves it implicit.
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
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
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.