Testing Sign Conditions on a Multivariate Polynomial and Applications

Let $f$ be a polynomial in $Q[X_{1},…,X_{n}]$ of degree $D$. We focus on testing the emptiness and computing at least one point in each connected component of the semi-algebraic set defined by $f > 0$ (or $f < 0$ or $f \neq 0$). To this end, the problem is reduced to computing at least one point in each connected component of a hypersurface defined by $f − e = 0$ for $e \in Q$ positive and small enough. We provide an algorithm allowing us to determine a positive rational number e which is small enough in this sense. This is based on the efficient computation of the set of generalized critical values of the mapping $f:y \in C^{n} \to f(y) \in C$ which is the union of the classical set of critical values of the mapping $f$ and the set of asymptotic critical values of the mapping $f$. Then, we show how to use the computation of generalized critical values in order to obtain an efficient algorithm deciding the emptiness of a semi-algebraic set defined by a single inequality or a single inequation. At last, we show how to apply our contribution to determining if a hypersurface contains real regular points. We provide complexity estimates for probabilistic versions of the latter algorithms which are within $O(n^{7}D^{4n})$ arithmetic operations in $Q$. The paper ends with practical experiments showing the efficiency of our approach on real-life applications.

[1]  J. Renegar,et al.  On the Computational Complexity and Geometry of the First-Order Theory of the Reals, Part I , 1989 .

[2]  B. Bank,et al.  Polar Varieties and Efficient Real Equation Solving: The Hypersurface Case , 1996 .

[3]  P. Zimmermann,et al.  Efficient isolation of polynomial's real roots , 2004 .

[4]  Fabrice Rouillier,et al.  Solving the Birkhoff Interpolation Problem via the Critical Point Method: An Experimental Study , 2000, Automated Deduction in Geometry.

[5]  Markus Schweighofer Global Optimization of Polynomials Using Gradient Tentacles and Sums of Squares , 2006, SIAM J. Optim..

[6]  Mohab Safey El Din,et al.  Properness defects of projections and computation of one point in each connected component of a real algebraic set , 2003 .

[7]  Grégoire Lecerf,et al.  Computing the equidimensional decomposition of an algebraic closed set by means of lifting fibers , 2003, J. Complex..

[8]  Marc Giusti,et al.  Generalized polar varieties and an efficient real elimination , 2004, Kybernetika.

[9]  J. E. Morais,et al.  When Polynomial Equation Systems Can Be "Solved" Fast? , 1995, AAECC.

[10]  K. Kurdyka,et al.  SEMIALGEBRAIC SARD THEOREM FOR GENERALIZED CRITICAL VALUES , 2000 .

[11]  J. E. Morais,et al.  Straight--Line Programs in Geometric Elimination Theory , 1996, alg-geom/9609005.

[12]  Jean Charles Faugère,et al.  A new efficient algorithm for computing Gröbner bases without reduction to zero (F5) , 2002, ISSAC '02.

[13]  C. Wampler,et al.  Basic Algebraic Geometry , 2005 .

[14]  Hazel Everett,et al.  The Voronoi diagram of three lines in R3 , 2006 .

[15]  Johan P. Hansen,et al.  INTERSECTION THEORY , 2011 .

[16]  Jean-Benoît Bost,et al.  Heights of projective varieties and positive Green forms , 1994 .

[17]  Zbigniew Jelonek,et al.  Testing sets for properness of polynomial mappings , 1999 .

[18]  J. Risler,et al.  Real algebraic and semi-algebraic sets , 1990 .

[19]  Y. N. Lakshman A Single Exponential Bound on the Complexity of Computing Gröbner Bases of Zero Dimensional Ideals , 1991 .

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

[21]  D. S. Arnon,et al.  Algorithms in real algebraic geometry , 1988 .

[22]  Daniel Lazard,et al.  Quantifier Elimination: Optimal Solution for Two Classical Examples , 1988, J. Symb. Comput..

[23]  Anneli Lax,et al.  On sums of squares , 1978 .

[24]  Fabrice Rouillier,et al.  Solving Zero-Dimensional Systems Through the Rational Univariate Representation , 1999, Applicable Algebra in Engineering, Communication and Computing.

[25]  Marc Giusti,et al.  Lower bounds for diophantine approximations , 1997 .

[26]  Zbigniew Jelonek,et al.  On asymptotic critical values of a complex polynomial , 2003 .

[27]  Joos Heintz,et al.  On the Theoretical and Practical Complexity of the Existential Theory of Reals , 1993, Comput. J..

[28]  Palaiseau Cedex,et al.  Computing Parametric Geometric Resolutions , 2001 .

[29]  Fabrice Rouillier,et al.  Real Solving for Positive Dimensional Systems , 2002, J. Symb. Comput..

[30]  Hazel Everett,et al.  The Voronoi Diagram of Three Lines , 2007, SCG '07.

[31]  Donal O'Shea,et al.  Ideals, varieties, and algorithms - an introduction to computational algebraic geometry and commutative algebra (2. ed.) , 1997, Undergraduate texts in mathematics.

[32]  J. Faugère A new efficient algorithm for computing Gröbner bases (F4) , 1999 .

[33]  S. Basu,et al.  A New Algorithm to Find a Point in Every Cell Defined by a Family of Polynomials , 1998 .

[34]  J. E. Morais,et al.  Lower Bounds for diophantine Approximation , 1996 .

[35]  Dima Grigoriev,et al.  Solving Systems of Polynomial Inequalities in Subexponential Time , 1988, J. Symb. Comput..

[36]  I. Shafarevich Basic algebraic geometry , 1974 .

[37]  Z. Jelonek Topological characterization of finite mappings , 2001 .

[38]  Éric Schost,et al.  Properness Defects of Projections and Computation of at Least One Point in Each Connected Component of a Real Algebraic Set , 2004, Discret. Comput. Geom..

[39]  Marc Giusti,et al.  Generalized polar varieties: geometry and algorithms , 2005, J. Complex..

[40]  Mohab Safey El Din Finding Sampling Points on Real Hypersurfaces in Easier in Singular Situations , 2005 .

[41]  Marc Giusti,et al.  A Gröbner Free Alternative for Polynomial System Solving , 2001, J. Complex..

[42]  Solen Corvez,et al.  Using Computer Algebra Tools to Classify Serial Manipulators , 2002, Automated Deduction in Geometry.

[43]  James Renegar,et al.  On the Computational Complexity and Geometry of the First-Order Theory of the Reals, Part III: Quantifier Elimination , 1992, J. Symb. Comput..

[44]  Marie-Françoise Roy,et al.  On the combinatorial and algebraic complexity of Quanti erEliminationS , 1994 .

[45]  Éric Schost,et al.  Polar varieties and computation of one point in each connected component of a smooth real algebraic set , 2003, ISSAC '03.

[46]  James Renegar On the computational complexity and geome-try of the first-order theory of the reals , 1992 .

[47]  Fabrice Rouillier,et al.  Finding at Least One Point in Each Connected Component of a Real Algebraic Set Defined by a Single Equation , 2000, J. Complex..