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 条
[1]  
[Anonymous], 1971, OPTIMIZING METHODS S
[2]  
[Anonymous], 2003, SPRINGER SERIES OPER, DOI DOI 10.1007/978-0-387-21815-16
[3]  
[Anonymous], SOC IND APPL MATH SI
[4]   Structure-texture image decomposition - Modeling, algorithms, and parameter selection [J].
Aujol, JF ;
Gilboa, G ;
Chan, T ;
Osher, S .
INTERNATIONAL JOURNAL OF COMPUTER VISION, 2006, 67 (01) :111-136
[5]  
Bauschke HH, 2011, CMS BOOKS MATH, P1, DOI 10.1007/978-1-4419-9467-7
[6]  
Bertsekas, 1982, COMPUTER SCI APPL MA
[7]  
Bertsekas D. P., 1996, NEURODYNAMIC PROGRAM
[8]  
Boyd S., 2011, FOUND TRENDS MACH LE, V3, P1, DOI DOI 10.1561/2200000016
[9]   The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent [J].
Chen, Caihua ;
He, Bingsheng ;
Ye, Yinyu ;
Yuan, Xiaoming .
MATHEMATICAL PROGRAMMING, 2016, 155 (1-2) :57-79
[10]   Proximal Splitting Methods in Signal Processing [J].
Combettes, Patrick L. ;
Pesquet, Jean-Christophe .
FIXED-POINT ALGORITHMS FOR INVERSE PROBLEMS IN SCIENCE AND ENGINEERING, 2011, 49 :185-+