Convergence of the Surrogate Lagrangian Relaxation Method

Studies have shown that the surrogate subgradient method, to optimize non-smooth dual functions within the Lagrangian relaxation framework, can lead to significant computational improvements as compared to the subgradient method. The key idea is to obtain surrogate subgradient directions that form acute angles toward the optimal multipliers without fully minimizing the relaxed problem. The major difficulty of the method is its convergence, since the convergence proof and the practical implementation require the knowledge of the optimal dual value. Adaptive estimations of the optimal dual value may lead to divergence and the loss of the lower bound property for surrogate dual values. The main contribution of this paper is on the development of the surrogate Lagrangian relaxation method and its convergence proof to the optimal multipliers, without the knowledge of the optimal dual value and without fully optimizing the relaxed problem. Moreover, for practical implementations, a stepsizing formula that guarantees convergence without requiring the optimal dual value has been constructively developed. The key idea is to select stepsizes in a way that distances between Lagrange multipliers at consecutive iterations decrease, and as a result, Lagrange multipliers converge to a unique limit. At the same time, stepsizes are kept sufficiently large so that the algorithm does not terminate prematurely. At convergence, the lower-bound property of the surrogate dual is guaranteed. Testing results demonstrate that non-smooth dual functions can be efficiently optimized, and the new method leads to faster convergence as compared to other methods available for optimizing non-smooth dual functions, namely, the simple subgradient method, the subgradient-level method, and the incremental subgradient method.

[1]  T. Koopmans,et al.  Assignment Problems and the Location of Economic Activities , 1957 .

[2]  Yu. M. Ermol’ev Methods of solution of nonlinear extremal problems , 1966 .

[3]  N. Z. Shor The rate of convergence of the generalized gradient descent method , 1968 .

[4]  Boris Polyak Minimization of unsmooth functionals , 1969 .

[5]  John W. Dickey,et al.  Campus building arrangement using topaz , 1972 .

[6]  Arthur M. Geoffrion,et al.  Scheduling Parallel Production Lines with Changeover Costs: Practical Application of a Quadratic Assignment/LP Approach , 1976, Oper. Res..

[7]  Rainer E. Burkard,et al.  Entwurf von Schreibmaschinentastaturen mittels quadratischer Zuordnungsprobleme , 1977, Math. Methods Oper. Res..

[8]  A. N. Elshafei,et al.  Hospital Layout as a Quadratic Assignment Problem , 1977 .

[9]  Jakob Krarup,et al.  Computer-aided layout design , 1978 .

[10]  Jeffery L. Kennington,et al.  A generalization of Polyak's convergence result for subgradient optimization , 1987, Math. Program..

[11]  Nicos Christofides,et al.  An Exact Algorithm for the Quadratic Assignment Problem on a Tree , 1989, Oper. Res..

[12]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[13]  J. P. Kelly,et al.  Tabu search for the multilevel generalized assignment problem , 1995 .

[14]  Yurii Nesterov,et al.  New variants of bundle methods , 1995, Math. Program..

[15]  John E. Beasley,et al.  A genetic algorithm for the generalised assignment problem , 1997, Comput. Oper. Res..

[16]  M. Caramanis,et al.  Efficient Lagrangian relaxation algorithms for industry size job-shop scheduling problems , 1998 .

[17]  T. Ibaraki,et al.  A variable depth search algorithm with branching search for the generalized assignment problem , 1998 .

[18]  Jean-Louis Goffin,et al.  Convergence of a simple subgradient level method , 1999, Math. Program..

[19]  X. Zhao,et al.  Surrogate Gradient Algorithm for Lagrangian Relaxation , 1999 .

[20]  D. Bertsekas,et al.  Convergen e Rate of In remental Subgradient Algorithms , 2000 .

[21]  Dimitri P. Bertsekas,et al.  Incremental Subgradient Methods for Nondifferentiable Optimization , 2001, SIAM J. Optim..

[22]  Masafumi Yamashita,et al.  Independentand cooperative parallel search methods for the generalized assignment problem , 2003, Optim. Methods Softw..

[23]  Yixin Chen,et al.  Subgoal Partitioning and Global Search for Solving Temporal Planning Problems in Mixed Space , 2004, Int. J. Artif. Intell. Tools.

[24]  Yurii Nesterov,et al.  Smooth minimization of non-smooth functions , 2005, Math. Program..

[25]  Fred W. Glover,et al.  A path relinking approach with ejection chains for the generalized assignment problem , 2006, Eur. J. Oper. Res..

[26]  Feng Zhao,et al.  Payment cost minimization auction for deregulated electricity markets using surrogate optimization , 2006, IEEE Transactions on Power Systems.

[27]  P. Luh,et al.  On the Surrogate Gradient Algorithm for Lagrangian Relaxation , 2007 .

[28]  T. Chang Comments on “Surrogate Gradient Algorithm for Lagrangian Relaxation” , 2008 .

[29]  Pasquale Avella,et al.  A computational study of exact knapsack separation for the generalized assignment problem , 2010, Comput. Optim. Appl..

[30]  Lale Özbakir,et al.  Bees algorithm for generalized assignment problem , 2010, Appl. Math. Comput..

[31]  Mikhail A. Bragin,et al.  Payment cost minimization using Lagrangian relaxation and modified surrogate optimization approach , 2011, 2011 IEEE Power and Energy Society General Meeting.

[32]  Dimitri P. Bertsekas,et al.  Incremental proximal methods for large scale convex optimization , 2011, Math. Program..

[33]  Jacques A. Ferland,et al.  An exact method with variable fixing for solving the generalized assignment problem , 2012, Comput. Optim. Appl..

[34]  Mikhail A. Bragin,et al.  An efficient surrogate subgradient method within Lagrangian relaxation for the Payment Cost Minimization problem , 2012, 2012 IEEE Power and Energy Society General Meeting.

[35]  Dimitri P. Bertsekas,et al.  Incremental Gradient, Subgradient, and Proximal Methods for Convex Optimization: A Survey , 2015, ArXiv.