Iteration complexity of inexact augmented Lagrangian methods for constrained convex programming
Yangyang Xu
DOI: https://doi.org/10.1007/s10107-019-01425-9
IF: 2.7
2019-08-21
Mathematical Programming
Abstract:Augmented Lagrangian method (ALM) has been popularly used for solving constrained optimization problems. Practically, subproblems for updating primal variables in the framework of ALM usually can only be solved inexactly. The convergence and local convergence speed of ALM have been extensively studied. However, the global convergence rate of the inexact ALM is still open for problems with nonlinear inequality constraints. In this paper, we work on general convex programs with both equality and inequality constraints. For these problems, we establish the global convergence rate of the inexact ALM and estimate its iteration complexity in terms of the number of gradient evaluations to produce a primal and/or primal-dual solution with a specified accuracy. We first establish an ergodic convergence rate result of the inexact ALM that uses constant penalty parameters or geometrically increasing penalty parameters. Based on the convergence rate result, we then apply Nesterov's optimal first-order method on each primal subproblem and estimate the iteration complexity of the inexact ALM. We show that if the objective is convex, then <span class="InlineEquation">\(O(\varepsilon ^{-1})\)</span> gradient evaluations are sufficient to guarantee a primal <span class="InlineEquation">\(\varepsilon \)</span>-solution in terms of both primal objective and feasibility violation. If the objective is strongly convex, the result can be improved to <span class="InlineEquation">\(O(\varepsilon ^{-\frac{1}{2}}|\log \varepsilon |)\)</span>. To produce a primal-dual <span class="InlineEquation">\(\varepsilon \)</span>-solution, more gradient evaluations are needed for convex case, and the number is <span class="InlineEquation">\(O(\varepsilon ^{-\frac{4}{3}})\)</span>, while for strongly convex case, the number is still <span class="InlineEquation">\(O(\varepsilon ^{-\frac{1}{2}}|\log \varepsilon |)\)</span>. Finally, we establish a nonergodic convergence rate result of the inexact ALM that uses geometrically increasing penalty parameters. This result is established only for the primal problem. We show that the nonergodic iteration complexity result is in the same order as that for the ergodic result. Numerical experiments on quadratically constrained quadratic programming are conducted to compare the performance of the inexact ALM with different settings.
mathematics, applied,computer science, software engineering,operations research & management science