BFGS can converge to non-Clarke-stationary points with positive probability
Lewis and Overton (2013) asked whether BFGS with Armijo–Wolfe line search, on locally Lipschitz semialgebraic functions with bounded sublevel sets, has only Clarke-stationary cluster points almost surely under random initialization.
A coercive, locally Lipschitz, semialgebraic function on R^2 is constructed that is positively homogeneous near the origin, which sustains a shrinking BFGS cycle. BFGS scaling and covariance identities let the cycle repeat under linear changes of coordinates. A stability argument extends the failure to nearby initializations. Explicit inequalities are certified by rigorous computation, and the main result is formalized in Lean.
BFGS converges to a non-Clarke-stationary point from a nonempty open set of initial points and inverse-Hessian estimates. This set has positive probability under Gaussian–Wishart initialization, which disproves the almost-sure stationarity assertion.
