Modified interior-point method for large-and-sparse low-rank semidefinite programs

Semidefinite programs (SDPs) are powerful theoretical tools that have been studied for over two decades, but their practical use remains limited due to computational difficulties in solving large-scale, realistic-sized problems. In this paper, we describe a modified interior-point method for the efficient solution of large-and-sparse low-rank SDPs, which finds applications in graph theory, approximation theory, control theory, sum-of-squares, etc. Given that the problem data is large-and-sparse, conjugate gradients (CG) can be used to avoid forming, storing, and factoring the large and fully-dense interior-point Hessian matrix, but the resulting convergence rate is usually slow due to ill-conditioning. Our central insight is that, for a rank-k, size-n SDP, the Hessian matrix is ill-conditioned only due to a rank-nk perturbation, which can be explicitly computed using a size-n eigendecomposition. We construct a preconditioner to “correct” the low-rank perturbation, thereby allowing preconditioned CG to solve the Hessian equation in a few tens of iterations. This modification is incorporated within SeDuMi, and used to reduce the solution time and memory requirements of large-scale matrix-completion problems by several orders of magnitude.

[1]  Anne Greenbaum,et al.  Iterative methods for solving linear systems , 1997, Frontiers in applied mathematics.

[2]  Martin S. Andersen,et al.  Chordal Graphs and Semidefinite Optimization , 2015, Found. Trends Optim..

[3]  Yurii Nesterov,et al.  Interior-point polynomial algorithms in convex programming , 1994, Siam studies in applied mathematics.

[4]  Javad Lavaei,et al.  ADMM for sparse semidefinite programming with applications to optimal power flow problem , 2015, 2015 54th IEEE Conference on Decision and Control (CDC).

[5]  Jos F. Sturm,et al.  A Matlab toolbox for optimization over symmetric cones , 1999 .

[6]  Iain S. Duff,et al.  MA57---a code for the solution of sparse symmetric definite and indefinite systems , 2004, TOMS.

[7]  Javad Lavaei,et al.  A fast distributed algorithm for decomposable semidefinite programs , 2015, 2015 54th IEEE Conference on Decision and Control (CDC).

[8]  Yong Gao,et al.  Treewidth of Erdős-Rényi random graphs, random intersection graphs, and scale-free random graphs , 2009, Discret. Appl. Math..

[9]  Stephen P. Boyd,et al.  Determinant Maximization with Linear Matrix Inequality Constraints , 1998, SIAM J. Matrix Anal. Appl..

[10]  Mohamadreza Ahmadi,et al.  Stability Analysis for a Class of Partial Differential Equations via Semidefinite Programming , 2016, IEEE Transactions on Automatic Control.

[11]  Yinyu Ye,et al.  An O(square-root[nL])-Iteration Homogeneous and Self-Dual Linear Programming Algorithm , 1992 .

[12]  Wotao Yin,et al.  Alternating direction augmented Lagrangian methods for semidefinite programming , 2010, Math. Program. Comput..

[13]  Kim-Chuan Toh,et al.  A Newton-CG Augmented Lagrangian Method for Semidefinite Programming , 2010, SIAM J. Optim..

[14]  Francis R. Bach,et al.  Low-Rank Optimization on the Cone of Positive Semidefinite Matrices , 2008, SIAM J. Optim..

[15]  Stephen P. Boyd,et al.  Conic Optimization via Operator Splitting and Homogeneous Self-Dual Embedding , 2013, Journal of Optimization Theory and Applications.

[16]  Michael L. Overton,et al.  Complementarity and nondegeneracy in semidefinite programming , 1997, Math. Program..

[17]  Richard Y. Zhang,et al.  Robust stability analysis for large-scale power systems , 2017 .

[18]  Jorge Nocedal,et al.  An Interior Point Algorithm for Large-Scale Nonlinear Programming , 1999, SIAM J. Optim..

[19]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[20]  Shinji Mizuno,et al.  An O(√nL)-Iteration Homogeneous and Self-Dual Linear Programming Algorithm , 1994, Math. Oper. Res..

[21]  Renato D. C. Monteiro,et al.  A nonlinear programming algorithm for solving semidefinite programs via low-rank factorization , 2003, Math. Program..

[22]  Emmanuel J. Candès,et al.  The Power of Convex Relaxation: Near-Optimal Matrix Completion , 2009, IEEE Transactions on Information Theory.

[23]  Kazuo Murota,et al.  Exploiting Sparsity in Semidefinite Programming via Matrix Completion I: General Framework , 2000, SIAM J. Optim..

[24]  Javad Lavaei,et al.  Exactness of Semidefinite Relaxations for Nonlinear Optimization Problems with Underlying Graph Structure , 2014, SIAM J. Optim..

[25]  Masakazu Kojima,et al.  Exploiting sparsity in linear and nonlinear matrix inequalities via positive semidefinite matrix completion , 2011, Math. Program..

[26]  Kim-Chuan Toh,et al.  Solving Some Large Scale Semidefinite Programs via the Conjugate Residual Method , 2002, SIAM J. Optim..

[27]  Pablo A. Parrilo,et al.  Semidefinite programming relaxations for semialgebraic problems , 2003, Math. Program..

[28]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[29]  Jacob K. White,et al.  GMRES-Accelerated ADMM for Quadratic Objectives , 2016, SIAM J. Optim..

[30]  T. Andô Concavity of certain maps on positive definite matrices and applications to Hadamard products , 1979 .

[31]  J. Lasserre Moments, Positive Polynomials And Their Applications , 2009 .

[32]  Gene H. Golub,et al.  Numerical solution of saddle point problems , 2005, Acta Numerica.

[33]  E. Yip A Note on the Stability of Solving a Rank-p Modification of a Linear System by the Sherman–Morrison–Woodbury Formula , 1986 .

[34]  Emmanuel J. Candès,et al.  Exact Matrix Completion via Convex Optimization , 2008, Found. Comput. Math..

[35]  P. Pardalos,et al.  Interior Point Methods for Combinatorial Optimization , 1998 .

[36]  Stephen J. Wright Primal-Dual Interior-Point Methods , 1997, Other Titles in Applied Mathematics.

[37]  Javad Lavaei,et al.  Convex Relaxation for Optimal Power Flow Problem: Mesh Networks , 2015, IEEE Transactions on Power Systems.