Computing rational solutions of linear matrix inequalities

  • Authors:
  • Qingdong Guo;Mohab Safey EI Din;Lihong Zhi

  • Affiliations:
  • Key Laboratory of Mathematics Mechanization, Chinese Academy of Sciences, Beijing, China;UPMC, Univ Paris 06/INRIA, Paris-Rocquencourt Center, POLSYS project-team/CNRS LIP6 UMR 7606/Institut Universitaire de France, Paris, France;Key Laboratory of Mathematics Mechanization, Chinese Academy of Sciences, Beijing, China

  • Venue:
  • Proceedings of the 38th international symposium on International symposium on symbolic and algebraic computation
  • Year:
  • 2013

Quantified Score

Hi-index 0.00

Visualization

Abstract

Consider a (D x D) symmetric matrix A whose entries are linear forms in Q[X1, ..., Xk] with coefficients of bit size ≤ τ. We provide an algorithm which decides the existence of rational solutions to the linear matrix inequality A ≥ 0 and outputs such a rational solution if it exists. This problem is of first importance: it can be used to compute algebraic certificates of positivity for multivariate polynomials. Our algorithm runs within (k≤)O(1)2O(min(k, D)D2)DO(D2) bit operations; the bit size of the output solution is dominated by τO(1)2O(\min(k, D)D2). These results are obtained by designing algorithmic variants of constructions introduced by Klep and Schweighofer. This leads to the best complexity bounds for deciding the existence of sums of squares with rational coefficients of a given polynomial. We have implemented the algorithm; it has been able to tackle Scheiderer's example of a multivariate polynomial that is a sum of squares over the reals but not over the rationals; providing the first computer validation of this counter-example to Sturmfels' conjecture.