Exact algorithms for linear matrix inequalities
|
|
- Daniella Lee
- 6 years ago
- Views:
Transcription
1 0 / 12 Exact algorithms for linear matrix inequalities Mohab Safey El Din (UPMC/CNRS/IUF/INRIA PolSys Team) Joint work with Didier Henrion, CNRS LAAS (Toulouse) Simone Naldi, Technische Universität Dortmund SMAI-MODE, 2016
2 1 / 12 Spectrahedra and LMI A 0, A 1,..., A n are m m real symmetric matrices Spectrahedron: S = {x R n : A(x) = A 0 + x 1 A x na n 0} It is basic semi-algebraic since, if det(a(x) + ti m) = f m(x) + f m 1 (x)t + + f 1 (x)t m 1 + t m then S = {x R n : f i (x) 0, i = 1,..., m}. A(x) 0 is called an LMI.
3 1 / 12 Spectrahedra and LMI A 0, A 1,..., A n are m m real symmetric matrices Spectrahedron: S = {x R n : A(x) = A 0 + x 1 A x na n 0} It is basic semi-algebraic since, if det(a(x) + ti m) = f m(x) + f m 1 (x)t + + f 1 (x)t m 1 + t m then S = {x R n : f i (x) 0, i = 1,..., m}. A(x) 0 is called an LMI. SDP : linear optimization over S (i.e. over LMI) 1 x 1 x 2 S = (x 1, x 2, x 3 ) R 3 : x 1 1 x 3 0 x 2 x 3 1 Figure: The Cayley spectrahedron
4 2 / 12 Why exact algorithms? 1. It is Hard to compute low-rank solutions to SDP Figure: Low-rank points : they minimize a cone of linear forms Figure: SEDUMI returns a floating point approximation of (0, 0) when maximizing x 2 2. The interior of S can be empty Interior point algorithms could fail 1 0 x 1 (1 x 2 4) x 1 x 2 x (1 x 2 4) x 3 x 4
5 2 / 12 Why exact algorithms? 1. It is Hard to compute low-rank solutions to SDP 2. The interior of S can be empty Interior point algorithms could fail x 1 (1 x 2 4) x 1 x 2 x (1 x 2 4) x 3 x Main motivations for the design of exact algorithms: 1. Can we manage algebraic constraints such as rank defects? 2. Can we handle degenerate non-full-dimensional examples? 3. Consequence: The output is a point whose coordinates may be real algebraic numbers {( q1 (t) (q, q 0, q 1,..., q n) Q[t] q 0 (t),..., qn(t) ) } : q(t) = 0 q 0 (t)
6 3 / 12 State of the art Decision / Sampling problem for real algebraic or semi-algebraic sets Cylindrical Algebraic Decomposition Tarski (1948), Seidenberg, Cohen,... Collins (1975) in O((2m) 22n+8 m 2n+6 ),... Critical Points Method local extrema of algebraic maps f on S Grigoriev, Vorobjov (1988) first singly exp: m O(n2 ) Renegar (1992), Heintz Roy Solernó (1989,1993), Basu Pollack Roy (1996,... ) linear exponent m O(n) Polar varieties local extrema of linear projections π on S Bank, Giusti, Heintz, Mbakop, Pardo (1997,... ) Safey El Din, Schost (2003,2004) regular in O(m 3n ), singular in O(m 4n ) The goal was: Better results for spectrahedra? How to take advantage of the structure?
7 4 / 12 Complexity of SDP Special case of SDP Khachiyan, Porkolab (1996) decide LMI-feasibility in time O(nm 4 ) + m O(min{n,m2 }) on (lm O(min{n,m2 }) )-bit numbers (l = input bit-size)
8 4 / 12 Complexity of SDP Special case of SDP Khachiyan, Porkolab (1996) decide LMI-feasibility in time O(nm 4 ) + m O(min{n,m2 }) on (lm O(min{n,m2 }) )-bit numbers (l = input bit-size) Positive aspects: 1. No assumptions, Deterministic 2. Binary complexity
9 4 / 12 Complexity of SDP Special case of SDP Khachiyan, Porkolab (1996) decide LMI-feasibility in time O(nm 4 ) + m O(min{n,m2 }) on (lm O(min{n,m2 }) )-bit numbers (l = input bit-size) Positive aspects: 1. No assumptions, Deterministic 2. Binary complexity Main drawbacks: 1. It relies on Quantifier Elimination 2. Too large constant in the exponent
10 5 / 12 Low rank positive semidefinite matrices Define: For any A(x) (not nec. symmetric): D r = {x C n : rank A(x) r} For A(x) symmetric, and S : r(a) = min{rank A(x) x S } So one has nested sequences D 0 D m 1 D 0 R n D m 1 R n
11 5 / 12 Low rank positive semidefinite matrices Define: For any A(x) (not nec. symmetric): D r = {x C n : rank A(x) r} For A(x) symmetric, and S : r(a) = min{rank A(x) x S } So one has nested sequences D 0 D m 1 D 0 R n D m 1 R n Smallest Rank Property Henrion-N.-Safey El Din 2015 A(x) symmetric, and S. Let C be a conn. comp. of D r(a) R n s.t. C S. Then C S. In particular C D r(a) \ D r(a) 1.
12 5 / 12 Low rank positive semidefinite matrices Define: For any A(x) (not nec. symmetric): D r = {x C n : rank A(x) r} For A(x) symmetric, and S : r(a) = min{rank A(x) x S } So one has nested sequences D 0 D m 1 D 0 R n D m 1 R n Smallest Rank Property Henrion-N.-Safey El Din 2015 A(x) symmetric, and S. Let C be a conn. comp. of D r(a) R n s.t. C S. Then C S. In particular C D r(a) \ D r(a) 1. Either the set S is empty
13 5 / 12 Low rank positive semidefinite matrices Define: For any A(x) (not nec. symmetric): D r = {x C n : rank A(x) r} For A(x) symmetric, and S : r(a) = min{rank A(x) x S } So one has nested sequences D 0 D m 1 D 0 R n D m 1 R n Smallest Rank Property Henrion-N.-Safey El Din 2015 A(x) symmetric, and S. Let C be a conn. comp. of D r(a) R n s.t. C S. Then C S. In particular C D r(a) \ D r(a) 1. Either the set S is empty Or r(a) is well-defined C D r(a) : C S
14 5 / 12 Low rank positive semidefinite matrices Smallest Rank Property Henrion-N.-Safey El Din 2015 A(x) symmetric, and S. Let C be a conn. comp. of D r(a) R n s.t. C S. Then C S. In particular C D r(a) \ D r(a) 1. Either the set S is empty Or r(a) is well-defined C D r(a) : C S
15 6 / 12 Problem statement Emptiness of spectrahedra Given A(x) symmetric, with entries in Q, compute a finite set meeting S = {x R n : A(x) 0}, or establish that S is empty. In other words: Decide the feasibility of an LMI A(x) 0. Particular instance of: Decide the emptiness of semi-algebraic sets.
16 6 / 12 Problem statement Emptiness of spectrahedra Given A(x) symmetric, with entries in Q, compute a finite set meeting S = {x R n : A(x) 0}, or establish that S is empty. In other words: Decide the feasibility of an LMI A(x) 0. Particular instance of: Decide the emptiness of semi-algebraic sets. Sample points on S + Smallest Rank Property = Sample points on D r(a) R n
17 Problem statement Emptiness of spectrahedra Given A(x) symmetric, with entries in Q, compute a finite set meeting S = {x R n : A(x) 0}, or establish that S is empty. In other words: Decide the feasibility of an LMI A(x) 0. Particular instance of: Decide the emptiness of semi-algebraic sets. Sample points on S + Smallest Rank Property = Sample points on D r(a) R n Real root finding on determinantal varieties Given any A(x) with entries in Q, compute a finite set meeting each connected component of D r R n = {x R n : rank A(x) r}. Particular instance of: Sampling real algebraic sets. 6 / 12
18 7 / 12 Strategy 1. The Smallest Rank Property ( C D r(a) : C S ) allows to reduce: Sampling/Optimization over One semi algebraic set Sampling/Optimization over Many algebraic sets This is somehow typical in PO. Ex. Polar Varieties for PO: Safey El Din, Greuet
19 7 / 12 Strategy 1. The Smallest Rank Property ( C D r(a) : C S ) allows to reduce: Sampling/Optimization over One semi algebraic set Sampling/Optimization over Many algebraic sets This is somehow typical in PO. Ex. Polar Varieties for PO: Safey El Din, Greuet But in our case the structure is preserved!
20 7 / 12 Strategy 1. The Smallest Rank Property ( C D r(a) : C S ) allows to reduce: Sampling/Optimization over One semi algebraic set Sampling/Optimization over Many algebraic sets This is somehow typical in PO. Ex. Polar Varieties for PO: Safey El Din, Greuet But in our case the structure is preserved! 2. For r = 1,..., m 1 compute sample points in D r R n 3. Output the minimum rank on S with a sample point.
21 7 / 12 Strategy 1. The Smallest Rank Property ( C D r(a) : C S ) allows to reduce: Sampling/Optimization over One semi algebraic set Sampling/Optimization over Many algebraic sets This is somehow typical in PO. Ex. Polar Varieties for PO: Safey El Din, Greuet But in our case the structure is preserved! 2. For r = 1,..., m 1 compute sample points in D r R n 3. Output the minimum rank on S with a sample point. Sampling determinantal varieties Either the empty list iff D r R n = Or (q, q 1,..., q n) Q[t] s.t. C D r R n t : x(t) C
22 7 / 12 Strategy 1. The Smallest Rank Property ( C D r(a) : C S ) allows to reduce: Sampling/Optimization over One semi algebraic set Sampling/Optimization over Many algebraic sets This is somehow typical in PO. Ex. Polar Varieties for PO: Safey El Din, Greuet But in our case the structure is preserved! 2. For r = 1,..., m 1 compute sample points in D r R n 3. Output the minimum rank on S with a sample point. Sampling determinantal varieties Either the empty list iff D r R n = Or (q, q 1,..., q n) Q[t] s.t. C D r R n t : x(t) C Emptiness of spectrahedra Either the empty list iff S = Or (q, q 1,..., q n) Q[t] s.t. t : x(t) S
23 8 / 12 Incidence varieties and critical points 1st step Lifting of the determinantal variety: y 1,1... y 1,m r A(x) Y (y) = A(x).. = 0. π1 (x2, y) C π1 y m,1... y m,m r x1 U Y (y) = I m r If A is generic, the lifted algebraic set V r is smooth and equidimensional 2nd step Compute critical points of the map π(x, y) = a 1 x a nx n on V r : When a 1.. a n are generic, there are finitely many critical points. 3rd step Intersect with any fiber of π and call point 2
24 8 / 12 Incidence varieties and critical points 1st step Lifting of the determinantal variety: y 1,1... y 1,m r A(x) Y (y) = A(x).. = 0. π1 (x2, y) C π 1 1 y m,1... y m,m r x1 U Y (y) = I m r If A is generic, the lifted algebraic set V r is smooth and equidimensional 2nd step Compute critical points of the map π(x, y) = a 1 x a nx n on V r : When a 1.. a n are generic, there are finitely many critical points. 3rd step Intersect with any fiber of π and call point 2
25 9 / 12 Computing critical points on incidence varieties A(X)Y = 0, UY = I m r F 1 (X, Y ) = = F m(m r) (Y ) = = F m(m r)+(m r) 2 = 0 Lagrange System L.jac(F, [X, Y]) = [a 0]
26 Computing critical points on incidence varieties Multi-Linear System F(X, A(X)Y Y) = = 0, 0, G(X, UY L) = = I m r 0, H(Y, L) = 0 All equations have multi-degree (1, 1, 0) or (1, 0, 1) or (0, 1, 1) F 1 (X, Y ) = Impact = F m(m r) on multi-linear/sparsity (Y ) = = F m(m r)+(m r) 2 = 0 structure on the number of solutions Multi-linear Bézout bounds Symbolic-homotopy algorithms Lagrange System cubicl.jac(f, complexity [X, in Y]) these = [abounds 0] Symbolic Newton iteration (Hensel lifting) Jeronimo/Matera/Solerno/Waissbein 9 / 12
27 Complexity bounds Complexity for Sampling determinantal varieties ( ) 6 O (n + m 2 r 2 ) 7 n + m(m r) n Complexity for Emptiness of spectrahedra O n ( ) ( ) 6 m (n + p r + r(m r)) 7 p r + n r n r r(a) O (k) = O(k log c k) c N with p r = (m r)(m + r + 1)/2. Remarkable aspects: Explicit constants in the exponent When m is fixed, polynomial in n Strictly depends on r(a) 10 / 12
28 11 / 12 SPECTRA: a library for real algebraic geometry and optimization What is SPECTRA? A MAPLE library, freely distributed Depends on Faugère s FGB for computations with Gröbner bases Addressed to researchers in Optimization, Convex alg. geom., Symb. comp. (m, r, n) RAGLIB SPECTRA deg (3, 2, 8) (3, 2, 9) (4, 2, 5) (4, 2, 6) (4, 2, 7) (4, 2, 8) (4, 2, 9) (4, 3, 10) (4, 3, 11) (5, 2, 9) (6, 5, 4) RAGLIB = Real algebraic geometry library - SPECTRA = new algorithms - deg = degree of Rational Parametrization - Time in seconds - = more than 2 days Download a beta version: homepages.laas.fr/snaldi/software.html
29 12 / 12 Scheiderer s spectrahedron f = u u 1 u u 4 2 3u 2 1u 2 u 3 4u 1 u 2 2u 3 + 2u 2 1u u 1 u u 2 u u 4 3 One can write f = v A(x)v with v = [u 2 1, u 1 u 2, u 2 2, u 1 u 3, u 2 u 3, u 2 3] 1 0 x 1 0 3/2 x 2 x 3 0 2x 1 1/2 x 2 2 x 4 x 5 x A(x) = 1 1/2 1 x 4 0 x 6 0 x 2 x 4 2x x 5 3/2 x 2 2 x 4 0 x 5 2x 6 1/2 x 3 x 5 x 6 1/2 1/2 1 What information can be extracted? No matrices of rank 1 s.t. A(x) 0 f g 2 Two matrices of rank 2 s.t. A(x) 0 f = g1 2 + g2 2 = g3 2 + g4 2 No matrices of rank 3 s.t. A(x) 0 f h1 2 + h2 2 + h3 2
30 12 / 12 Perspectives 1. Remove genericity assumptions on the input linear matrix A 2. Use of numerical homotopy for studying incidence varieties 3. Theoretical toolbox for analyzing singularities of determinantal varieties Surprising applications in optimal control techniques for the contrast imaging problem in medical imagery joint work with B. Bonnard, J.-C. Faugère, A. Jacquemard, T. Verron.
SPECTRA - a Maple library for solving linear matrix inequalities in exact arithmetic
SPECTRA - a Maple library for solving linear matrix inequalities in exact arithmetic Didier Henrion Simone Naldi Mohab Safey El Din Version 1.0 of November 5, 2016 Abstract This document briefly describes
More informationSolving rank-constrained semidefinite programs in exact arithmetic
Solving rank-constrained semidefinite programs in exact arithmetic Simone Naldi Technische Universität Dortmund Fakultät für Mathematik Vogelpothsweg 87, 44227 Dortmund Abstract We consider the problem
More informationSymbolic computation in hyperbolic programming
Symbolic computation in hyperbolic programming Simone Naldi and Daniel Plaumann MEGA 2017 Nice, 14 June 2017 1 Polyhedra P = {x R n : i, l i (x) 0} Finite # linear inequalities : l i (x) 0, i = 1,...,
More informationExact algorithms for linear matrix inequalities
Exact algorithms for linear matrix inequalities Didier Henrion Simone Naldi Mohab Safey El Din arxiv:1508.03715v3 [math.oc] 19 Sep 2016 September 20, 2016 Abstract Let A(x) = A 0 + x 1 A 1 + + x n A n
More informationExact algorithms for semidefinite programs with degenerate feasible set
Exact algorithms for semidefinite programs with degenerate feasible set Didier Henrion Simone Naldi Mohab Safey El Din 2018 Abstract Let A 0,..., A n be m m symmetric matrices with entries in Q, and let
More informationExact algorithms: from Semidefinite to Hyperbolic programming
Exact algorithms: from Semidefinite to Hyperbolic programming Simone Naldi and Daniel Plaumann PGMO Days 2017 Nov 14 2017 - EDF Lab Paris Saclay 1 Polyhedra P = {x R n : i, l i (x) 0} Finite # linear inequalities
More information6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC
6-1 The Positivstellensatz P. Parrilo and S. Lall, ECC 2003 2003.09.02.10 6. The Positivstellensatz Basic semialgebraic sets Semialgebraic sets Tarski-Seidenberg and quantifier elimination Feasibility
More informationComputing the Global Optimum of a Multivariate Polynomial over the Reals
Computing the Global Optimum of a Multivariate Polynomial over the Reals Mohab Safey El Din INRIA, Paris-Rocquencourt Center, SALSA Project UPMC, Univ Paris 06, LIP6 CNRS, UMR 7606, LIP6 UFR Ingéniérie
More informationLMI MODELLING 4. CONVEX LMI MODELLING. Didier HENRION. LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ. Universidad de Valladolid, SP March 2009
LMI MODELLING 4. CONVEX LMI MODELLING Didier HENRION LAAS-CNRS Toulouse, FR Czech Tech Univ Prague, CZ Universidad de Valladolid, SP March 2009 Minors A minor of a matrix F is the determinant of a submatrix
More informationGeneralized critical values and testing sign conditions
Generalized critical values and testing sign conditions on a polynomial Mohab Safey El Din. Abstract. Let f be a polynomial in Q[X 1,..., X n] of degree D. We focus on testing the emptiness of the semi-algebraic
More informationGlobal optimization, polynomial optimization, polynomial system solving, real
PROBABILISTIC ALGORITHM FOR POLYNOMIAL OPTIMIZATION OVER A REAL ALGEBRAIC SET AURÉLIEN GREUET AND MOHAB SAFEY EL DIN Abstract. Let f, f 1,..., f s be n-variate polynomials with rational coefficients of
More informationInformatique Mathématique une photographie en X et Y et Z et... (Eds.)
Informatique Mathématique une photographie en 2015 X et Y et Z et... (Eds.) 9 mars 2015 Table des matières 1 Algorithms in real algebraic geometry 3 1.1 Introduction............................ 4 1.1.1
More informationComputing Rational Points in Convex Semi-algebraic Sets and Sums-of-Squares Decompositions
Computing Rational Points in Convex Semi-algebraic Sets and Sums-of-Squares Decompositions Mohab Safey El Din 1 Lihong Zhi 2 1 University Pierre et Marie Curie, Paris 6, France INRIA Paris-Rocquencourt,
More informationOn the intrinsic complexity of point finding in. real singular hypersurfaces
On the intrinsic complexity of point finding in real singular hypersurfaces Bernd Bank 2, Marc Giusti 3, Joos Heintz 4, Luis Miguel Pardo 5 May 20, 2009 Abstract In previous work we designed an efficient
More informationAlgebraic classification methods for an optimal control problem in medical imagery
Algebraic classification methods for an optimal control problem in medical imagery Bernard Bonnard, Olivier Cots 5 Jean-Charles Faugère 3 Alain Jacquemard Jérémy Rouot 4 Mohab Safey El Din 3 Thibaut Verron
More informationStrong Bi-homogeneous Bézout s Theorem and Degree Bounds for Algebraic Optimization
Strong Bi-homogeneous Bézout s Theorem and Degree Bounds for Algebraic Optimization Mohab Safey El Din LIP6, Université Paris 6 75005 Paris, France MohabSafey@lip6fr Philippe Trébuchet LIP6, Université
More informationSolving the Birkhoff Interpolation Problem via the Critical Point Method: an Experimental Study
Solving the Birkhoff Interpolation Problem via the Critical Point Method: an Experimental Study F. Rouillier 1, M. Safey El Din 2, É. Schost3 1 LORIA, INRIA-Lorraine, Nancy, France Fabrice.Rouillier@loria.fr
More informationCOURSE ON LMI PART I.2 GEOMETRY OF LMI SETS. Didier HENRION henrion
COURSE ON LMI PART I.2 GEOMETRY OF LMI SETS Didier HENRION www.laas.fr/ henrion October 2006 Geometry of LMI sets Given symmetric matrices F i we want to characterize the shape in R n of the LMI set F
More informationALGORITHMS IN REAL ALGEBRAIC GEOMETRY: A SURVEY
ALGORITHMS IN REAL ALGEBRAIC GEOMETRY: A SURVEY SAUGATA BASU Abstract. We survey both old and new developments in the theory of algorithms in real algebraic geometry starting from effective quantifier
More informationALGORITHMS IN REAL ALGEBRAIC GEOMETRY: A SURVEY
ALGORITHMS IN REAL ALGEBRAIC GEOMETRY: A SURVEY SAUGATA BASU Abstract. We survey both old and new developments in the theory of algorithms in real algebraic geometry starting from effective quantifier
More informationCONVEX ALGEBRAIC GEOMETRY. Bernd Sturmfels UC Berkeley
CONVEX ALGEBRAIC GEOMETRY Bernd Sturmfels UC Berkeley Multifocal Ellipses Given m points (u 1, v 1 ),...,(u m, v m )intheplaner 2,and aradiusd > 0, their m-ellipse is the convex algebraic curve { (x, y)
More informationLecture Note 5: Semidefinite Programming for Stability Analysis
ECE7850: Hybrid Systems:Theory and Applications Lecture Note 5: Semidefinite Programming for Stability Analysis Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio State
More informationSPECTRAHEDRA. Bernd Sturmfels UC Berkeley
SPECTRAHEDRA Bernd Sturmfels UC Berkeley GAeL Lecture I on Convex Algebraic Geometry Coimbra, Portugal, Monday, June 7, 2010 Positive Semidefinite Matrices For a real symmetric n n-matrix A the following
More informationComputer algebra methods for testing the structural stability of multidimensional systems
Computer algebra methods for testing the structural stability of multidimensional systems Yacine Bouzidi INRIA Saclay L2S, Supelec, Gif-sur-Yvette, France Email: yacinebouzidi@inriafr Alban Quadrat INRIA
More informationA nearly optimal algorithm for deciding connectivity queries in smooth and bounded real algebraic sets
A nearly optimal algorithm for deciding connectivity queries in smooth and bounded real algebraic sets Mohab Safey El Din, Éric Schost To cite this version: Mohab Safey El Din, Éric Schost. A nearly optimal
More informationThe moment-lp and moment-sos approaches
The moment-lp and moment-sos approaches LAAS-CNRS and Institute of Mathematics, Toulouse, France CIRM, November 2013 Semidefinite Programming Why polynomial optimization? LP- and SDP- CERTIFICATES of POSITIVITY
More informationQuartic Spectrahedra
Quartic Spectrahedra arxiv:1311.3675v2 [math.oc] 7 Nov 2014 John Christian Ottem, Kristian Ranestad, Bernd Sturmfels and Cynthia Vinzant Abstract Quartic spectrahedra in 3-space form a semialgebraic set
More informationThe Geometry of Semidefinite Programming. Bernd Sturmfels UC Berkeley
The Geometry of Semidefinite Programming Bernd Sturmfels UC Berkeley Positive Semidefinite Matrices For a real symmetric n n-matrix A the following are equivalent: All n eigenvalues of A are positive real
More informationSolving Polynomial Equations in Smoothed Polynomial Time
Solving Polynomial Equations in Smoothed Polynomial Time Peter Bürgisser Technische Universität Berlin joint work with Felipe Cucker (CityU Hong Kong) Workshop on Polynomial Equation Solving Simons Institute
More informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) G. Sum of Squares (SOS) G.1 SOS Program: SOS/PSD and SDP G.2 Duality, valid ineqalities and Cone G.3 Feasibility/Optimization
More informationEE 227A: Convex Optimization and Applications October 14, 2008
EE 227A: Convex Optimization and Applications October 14, 2008 Lecture 13: SDP Duality Lecturer: Laurent El Ghaoui Reading assignment: Chapter 5 of BV. 13.1 Direct approach 13.1.1 Primal problem Consider
More informationRepresentations of Positive Polynomials: Theory, Practice, and
Representations of Positive Polynomials: Theory, Practice, and Applications Dept. of Mathematics and Computer Science Emory University, Atlanta, GA Currently: National Science Foundation Temple University
More informationBipolar varieties and real solving of a singular polynomial equation
Logo ISSN: 1889-3066 c 20xx Universidad de Jaén Web site: jja.ujaen.es Jaen J. Approx. x(x) (20xx), 1 xx Logo Bipolar varieties and real solving of a singular polynomial equation Bernd Bank, Marc Giusti,
More informationHYPERBOLIC POLYNOMIALS, INTERLACERS AND SUMS OF SQUARES
HYPERBOLIC POLYNOMIALS, INTERLACERS AND SUMS OF SQUARES DANIEL PLAUMANN Universität Konstanz joint work with Mario Kummer (U. Konstanz) Cynthia Vinzant (U. of Michigan) Out[121]= POLYNOMIAL OPTIMISATION
More informationDescribing convex semialgebraic sets by linear matrix inequalities. Markus Schweighofer. Universität Konstanz
Describing convex semialgebraic sets by linear matrix inequalities Markus Schweighofer Universität Konstanz Convex algebraic geometry, optimization and applications American Institute of Mathematics, Palo
More informationGram Spectrahedra. Lynn Chua UC Berkeley. joint work with Daniel Plaumann, Rainer Sinn, and Cynthia Vinzant
Gram Spectrahedra Lynn Chua UC Berkeley joint work with Daniel Plaumann, Rainer Sinn, and Cynthia Vinzant SIAM Conference on Applied Algebraic Geometry 017 1/17 Gram Matrices and Sums of Squares Let P
More informationEstimating the Region of Attraction of Ordinary Differential Equations by Quantified Constraint Solving
Estimating the Region of Attraction of Ordinary Differential Equations by Quantified Constraint Solving Henning Burchardt and Stefan Ratschan October 31, 2007 Abstract We formulate the problem of estimating
More informationUnbounded Convex Semialgebraic Sets as Spectrahedral Shadows
Unbounded Convex Semialgebraic Sets as Spectrahedral Shadows Shaowei Lin 9 Dec 2010 Abstract Recently, Helton and Nie [3] showed that a compact convex semialgebraic set S is a spectrahedral shadow if the
More informationQuartic Spectrahedra
Quartic Spectrahedra arxiv:1311.3675v1 [math.oc] 14 Nov 2013 John Christian Ottem, Kristian Ranestad, Bernd Sturmfels and Cynthia Vinzant Abstract Quartic spectrahedra in 3-space form a semialgebraic set
More informationd A 0 + m t k A k 0 whenever λ min (B k (x)) t k λ max (B k (x)) for k = 1, 2,..., m x n B n (k).
MATRIX CUBES PARAMETERIZED BY EIGENVALUES JIAWANG NIE AND BERND STURMFELS Abstract. An elimination problem in semidefinite programming is solved by means of tensor algebra. It concerns families of matrix
More informationConvergence rates of moment-sum-of-squares hierarchies for volume approximation of semialgebraic sets
Convergence rates of moment-sum-of-squares hierarchies for volume approximation of semialgebraic sets Milan Korda 1, Didier Henrion,3,4 Draft of December 1, 016 Abstract Moment-sum-of-squares hierarchies
More informationSemidefinite Programming
Semidefinite Programming Notes by Bernd Sturmfels for the lecture on June 26, 208, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra The transition from linear algebra to nonlinear algebra has
More informationIntrinsic complexity estimates in polynomial optimization 1
Intrinsic complexity estimates in polynomial optimization 1 Bernd Bank 2 Marc Giusti 3 Joos Heintz 4 Mohab Safey El Din 5 January 27, 2014 Abstract It is known that point searching in basic semialgebraic
More informationCOMPUTING RATIONAL POINTS IN CONVEX SEMI-ALGEBRAIC SETS AND SOS DECOMPOSITIONS
COMPUTING RATIONAL POINTS IN CONVEX SEMI-ALGEBRAIC SETS AND SOS DECOMPOSITIONS MOHAB SAFEY EL DIN AND LIHONG ZHI Abstract. Let P = {h 1,..., h s} Z[Y 1,..., Y k ], D deg(h i ) for 1 i s, σ bounding the
More informationReal solving polynomial equations with semidefinite programming
.5cm. Real solving polynomial equations with semidefinite programming Jean Bernard Lasserre - Monique Laurent - Philipp Rostalski LAAS, Toulouse - CWI, Amsterdam - ETH, Zürich LAW 2008 Real solving polynomial
More informationVariant Quantifier Elimination
Variant Quantifier Elimination Hoon Hong Department of Mathematics North Carolina State University Raleigh NC 27695, USA Mohab Safey El Din INRIA, Paris-Rocquencourt, SALSA Project UPMC, Univ Paris 06,
More informationSemi-definite representibility. For fun and profit
Semi-definite representation For fun and profit March 9, 2010 Introduction Any convex quadratic constraint x T C T Cx a + b T x Can be recast as the linear matrix inequality (LMI): ( ) In Cx (Cx) T a +
More informationELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications
ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications Professor M. Chiang Electrical Engineering Department, Princeton University March
More informationPOSITIVE SEMIDEFINITE INTERVALS FOR MATRIX PENCILS
POSITIVE SEMIDEFINITE INTERVALS FOR MATRIX PENCILS RICHARD J. CARON, HUIMING SONG, AND TIM TRAYNOR Abstract. Let A and E be real symmetric matrices. In this paper we are concerned with the determination
More informationCSCI 1951-G Optimization Methods in Finance Part 01: Linear Programming
CSCI 1951-G Optimization Methods in Finance Part 01: Linear Programming January 26, 2018 1 / 38 Liability/asset cash-flow matching problem Recall the formulation of the problem: max w c 1 + p 1 e 1 = 150
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationProving Unsatisfiability in Non-linear Arithmetic by Duality
Proving Unsatisfiability in Non-linear Arithmetic by Duality [work in progress] Daniel Larraz, Albert Oliveras, Enric Rodríguez-Carbonell and Albert Rubio Universitat Politècnica de Catalunya, Barcelona,
More informationChange of Ordering for Regular Chains in Positive Dimension
Change of Ordering for Regular Chains in Positive Dimension X. Dahan, X. Jin, M. Moreno Maza, É. Schost University of Western Ontario, London, Ontario, Canada. École polytechnique, 91128 Palaiseau, France.
More informationReal Interactive Proofs for VPSPACE
Brandenburgische Technische Universität, Cottbus-Senftenberg, Germany Colloquium Logicum Hamburg, September 2016 joint work with M. Baartse 1. Introduction Blum-Shub-Smale model of computability and complexity
More informationI.3. LMI DUALITY. Didier HENRION EECI Graduate School on Control Supélec - Spring 2010
I.3. LMI DUALITY Didier HENRION henrion@laas.fr EECI Graduate School on Control Supélec - Spring 2010 Primal and dual For primal problem p = inf x g 0 (x) s.t. g i (x) 0 define Lagrangian L(x, z) = g 0
More informationALGEBRAIC DEGREE OF POLYNOMIAL OPTIMIZATION. 1. Introduction. f 0 (x)
ALGEBRAIC DEGREE OF POLYNOMIAL OPTIMIZATION JIAWANG NIE AND KRISTIAN RANESTAD Abstract. Consider the polynomial optimization problem whose objective and constraints are all described by multivariate polynomials.
More informationSemidefinite Programming Basics and Applications
Semidefinite Programming Basics and Applications Ray Pörn, principal lecturer Åbo Akademi University Novia University of Applied Sciences Content What is semidefinite programming (SDP)? How to represent
More information4. Algebra and Duality
4-1 Algebra and Duality P. Parrilo and S. Lall, CDC 2003 2003.12.07.01 4. Algebra and Duality Example: non-convex polynomial optimization Weak duality and duality gap The dual is not intrinsic The cone
More informationConvex Optimization & Parsimony of L p-balls representation
Convex Optimization & Parsimony of L p -balls representation LAAS-CNRS and Institute of Mathematics, Toulouse, France IMA, January 2016 Motivation Unit balls associated with nonnegative homogeneous polynomials
More informationEffectivity Issues and Results for Hilbert 17 th Problem
Effectivity Issues and Results for Hilbert 17 th Problem Marie-Françoise Roy Université de Rennes 1, France joint work with Henri Lombardi Université de Franche-Comté, France and Daniel Perrucci Universidad
More informationSPECTRAHEDRA. Bernd Sturmfels UC Berkeley
SPECTRAHEDRA Bernd Sturmfels UC Berkeley Mathematics Colloquium, North Carolina State University February 5, 2010 Positive Semidefinite Matrices For a real symmetric n n-matrix A the following are equivalent:
More informationA baby steps/giant steps probabilistic algorithm for computing roadmaps in smooth bounded real hypersurface
A baby steps/giant steps probabilistic algorithm for computing roadmaps in smooth bounded real hypersurface Mohab Safey el Din Université Paris 6 and INRIA Paris-Rocquencourt Mohab.Safey@lip6.fr Éric Schost
More informationBounding the radii of balls meeting every connected component of semi-algebraic sets
Bounding the radii of balls meeting every connected component of semi-algebraic sets Saugata Basu Department of Mathematics Purdue University West Lafayette, IN 47907 USA Marie-Françoise Roy IRMAR (URA
More informationQuantifier Elimination
J.JSSAC (2003) Vol. 10, No. 1, pp. 3-12 Quantifier Elimination Quantifier Elimination ( ) IT 1 (quantifier elimination problem) (formal theory) ( ) (quantifier elimination :QE) (elementary theory of real
More informationGlobal Optimization of Polynomials
Semidefinite Programming Lecture 9 OR 637 Spring 2008 April 9, 2008 Scribe: Dennis Leventhal Global Optimization of Polynomials Recall we were considering the problem min z R n p(z) where p(z) is a degree
More informationSOLVING ALGEBRAIC INEQUALITIES
SOLVING ALGEBRAIC INEQUALITIES ADAM STRZEBOŃSKI ABSTRACT. We study the problem of solving, possibly quantified, systems of real algebraic equations and inequalities. We propose a way of representing solution
More informationVICTORIA POWERS AND THORSTEN W ORMANN Suppose there exists a real, symmetric, psd matrix B such that f = x B x T and rank t. Since B is real symmetric
AN ALGORITHM FOR SUMS OF SQUARES OF REAL POLYNOMIALS Victoria Powers and Thorsten Wormann Introduction We present an algorithm to determine if a real polynomial is a sum of squares (of polynomials), and
More informationOptimization on linear matrix inequalities for polynomial systems control
arxiv:1309.3112v1 [math.oc] 12 Sep 2013 Optimization on linear matrix inequalities for polynomial systems control Didier Henrion 1,2 Draft lecture notes of September 13, 2013 1 LAAS-CNRS, Univ. Toulouse,
More informationIntroduction to Semidefinite Programming I: Basic properties a
Introduction to Semidefinite Programming I: Basic properties and variations on the Goemans-Williamson approximation algorithm for max-cut MFO seminar on Semidefinite Programming May 30, 2010 Semidefinite
More informationLecture 7: Positive Semidefinite Matrices
Lecture 7: Positive Semidefinite Matrices Rajat Mittal IIT Kanpur The main aim of this lecture note is to prepare your background for semidefinite programming. We have already seen some linear algebra.
More informationDetecting global optimality and extracting solutions in GloptiPoly
Detecting global optimality and extracting solutions in GloptiPoly Didier HENRION 1,2 Jean-Bernard LASSERRE 1 1 LAAS-CNRS Toulouse 2 ÚTIA-AVČR Prague Part 1 Description of GloptiPoly Brief description
More informationConvex Optimization 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Convex Optimization 1 Many optimization objectives generated
More informationMATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.
MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If
More information15. Conic optimization
L. Vandenberghe EE236C (Spring 216) 15. Conic optimization conic linear program examples modeling duality 15-1 Generalized (conic) inequalities Conic inequality: a constraint x K where K is a convex cone
More informationAnalysis and synthesis: a complexity perspective
Analysis and synthesis: a complexity perspective Pablo A. Parrilo ETH ZürichZ control.ee.ethz.ch/~parrilo Outline System analysis/design Formal and informal methods SOS/SDP techniques and applications
More informationResearch overview. Seminar September 4, Lehigh University Department of Industrial & Systems Engineering. Research overview.
Research overview Lehigh University Department of Industrial & Systems Engineering COR@L Seminar September 4, 2008 1 Duality without regularity condition Duality in non-exact arithmetic 2 interior point
More informationConvex algebraic geometry, optimization and applications
Convex algebraic geometry, optimization and applications organized by William Helton and Jiawang Nie Workshop Summary We start with a little bit of terminology. A linear pencil is a function L : R g d
More informationSemidefinite representation of convex hulls of rational varieties
Semidefinite representation of convex hulls of rational varieties Didier Henrion 1,2 June 1, 2011 Abstract Using elementary duality properties of positive semidefinite moment matrices and polynomial sum-of-squares
More informationNonsmooth optimization: conditioning, convergence, and semi-algebraic models
Nonsmooth optimization: conditioning, convergence, and semi-algebraic models Adrian Lewis ORIE Cornell International Congress of Mathematicians Seoul, August 2014 1/16 Outline I Optimization and inverse
More informationDeterminant maximization with linear. S. Boyd, L. Vandenberghe, S.-P. Wu. Information Systems Laboratory. Stanford University
Determinant maximization with linear matrix inequality constraints S. Boyd, L. Vandenberghe, S.-P. Wu Information Systems Laboratory Stanford University SCCM Seminar 5 February 1996 1 MAXDET problem denition
More informationQUARTIC SPECTRAHEDRA. Bernd Sturmfels UC Berkeley and MPI Bonn. Joint work with John Christian Ottem, Kristian Ranestad and Cynthia Vinzant
QUARTIC SPECTRAHEDRA Bernd Sturmfels UC Berkeley and MPI Bonn Joint work with John Christian Ottem, Kristian Ranestad and Cynthia Vinzant 1 / 20 Definition A spectrahedron of degree n in R 3 is a convex
More informationIdentifying Redundant Linear Constraints in Systems of Linear Matrix. Inequality Constraints. Shafiu Jibrin
Identifying Redundant Linear Constraints in Systems of Linear Matrix Inequality Constraints Shafiu Jibrin (shafiu.jibrin@nau.edu) Department of Mathematics and Statistics Northern Arizona University, Flagstaff
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More information1 Strict local optimality in unconstrained optimization
ORF 53 Lecture 14 Spring 016, Princeton University Instructor: A.A. Ahmadi Scribe: G. Hall Thursday, April 14, 016 When in doubt on the accuracy of these notes, please cross check with the instructor s
More informationRésolution de systèmes polynomiaux structurés et applications en Cryptologie
Résolution de systèmes polynomiaux structurés et applications en Cryptologie Pierre-Jean Spaenlehauer University of Western Ontario Ontario Research Center for Computer Algebra Magali Bardet, Jean-Charles
More informationOverview of Computer Algebra
Overview of Computer Algebra http://cocoa.dima.unige.it/ J. Abbott Universität Kassel J. Abbott Computer Algebra Basics Manchester, July 2018 1 / 12 Intro Characteristics of Computer Algebra or Symbolic
More information4. Convex optimization problems
Convex Optimization Boyd & Vandenberghe 4. Convex optimization problems optimization problem in standard form convex optimization problems quasiconvex optimization linear optimization quadratic optimization
More informationarxiv: v2 [math.ag] 24 Jun 2015
TRIANGULATIONS OF MONOTONE FAMILIES I: TWO-DIMENSIONAL FAMILIES arxiv:1402.0460v2 [math.ag] 24 Jun 2015 SAUGATA BASU, ANDREI GABRIELOV, AND NICOLAI VOROBJOV Abstract. Let K R n be a compact definable set
More informationRank-one Generated Spectral Cones Defined by Two Homogeneous Linear Matrix Inequalities
Rank-one Generated Spectral Cones Defined by Two Homogeneous Linear Matrix Inequalities C.J. Argue Joint work with Fatma Kılınç-Karzan October 22, 2017 INFORMS Annual Meeting Houston, Texas C. Argue, F.
More informationPolytopes and Algebraic Geometry. Jesús A. De Loera University of California, Davis
Polytopes and Algebraic Geometry Jesús A. De Loera University of California, Davis Outline of the talk 1. Four classic results relating polytopes and algebraic geometry: (A) Toric Geometry (B) Viro s Theorem
More informationLecture: Convex Optimization Problems
1/36 Lecture: Convex Optimization Problems http://bicmr.pku.edu.cn/~wenzw/opt-2015-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghe s lecture notes Introduction 2/36 optimization
More informationHandout 8: Dealing with Data Uncertainty
MFE 5100: Optimization 2015 16 First Term Handout 8: Dealing with Data Uncertainty Instructor: Anthony Man Cho So December 1, 2015 1 Introduction Conic linear programming CLP, and in particular, semidefinite
More informationALGEBRA: From Linear to Non-Linear. Bernd Sturmfels University of California at Berkeley
ALGEBRA: From Linear to Non-Linear Bernd Sturmfels University of California at Berkeley John von Neumann Lecture, SIAM Annual Meeting, Pittsburgh, July 13, 2010 Undergraduate Linear Algebra All undergraduate
More informationSEMIDEFINITE PROGRAM BASICS. Contents
SEMIDEFINITE PROGRAM BASICS BRIAN AXELROD Abstract. A introduction to the basics of Semidefinite programs. Contents 1. Definitions and Preliminaries 1 1.1. Linear Algebra 1 1.2. Convex Analysis (on R n
More informationThe Algebraic Degree of Semidefinite Programming
The Algebraic Degree ofsemidefinite Programming p. The Algebraic Degree of Semidefinite Programming BERND STURMFELS UNIVERSITY OF CALIFORNIA, BERKELEY joint work with and Jiawang Nie (IMA Minneapolis)
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationOptimal Rate and Maximum Erasure Probability LDPC Codes in Binary Erasure Channel
Optimal Rate and Maximum Erasure Probability LDPC Codes in Binary Erasure Channel H. Tavakoli Electrical Engineering Department K.N. Toosi University of Technology, Tehran, Iran tavakoli@ee.kntu.ac.ir
More informationLecture 14: Optimality Conditions for Conic Problems
EE 227A: Conve Optimization and Applications March 6, 2012 Lecture 14: Optimality Conditions for Conic Problems Lecturer: Laurent El Ghaoui Reading assignment: 5.5 of BV. 14.1 Optimality for Conic Problems
More informationConvex optimization. Javier Peña Carnegie Mellon University. Universidad de los Andes Bogotá, Colombia September 2014
Convex optimization Javier Peña Carnegie Mellon University Universidad de los Andes Bogotá, Colombia September 2014 1 / 41 Convex optimization Problem of the form where Q R n convex set: min x f(x) x Q,
More informationRank minimization via the γ 2 norm
Rank minimization via the γ 2 norm Troy Lee Columbia University Adi Shraibman Weizmann Institute Rank Minimization Problem Consider the following problem min X rank(x) A i, X b i for i = 1,..., k Arises
More information