English

Probabilistic algorithm for computing all local minimizers of Morse functions on a compact domain

Symbolic Computation 2026-04-06 v2 Optimization and Control

Abstract

Let K be the unit-cube in Rn and f\,: K \rightarrow R^n be a Morse function. We assume that the function f is given by an evaluation program Γ\Gamma in the noisy model, i.e., the evaluation program Γ\Gamma takes an extra parameter η\eta as input and returns an approximation that is η\eta-close to the true value of f . In this article, we design an algorithm able to compute all local minimizers of f on K . Our algorithm takes as input Γ\Gamma, η\eta, a numerical accuracy parameter ϵ\epsilon as well as some extra regularity parameters which are made explicit. Under assumptions of probabilistic nature -- related to the choice of the evaluation points used to feed Γ\Gamma --, it returns finitely many rational points of K , such that the set of balls of radius ϵ\epsilon centered at these points contains and separates the set of all local minimizers of f . Our method is based on approximation theory, yielding polynomial approximants for f , combined with computer algebra techniques for solving systems of polynomial equations. We provide bit complexity estimates for our algorithm when all regularity parameters are known. Practical experiments show that our implementation of this algorithm in the Julia package Globtim can tackle examples that were not reachable until now.

Cite

@article{arxiv.2507.13071,
  title  = {Probabilistic algorithm for computing all local minimizers of Morse functions on a compact domain},
  author = {Mohab Safey El Din and Georgy Scholten and Emmanuel Trélat},
  journal= {arXiv preprint arXiv:2507.13071},
  year   = {2026}
}
R2 v1 2026-07-01T04:05:59.204Z