p. 87 — the localizing ellipsoids shrink in volume by the ratio per iteration
OpenShorNonsmooth.Ellipsoid.ellipsoid_volume_ratioLet , , any vector field and , and run the algorithm (3.57)–(3.60). Suppose the -st iteration is performed, i.e. . With and for arbitrary centers , the volume is positive and finite, and
Thus the volume of the ellipsoid that localizes the solution according to (3.61) decreases geometrically with ratio , which depends only on the dimension ( for large ).
Formalization Note The chain of equalities printed on p. 87 carries the exponents and on where is meant, and writes where is meant; the statement uses the value of printed on p. 88, which is what the book's own derivation gives (). The volumes are compared in ℝ≥0∞; positivity and finiteness of are stated so that the identity is the book's ratio.
import Mathlib import Definitions.Def_ShorNonsmooth_Ellipsoid_EllipsoidMethod open MeasureTheory
namespace ShorNonsmooth.Ellipsoid
/-- Shor (1985), pp. 87–88, the volume ratio after the proof of Theorem 3.14, corrected. Let
`n > 1`, `R > 0`, and let the `(k+1)`-st iteration of (3.57)–(3.60) be performed (`g(x_k) ≠ 0`).
With `Φ_k = {x : ‖A_k (x - c)‖ ≤ (n + 1) h_k}`, `A_k = B_k⁻¹`, the volume of `Φ_k` is positive and
finite and, for any centers `c, c'`,
`v(Φ_{k+1}) = q_n v(Φ_k)` with `q_n = √((n - 1)/(n + 1)) (n/√(n² - 1))ⁿ < 1` (p. 88).
The printed chain on p. 87 carries the exponents `2` and `1` on `n/√(n² - 1)` and drops the
square root on `(n - 1)/(n + 1)`; the statement follows the value of `q_n` on p. 88. -/
theorem ellipsoid_volume_ratio {n : ℕ} (hn : 1 < n)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n)) (R : ℝ) (hR : 0 < R)
(x₀ : EuclideanSpace ℝ (Fin n)) (k : ℕ)
(hstep : g (ellipsoidMethod g R x₀ k).x ≠ 0) :
qRatio n < 1 ∧
∀ c c' : EuclideanSpace ℝ (Fin n),
0 < volume (ellipsoid (ellipsoidMethod g R x₀ k).B⁻¹ c
(((n : ℝ) + 1) * (ellipsoidMethod g R x₀ k).h)) ∧
volume (ellipsoid (ellipsoidMethod g R x₀ k).B⁻¹ c
(((n : ℝ) + 1) * (ellipsoidMethod g R x₀ k).h)) < ⊤ ∧
volume (ellipsoid (ellipsoidMethod g R x₀ (k + 1)).B⁻¹ c'
(((n : ℝ) + 1) * (ellipsoidMethod g R x₀ (k + 1)).h)) =
ENNReal.ofReal (qRatio n) *
volume (ellipsoid (ellipsoidMethod g R x₀ k).B⁻¹ c
(((n : ℝ) + 1) * (ellipsoidMethod g R x₀ k).h)) := by sorry
end ShorNonsmooth.Ellipsoid
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.