Theorem 3.2 (second part) — convergence of order
ProvedNonsmoothNewton.Local.newton_local_orderLet be locally Lipschitz, let satisfy , let be semismooth at , and let every be nonsingular. Suppose in addition that and is -order semismooth at . Then there are and a constant such that every run of the nonsmooth Newton method
started at a point with (for every choice of the ) satisfies
For this is local quadratic convergence of the generalized-Jacobian Newton method.
Formalization Note The paper's "the convergence of (3.2) is of order " is pinned as the displayed inequality, with and chosen before the run (uniform over all starting points and all choices of ), which is what the proof gives. A run is a sequence with and .
import Mathlib import Definitions.Def_NonsmoothNewton_Shared_clarkeJac import Definitions.Def_NonsmoothNewton_Local_SemismoothAt import Definitions.Def_NonsmoothNewton_Local_IsNewtonRun open Filter Topology
namespace NonsmoothNewton.Local
/-- Qi–Sun (1993), Theorem 3.2 (second sentence), p. 359. Under the hypotheses of Theorem 3.2,
if in addition `F` is `p`-order semismooth at `x*` (`0 < p ≤ 1`), the convergence of (3.2) is
of order `1 + p`: there are `δ > 0` and `C` such that every run started within `δ` of `x*`
satisfies `‖x^{k+1} - x*‖ ≤ C ‖x^k - x*‖^{1+p}` for all `k`. -/
theorem newton_local_order {n : ℕ}
(F : EuclideanSpace ℝ (Fin n) → EuclideanSpace ℝ (Fin n)) (hF : LocallyLipschitz F)
(xstar : EuclideanSpace ℝ (Fin n)) (hroot : F xstar = 0) (hsemi : SemismoothAt F xstar)
(hns : ∀ V ∈ NonsmoothNewton.Shared.clarkeJac F xstar, IsUnit V)
(p : ℝ) (hp0 : 0 < p) (hp1 : p ≤ 1) (hpord : POrderSemismoothAt p F xstar) :
∃ δ > 0, ∃ C : ℝ, ∀ (x : ℕ → EuclideanSpace ℝ (Fin n))
(V : ℕ → EuclideanSpace ℝ (Fin n) →L[ℝ] EuclideanSpace ℝ (Fin n)),
‖x 0 - xstar‖ < δ → IsNewtonRun F x V →
∀ k, ‖x (k + 1) - xstar‖ ≤ C * ‖x k - xstar‖ ^ (1 + p) := by sorry
end NonsmoothNewton.Local
Read-back
What the Lean code literally says, in plain math · claude-opus-5-5
Setting. may be . carries the Euclidean norm. is locally Lipschitz on all of . is the set of linear maps given by the external definition clarkeJac, which is not shown here.
Hypotheses. A point and a real number satisfy all of the following.
- Root: .
- Semismoothness at : is -Lipschitz on a neighbourhood of for some . Also, for every there is such that for every there is a with whenever , and .
- Invertibility: every is invertible.
- Range of : .
- -order semismoothness at : is semismooth at , and there are real and with
for all , all , and every one-sided directional derivative .
Conclusion. There exist and a real constant with the following property. Take any sequence in and any sequence of linear maps such that:
- , and
- for every , and .
For every such pair and every ,
The theorem does not assert that any such sequence exists. It also does not state separately that the iterates remain within of or converge; only the displayed inequality for all is asserted.
Degenerate cases.
- Landing on . If some , the inequality forces .
- No runs. If no sequences satisfying the run conditions start within of , the conclusion is vacuous. This includes the case where is empty at the relevant points.
- Empty . The invertibility hypothesis is vacuous.
- . Every point equals , and the statement holds trivially.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.