ZIVR: An Incremental Variance Reduction Technique For Zeroth-Order Composite Problems
Silan Zhang, Yujie Tang
TL;DR
This work tackles zeroth-order finite-sum composite optimization where gradient information is unavailable and standard first-order variance reduction cannot be directly applied. It introduces ZIVR, a general framework that maintains an incremental Jacobian estimator $J_k$ to construct a variance-reduced gradient $g_k$ by blending historical gradient information with current two-point zeroth-order estimates, enabling pure $2$-point updates and memory-efficient variants. The authors prove comprehensive convergence guarantees for strongly convex, convex, and non-convex regimes, achieving rates that match first-order methods up to a factor of $O(d)$ and providing oracle-complexity bounds that scale favorably with problem dimensions. They validate the approach through numerical experiments on logistic and Cox regression, showing that ZIVR outperforms competing zeroth-order methods while reducing batch-size requirements. Overall, ZIVR offers a scalable and principled way to perform constrained or composite zeroth-order optimization with provable efficiency guarantees.
Abstract
This paper investigates zeroth-order (ZO) finite-sum composite optimization. Recently, variance reduction techniques have been applied to ZO methods to mitigate the non-vanishing variance of 2-point estimators in constrained/composite optimization, yielding improved convergence rates. However, existing ZO variance reduction methods typically involve batch sampling of size at least $Θ(n)$ or $Θ(d)$, which can be computationally prohibitive for large-scale problems. In this work, we propose a general variance reduction framework, Zeroth-Order Incremental Variance Reduction (ZIVR), which supports flexible implementations$\unicode{x2014}$including a pure 2-point zeroth-order algorithm that eliminates the need for large batch sampling. Furthermore, we establish comprehensive convergence guarantees for ZIVR across strongly-convex, convex, and non-convex settings that match their first-order counterparts. Numerical experiments validate the effectiveness of our proposed algorithm.
