A SEQUENTIAL UPDATING SCHEME OF THE LAGRANGE MULTIPLIER FOR SEPARABLE CONVEX PROGRAMMING

被引:31
作者
Dai, Yu-Hong [1 ]
Han, Deren [2 ]
Yuan, Xiaoming [3 ]
Zhang, Wenxing [4 ]
机构
[1] Chinese Acad Sci, Inst Computat Math, LSEC, POB 2719, Beijing 100190, Peoples R China
[2] Nanjing Normal Univ, Jiangsu Key Lab NSLSCS, Sch Math Sci, Nanjing 210023, Jiangsu, Peoples R China
[3] Hong Kong Baptist Univ, Dept Math, Kowloon, Hong Kong, Peoples R China
[4] Univ Elect Sci & Technol China, Sch Math Sci, Chengdu 611731, Peoples R China
关键词
Convex programming; augmented Lagrangian method; splitting method; method of multipliers; image processing; ALTERNATING DIRECTION METHOD; TOTAL VARIATION MINIMIZATION; IMAGE DECOMPOSITION; CONVERGENCE RATE; ALGORITHMS;
D O I
10.1090/mcom/3104
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The augmented Lagrangian method (ALM) is a benchmark for solving convex minimization problems with linear constraints. Solving the augmented subproblems over the primal variables can be regarded as sequentially providing inputs for updating the Lagrange multiplier (i.e., the dual variable). We consider the separable case of a convex minimization problem where its objective function is the sum of more than two functions without coupled variables. When applying the ALM to this case, at each iteration we can (sometimes must) split the resulting augmented subproblem in order to generate decomposed subproblems which are often easy enough to have closed-form solutions. But the decomposition of primal variables only provides less accurate inputs for updating the Lagrange multiplier, and it points out the lack of convergence for such a decomposition scheme. To remedy this difficulty, we propose to update the Lagrange multiplier sequentially once each decomposed subproblem over the primal variables is solved. This scheme updates both the primal and dual variables in Gauss- Seidel fashion. In addition to the exact version which is useful enough for the case where the functions in the objective are all simple such that the decomposed subproblems all have closed- form solutions, we investigate an inexact version of this scheme which allows the decomposed subproblems to be solved approximately subject to certain inexactness criteria. Some preliminary numerical results when the proposed scheme is respectively applied to an image decomposition problem and an allocation problem are reported.
引用
收藏
页码:315 / 343
页数:29
相关论文
共 46 条
[11]   A GENERALIZED PROXIMAL POINT ALGORITHM AND ITS CONVERGENCE RATE [J].
Corman, Etienne ;
Yuan, Xiaoming .
SIAM JOURNAL ON OPTIMIZATION, 2014, 24 (04) :1614-1638
[12]  
Douglas J., 1956, T AM MATH SOC, V82, P421, DOI [10.1090/S0002-9947-1956-0084194-4, DOI 10.1090/S0002-9947-1956-0084194-4, 10.2307/1993056]
[13]  
Eckstein J, 2015, PAC J OPTIM, V11, P619
[14]  
Gabay D., 1976, Computers & Mathematics with Applications, V2, P17, DOI 10.1016/0898-1221(76)90003-1
[15]  
Gabay D., 1983, Studies in Mathematics and its Applications
[16]  
GLOWINSKI R, 1975, REV FR AUTOMAT INFOR, V9, P41
[17]  
Glowinski R., 1984, NUMERICAL METHODS NO
[18]  
GLOWINSKI R, 1989, SIAM STUDIES APPL MA, V0009
[19]  
Glowinski R., 2003, Numerical Methods for Scientific Computing. Variational Problems and Applications, P67
[20]  
GLOWINSKI R., 2014, COMPUTATIONAL METHOD, P59, DOI DOI 10.1007/978-94-017-9054-3_4