The whole method assumes a coarsely (or finely) central x̃(λ₀) to start from. This closing section sketches how to manufacture one — precision calculations are deliberately omitted and deferred to the notes 1:13:18, 1:17:13.
Step 1 — Modify the LP so an interior point is obvious
-
minimize cᵀx + N·z subject to Ax + 𝟙·z ≥ b 0 ≤ z ≤ 2^{L+1} −2^{L+1}·𝟙 ≤ x ≤ 2^{L+1}·𝟙 -
Claim. If
Nis exponential in the precision (e.g.10^Lor100^L), then whenever the original LP is bounded and feasible, every optimum of LP′ hasz = 0— so solving LP′ solves the original LP1:18:40,1:19:20. -
Proved via precision arguments; a reference is promised in the notes. Separately, any solution of the original LP has bounded precision anyway
1:19:20,1:19:36. -
LP′ has an obvious interior point:
x = 0,z = ‖b‖_∞1:19:36,1:20:02. -
But an interior point is not enough — it must be central for some known
λ, and this one isn't known to be1:20:02,1:20:35.
Step 2 — The trick: change the cost function instead
- Set
λ = 1and look at∇f_λ(x′) = λc − AᵀS(x′)⁻¹𝟙. Optimality means zero gradient1:20:35,1:21:25,1:21:55. - The trick:
x′is perfectly central for a different cost function, namelyc := AᵀS(x′)⁻¹𝟙, which makes the gradient identically zero1:21:55,1:22:23.
Step 3 — Run the algorithm backwards, then swap the cost
- With a perfectly central point for some cost function in hand, run all the update steps in reverse: decrement
λa little at a time and run Newton's method, untilλis exponentially small in the precision1:22:23,1:22:48. - The result is a point that is finely central near the analytic center
1:23:04. - Key property: once
λis that tiny, the cost function barely matters — so you can swap in the real costcand remain finely central for it (for any cost function of bounded precision)1:23:04,1:23:27. - Then run everything forward with the true cost function
1:23:27,1:24:18.
graph TD
A["LP′ with slack variable z<br/>and box constraints"] --> B["obvious interior point<br/>x = 0, z = ‖b‖∞"]
B --> C["set λ = 1; declare cost<br/>c := AᵀS(x′)⁻¹𝟙"]
C --> D["perfectly central for<br/>that fake cost"]
D --> E["run IPM in REVERSE:<br/>decrease λ + Newton steps"]
E --> F["finely central near<br/>the analytic center"]
F --> G["swap in the real cost c<br/>(λ tiny ⇒ cost irrelevant)"]
G --> H["run IPM FORWARD:<br/>increase λ by 1+1/√m"]