Source-linked AI summary
A Single-Loop Smoothed Gradient Descent-Ascent Algorithm for Nonconvex-Concave Min-Max Problems
Jiawei Zhang, Peijun Xiao, Ruoyu Sun, Zhi-Quan Luo
TL;DR
Unlike convex-concave settings, where a global stationary solution can be computed efficiently, finding a global optimum when f is nonconvex in x is difficult. SmoothedGDA uses an exponentially weighted primal sequence to center a quadratic proximal term, then alternates gradient descent and ascent on the resulting objective. O(1/ε^2) iterations yield an ε-stationary solution for the finite pointwise-maximum problem, while general nonconvex-concave problems have O(1/ε^4) complexity.
Problem
Unlike convex-concave settings, where a global stationary solution can be computed efficiently, finding a global optimum when f is nonconvex in x is difficult.
Source: §1 · p. 2Method
SmoothedGDA uses an exponentially weighted primal sequence to center a quadratic proximal term, then alternates gradient descent and ascent on the resulting objective.
Source: §1 · p. 3, §3.1 · p. 5Results
O(1/ε^2) iterations yield an ε-stationary solution for the finite pointwise-maximum problem, while general nonconvex-concave problems have O(1/ε^4) complexity.
Source: §1 · p. 3Takeaways & Limitations
The single-loop method is presented as suitable for large-scale and distributed settings, and its framework extends to multi-block problems.
Source: §6 · p. 12For general nonconvex-concave problems, the method requires β proportional to ε^2, giving higher complexity than the finite pointwise-maximum case.
Source: §4.1 · p. 10
Abstract
from arXiv · showhide
Nonconvex-concave min-max problem arises in many machine learning applications including minimizing a pointwise maximum of a set of nonconvex functions and robust adversarial training of neural networks. A popular approach to solve this problem is the gradient descent-ascent (GDA) algorithm which unfortunately can exhibit oscillation in case of nonconvexity. In this paper, we introduce a "smoothing" scheme which can be combined with GDA to stabilize the oscillation and ensure convergence to a stationary solution. We prove that the stabilized GDA algorithm can achieve an $O(1/\epsilon^2)$ iteration complexity for minimizing the pointwise maximum of a finite collection of nonconvex functions. Moreover, the smoothed GDA algorithm achieves an $O(1/\epsilon^4)$ iteration complexity for general nonconvex-concave problems. Extensions of this stabilized GDA algorithm to multi-block cases are presented. To the best of our knowledge, this is the first algorithm to achieve $O(1/\epsilon^2)$ for a class of nonconvex-concave problem. We illustrate the practical efficiency of the stabilized GDA algorithm on robust training.
1 Introduction
The paper studies min-max problems with objectives nonconvex in the minimizing variable and concave in the maximizing variable, including finite pointwise maxima. It introduces SmoothedGDA, a single-loop method intended to stabilize GDA and improve convergence complexity.
1 Introduction
- The formulation covers objectives nonconvex in x and concave in y, including minimizing the pointwise maximum of finitely many nonconvex functions.The paper connects this formulation to adversarial training, fairness training, and distribution-agnostic meta-learning.Source: §1 · p. 2
- Prior multi-loop methods have a best-known O(1/ε^2.5) rate for problem (1.2), while their inner-loop costs depend on ε and their acceleration is difficult to extend to multi-block problems.The paper identifies multi-block structure as important for distributed training.Source: §1 · pp. 2–3
- With constant step sizes, GDA can oscillate when the maximization problem is not strongly concave; for nonconvex-concave problems, its stated complexity is O(1/ε^5).The cited comparison also notes a single-loop method with O(1/ε^4) complexity that requires diminishing step sizes.Source: §1 · p. 3
- SmoothedGDA averages the primal iterates into an auxiliary sequence and performs GDA on an objective with a quadratic proximal term centered there.The paper presents this smoothing scheme as a way to address GDA oscillation.Source: §1 · p. 3
- O(1/ε^2) iterations suffice for an ε-stationary solution on the finite pointwise-maximum problem under the stated assumptions.For general nonconvex-concave problems, the reported complexity is O(1/ε^4), without requiring a compact domain X.Source: §1 · p. 3
- SmoothedGDA extends to multi-block problems with the same convergence guarantee as in the one-block case.The paper also describes the single-loop framework as suitable for large-scale and distributed settings.Source: §1 · p. 3
2 Representative Applications
The paper illustrates the min-max formulation through robust learning across distributions, fair models, and distribution-agnostic meta-learning. Each example frames training around losses across multiple distributions, categories, or tasks.
Robust learning from multiple distributions
- Robust learning across perturbed data distributions minimizes the maximum expected loss over those distributions.In adversarial learning, the distributions can represent those used to generate adversarial examples.Source: §2 · p. 4
Fair models
- Fair-model training can minimize the maximum loss across categories, with each category having a corresponding loss function.The passage motivates this formulation by the possibility that models discriminate against protected groups.Source: §2 · p. 4
Distribution-agnostic meta-learning
- Distribution-agnostic MAML chooses an initialization that minimizes the post-update objective across task losses, and is reported to improve worst-case performance over original MAML.The formulation evaluates each task loss after one gradient step from the initialization.Source: §2 · p. 4
3 Smoothed GDA Algorithm and Its Convergence
Smoothed-GDA modifies GDA with an auxiliary sequence and proximal term, then establishes convergence rates for general nonconvex-concave problems and a faster rate for minimizing a finite pointwise maximum. The section also extends the method to multi-block problems and gives conditions for its guarantees.
3 Smoothed GDA Algorithm and Its Convergence
- An ε-solution bounds the norms of primal and dual stationarity residuals by ε; a game stationary point has both residuals equal to zero.The residuals include subgradients of the indicator functions for the feasible sets.Source: §3 · p. 5
3.1 Smoothed Gradient Descent-Ascent (Smoothed-GDA)
- Smoothed-GDA applies alternating descent and ascent to a proximal objective centered at an auxiliary primal sequence, then averages that sequence toward the new primal iterate.The proximal coefficient p is positive, and β = 1 reduces the method to standard GDA.Source: §3.1 · p. 5
- Smoothed-BGDA alternates updates across primal variable blocks while retaining the single-block strategy for the dual and auxiliary variables.The method applies when the primal variable has a block structure.Source: §3.1 · pp. 5–6
3.2 Iteration Complexity for Nonconvex-Concave Problems
- O(ε^-4) iterations suffice for an ε-stationary solution in general nonconvex-concave problems for both one-block and multi-block algorithms.The one-block result uses β < 1/√T to obtain an O(T^-1/4)-stationary iterate; the multi-block result requires β ≤ ε².Source: §3.2 · p. 6
3.3 The Iteration Complexity of Achieving ǫ-Optimization-Stationary Point
- The same O(ε^-4) complexity holds for the paper’s optimization-stationary point, and the authors identify this as the only single-loop algorithm with that guarantee.The theorem uses the same parameter choices as Theorem 3.4.Source: §3.3 · p. 7
3.4 Convergence Results for Minimizing the Point-Wise Maximum of Finite Functions
- Every limit point of the one-block iterates is a solution of the finite pointwise-maximum problem under the stated assumptions.The theorem requires Assumptions 3.3 and 3.7, plus a bounded level set or bounded auxiliary iterates.Source: §3.4 · p. 8
- O(1/ε²) iterations suffice for an ε-stationary solution of the finite pointwise-maximum problem in both one-block and multi-block cases.The guarantee assumes Assumptions 3.3 and 3.7, plus a bounded level set or bounded auxiliary iterates; the multi-block case also requires β < ε².Source: §3.4 · p. 8
4 Proof Sketch
The proof uses a potential function and error bounds to establish convergence rates for the structured problem (1.2) and the general problem (1.1). The structured case has an O(1/ε²) rate, while the general case requires O(1/ε⁴) iterations.
4.1 The Potential Function and Basic Estimates
- The potential φt combines the primal, dual, and proximal functions and is bounded below; the proof seeks sufficient decrease at each iteration.φt = K(xt, zt; yt) − 2d(yt, zt) + 2P(zt).Source: §4.1 · p. 9
- For (1.2), the dual error bound supports sufficient potential decrease, yielding O(1/ε²) iterations and min-max limit points.The proof applies when the sufficient-decrease condition holds at every iteration.Source: §4.1 · pp. 9–10
- For (1.2), small enough β lets compactness of Y handle large residuals, while the dual error bound handles small residuals.Together these bounds establish sufficient decrease and yield Theorem 3.10.Source: §4.1 · pp. 9–10
- For general (1.1), a weaker nonhomogeneous error bound requires β proportional to ε², giving O(1/ε⁴) iterations.The paper notes that stronger bounds for specially structured problems may improve this complexity.Source: §4.1 · p. 10
4 Proof Sketch
- The proof bounds a negative term in the potential difference to establish per-iteration decrease, then uses that decrease to count iterations.The negative term is identified as the main proof difficulty.Source: §4.1 · p. 10
5 Numerical Results on Robust Neural Network Training
The robust-training experiment reports comparable or slightly better accuracy for formulation (1.2) than the other algorithms. Smoothed-GDA also reaches a loss below 0.2 in fewer epochs than the algorithm in [20].
5 Numerical Results on Robust Neural Network Training
- The experiment applies Smoothed-GDA to robust neural-network training on MNIST against adversarial attacks, using a nonconvex-concave approximation of the formulation.The approximation is used because the original formulation is difficult to solve directly.Source: §5 · p. 11
- Formulation (1.2) achieves comparable or slightly better performance than the other algorithms on test accuracy under FGSM and PGD attacks.Table 2 summarizes the accuracy comparison.Source: Table 2 · p. 11, §5 · p. 12
- 5 epochs bring Smoothed-GDA loss below 0.2, versus more than 14 epochs for the algorithm in [20], and its loss has smaller variance.Source: §5 · p. 12
6 Conclusion
The conclusion presents Smoothed-GDA as a single-loop method for nonconvex min-max problems, especially efficient for a family with a dual error bound. It also emphasizes the framework's potential for large-scale and distributed settings.
6 Conclusion
- The paper proposes a simple single-loop algorithm for nonconvex min-max problems.Source: §6 · p. 12
- A dual error bound makes the method more efficient for an important family of problems, and the paper describes it as well-suited to large dimensions and distributed settings.Source: §6 · p. 12
- The authors identify extending the framework to practical problems and deriving stronger error bounds for lower iteration complexity as future work.Source: §6 · p. 12
Broader Impact
The broader-impact discussion proposes possible benefits of applying the algorithm to training, while the appendices introduce notation and outline the one-block proof approach.
Broader Impact
- The authors suggest that faster training for minimax problems could reduce energy use, including in robust training, GANs, and AI for games.Source: Broader Impact · p. 13
- The paper argues that min-max training can promote fairness by preventing some objectives from dominating average-loss minimization when data are biased.Source: Broader Impact · p. 13
- The authors describe the multi-block framework as extendable to distributed data, decentralized control, and privacy-constrained settings.Source: Broader Impact · p. 13
A Notations
- The notation appendix defines index sets, solution sets, Euclidean balls and distances, vector components, and matrix row-submatrices.Source: A Notations · p. 18
- It also defines the smallest singular value of a matrix and projection of a point onto a set.Source: A Notations · p. 18
B Proof of the two main theorems: one-block case
- The one-block proof is presented first, with the multi-block proof deferred to the next section.Source: B Proof of the two main theorems:… · p. 18
Proof Sketch.
- The proof sketch establishes a lower-bounded potential, analyzes a negative term in its decrease, and uses sufficient decrease to derive the iteration count.Source: Proof Sketch · pp. 18–19
B.1 The potential function and basic estimate
This section develops the potential and error-bound estimates used in the convergence analysis, including complexity guarantees for the general nonconvex-concave case and potential decrease for the discrete max problem under strict complementarity. It also completes the discrete-max proof, outlines proof replacements for the multi-block case, and derives an estimate used in Lemma B.10.
B.1 The potential function and basic estimate
- The potential φ_t is bounded below, and Lemma B.2 provides error bounds with constants independent of y.The potential is defined as φ_t = Φ(x_t, z_t; y_t) = K(x_t, z_t; y_t) − 2d(y_t, z_t) + 2P(z_t).Source: §B.1 · pp. 19–20
- The dual function d(·, z) is differentiable in y with Lipschitz-continuous gradient, with Lipschitz constant L_d = L + Lσ_2.Source: §B.1 · p. 21
B.2 General nonconvex-concave case
- When (B.47) holds at every iteration, some iterate is a C/√(Tβ)-solution, and the iteration complexity for an ε-solution is O(1/ε^2).Under the same condition, every limit point of (x_t, y_t) is a solution of (1.2).Source: §B.2 · pp. 28–29
B.3 The max problem is over a discrete set
- Under strict complementarity and bounded level sets, the potential decreases sufficiently after every iteration in the discrete-set max problem.The appendix states this decrease as the basis for proving Theorem 3.10.Source: §B.3 · p. 31
- For the discrete-max case, every iterate z_t remains bounded by R(x_0, y_0, z_0), and (B.47) holds for all t, yielding the theorem via Lemma B.12.Source: §B.3 · p. 32
C The multi-block cases
- The multi-block proof replaces one-block ingredients with primal descent, an error bound, and Proposition C.3, using arguments similar to Lemmas 5.2 and 5.3 and Proposition B.4 in.The passage identifies Lemma B.5, (B.5), and Proposition 4.1 as the differences from the one-block proof.Source: C The multi-block cases · pp. 32–33
D.1 Proof of lemma B.10
- Lemma B.10's proof derives an estimate bounding α(p − L)∥x*(z) − x(y+(z), z)∥² by (1 + αL + αLσ²)∥y − y+(z)∥·dist(y+(z), Y(z)).The derivation uses strong convexity and the fact that y+(z) maximizes the stated problem.Source: §D.1 · pp. 33–34
D.2 Proof of Lemma B.13
The proof uses strict complementarity to establish convergence and uniqueness properties, then discusses a weaker regularity condition that preserves the theorem’s result and holds almost surely in a specified robust regression setting.
D.2 Proof of Lemma B.13
- Under strict complementarity, vanishing residuals imply convergence of iteration points toward a solution, with a convergent subsequence under the stated boundedness condition.The proof also derives local regularity and active-set properties used to establish its error bounds.Source: §D.2 · p. 35
- For any solution x* under strict complementarity, the associated y and the multipliers satisfying the KKT conditions are unique.The proof establishes full row rank of M(x) at the solution before deriving uniqueness.Source: §D.2 · p. 36
- When residuals are sufficiently small, the active sets of y+(z) and y(z) coincide, supporting the local comparison used in the proof.The proof also shows the top-coordinate sets for x*(z) and x(y+(z), z) agree under the stated conditions.Source: §D.2 · pp. 37–38
E Discussion of the strict complementarity condition
- Assumption E.1 is weaker than strict complementarity, preserves Theorem 3.10’s result, and holds with probability 1 for robust regression with jointly continuous data outputs.The robust-regression claim is for square loss.Source: E Discussion of the strict… · pp. 39–40
E.1 Replacing Assumption D.4 by Assumption E.1 in Theorem 3.10
The appendix shows that the results of Theorem 3.10 hold under the weaker Assumption E.1. It also establishes that this regularity condition holds almost surely for robust regression under specified data assumptions.
E.1 Replacing Assumption D.4 by Assumption E.1 in Theorem 3.10
- Under Assumption E.1, Algorithm 2 achieves O(1/ε²) iteration complexity and every limit point is a solution; the multi-block case has the same results under its stated parameter conditions.The multi-block result requires β < ε² and β < β′′.Source: §E.1 · p. 40
E.2 The rationality of Assumption E.1
- For robust regression with square loss, Assumption E.1 holds with probability 1 when the outputs are jointly drawn from a continuous distribution.The appendix also states that, under Assumption E.4, Assumption E.5 holds with probability 1 when m > n.Source: §E.2 · pp. 40–41
E.3 Proof of Lemma E.2 and Theorem E.3
This appendix develops the proof of Theorem E.3 under Assumption E.1 and proves the probability-one regularity claim for continuous-output robust regression using a measure-zero argument.
E.3 Proof of Lemma E.2 and Theorem E.3
- The proof bounds the relevant matrix singular values away from zero under Assumption E.1, then uses this regularity and Lipschitz continuity to establish the needed error bound.The resulting lemma is used with Lemma 4.3 to obtain Theorem E.3.Source: §E.3 · pp. 41, 43
E.4 Proof of Proposition E.7
- The probability-one argument partitions the possible active-index sets into finitely many cases and represents the corresponding data values through a smooth map.The map is defined in the proof as G(x1, …, xn+1).Source: §E.4 · p. 43
- The proof reduces failure of the regularity condition to a critical-value set of measure zero, then uses continuity of the output distribution to conclude probability zero.The measure-zero step invokes Sard’s Theorem [53].Source: §E.4 · p. 44
F Proof of Theorem 3.6
The proof combines cited inequalities and lemmas, sets β = ε², and then applies the triangle inequality.
- The proof invokes inequalities B.54 and B.62 together with weaker bounds from Lemma 4.3 and Lemma B.2.Source: F Proof of Theorem 3.6 · p. 44
- The proof takes β = ε² when applying those bounds.Source: F Proof of Theorem 3.6 · p. 44
- It then uses the triangle inequality.Source: F Proof of Theorem 3.6 · p. 44
G Details in Experiments
The experiments detail a robust-training formulation using targeted attacks and report Smoothed-GDA's CIFAR10 convergence-speed advantage over.
- The robust-training approximation replaces the inner maximization with a finite maximum over targeted attacks, then a concave simplex formulation.Source: G Details in Experiments · p. 45
- Targeted attacks use gradient ascent; the stated settings are K = 40, step size 0.01, and perturbation levels from 0.0 to 0.4.Source: G Details in Experiments · p. 45
- The MNIST experiments use the CNN architecture detailed in Table 3, with the setting matching.Source: G Details in Experiments · p. 45, Table 3 · p. 46
- For the Table 2 comparison, Algorithm 2 runs for 100 iterations with p = 0.2, β = 0.8, and α = 0.5; simplex projection costs O(d log d).Source: G Details in Experiments · p. 45
- 38.5% robust accuracy and 82.6% testing accuracy; after 30 epochs, Smoothed-GDA outperforms in CIFAR10 convergence speed.Source: G Details in Experiments · p. 45, Figure 2 · p. 46