Certification of inequalities involving transcendental functions: Combining SDP and max-plus approximation

We consider the problem of certifying an inequality of the form f(x) ≥ 0, ∀x ∈ K, where f is a multivariate transcendental function, and K is a compact semialgebraic set. We introduce a certification method, combining semialgebraic optimization and max-plus approximation. We assume that f is given by a syntaxic tree, the constituents of which involve semialgebraic operations as well as some transcendental functions like cos, sin, exp, etc. We bound some of these constituents by suprema or infima of quadratic forms (max-plus approximation method, initially introduced in optimal control), leading to semialgebraic optimization problems which we solve by semidefinite relaxations. The max-plus approximation is iteratively refined and combined with branch and bound techniques to reduce the relaxation gap. Illustrative examples of application of this algorithm are provided, explaining how we solved tight inequalities issued from the Flyspeck project (one of the main purposes of which is to certify numerical inequalities used in the proof of the Kepler conjecture by Thomas Hales).

[1]  Makoto Yamashita,et al.  A high-performance software package for semidefinite programs: SDPA 7 , 2010 .

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

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

[4]  William M. McEneaney,et al.  A Max-Plus-Based Algorithm for a Hamilton--Jacobi--Bellman Equation of Nonlinear Filtering , 2000, SIAM J. Control. Optim..

[5]  Lawrence C. Paulson,et al.  MetiTarski: An Automatic Theorem Prover for Real-Valued Special Functions , 2010, Journal of Automated Reasoning.

[6]  Bin Li,et al.  Exact certification in global polynomial optimization via sums-of-squares of rational functions with rational coefficients , 2012, J. Symb. Comput..

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

[8]  William M. McEneaney,et al.  Max-plus methods for nonlinear control and estimation , 2005 .

[9]  G. Calafiore,et al.  Reduced Vertex Set Result for Interval Semidefinite Optimization Problems , 2008 .

[10]  John Harrison,et al.  Verifying Nonlinear Real Formulas Via Sums of Squares , 2007, TPHOLs.

[11]  Didier Henrion,et al.  GloptiPoly 3: moments, optimization and semidefinite programming , 2007, Optim. Methods Softw..

[12]  Jean B. Lasserre,et al.  Positivity and Optimization for Semi-Algebraic Functions , 2009, SIAM J. Optim..

[13]  Eldon R. Hansen Sharpening Interval Computations , 2006, Reliab. Comput..

[14]  Thomas C. Hales,et al.  Introduction to the Flyspeck Project , 2005, Mathematics, Algorithms, Proofs.

[15]  B. Borchers CSDP, A C library for semidefinite programming , 1999 .

[16]  W.M. McEneaney,et al.  Curse-of-complexity attenuation in the curse-of-dimensionality-free method for HJB PDEs , 2008, 2008 American Control Conference.

[17]  Nicholas I. M. Gould,et al.  Adaptive cubic regularisation methods for unconstrained optimization. Part I: motivation, convergence and numerical results , 2011, Math. Program..

[18]  Stéphane Gaubert,et al.  The Max-Plus Finite Element Method for Solving Deterministic Optimal Control Problems: Basic Properties and Convergence Analysis , 2008, SIAM J. Control. Optim..

[19]  Tobias Nipkow,et al.  A Revision of the Proof of the Kepler Conjecture , 2009, Discret. Comput. Geom..

[20]  E. Hansen,et al.  An interval Newton method , 1983 .

[21]  Frédéric Besson,et al.  Fast Reflexive Arithmetic Tactics the Linear Case and Beyond , 2006, TYPES.

[22]  William M. McEneaney,et al.  Curse of dimensionality reduction in max-plus based approximation methods: Theoretical estimates and improved pruning algorithms , 2011, IEEE Conference on Decision and Control and European Control Conference.

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

[24]  Pierre Corbineau,et al.  On the Generation of Positivstellensatz Witnesses in Degenerate Cases , 2011, ITP.

[25]  Thorsten Theobald,et al.  Exploiting Symmetries in SDP-Relaxations for Polynomial Optimization , 2011, Math. Oper. Res..

[26]  Victor Magron,et al.  Certification of Bounds of Non-linear Functions: The Templates Method , 2013, MKM/Calculemus/DML.

[27]  William M. McEneaney,et al.  Reduced-complexity numerical method for optimal gate synthesis , 2010, 1011.6013.

[28]  S. Gaubert,et al.  Set coverings and invertibility of functional Galois connections , 2004, math/0403441.

[29]  William M. McEneaney,et al.  A Curse-of-Dimensionality-Free Numerical Method for Solution of Certain HJB PDEs , 2007, SIAM J. Control. Optim..

[30]  Martin Berz,et al.  Rigorous global search using taylor models , 2009, SNC '09.

[31]  William M. McEneaney,et al.  Convergence Rate for a Curse-of-Dimensionality-Free Method for a Class of HJB PDEs , 2009, SIAM J. Control. Optim..