A Fast Newton Method Under Local Lipschitz Smoothness
Serge Gratton, Sadok Jerad, Philippe L. Toint
TL;DR
The paper tackles unconstrained nonconvex optimization under a local Lipschitz smoothness condition instead of the global Hessian Lipschitz assumption. It introduces AN2CLS, a fast second-order method that alternates Newton-like steps with negative-curvature steps and uses a Stepcomp procedure to adaptively regularize and choose steps, achieving the optimal first-order rate $O(|\log(\epsilon)|\epsilon^{-3/2})$ under local smoothness. It extends this framework to second-order critical points via SOAN2CLS, providing a corresponding complexity bound and preserving the same adaptive, geometry-aware philosophy. Two practical implementations are explored—AN2CLSE (exact solves) and AN2CLSK (Krylov-based)—with preconditioning discussed—demonstrating competitiveness against standard trust-region and cubic-regularization baselines on CUTEst problems and illustrating the method’s practical potential for large-scale nonconvex optimization.
Abstract
A new, fast second-order method is proposed that achieves the optimal $\mathcal{O}\left(|\log(ε)|ε^{-3/2}\right)$ complexity to obtain first-order $ε$-stationary points. Crucially, this is deduced without assuming the standard global Lipschitz Hessian continuity condition, but only using an appropriate local smoothness requirement. The algorithm exploits Hessian information to compute a Newton step and a negative curvature step when needed, in an approach similar to that of the AN2C method. Inexact versions of the Newton step and negative curvature are proposed in order to reduce the cost of evaluating second-order information. Details are given of such an iterative implementation using Krylov subspaces. An extended algorithm for finding second-order critical points is also developed and its complexity is again shown to be within a log factor of the optimal one. Initial numerical experiments are discussed for both factorised and Krylov variants, which demonstrate the competitiveness of the proposed algorithm.
