Theorem 3.3 — the SDG iterates satisfy
OpenShorNonsmooth.SpaceDilation.dilated_distance_invariantLet , let and , and let . Let satisfy, for every ,
where . Run the SDG method with and
- ;
- , (3.19)
- , (3.20)
Then
The distance from the iterate to , measured in the transformed space, never exceeds its initial value; in particular all iterates stay in . This invariant is what lets Theorems 3.1 and 3.2 be applied to function values in Theorem 3.4.
Formalization Note The book takes almost differentiable on , its almost-gradient and a local minimum point; the proof uses only (3.18), so the Lean statement holds for every on and every satisfying (3.18) on (a generalization). Condition (3.18) with already forces on . The proof's is assumed (); for another the claim fails at . The stepsize (3.19) is used only while , when ; if the method stops and the state is repeated.
import Mathlib import Definitions.Def_ShorNonsmooth_SpaceDilation_SDGMethod
namespace ShorNonsmooth.SpaceDilation
/-- Shor (1985), pp. 56–57, Theorem 3.3. Let `x*` be a point, `d > 0`, `S_d = {x : ‖x - x*‖ ≤ d}`,
and let the selection `g` satisfy on `S_d`
`N [f(x) - f(x*)] ≤ (g(x), x - x*) ≤ M [f(x) - f(x*)]` (3.18) with `M > N > 0`.
Run the SDG method with `B₀ = I`, `x₀ ∈ S_d`, stepsizes
`h_{k+1} = (2MN/(M+N)) (f(x_k) - f(x*)) / ‖g̃_k‖` (3.19) and coefficients
`1 < α_{k+1} = α ≤ (M+N)/(M-N)` (3.20). Then `‖A_k (x_k - x*)‖ ≤ d` for `k = 0, 1, 2, …`.
(The book takes `f` almost differentiable, `g` its almost-gradient and `x*` a local minimum;
the proof uses only (3.18), which already forces `f ≥ f(x*)` on `S_d`.) -/
theorem dilated_distance_invariant {n : ℕ}
(f : EuclideanSpace ℝ (Fin n) → ℝ)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(xstar x₀ : EuclideanSpace ℝ (Fin n)) (d M N α : ℝ)
(hd : 0 < d) (hN : 0 < N) (hNM : N < M)
(h318 : ∀ x ∈ Metric.closedBall xstar d,
N * (f x - f xstar) ≤ inner ℝ (g x) (x - xstar) ∧
inner ℝ (g x) (x - xstar) ≤ M * (f x - f xstar))
(hx₀ : x₀ ∈ Metric.closedBall xstar d)
(hα : 1 < α) (hαMN : α ≤ (M + N) / (M - N)) (k : ℕ) :
‖(sdg g (fun _ x gt => 2 * M * N / (M + N) * (f x - f xstar) / ‖gt‖) (fun _ => α) x₀
(ContinuousLinearEquiv.refl ℝ _) k).A
((sdg g (fun _ x gt => 2 * M * N / (M + N) * (f x - f xstar) / ‖gt‖) (fun _ => α) x₀
(ContinuousLinearEquiv.refl ℝ _) k).x - xstar)‖ ≤ d := by sorry
end ShorNonsmooth.SpaceDilation
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.