Lemma 5.1 — the while-loop of Algorithms MR and TL terminates, with
ProvedGlobalInexactNewton.Backtracking.lemma_5_1Let be a finite-dimensional real normed space and continuously differentiable. Fix and , and consider the th iteration of Algorithm MR or Algorithm TL: a current iterate , an initial level and a curve satisfying (5.1) on . The while-loop starts at and, as long as , replaces by for some .
Assume and that there is with
Then, whatever factors the loop chooses:
- the while-loop terminates;
- on termination,
for every small enough that
The lemma says when neither algorithm breaks down inside the while-loop, and the lower bound (5.4) on is what makes diverge in Theorem 5.2.
Formalization Note The while-loop is identical in MR and TL, so the lemma is stated once for a single iteration, which covers both. The factors are an arbitrary sequence in , and are the levels the loop visits. "Terminates" is: some index passes (2.2) while every fails it. Part 2 is stated for every such exit index . The hypotheses force (at , (5.1) and (5.3) are incompatible with when ), so the divisions by are genuine.
import Mathlib import Definitions.Def_GlobalInexactNewton_Backtracking_Method open Filter Topology
namespace GlobalInexactNewton.Backtracking
/-- Eisenstat–Walker (1994), p. 408, Lemma 5.1, for one iteration of the while-loop shared by
Algorithms MR and TL. Data: `t ∈ (0,1)`, `0 < θ_min < θ_max < 1`, the current iterate `x_k`, an
initial level `η̄_k ∈ [0,1)`, a curve `σ_k` with (5.1) on `[η̄_k, 1]`, and an arbitrary choice of
factors `θ_j ∈ [θ_min, θ_max]`. If `F(x_k) ≠ 0` and `Γ` satisfies (5.3)
`‖σ_k(η)‖ ≤ Γ (1 - η) ‖F(x_k)‖` on `[η̄_k, 1]`, then the while-loop terminates (some trial passes (2.2)
after all earlier trials failed), and at termination (5.4)
`1 - η_k ≥ min {1 - η̄_k, θ_min δ / (Γ ‖F(x_k)‖)}` for every `δ > 0` with (5.5)
`‖F(x) - F(x_k) - F'(x_k)(x - x_k)‖ ≤ ((1 - t)/Γ) ‖x - x_k‖` on `N_δ(x_k)`. -/
theorem lemma_5_1 {E : Type*} [NormedAddCommGroup E] [NormedSpace ℝ E] [FiniteDimensional ℝ E]
(F : E → E) (hF : ContDiff ℝ 1 F) (t θmin θmax : ℝ) (ht : t ∈ Set.Ioo (0 : ℝ) 1)
(hθ : 0 < θmin ∧ θmin < θmax ∧ θmax < 1)
(xk : E) (ηbar : ℝ) (hηbar : ηbar ∈ Set.Ico (0 : ℝ) 1) (σ : ℝ → E)
(hσ : ∀ η ∈ Set.Icc ηbar 1, InexactNewtonCond F xk (σ η) η)
(hFxk : F xk ≠ 0) (Γ : ℝ) (hΓ : ∀ η ∈ Set.Icc ηbar 1, ‖σ η‖ ≤ Γ * (1 - η) * ‖F xk‖)
(θ : ℕ → ℝ) (hθmem : ∀ j, θ j ∈ Set.Icc θmin θmax) :
(∃ m : ℕ, (∀ j < m, ¬ SuffDecrease F t xk (σ (trialLevel ηbar θ j)) (trialLevel ηbar θ j)) ∧
SuffDecrease F t xk (σ (trialLevel ηbar θ m)) (trialLevel ηbar θ m)) ∧
∀ m : ℕ, (∀ j < m, ¬ SuffDecrease F t xk (σ (trialLevel ηbar θ j)) (trialLevel ηbar θ j)) →
SuffDecrease F t xk (σ (trialLevel ηbar θ m)) (trialLevel ηbar θ m) →
∀ δ > 0, (∀ y ∈ Metric.ball xk δ,
‖F y - F xk - fderiv ℝ F xk (y - xk)‖ ≤ (1 - t) / Γ * ‖y - xk‖) →
min (1 - ηbar) (θmin * δ / (Γ * ‖F xk‖)) ≤ 1 - trialLevel ηbar θ m := by sorry
end GlobalInexactNewton.Backtracking
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.