Eq. (2.5) —
ProvedCalamaiMore.QP.projection_step_descentdescentgradient-projectionp2o-batch-p200ap2o-gran-per-chapterp2o-plan-paperp2o-v1
Let be a nonempty closed convex subset of a finite-dimensional real inner product space , and let be continuously differentiable on . Let , let , and let be the gradient projection step. Then
The estimate says that the step direction is a descent direction whenever it is nonzero; combined with the sufficient decrease condition (2.1) it shows that a gradient projection step strictly decreases unless .
Formalization Note "Continuously differentiable on " is stated as differentiability at every point of together with continuity of gradient f on , the standing assumption of Section 2 of the paper.
Preamble
import Mathlib import Definitions.Def_CalamaiMore_QP_proj
Formal statement
namespace CalamaiMore.QP
/-- Calamai–Moré, Eq. (2.5) (p. 98): for a nonempty closed convex `Ω`, an iterate `x_k ∈ Ω`, a
step `α_k > 0` and `x_{k+1} = P(x_k - α_k ∇f(x_k))`,
`⟨∇f(x_k), x_k - x_{k+1}⟩ ≥ ‖x_{k+1} - x_k‖² / α_k`. -/
theorem projection_step_descent {E : Type*} [NormedAddCommGroup E] [InnerProductSpace ℝ E] [FiniteDimensional ℝ E]
(Ω : Set E) (hΩne : Ω.Nonempty) (hΩc : IsClosed Ω) (hΩcv : Convex ℝ Ω)
(f : E → ℝ) (hfd : ∀ y ∈ Ω, DifferentiableAt ℝ f y) (hfc : ContinuousOn (gradient f) Ω)
(xk xnext : E) (hxk : xk ∈ Ω) (αk : ℝ) (hαk : 0 < αk)
(hstep : xnext = projPath f Ω xk αk) :
‖xnext - xk‖ ^ 2 / αk ≤ inner ℝ (gradient f xk) (xk - xnext) := by sorry
end CalamaiMore.QP
Source
Calamai & Moré, Projected gradient methods for linearly constrained problems, Math. Programming 39 (1987), p. 98, Eq. (2.5)
Human review
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.