Theorem 3.1 — along a subsequence
OpenShorNonsmooth.SpaceDilation.gTilde_subsequence_boundRun the SDG method in () from a starting point with a nonsingular initial operator , an arbitrary stepsize rule and space-dilation coefficients ; write for the transformed gradients. Suppose there are positive numbers , , with
- for ;
- for .
Then there exist a constant and indices such that
Since , the transformed gradients become small at a geometric rate along a subsequence. This is the first step towards the convergence rate of function values in Theorem 3.4.
Formalization Note The book states the theorem for an almost differentiable with its almost-gradient; its proof uses only the bound , so the Lean statement holds for every map and every stepsize rule, which is a generalization, not a weakening. The constant and the subsequence depend on the run and are chosen after all the data. If the method stops (), the state is repeated and from then on.
import Mathlib import Definitions.Def_ShorNonsmooth_SpaceDilation_SDGMethod
namespace ShorNonsmooth.SpaceDilation
/-- Shor (1985), p. 53, Theorem 3.1. Let the SDG method be run with a nonsingular initial operator
`B₀`, any stepsize rule `h` and any generalized-gradient selection `g`. If for positive `d`, `α*`,
`δ` one has `‖g(x_k)‖ ≤ d` for `k = 0, 1, …` and `1 + δ ≤ α_k ≤ α*` for `k = 1, 2, …`, then there
are a constant `c > 0` and a strictly increasing sequence of indices `k_p` with
`‖g̃_{k_p}‖ < c (∏_{j=1}^{k_p} α_j)^{-1/n}` for every `p`.
(The book states it for an almost differentiable `f` with `g` an almost-gradient; the proof uses
only the bound `‖g(x_k)‖ ≤ d`, so the theorem is stated for every selection `g`.) -/
theorem gTilde_subsequence_bound {n : ℕ} (hn : 0 < n)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n))
(h : ℕ → EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n) → ℝ) (α : ℕ → ℝ)
(x₀ : EuclideanSpace ℝ (Fin n))
(B₀ : EuclideanSpace ℝ (Fin n) ≃L[ℝ] EuclideanSpace ℝ (Fin n))
(d αstar δ : ℝ) (hd : 0 < d) (hαstar : 0 < αstar) (hδ : 0 < δ)
(hg : ∀ k : ℕ, ‖g (sdg g h α x₀ B₀ k).x‖ ≤ d)
(hα : ∀ k : ℕ, 1 ≤ k → 1 + δ ≤ α k ∧ α k ≤ αstar) :
∃ c : ℝ, 0 < c ∧ ∃ kp : ℕ → ℕ, StrictMono kp ∧
∀ p : ℕ, ‖gTilde g h α x₀ B₀ (kp p)‖ <
c * (∏ j ∈ Finset.Icc 1 (kp p), α j) ^ (-(1 : ℝ) / n) := by sorry
end ShorNonsmooth.SpaceDilation
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.