Theorem 6.1 — global convergence of Algorithm INB
ProvedGlobalInexactNewton.Backtracking.theorem_6_1Let be a finite-dimensional real normed space and continuously differentiable. Assume that Algorithm INB, with parameters , and , does not break down. Iteration chooses and an inexact Newton step with , and backtracks , until
then .
If is a limit point of such that is invertible, then
Furthermore, and for all sufficiently large : eventually the initial inexact Newton step is accepted without backtracking.
The last clause means that, near a solution, the local convergence rate of the method is governed by the choice of initial levels , as in the local theory of inexact Newton methods.
Formalization Note "Does not break down" is the run predicate IsINBRun: the run is an infinite sequence, and every while-loop exits after finitely many passes . and are the values on exit from the loop. A limit point is MapClusterPt, and invertibility of is ContinuousLinearMap.IsInvertible.
import Mathlib import Definitions.Def_GlobalInexactNewton_Backtracking_Method open Filter Topology
namespace GlobalInexactNewton.Backtracking
/-- Eisenstat–Walker (1994), p. 410, Theorem 6.1 (global convergence of Algorithm INB): assume
Algorithm INB does not break down. If `x_*` is a limit point of `{x_k}` such that `F'(x_*)` is
invertible, then `F(x_*) = 0` and `x_k → x_*`. Furthermore, `s_k = s̄_k` and `η_k = η̄_k` for all
sufficiently large `k`, where `s_k = trialStep (sbar k) (θ k) (m k)` and
`η_k = trialLevel (ηbar k) (θ k) (m k)` are the values on exit from the while-loop. -/
theorem theorem_6_1 {E : Type*} [NormedAddCommGroup E] [NormedSpace ℝ E] [FiniteDimensional ℝ E]
(F : E → E) (hF : ContDiff ℝ 1 F) (ηmax t θmin θmax : ℝ) (x : ℕ → E) (ηbar : ℕ → ℝ)
(sbar : ℕ → E) (θ : ℕ → ℕ → ℝ) (m : ℕ → ℕ)
(hrun : IsINBRun F ηmax t θmin θmax x ηbar sbar θ m)
(xstar : E) (hlim : MapClusterPt xstar atTop x) (hinv : (fderiv ℝ F xstar).IsInvertible) :
F xstar = 0 ∧ Tendsto x atTop (𝓝 xstar) ∧
∀ᶠ k in atTop, trialStep (sbar k) (θ k) (m k) = sbar k ∧
trialLevel (ηbar k) (θ k) (m k) = ηbar k := by sorry
end GlobalInexactNewton.Backtracking
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.