Complexity guarantees and polling strategies for Riemannian direct-search methods
Bastien Cavarretta, Florentin Goyens, Clément W. Royer, Florian Yger
TL;DR
The paper develops a complexity framework for direct-search methods on smooth Riemannian manifolds, extending the Euclidean cosine measure to tangent spaces and establishing an $\mathcal{O}(\varepsilon^{-2})$ evaluation bound under appropriate PSS quality. It introduces two PSS construction paradigms—intrinsic PSSs built in tangent spaces and projected PSSs from ambient space—and analyzes their impact on the complexity measure, with a detailed sphere-manifold study showing intrinsic PSSs offer provable advantages. Empirical experiments on sphere and subspace manifolds validate the theoretical results, showing intrinsic directions yield better performance, particularly at higher codimensions, and that random rotations can enhance efficiency. The work provides practical guidance for constructing high-quality polling sets on manifolds and points to future extensions to other manifolds and probabilistic direct-search variants. Code and experiments are available on GitHub, enabling reproducibility and broader adoption.
Abstract
Direct-search algorithms are derivative-free optimization techniques that operate by polling the variable space along specific directions forming positive spanning sets (PSSs). When the problem variables are constrained to lie on a Riemannian manifold, polling must be performed along tangent directions. Although Riemannian variants of direct search have already been proposed and endowed with asymptotic guarantees, a proper generalization of PSSs on manifolds remains to be investigated. In particular, a measure of quality for those PSSs is required to obtain complexity bounds for direct search. In this paper, we derive complexity guarantees for a class of Riemannian direct-search techniques, and study two ways of generating positive spanning sets in tangent spaces. We pay particular attention the unit hypersphere case, for which we establish that generating directions directly within the tangent space leads to better complexity properties than projecting PSSs from the ambient space onto the tangent space. Our numerical experiments highlight the impact of dimension and codimension in more general settings.
