A New Sparse SOS Decomposition Algorithm Based on Term Sparsity

A new sparse SOS decomposition algorithm is proposed based on a new sparsity pattern, called cross sparsity patterns. The new sparsity pattern focuses on the sparsity of terms and thus is different from the well-known correlative sparsity pattern which focuses on the sparsity of variables though the sparse SOS decomposition algorithms based on these two sparsity patterns both take use of chordal extensions/chordal decompositions. Moreover, it is proved that the SOS decomposition obtained by the new sparsity pattern is always a refinement of the block-diagonalization obtained by the sign-symmetry method. %Because the new sparsity pattern covers more sparse polynomials than correlative sparsity pattern, Various experiments show that the new algorithm dramatically saves the computational cost compared to existing tools and can handle some really huge polynomials.

[1]  P. Parrilo Structured semidefinite programs and semialgebraic geometry methods in robustness and optimization , 2000 .

[2]  Mohab Safey El Din,et al.  On Exact Polya, Hilbert-Artin and Putinar's Representations , 2018, ArXiv.

[3]  Pinar Heggernes,et al.  Minimal triangulations of graphs: A survey , 2006, Discret. Math..

[4]  Etienne de Klerk,et al.  Solving sparse polynomial optimization problems with chordal structure using the sparse bounded-degree sum-of-squares hierarchy , 2020, Discret. Appl. Math..

[5]  Pablo A. Parrilo,et al.  Computing sum of squares decompositions with rational coefficients , 2008 .

[6]  B. Reznick Extremal PSD forms with few terms , 1978 .

[7]  J. Malick,et al.  Projection methods in conic optimization , 2011, 1103.1511.

[8]  J. Lofberg,et al.  YALMIP : a toolbox for modeling and optimization in MATLAB , 2004, 2004 IEEE International Conference on Robotics and Automation (IEEE Cat. No.04CH37508).

[9]  Jean B. Lasserre,et al.  Sparse-BSOS: a bounded degree SOS hierarchy for large scale polynomial optimization with sparsity , 2016, Mathematical Programming Computation.

[10]  B. Peyton,et al.  An Introduction to Chordal Graphs and Clique Trees , 1993 .

[11]  Lihong Zhi,et al.  Global optimization of polynomials using generalized critical values and sums of squares , 2010, ISSAC.

[12]  Juan Campos Salazar,et al.  A multigrid approach to SDP relaxations of sparse polynomial optimization problems , 2017 .

[13]  Barry W. Peyton,et al.  Maximum Cardinality Search for Computing Minimal Triangulations of Graphs , 2004, Algorithmica.

[14]  Johan Löfberg,et al.  Pre- and Post-Processing Sum-of-Squares Programs in Practice , 2009, IEEE Transactions on Automatic Control.

[15]  Yang Zheng,et al.  Exploiting Sparsity in the Coefficient Matching Conditions in Sum-of-Squares Programming Using ADMM , 2017, IEEE Control Systems Letters.

[16]  M. Golumbic Algorithmic graph theory and perfect graphs , 1980 .

[17]  Pablo A. Parrilo,et al.  Finding sparse, equivalent SDPs using minimal coordinate projections , 2015, 2015 54th IEEE Conference on Decision and Control (CDC).

[18]  Pablo A. Parrilo,et al.  Basis selection for SOS programs via facial reduction and polyhedral approximations , 2014, 53rd IEEE Conference on Decision and Control.

[19]  Liyun Dai,et al.  Smaller SDP for SOS decomposition , 2014, J. Glob. Optim..

[20]  Masakazu Muramatsu,et al.  SparsePOP: a Sparse Semidefinite Programming Relaxation of Polynomial Optimization Problems , 2005 .

[21]  Ken Ono,et al.  THEOREMS AND CONJECTURES INVOLVING ROOK POLYNOMIALS WITH REAL ROOTS , 1997 .

[22]  Martin Charles Golumbic Not So Perfect Graphs , 2004 .

[23]  P. Parrilo,et al.  Symmetry groups, semidefinite programs, and sums of squares , 2002, math/0211450.

[24]  Dimitris Bertsimas,et al.  An accelerated first-order method for solving SOS relaxations of unconstrained polynomial optimization problems , 2013, Optim. Methods Softw..

[25]  B. Reznick,et al.  Sums of squares of real polynomials , 1995 .

[26]  George E. Collins,et al.  Quantifier elimination for real closed fields by cylindrical algebraic decomposition , 1975 .

[27]  Masakazu Kojima,et al.  Sparsity in sums of squares of polynomials , 2005, Math. Program..

[28]  Masakazu Muramatsu,et al.  Sums of Squares and Semidefinite Programming Relaxations for Polynomial Optimization Problems with Structured Sparsity , 2004 .

[29]  Claus Scheiderer,et al.  Sums of squares of polynomials with rational coefficients , 2012, 1209.2976.

[30]  Jean B. Lasserre,et al.  Global Optimization with Polynomials and the Problem of Moments , 2000, SIAM J. Optim..

[31]  Masakazu Muramatsu,et al.  A facial reduction algorithm for finding sparse SOS representations , 2010, Oper. Res. Lett..

[32]  Peter J Seiler,et al.  SOSTOOLS: Sum of squares optimization toolbox for MATLAB , 2002 .

[33]  L. Rodman,et al.  Positive semidefinite matrices with a given sparsity pattern , 1988 .

[34]  James Demmel,et al.  Sparse SOS Relaxations for Minimizing Functions that are Summations of Small Polynomials , 2008, SIAM J. Optim..

[35]  Erich Kaltofen,et al.  A proof of the monotone column permanent (MCP) conjecture for dimension 4 via sums-of-squares of rational functions , 2009, SNC '09.

[36]  D. R. Fulkerson,et al.  Incidence matrices and interval graphs , 1965 .

[37]  Yang Zheng,et al.  Sparse sum-of-squares (SOS) optimization: A bridge between DSOS/SDSOS and SOS optimization for sparse polynomials , 2018, 2019 American Control Conference (ACC).

[38]  Timothy A. Davis,et al.  Algorithm 837: AMD, an approximate minimum degree ordering algorithm , 2004, TOMS.

[39]  Mohab Safey El Din,et al.  On Exact Polya and Putinar's Representations , 2018, ISSAC.