Optimality relations (11) of the proximal ADMM iterates
ProvedNonconvexSplitting.ProxADMM.eq11_optimalityadmmnonconvex-optimizationp2o-batch-p200bp2o-gran-per-chapterp2o-plan-paperp2o-v1
Assume is with bounded Hessian, is proper and lower semicontinuous, is and convex, and . Let be generated by the proximal ADMM. Then for every
where is the limiting subdifferential.
These relations rewrite the first-order optimality conditions of the two subproblems in terms of the updated multiplier; Assumption 1 is not needed.
Preamble
import Mathlib import Definitions.Def_NonconvexSplitting_Shared_LimitingSubdiff import Definitions.Def_NonconvexSplitting_Shared_StandingAssumptions import Definitions.Def_NonconvexSplitting_Shared_IsProxADMMSeq open NonconvexSplitting.Shared open Filter Topology open scoped InnerProductSpace
Formal statement
namespace NonconvexSplitting.ProxADMM
/-- Eq. (11) of Li–Pong (p. 6): the relations satisfied by every proximal ADMM iterate. -/
theorem eq11_optimality {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 β)
(x : ℕ → EuclideanSpace ℝ (Fin n)) (y z : ℕ → EuclideanSpace ℝ (Fin m))
(hseq : IsProxADMMSeq h P phi M β x y z) (t : ℕ) :
-(z (t + 1)) - β • M (x (t + 1) - x t) ∈ LimitingSubdiff P (y (t + 1)) ∧
gradient h (x (t + 1)) - ContinuousLinearMap.adjoint M (z (t + 1)) =
-(gradient phi (x (t + 1))) + gradient phi (x t) ∧
M (x (t + 1)) - y (t + 1) = (1 / β) • (z t - z (t + 1)) := by sorry
end NonconvexSplitting.ProxADMM
Source
Li & Pong, Global Convergence of Splitting Methods for Nonconvex Composite Optimization, arXiv:1407.0753v6, p. 6, Eq. (11) (derived from (8))
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.