Two-stage stochastic algorithm for solving large-scale (non)-convex separable optimization problems under affine constraints

TL;DR

This paper introduces a two-stage stochastic algorithm for large-scale (non)-convex separable optimization problems with affine constraints, significantly reducing the number of Fenchel conjugate computations needed for primal solutions.

Contribution

It proposes a novel two-stage method combining stochastic dual subgradient and block-coordinate Frank-Wolfe algorithms, improving efficiency for large-scale problems.

Findings

Requires fewer Fenchel conjugate calls: O(1/ε² + N/ε^{2/3})

Applicable to nonconvex functions with similar convergence rates

Achieves ε-optimal primal solutions with high probability

Abstract

We consider nonsmooth optimization problems under affine constraints, where the objective consists of the average of the component functions of a large number $N$ of agents, and we only assume access to the Fenchel conjugate of the component functions. The algorithm of choice for solving such problems is the dual subgradient method, also known as dual decomposition, which requires $O(\frac{1}{\epsilon^2})$ iterations to reach $\epsilon$-optimality in the convex case. However, each iteration requires computing the Fenchel conjugate of each of the $N$ agents, leading to a complexity $O(\frac{N}{\epsilon^2})$ which might be prohibitive in practical applications. To overcome this, we propose a two-stage algorithm, combining a stochastic subgradient algorithm on the dual problem, followed by a block-coordinate Frank-Wolfe algorithm to obtain primal solutions. The resulting algorithm requires only $O(\frac{1}{\epsilon^2} + \frac{N}{\epsilon^{2/3}})$ calls to Fenchel conjugates to obtain an $\epsilon$-optimal primal solution in expectation in the convex case. We extend our results to nonconvex component functions and show that our method still applies and gets (almost) the same convergence rate, this time only to an approximate primal solution recovering the classical duality gap bounds usually obtained using the Shapley-Folkman theorem.