Theorem 2: the proximal ADMM sequence is bounded under coercivity
ProvedNonconvexSplitting.ADMMBounded.theorem2Let be twice continuously differentiable with bounded Hessian, proper and lower semicontinuous, twice continuously differentiable and convex, linear with adjoint , and . Let be a sequence generated by the proximal ADMM with penalty and proximal function (initial points arbitrary, unconstrained). Suppose Assumption 1 holds with witnesses and . For a self-map write . Suppose moreover that condition (29) holds with the same and : there are with and such that
Suppose that either
- is invertible and ; or
- and .
Then the sequence is bounded: there is with
Boundedness guarantees that the sequence has a cluster point, which is the hypothesis under which the convergence results for the proximal ADMM (stationarity of cluster points, and whole-sequence convergence under the Kurdyka–Łojasiewicz property) apply.
Formalization Note The paper defines as the infimum and requires it to be ; here is any real lower bound. The infimum is finite exactly when a real lower bound exists, and the proof uses only as a lower bound, so the two readings give the same result. "Invertible" is bijectivity of as a map (the dimensions are kept separate). is encoded as: for every there is with whenever ; likewise for . is: some real satisfies for all . The two cases are one disjunctive hypothesis of a single statement. The bound covers all three blocks and all ; is a single (arbitrary) point and does not affect boundedness.
import Mathlib import Definitions.Def_NonconvexSplitting_Shared_StandingAssumptions import Definitions.Def_NonconvexSplitting_Shared_IsProxADMMSeq import Definitions.Def_NonconvexSplitting_Shared_Assumption1 open NonconvexSplitting.Shared open Filter Topology open scoped InnerProductSpace
namespace NonconvexSplitting.ADMMBounded
/-- Theorem 2 of Li–Pong (p. 11): boundedness of the proximal ADMM sequence.
Under the standing assumptions, Assumption 1 and condition (29) (with the same `σ`, `γ`),
if either (i) `M` is invertible and `P` is coercive, or (ii) `h` is coercive and `P` is
bounded below, then every proximal ADMM sequence `(x^t, y^t, z^t)` is bounded. -/
theorem theorem2 {n m : ℕ} (h : EuclideanSpace ℝ (Fin n) → ℝ) (P : EuclideanSpace ℝ (Fin m) → EReal)
(phi : EuclideanSpace ℝ (Fin n) → ℝ) (M : EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin m)) (β : ℝ)
(hstd : StandingAssumptions h P phi β)
(σ δ γ : ℝ) (Q1 Q2 T1 T2 Q3 : EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n))
(hA : Assumption1 h phi M β σ Q1 Q2 T1 T2 δ Q3 γ)
(ζ h₀ : ℝ) (hζ : 0 < ζ ∧ ζ < 2 * β * γ)
(h29 : ∀ x, h₀ ≤ h x - 1 / (σ * ζ) * ‖gradient h x‖ ^ 2)
(hcase :
(Function.Bijective M ∧
∀ r : ℝ, ∃ R : ℝ, ∀ w : EuclideanSpace ℝ (Fin m), R < ‖w‖ → (r : EReal) < P w) ∨
((∀ r : ℝ, ∃ R : ℝ, ∀ v : EuclideanSpace ℝ (Fin n), R < ‖v‖ → r < h v) ∧
∃ c : ℝ, ∀ w : EuclideanSpace ℝ (Fin m), (c : EReal) ≤ P w))
(x : ℕ → EuclideanSpace ℝ (Fin n)) (y z : ℕ → EuclideanSpace ℝ (Fin m))
(hseq : IsProxADMMSeq h P phi M β x y z) :
∃ R : ℝ, ∀ t : ℕ, ‖x t‖ ≤ R ∧ ‖y t‖ ≤ R ∧ ‖z t‖ ≤ R := by sorry
end NonconvexSplitting.ADMMBounded
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.