Theorem 4.1 — steepest descent with relative gradient error contracts by per step
ProvedDynSampleSize.Det.theorem_4_1Let be twice continuously differentiable and uniformly convex, satisfying (4.2):
with constants . Let be its minimizer and assume . Let , let be arbitrary vectors (approximations to ), and let be generated by iteration (4.1) with the steplength (4.7),
Put . Then:
- (4.9) If at an iteration the approximate gradient satisfies (4.8), , then
If moreover (4.8) holds at every iteration, then:
- (4.10) .
- (4.11) For every , every iteration count with
satisfies ; that is, the number of iterations needed to reach is at most . 4. (4.12) For every ,
The theorem shows that a gradient approximation whose error is at most a fixed fraction of its own norm preserves the linear convergence rate of steepest descent, at the price of the factor . It motivates the sample-size rule of the stochastic analysis in §4.2.
Formalization Note The space is EuclideanSpace ℝ (Fin m), lam is the paper's , and assumption (4.2) is HessianBounds J lam L. The minimizer is a point with for all , and is a hypothesis. (4.11) is stated in the form the proof establishes ((4.18), p. 10): every beyond the bound is -accurate; is the natural logarithm. If Lean's convention applies, and the statement remains true because then every iterate is optimal.
import Mathlib import Definitions.Def_DynSampleSize_Det_Setting open Filter Topology
namespace DynSampleSize.Det
theorem theorem_4_1 {m : ℕ} (J : EuclideanSpace ℝ (Fin m) → ℝ) (lam L : ℝ)
(hJ : HessianBounds J lam L) (wstar : EuclideanSpace ℝ (Fin m))
(hmin : ∀ w, J wstar ≤ J w) (hJstar : J wstar = 0)
(θ : ℝ) (hθ0 : 0 < θ) (hθ1 : θ < 1) (g w : ℕ → EuclideanSpace ℝ (Fin m))
(hw : ∀ k, w (k + 1) = w k - ((1 - θ) / L) • g k) :
-- (4.9): one step, under (4.8) at that iteration only
(∀ k, ‖g k - gradient J (w k)‖ ≤ θ * ‖g k‖ →
J (w (k + 1)) ≤ (1 - beta θ * lam / L) * J (w k)) ∧
-- (4.10)–(4.12): under (4.8) at every iteration
((∀ k, ‖g k - gradient J (w k)‖ ≤ θ * ‖g k‖) →
Tendsto w atTop (𝓝 wstar) ∧
(∀ ε : ℝ, 0 < ε → ∀ k : ℕ,
L / (beta θ * lam) * (Real.log (1 / ε) + Real.log (J (w 0))) < (k : ℝ) →
J (w k) < ε) ∧
(∀ k : ℕ, ‖g k‖ ^ 2 ≤
1 / (1 - θ) ^ 2 * (2 * L ^ 2 * J (w 0) / lam) * (1 - beta θ * lam / L) ^ k)) := by sorry
end DynSampleSize.Det
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.