Source-linked AI summary

The Complexity of Large-scale Convex Programming under a Linear Optimization Oracle

Guanghui Lan

arXiv:1309.5550v2math.OC

TL;DR

Large-scale convex programming motivates methods with cheaper iterations, but complexity under a linear optimization oracle remains limited. The paper establishes lower bounds and develops LCP methods, showing optimal or nearly optimal complexity for conditional-gradient variants and new methods in supported problem classes.

  • Problem

    Large-scale convex programs can be beyond second-order interior-point methods, while the complexity of solving them under a linear optimization oracle remains limited.

    Source: §1 · pp. 1–2
  • Method

    The paper analyzes LCP methods, which use linear optimization subproblems, and develops conditional-gradient variants and PA-CndG and PDA-CndG algorithms.

    Source: §1 · p. 3, §5 · p. 24
  • Results

    Conditional gradient and its variants attain optimal or nearly optimal complexity for several smooth, nonsmooth, and saddle-point classes, while PA-CndG and PDA-CndG also exhibit optimal convergence rates for smooth problems under an LO oracle.

    Source: §1 · p. 3, §5 · p. 24
  • Takeaways & Limitations

    The results characterize when linear-optimization-based methods can match lower complexity bounds across several classes of large-scale convex programs.

    Source: §1 · p. 3, §5 · p. 24

    For smooth problems, the conditional-gradient convergence rate depends on summed distances that need not vanish unless the optimum is an extreme point of the feasible set.

    Source: §3.1 · p. 11

Abstract

from arXiv · show

This paper considers a general class of iterative optimization algorithms, referred to as linear-optimization-based convex programming (LCP) methods, for solving large-scale convex programming (CP) problems. The LCP methods, covering the classic conditional gradient (CG) method (a.k.a., Frank-Wolfe method) as a special case, can only solve a linear optimization subproblem at each iteration. In this paper, we first establish a series of lower complexity bounds for the LCP methods to solve different classes of CP problems, including smooth, nonsmooth and certain saddle-point problems. We then formally establish the theoretical optimality or nearly optimality, in the large-scale case, for the CG method and its variants to solve different classes of CP problems. We also introduce several new optimal LCP methods, obtained by properly modifying Nesterov's accelerated gradient method, and demonstrate their possible advantages over the classic CG for solving certain classes of large-scale CP problems.

1 Introduction

Large-scale convex programs can exceed the capability of second-order interior-point methods, while complexity under a linear optimization oracle remains comparatively underexplored. The paper develops lower bounds and analyzes LCP methods, including conditional gradient and new variants.

1 Introduction

  • LCP methods solve large-scale convex programs by iteratively calling a linear optimization oracle, and a generic algorithm can form outputs as convex combinations of earlier iterates.The framework permits oracle input vectors based on gradients, approximations, random noise, or second-order information.Source: §1 · p. 1
  • The paper establishes LO-oracle lower bounds for smooth, nonsmooth, and saddle-point problems, including bounds that apply to LCP methods using higher-order information.For smooth problems, the bounds make dependence on dimension and problem parameters explicit; the supplied passages do not provide the bound expressions.Source: §1 · p. 3
  • Conditional gradient is optimal for smooth problems in sufficiently large dimensions and can attain optimal or nearly optimal complexity for several nonsmooth and saddle-point settings.For smooth problems, the stated sufficient condition is n ≥ LD_X^2/ε; variants use smoothing, including randomized smoothing for general nonsmooth problems.Source: §1 · p. 3
  • PA-CndG and PDA-CndG replace projection subproblems in Nesterov’s accelerated methods with linear optimization subproblems.Their convergence rates depend on sums of distances among solutions; under specified conditions on the LO oracle they achieve O(1/√ε) iteration complexity for smooth problems. Preliminary experiments report that PDA-CndG can significantly outperform conditional gradient on some box-constrained problems.Source: §1 · pp. 3–4

1.1 Notation and terminology

  • The notation section defines set diameters, conjugate norms, operator norms, and linear approximations used in the paper.It also notes that the smoothness constant L depends on the chosen norm.Source: §1.1 · p. 4

2 Lower Complexity Bounds for CP under an LO oracle

The section develops lower complexity bounds for linear-optimization-based convex programming methods across smooth, nonsmooth, and certain saddle-point problems. It introduces a flexible LCP framework and shows that the bounds apply beyond first-order methods.

2.1 A generic LCP algorithm

  • The generic LCP framework calls a linear optimization oracle each iteration and allows flexible search directions and convex-combination outputs.It includes the classic conditional gradient method and several new LCP methods as special cases.Source: §2.1 · pp. 4–5

2 Lower Complexity Bounds for CP under an LO oracle

  • The section extends earlier LCP lower-bound work by making problem-dimension and parameter dependence explicit and treating nonsmooth and saddle-point problems as well as smooth optimization.The analysis assumes a resisting LO oracle that provides the least information needed by the algorithm in the worst case.Source: §2.2 · p. 5

2.2 Lower complexity bounds for smooth minimization

  • For smooth problems, any LCP method faces a dimension-dependent lower complexity bound; assuming strong convexity does not improve it when n is sufficiently large.The smooth lower-bound analysis is developed for a class of convex problems with Lipschitz gradients.Source: §2.2 · pp. 5, 7

2.3 Lower complexity bounds for nonsmooth minimization

  • For nonsmooth problems, the authors derive bounds for general convex minimization and bilinear saddle-point classes; the general bound applies even to LCP methods using higher-order information.The general nonsmooth bound is reported to be of the same order as the bound for general first-order methods.Source: §2.3 · pp. 7–8

3 The Optimality of CndG Methods for CP under an LO oracle

This section establishes optimality or near-optimality results for conditional-gradient methods and variants under a linear optimization oracle across smooth, saddle-point, strongly convex, and general nonsmooth convex problems. For strongly convex problems, it analyzes a shrinking method with an enhanced oracle; for general nonsmooth problems, the stated guarantee is in expectation, with a possible high-probability result omitted.

3.1 Optimal CndG methods for smooth problems under an LO oracle

  • With either standard stepsize policy, classic CndG achieves O(1/k) convergence and is an optimal LCP method for smooth problems when the dimension is sufficiently large.Its convergence analysis makes the rate depend explicitly on the distance between the previous iterate and the linear-oracle output; this distance need not vanish and can limit the bound.Source: §3.1 · pp. 9–11

3.2 Optimal CndG methods for saddle point problems under an LO oracle

  • For saddle-point problems, a CndG variant using variable smoothing parameters achieves optimal complexity when the dimension is sufficiently large.The method need not know the target accuracy in advance, and inaccurate estimates of problem parameters affect its convergence rate only by a constant factor.Source: §3.2 · pp. 11–13

3.3 Nearly optimal CndG methods for general nonsmooth problems under an LO oracle

  • For saddle-point problems, randomized CndG avoids the smoothing subproblems used by the method in Subsection 3.2 and can reduce them to linear optimization over Y when the stated term is zero.The authors describe this as the first proposal of this type of optimization algorithm in the literature.Source: §3.3 · pp. 15–16
  • A randomized CndG method for general nonsmooth problems achieves a nearly optimal complexity bound relative to the lower bound.It uses convolution-based smoothing and sampled gradients; the analysis sets the sample count at iteration k to k, and generated samples can be reused in later iterations.Source: §3.3 · pp. 14–16
  • The analysis gives relation (3.29) in expectation; an overwhelming-probability result via large-deviation analysis is possible but omitted for simplicity, with similar developments noted in [24].Source: §3.3 · p. 16

3.4 CndG methods for strongly convex problems under an enhanced LO oracle

  • The shrinking CndG method uses an enhanced LO oracle and contracts the feasible region around the current iterate at each outer iteration.For a polytope, the passage suggests choosing the norm so the enhanced oracle’s local constrained problem is comparable in complexity to the original LO problem, though this may increase L/μ.Source: §3.4 · pp. 16–17
  • The shrinking method performs 8L/μ inner iterations per outer round and at most ⌈max(log(μR0/ε), 1)⌉ outer rounds when finding an ε-solution.The theorem bounds the enhanced-oracle calls; the iteration counts follow from the proof and algorithm description.Source: §3.4 · p. 17

4 New variants of for LCP methods

Section 4 introduces LCP variants that replace projection or prox-mapping steps with linear optimization, then analyzes primal- and primal-dual-averaging CndG methods. In the numerical illustration, PDA-CndG is particularly advantageous on box-type constraints despite sharing similar worst-case complexity bounds with the other methods.

4 New variants of for LCP methods

  • The new methods are obtained by replacing accelerated gradient methods' projection or prox-mapping subproblems with linear optimization; the algorithms can also be modified for saddle-point, nonsmooth, and strongly convex problems.Source: §4 · p. 17

4.1 Primal averaging CndG method

  • PA-CndG averages the primal iterates before forming the gradient direction sent to the LO oracle, unlike classic CndG, which uses the previous iterate.The paper characterizes PA-CndG as an optimal LCP method for the stated smooth problem class when the dimension is sufficiently large.Source: §4.1 · pp. 18–19
  • PA-CndG's convergence bound depends on the distance between consecutive LO-oracle outputs, and a faster O(1/k^2) rate follows under the stated Hölder continuity condition with ρ in (0.5,1].The paper notes that the continuity assumptions needed for this improvement are strong and difficult to satisfy globally.Source: §4.1 · p. 19

4.2 Primal-dual averaging CndG methods

  • PDA-CndG uses a weighted average of past gradients as its LO-oracle direction and achieves the same stated optimal LCP status as PA-CndG for sufficiently large dimension.It also provides an online lower bound on the optimal value, enabling termination based on the gap between that bound and the current objective value.Source: §4.2 · pp. 20–21
  • PDA-CndG's faster-rate connection to general optimal first-order methods is conditional: the LO oracle need not satisfy any positive-ρ Hölder condition.The paper says the oracle may satisfy only the ρ = 0 case with Q = D_X.Source: §4.2 · p. 22

4.3 Numerical Illustration

  • On simplex/spectrahedron QPs, the three methods perform similarly, with CndG slightly ahead; PDA-CndG outperforms the others on hypercubes and by up to 2 orders of magnitude on hypercube–simplex intersections.For hypercube QPs, PDA-CndG's function values at 100 iterations are already comparable to those of CndG and PA-CndG at 1,000 iterations.Source: Table 3 · p. 23, Table 2 · p. 23, Table 4 · p. 24

5 Concluding remarks

The paper establishes complexity bounds and optimality results for LCP methods across several classes of convex programs. It also introduces new methods with optimal smooth-problem convergence rates and reports preliminary gains for PDA-CndG on certain large-scale problems.

  • The paper establishes lower complexity bounds for LCP methods across different classes of convex programming problems.Source: §5 · p. 24
  • The classic CndG method is optimal among LCP methods for solving smooth convex programming problems.Source: §5 · p. 24
  • New CndG variants are optimal or nearly optimal for certain saddle-point and general nonsmooth problems under a linear optimization oracle.Source: §5 · p. 24
  • PA-CndG and PDA-CndG, developed by modifying Nesterov’s accelerated gradient method, have optimal smooth-problem convergence rates under a linear optimization oracle.Source: §5 · p. 24
  • Preliminary numerical experiments show that PDA-CndG can significantly outperform classic CndG on certain classes of large-scale problems.Source: §5 · p. 24
Loading 1309.5550v2…