We consider a separable convex minimization model whose variables are coupled by linear constraints and they are subject to the positive orthant constraints, and its objective function is in form of $m$ functions without coupled variables. It is well recognized that when the augmented Lagrangian method (ALM) is applied to solve some concrete applications, the resulting subproblem at each iteration should be decomposed to generate solvable subproblems. When the Gauss-Seidel decomposition is implemented, this idea has inspired the alternating direction method of multiplier (for $m=2$) and its variants (for $m\ge 3$). When the Jacobian decomposition is considered, it has been shown that the ALM with Jacobian decomposition in its subproblem is not necessarily convergent even when $m=2$ and it was suggested to regularize the decomposed subproblems with quadratic proximal terms to ensure the convergence. In this paper, we focus on the multiple-block case with $m\ge 3$. We consider implementing the full Jacobian decomposition to ALM’s subproblems and using the logarithmic-quadratic proximal (LQP) terms to regularize the decomposed subproblems. The resulting subproblems are all unconstrained minimization problems because the positive orthant constraints are all inactive; and they are fully eligible for parallel computation. Accordingly, the ALM with full Jacobian decomposition and LQP regularization is proposed. We also consider its inexact version which allows the subproblems to be solved inexactly. For both the exact and inexact versions, we comprehensively discuss their convergence, including their global convergence, worst-case convergence rates measured by the iteration-complexity in both the ergodic and nonergodic senses, and linear convergence rates under additional assumptions. Some preliminary numerical results are reported to demonstrate the efficiency of the ALM with full Jacobian decomposition and LQP regularization.

DOI: https://doi.org/10.5802/smai-jcm.30

Classification: 90C25, 90C33, 65K05

Keywords: convex programming, splitting methods, augmented Lagrangian method, logarithmic-quadratic proximal, parallel computation, convergence rate

@article{SMAI-JCM_2018__4__81_0, author = {Min Li and Xiaoming Yuan}, title = {The augmented {Lagrangian} method with full {Jacobian} decomposition and logarithmic-quadratic proximal regularization for multiple-block separable convex programming}, journal = {The SMAI journal of computational mathematics}, pages = {81--120}, publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles}, volume = {4}, year = {2018}, doi = {10.5802/smai-jcm.30}, mrnumber = {3774649}, zbl = {1418.90197}, language = {en}, url = {https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.30/} }

Min Li; Xiaoming Yuan. The augmented Lagrangian method with full Jacobian decomposition and logarithmic-quadratic proximal regularization for multiple-block separable convex programming. The SMAI journal of computational mathematics, Volume 4 (2018) , pp. 81-120. doi : 10.5802/smai-jcm.30. https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.30/

[1] Entropic proximal decomposition methods for convex programs and variational inequalities, Math. Program., Volume 91 (2001), pp. 33-47 | Article | MR 1865260 | Zbl 1051.90017

[2] Interior projection-like methods for monotone variational inequalities, Math. Program., Volume 104 (2005), pp. 39-68 | Article | MR 2167425 | Zbl 1159.90517

[3] A logarithmic-quadratic proximal method for variational inequalities, Comput. Optim. Appl., Volume 12 (1999), pp. 31-40 | Article | MR 1704099 | Zbl 1039.90529

[4] Constrained Optimization and Lagrange Multiplier Methods, Academic Press, New York, 1982 | Zbl 0572.90067

[5] Distributed optimization and statistical learning via the alternating direction method of multipliers, Foundations and Trends in Machine Learning, Volume 3 (2010), pp. 1-122 | Article | Zbl 1229.90122

[6] The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent, Math. Program., Ser A, Volume 155 (2016), pp. 57-79 | Article | MR 3439797 | Zbl 1332.90193

[7] Further study on the convergence rate of alternating direction method of multipliers with logarithmic-quadratic proximal regularization, J. Optim. Theory Appli., Volume 166 (2015), pp. 906-929 | Article | MR 3375607 | Zbl 1327.90199

[8] Proximal splitting methods in signal processing, Fixed-Point Algorithms for Inverse Problems in Science and Engineering (H. H. Bauschke; R. S. Burachik; P. L. Combettes; V. Elser; D. R. Luke; H. Wolkowicz, eds.), Springer, 2011, pp. 185-212 | Article | Zbl 1242.90160

[9] Convergence rate analysis of several splitting schemes, Splitting Methods in Communication, Imaging, Science, and Engineering (R. Glowinski; S. J. Osher; W. T. Yin, eds.), Springer, 2017, pp. 115-163

[10] Parallel multi-block ADMM with $o(1/k)$ convergence, J. Sci. Comput., Volume 71 (2017), pp. 712-736 | Article | MR 3627538 | Zbl 1398.65121

[11] On the basic theorem of complementarity, Math. Program., Volume 1 (1971), pp. 68-75 | Article | MR 287901 | Zbl 0227.90044

[12] On the Douglas-Rachford splitting method and the proximal point algorithm for maximal monotone operators, Math. Program., Volume 55 (1992), pp. 293-318 | Article | MR 1168183 | Zbl 0765.90073

[13] Augmented Lagrangian and alternating direction methods for convex optimization: A tutorial and some illustrative computational results (2012) (RUTCOR Research Report RRR 32-2012)

[14] Finite-Dimensional Variational Inequalities and Complementarity Problems. Vol. I. Springer Series in Operations Research, Springer-Verlag, New York, 2003

[15] Numerical Methods for Nonlinear Variational Problems, Springer-Verlag, New York, Berlin, Heidelberg, Tokyo, 1984 | Article

[16] On alternating direction methods of multipliers: A historical perspective, Modeling, Simulation and Optimization for Science and Technology, Volume 34 of the series Computational Methods in Applied Sciences, Springer, 2014, pp. 59-82 | Zbl 1320.65098

[17] Approximation par éléments finis d’ordre un et résolution par pénalisation-dualité d’une classe de problèmes non linéaires, R.A.I.R.O., Volume R2 (1975), pp. 41-76 | Numdam | Zbl 0368.65053

[18] Augmented Lagrangian and Operator-Splitting Methods in Nonlinear Mechanics, SIAM Studies in Applied Mathematics, Philadelphia, PA, 1989 | Article | Zbl 0698.73001

[19] Modified Lagrangian in convex programming and their generalizations, Math. Program. Studies, Volume 10 (1979), pp. 86-97 | Article | MR 527059 | Zbl 0404.90069

[20] An augmented-Lagrangian-based parallel splitting method for separable convex programming with applications to image processing, Math. Comput., Volume 83 (2014), pp. 2263-2291 | Article

[21] Deblurring Images: Matrices, Spectra, and Filtering, SIAM, Philadelphia, 2006 | Zbl 1112.68127

[22] On full Jacobian decomposition of the augmented Lagrangian method for separable convex programming, SIAM J. Optim., Volume 25 (2015), pp. 2274-2312 | MR 3424071 | Zbl 1327.90209

[23] A strictly contractive Peaceman-Rachford splitting method for convex programming, SIAM J. Optim., Volume 24 (2014), pp. 1101-1140 | MR 3231988 | Zbl 1327.90210

[24] Alternating direction method with Gaussian back substitution for separable convex programming, SIAM J. Optim., Volume 22 (2012), pp. 313-340 | MR 2968856 | Zbl 1273.90152

[25] On the proximal Jacobian decomposition of ALM for multiple-block separable convex minimization problems and its relationship to ADMM, J. Sci. Comput., Volume 66 (2016), pp. 1204-1217 | MR 3456970 | Zbl 1371.65052

[26] On the O($1/n$) convergence rate of Douglas-Rachford alternating direction method, SIAM J. Numer. Anal., Volume 50 (2012), pp. 700-709 | MR 2914282 | Zbl 1245.90084

[27] On nonergodic convergence rate of Douglas-Rachford alternating direction method of multipliers, Numer. Math., Volume 130 (2015), pp. 567-577

[28] On the linear convergence of the alternating direction method of multipliers, Math. Program., Volume 162 (2017), pp. 165-199 | Article | MR 3612937

[29] Inexact alternating direction method of multipliers with logarithmic-quadratic proximal regularization, J. Optim. Theory Appli., Volume 159 (2013), pp. 412-436 | Article | MR 3118276

[30] A majorized ADMM with indefinite proximal terms for linearly constrained convex composite optimization, SIAM J. Optim., Volume 26 (2016), pp. 922-950 | Article | MR 3484404 | Zbl 1338.90305

[31] A strictly contractive Peaceman-Rachford splitting method with logarithmic-quadratic proximal regularization for convex programming, Math. Oper. Res., Volume 40 (2015), pp. 842-858 | Article | MR 3423741 | Zbl 1329.90106

[32] On the global linear convergence of the ADMM with multiblock variables, SIAM J. Optim., Volume 25 (2015), pp. 1478-1497 | Article | MR 3374652

[33] On the sublinear convergence rate of multi-block ADMM, J. Oper. Res. Soc. China, Volume 3 (2015), pp. 251-274 | Article | MR 3394471 | Zbl 1323.90052

[34] Regularization d’inequations variationelles par approximations successives, Revue Francaise d’Informatique et de Recherche Opérationelle, Volume 4 (1970), pp. 154-159 | Numdam | MR 298899 | Zbl 0215.21103

[35] Iteration-complexity of a Jacobi-type non-Euclidean ADMM for multi-block linearly constrained nonconvex programs (2017) (arXiv:1705.07229v1)

[36] Problem Complexity and Method Efficiency in Optimization, Wiley-Interscience Series in Discrete Mathematics, John Wiley & Sons, New York, 1983

[37] A method for unconstrained convex minimization problem with the rate of convergence O($1/{k}^{2}$), Doklady AN SSSR, Volume 269 (1983), pp. 543-547

[38] Gradient methods for minimizing composite objective function, Math. Program., Ser. B, Volume 140 (2013), pp. 125-161 | Article

[39] A survey on the continuous nonlinear resource allocation problem, European J. Oper. Res., Volume 185 (2008), pp. 1-46 | Article | MR 2358684 | Zbl 1146.90493

[40] Robust alignment by sparse and low-rank decomposition for linearly correlated images, IEEE Tran. Pattern Anal. Mach. Intel., Volume 34 (2012), pp. 2233-2246 | Article

[41] A method for nonlinear constraints in minimization problems, Optimization (R. Fletcher, ed.), Academic Press, New York, 1969, pp. 283-298 | Zbl 0194.47701

[42] Augmented Lagrangians and applications of the proximal point algorithm in convex programming, Math. Oper. Res., Volume 1 (1976), pp. 97-116 | Article | MR 418919 | Zbl 0402.90076

[43] Recovering low-rank and sparse components of matrices from incomplete and noisy observations, SIAM J. Optim., Volume 21 (2011), pp. 57-81 | Article | MR 2765489 | Zbl 1218.90115

[44] On the O$(1/t)$ convergence rate of alternating direction method with logarithmic-quadratic proximal regularization, SIAM J. Optim., Volume 22 (2012), pp. 1431-1448 | Article | MR 3023778 | Zbl 1263.90103

[45] Market mechanisms and mathematical programming, Econometrica, Volume 28 (1960), pp. 872-881 | Article | MR 136447 | Zbl 0098.33602

[46] An LQP-based decomposition method for solving a class of variational inequalities, SIAM J. Optim, Volume 21 (2011), pp. 1309-1318 | Article | MR 2854585