p. 87 — the localizing ellipsoid has volume
ProvedShorNonsmooth.Ellipsoid.ellipsoid_volume_formulaLet , , let be any vector field and , and run the algorithm (3.57)–(3.60) from , , . Suppose that the first iterations are performed, i.e. for , and put . Then
and for every center the ellipsoid has Lebesgue volume
where is the volume of the closed unit ball of .
This is the volume computation that turns the localization inequality (3.61) of Theorem 3.14 into a rate: the solution lies in an ellipsoid whose volume is known explicitly.
Formalization Note The book centers at ; since Lebesgue measure is translation invariant, the statement is given for an arbitrary center, which covers both and the iterate . The volume is an element of ℝ≥0∞ and the right-hand side is volume (closedBall 0 1) * ENNReal.ofReal (…).
import Mathlib import Definitions.Def_ShorNonsmooth_Ellipsoid_EllipsoidMethod open MeasureTheory
namespace ShorNonsmooth.Ellipsoid
/-- Shor (1985), p. 87, display after the proof of Theorem 3.14. Let `n > 1`, `R > 0`, and let
the algorithm (3.57)–(3.60) run for `k` iterations without stopping (`g(x_j) ≠ 0` for `j < k`).
Then `(n + 1) h_k = R (n/√(n² - 1))^k`, and for every center `c` the ellipsoid
`Φ_k = {x : ‖A_k (x - c)‖ ≤ (n + 1) h_k}`, `A_k = B_k⁻¹`, has volume
`v₀ Rⁿ (n/√(n² - 1))^{nk} / det A_k`, where `v₀` is the volume of the unit ball.
(The book centers `Φ_k` at `x*`; the volume does not depend on the center.) -/
theorem ellipsoid_volume_formula {n : ℕ} (hn : 1 < n)
(g : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n)) (R : ℝ) (hR : 0 < R)
(x₀ : EuclideanSpace ℝ (Fin n)) (k : ℕ)
(hrun : ∀ j < k, g (ellipsoidMethod g R x₀ j).x ≠ 0) :
((n : ℝ) + 1) * (ellipsoidMethod g R x₀ k).h = R * ratio n ^ k ∧
∀ c : EuclideanSpace ℝ (Fin n),
volume (ellipsoid (ellipsoidMethod g R x₀ k).B⁻¹ c
(((n : ℝ) + 1) * (ellipsoidMethod g R x₀ k).h)) =
volume (Metric.closedBall (0 : EuclideanSpace ℝ (Fin n)) 1) *
ENNReal.ofReal (R ^ n * ratio n ^ (n * k) / ((ellipsoidMethod g R x₀ k).B⁻¹).det) := by sorry
end ShorNonsmooth.Ellipsoid
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.