One-Shot Klein Cutting Planes for Lipschitz Geodesically Convex Optimization in Hyperbolic Space

TL;DR

This paper introduces a one-shot Klein cutting-plane method for Lipschitz geodesically convex optimization in hyperbolic space, achieving near-optimal oracle complexity and localization without convex coordinate pullback.

Contribution

It provides the first deterministic first-order method for negative curvature hyperbolic spaces with provable guarantees, extending convex optimization techniques to quasiconvex settings.

Findings

Achieves oracle complexity of O(d^2 ζ_s log(e/ε)) for hyperbolic convex optimization.

Localization update costs O(d^2) operations for dimensions d≥2.

Method works in quasiconvex setting without convex coordinate pullback.

Abstract

We solve the negative constant-curvature case of the COLT 2023 open problem of Criscitiello, Mart\'inez-Rubio, and Boumal on deterministic first-order methods for Lipschitz geodesically convex optimization. Let \[ \HH^d_{-\kappaC^2}=\{X\in\R^{d+1}:\ipL{X}{X}=-1,\ X_0>0\}, \qquad \ip{U}{V}_{X}=\kappaC^{-2}\ipL{U}{V}, \] so the sectional curvature is $-\kappaC^2$. If \[ f:\bar B_{\HH}(x_0,r)\to\R \] is geodesically convex and $M$-Lipschitz, and $s=\kappaC r$, our one-shot Klein cutting-plane method returns a queried point $\hat x$ with \[ f(\hat x)-\min_{\bar B_{\HH}(x_0,r)}f\le \eps Mr \] using at most \[ \left\lceil 2d(d+1) \log\!\left(\frac{16\sinh s\cosh s}{s\eps}\right)\right\rceil \] oracle calls. For $d\ge2$ each localization update costs $O(d^2)$ arithmetic operations; for $d=1$ an interval variant satisfies the same bound. Consequently \[ N=O\bigl(d^2(s+\log(e/\eps))\bigr) =O\bigl(d^2\zeta_s\log(e/\eps)\bigr), \qquad \zeta_s=s/\tanh s . \] The argument is not a convex coordinate pullback: in the Beltrami--Klein chart the objective is generally only quasiconvex. The key point is that every Riemannian subgradient halfspace becomes an exact Euclidean central cut. For \[ \theta=\kappaC\dist(X,Y), \] \[ \ip{g}{\log_XY}_{X} =\frac{\theta}{\kappaC^2\sinh\theta}\ipL{g}{Y}, \] and tangency at $X$ turns $\ipL{g}{Y}\le0$ into \[ \gbar^{\mathsf T}(u-c)\le0, \qquad u=\Phi(Y),\quad c=\Phi(X). \] Thus a fixed Euclidean ellipsoid localizes the whole hyperbolic ball. The only curvature payment is the Klein distortion factor \[ \log\left(\frac{\sinh s\cosh s}{s\eps}\right) =\log(1/\eps)+2s-\log(4s)+O(e^{-4s}). \]