math.OCOct 9, 2026

Optimal Algorithms for Nonsmooth Nonconvex Stochastic Optimization with Constraints

Authors: Gon Buzaglo, Guy Kornowski

Organizations: Princeton University · Institute for Advanced Study

Abstract

We present new algorithms for stochastic optimization problems in which the objective is neither smooth nor convex, subject to a convex constraint set. Our algorithms return a (δ,ε)(δ,\varepsilon)-Goldstein stationary point, as recently defined for constrained problems by Liu et al. (2024), with improved complexities, and are the first to achieve optimal convergence rates for this task in several settings. In particular, we design projection-based and projection-free algorithms that achieve a O~(δ−1ε−3)\widetilde O(δ^{-1}\varepsilon^{-3}) rate using stochastic gradients, which we show to be optimal, and a O~(dδ−1ε−3)\widetilde O(dδ^{-1}\varepsilon^{-3}) rate using zeroth-order queries, which is also optimal in its dependence on the dimension dd and improves over prior work. Notably, our lower bound is given in terms of the Frank-Wolfe gap and holds even for smooth problems, which is of independent interest. These results show that in nonconvex stochastic optimization, the convergence rate is degraded neither by the presence of constraints nor by nonsmoothness. Our analysis is based on a reduction to online linear optimization, extending a related technique from unconstrained optimization.

Explore similar work

May 23, 2026cs.LG

Zeroth-Order Nonconvex Nonsmooth Optimization with Heavy-Tailed Noise

This paper considers the nonconvex nonsmooth problem in which the objective function is Lipschitz continuous. We focus on the stochastic setting where the algorithm can access stochastic function value evaluations with heavy-tailed noise, which is prevalent in many popular machine learning applications. We propose a stochastic zeroth-order algorithm that refines the framework of online-to-nonconvex conversion by clipping the two-point gradient estimator. The theoretical analysis shows that our algorithm can find a (δ,ε)(δ, ε)-Goldstein stationary point with zeroth-order oracle complexity of O(dp2(p−1)δ−1ε−2p−1p−1){\mathcal O}(d^{\frac{p}{2(p-1)}}δ^{-1}ε^{-\frac{2p-1}{p-1}}), where dd is the problem dimension and p∈(1,2]p\in(1,2] is the order of bounded moments. Note that our dependence on dimension dd matches the best-known results of stochastic zeroth-order optimization for finding the sub-optimal solution of a stochastic convex nonsmooth problem. In addition, our dependence on accuracy parameters δδ and εε is consistent with that of the best-known stochastic first-order algorithms for stochastic nonconvex nonsmooth problems. Finally, we conduct numerical experiments to demonstrate the effectiveness of the proposed method.
Oct 5, 2026math.OC

Dimension-Free Decentralized Nonsmooth Nonconvex Stochastic Optimization

We investigate decentralized nonsmooth nonconvex stochastic optimization over a network of nn nodes, with the goal of finding an (δ,ε)(δ,ε)-Goldstein stationary point. The best existing algorithm achieves O(δ−1(ε−3+dε−1))O(δ^{-1}(ε^{-3}+dε^{-1})) sample complexity and O~(γ−1/2δ−1(ε−3+dε−1))\widetilde{O}(γ^{-1/2}δ^{-1}(ε^{-3}+dε^{-1})) communication complexity, where dd is the problem dimension and γγ is the spectral gap of the communication matrix. However, the polynomial dependence on dd can be a major bottleneck in high-dimensional regimes. In this paper, we propose a novel algorithm that achieves O(δ−1ε−3)O(δ^{-1}ε^{-3}) sample complexity and O~(γ−1/2δ−1ε−3)\widetilde{O}(γ^{-1/2}δ^{-1}ε^{-3}) communication complexity. The primary technique is an elegant decentralized online-to-nonconvex conversion that reduces the original problem to a decentralized online convex optimization (D-OCO) problem. A key property of our conversion is that its consensus requirements can be inherited directly from the consensus of the underlying D-OCO decisions. In particular, this property enables us to establish an explicit connection between the dimension dependence and the consensus error, which in turn shows that the polynomial dependence on dd can be removed with only logarithmic additional communication.
Sep 20, 2026math.OC

The Exponential Price of Determinism in Nonsmooth Nonconvex Optimization

We study the complexity of finding (δ,ε)(δ,ε)-Goldstein stationary points of nonsmooth nonconvex Lipschitz functions. By now, it is known that randomized first-order algorithms can solve this task with a dimension-free oracle complexity [Zhang et al., 2020], whereas deterministic algorithms cannot, as their complexity must scale at least linearly with the dimension dd [Jordan et al., 2023, Tian and So, 2024]. This leaves open whether deterministic algorithms can nevertheless solve the problem with oracle complexity polynomial in dd. We answer this question negatively by proving a lower bound of order (1/ε)Ω(d)(1/ε)^{Ω(d)} for deterministic algorithm, closing the exponential gap between the previously known lower and upper bounds and resolving an open problem posed by Jordan et al. [2023]. We further discuss several extensions and implications of this result to weaker stationarity notions, finding a descent direction and deterministic smoothing. Overall, our results establish an exponential computational advantage in nonsmooth nonconvex optimization offered by randomization.