Theorem 3.11 — along the -algorithm, is infinitely often at least
OpenShorNonsmooth.RAlgorithm.pRatio_frequently_geLet , let satisfy
let and , and let be a sequence constructed by the -algorithm applied to (with any admissible choices of almost-gradients and stepsizes) such that
Then for every fixed with , every , and every positive integer there exists such that
In words: infinitely often, the local set of almost-gradients around the iterate cannot be both thin and far from the origin. This is the quantitative core from which the convergence results of Section 3.7 are derived.
Formalization Note Condition (3.50) is the standing assumption of the section ("from now on we shall assume", p. 79) and the proof uses the boundedness of it implies, so it is a hypothesis here although the theorem's sentence does not repeat it. takes values in ; the right-hand side is embedded with ENNReal.ofReal, and is the real power .
import Mathlib import Definitions.Def_ShorNonsmooth_RAlgorithm_RAlgorithm open scoped InnerProductSpace ENNReal open Filter Topology
namespace ShorNonsmooth.RAlgorithm
/-- Shor (1985), p. 82, **Theorem 3.11**. Let `f ∈ K` (formed from the data `P`), satisfying the
standing condition (3.50) `f(x) → +∞` as `‖x‖ → ∞` (assumed "from now on", p. 79), let `α > 1`,
and let `{x_k}` be constructed by the `r(α)`-algorithm (the `r_μ(α)`-algorithm with `μ = 0`)
applied to `f`, with `lim_{k→∞} ‖x_{k+1} - x_k‖ = 0` (3.52). Then for each fixed `v` with
`ⁿ√β < v < 1` (`β = 1/α`), `ε > 0`, `δ > 0` and positive integer `r` there is `k̄ > r` with
`p(P̄_{δ,ε}(x_k̄)) ≥ √((v² ⁿ√(α²) - 1)/(α² - 1))`. -/
theorem pRatio_frequently_ge {n : ℕ} (hn : 0 < n) (P : KRep n)
(f : EuclideanSpace ℝ (Fin n) → ℝ) (hf : P.Forms f)
(hf_coercive : Tendsto f (cocompact (EuclideanSpace ℝ (Fin n))) atTop)
(α : ℝ) (hα : 1 < α)
(x gt g : ℕ → EuclideanSpace ℝ (Fin n))
(B : ℕ → EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n)) (h : ℕ → ℝ)
(hrun : IsRun P f α 0 x gt g B h)
(hstep : Tendsto (fun k => ‖x (k + 1) - x k‖) atTop (𝓝 0))
(v : ℝ) (hv_gt : (1 / α) ^ (1 / (n : ℝ)) < v) (hv_lt : v < 1)
(ε : ℝ) (hε : 0 < ε) (δ : ℝ) (hδ : 0 < δ) (r : ℕ) (hr : 0 < r) :
∃ kbar : ℕ, r < kbar ∧
ENNReal.ofReal (Real.sqrt ((v ^ 2 * (α ^ 2) ^ (1 / (n : ℝ)) - 1) / (α ^ 2 - 1))) ≤
pRatio (P.Pbar δ ε (x kbar)) := by sorry
end ShorNonsmooth.RAlgorithm
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.