The lecturer's own summary: "conceptually there's not really a whole lot going on here, but it's just a lot of symbols" 25:32.
Step 1 — express the new gradient exactly
-
By the fundamental theorem of calculus (equivalently, a zeroth-order Taylor expansion with explicit error)
18:11,18:46:∇f(x₁) = ∇f(x₀) + ∫₀¹ ∇²f(x_t)(x₁ − x₀) dt -
Substituting the Newton step
x₁ − x₀ = −(∇²f(x₀))⁻¹∇f(x₀), the leading∇f(x₀)cancels and one is left with an integral of the Hessian difference19:03,20:06:∇f(x₁) = ∫₀¹ [∇²f(x₀) − ∇²f(x_t)] (∇²f(x₀))⁻¹ ∇f(x₀) dt -
A student catches a sign error in the expansion; the lecturer confirms the middle term is subtracted and that the remaining piece is by definition
x₁ − x₀22:41,23:12,23:51.
Step 2 — swap the norm to the x₀ norm, at cost 1/(1−ε)
- Because the Hessians are spectrally similar along the whole segment, the
∇²f(x₁)⁻¹-norm can be replaced by the∇²f(x₀)⁻¹-norm at a factor of at most1/(1−ε)21:04,22:24. - Reason: inverting the sandwich flips the constants — the inverses satisfy the same inequality with
1/(1+ε)on one side and1/(1−ε)on the other22:41. - Justification given when questioned: apply the spectral theorem,
(UΛUᵀ)⁻¹ = UΛ⁻¹Uᵀ, and the eigenvalues agree up to1 ± ε. Worst case the exponent is off by a square root, which only changesεby a constant factor30:25,30:42,31:38,31:51. - Pulling the norm inside the integral only increases the bound (triangle inequality for integrals)
24:27,24:51.
Step 3 — rewrite as a whitened quadratic form
-
Expanding both sides using
‖x‖_A = √(xᵀAx)shows the integrand equals a norm of∇²f(x₀)^{−1/2}∇f(x₀)conjugated by the matrixM := ∇²f(x₀)^{−1/2}[∇²f(x₀) − ∇²f(x_t)]∇²f(x₀)^{−1/2} -
This identity is pure symbol-pushing; it is left to the notes and to the reader to verify, and is "not something deep"
27:25,27:56.
Step 4 — bound M²
-
By the whitened sandwich from topic 4:
−εI ⪯ M ⪯ εI27:56,28:34. -
The quadratic form involves
M², andM²is PSD, so28:34,29:22:0 ⪯ M² ⪯ ε²·I
Step 5 — collapse to the conclusion
- With
ε²Iin the middle, the quantity becomes a scaled L2 norm: the plain L2 norm scaled byε(from the square root)32:10,32:50. - Computing
zᵀzonz = ∇²f(x₀)^{−1/2}∇f(x₀)turns the two−1/2powers into∇²f(x₀)⁻¹, recovering exactly the original centrality-style norm32:50,33:20. - Final factor:
ε/(1−ε), the extra1/(1−ε)coming from the norm swap in step 233:20.
graph TD
A["∇f(x₁) via FTC"] --> B["substitute Newton step<br/>→ integral of Hessian difference"]
B --> C["swap to x₀-norm<br/>cost 1/(1−ε)"]
C --> D["push norm inside integral"]
D --> E["whiten: matrix M"]
E --> F["−εI ⪯ M ⪯ εI ⇒ 0 ⪯ M² ⪯ ε²I"]
F --> G["‖∇f(x₁)‖ ≤ ε/(1−ε) · ‖∇f(x₀)‖"]