Corollary of Theorem 4.1 — formula (4.7) for functions differentiable in
OpenShorNonsmooth.Decomposition.subgradient_formula_differentiableAssume the hypotheses of Theorem 4.1: and , , are jointly convex; is a convex set of -values at which the minimum in (4.5) is attained; and the Slater condition holds for (4.4) at ; is optimal and are Kuhn–Tucker multipliers. Suppose in addition that each , , is continuously differentiable in . Let be arbitrary subgradients of at . Then
is a subgradient of at on : for all .
Formula (4.7) lets one compute a subgradient of from subgradients of the problem functions, without choosing a special subgradient of ; it underlies the decomposition algorithm of p. 96 for linear and quadratic subproblems.
Formalization Note "Continuously differentiable with respect to " is ContDiff ℝ 1 of for every .
import Mathlib import Definitions.Def_ShorNonsmooth_Decomposition_ValueFunction
namespace ShorNonsmooth.Decomposition
/-- Shor (1985), **Corollary** of Theorem 4.1 (p. 95), formula (4.7): under the hypotheses of Theorem 4.1
(jointly convex `f_α`, a convex set `W` on which the minimum (4.5) is attained, `xbar ∈ W` with the Slater
condition, an optimal `ybar` and Kuhn–Tucker multipliers `U`), if in addition every `f_α(x, ·)`,
`α = 0, 1, …, n`, is continuously differentiable in `y`, then for **arbitrary** subgradients
`(gx₀, gy₀)` of `f₀` and `(gx_i, gy_i)` of `f_i` at `(xbar, ybar)`, the vector
`gx₀ + Σ_i U_i gx_i` is a subgradient of `Φ` at `xbar` on `W`. -/
theorem subgradient_formula_differentiable {l m n : ℕ}
(f₀ : EuclideanSpace ℝ (Fin l) → EuclideanSpace ℝ (Fin m) → ℝ)
(f : Fin n → EuclideanSpace ℝ (Fin l) → EuclideanSpace ℝ (Fin m) → ℝ)
(hf₀ : JointlyConvex f₀) (hf : ∀ i, JointlyConvex (f i))
(hd₀ : ∀ x, ContDiff ℝ 1 (f₀ x)) (hd : ∀ i x, ContDiff ℝ 1 (f i x))
(W : Set (EuclideanSpace ℝ (Fin l))) (hW : Convex ℝ W)
(hWmin : ∀ x ∈ W, MinAttained f₀ f x)
(xbar : EuclideanSpace ℝ (Fin l)) (hxbar : xbar ∈ W) (hslater : SlaterAt f xbar)
(ybar : EuclideanSpace ℝ (Fin m)) (hybar : IsOptimalY f₀ f xbar ybar)
(U : Fin n → ℝ) (hU : IsKuhnTuckerMultiplier f₀ f xbar ybar U)
(gx₀ : EuclideanSpace ℝ (Fin l)) (gy₀ : EuclideanSpace ℝ (Fin m))
(hg₀ : IsJointSubgradient f₀ xbar ybar gx₀ gy₀)
(gx : Fin n → EuclideanSpace ℝ (Fin l)) (gy : Fin n → EuclideanSpace ℝ (Fin m))
(hg : ∀ i, IsJointSubgradient (f i) xbar ybar (gx i) (gy i)) :
IsSubgradientOn (valueFn f₀ f) W xbar (gx₀ + ∑ i, U i • gx i) := by sorry
end ShorNonsmooth.Decomposition
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.