Löwner–John rounding for polytopes
ProvedConvexOptimization.lowner_john_polytope_roundingLöwner–John rounding for polytopes. Shrinking the minimum-volume covering ellipsoid of a polytope about its centre by the factor lands inside the polytope.
Let , let and let . Let be a Löwner–John pair for the points: symmetric positive definite, , and maximal among all such covering pairs. Write for the centre of . Then
that is, every with — equivalently — belongs to . Combined with the covering property that is part of the hypothesis, the polytope is sandwiched between two concentric homothetic ellipsoids whose ratio is the dimension .
This is the quantitative statement that a polytope can be rounded: after the affine change of coordinates taking to the unit ball, lies between the balls of radius and . Such a sandwich is what makes ellipsoid-method volume arguments and Banach–Mazur distance estimates work. The factor cannot be improved — it is attained by the simplex (Boyd & Vandenberghe, exercise 8.13) — and for general convex bodies, rather than polytopes, the same bound holds by an approximation argument.
Formalization Note The conclusion is stated in the scaled-preimage form , using the dot product ⬝ᵥ on Fin n → ℝ, so that no inverse or explicit centre needs to appear; the covering hypothesis is imposed on the range of the family x : Fin m → (Fin n → ℝ), which is equivalent to imposing it on because ellipsoids are convex. The dimension nn is assumed positive. Source: Boyd & Vandenberghe §8.4.1, p. 412.
import Mathlib import Definitions.Def_ConvexOptimization_ellipsoidBody import Definitions.Def_ConvexOptimization_IsLownerJohn open scoped RealInnerProductSpace ENNReal open MeasureTheory
theorem ConvexOptimization.lowner_john_polytope_rounding {nn m : ℕ} (hnn : 0 < nn)
(x : Fin m → Fin nn → ℝ) (A : Matrix (Fin nn) (Fin nn) ℝ) (b : Fin nn → ℝ)
(hopt : IsLownerJohn A b (Set.range x)) (v : Fin nn → ℝ)
(hv : (A.mulVec v + b) ⬝ᵥ (A.mulVec v + b) ≤ 1 / (nn : ℝ) ^ 2) :
v ∈ convexHull ℝ (Set.range x) := by
sorry
Read-back
What the Lean code literally says, in plain math · claude-fable-5
Theorem. Fix with ( excluded, so is a genuine positive real), points , a matrix and a vector . Assume satisfies the Löwner–John predicate for — unfolded: is symmetric; is positive definite (Mathlib's notion, which includes symmetry; in particular is invertible); every satisfies ; and for every symmetric positive-definite with any offset whose body contains all the . Then every vector with — the same quadratic form as the ellipsoid body but at level instead of , i.e. , which (since is invertible under the hypothesis) is the body concentrically shrunk by the factor about its center — belongs to the convex hull over of . Edge cases: the shrunk region always contains the center itself; there is no separate hypothesis that the hull has nonempty interior or that — for point families for which the Löwner–John predicate is unsatisfiable (e.g. , or all points in a proper affine subspace, with ), the statement is vacuous.
Confirmed by the mission captain (proposal self-audit).