Finding points on real solution components and applications to differential polynomial systems

In this paper we extend complex homotopy methods to finding witness points on the irreducible components of real varieties. In particular we construct such witness points as the isolated real solutions of a constrained optimization problem. First a random hyperplane characterized by its random normal vector is chosen. Witness points are computed by a polyhedral homotopy method. Some of them are at the intersection of this hyperplane with the components. Other witness points are the local critical points of the distance from the plane to components. A method is also given for constructing regular witness points on components, when the critical points are singular. The method is applicable to systems satisfying certain regularity conditions. Illustrative examples are given. We show that the method can be used in the consistent initialization phase of a popular method due to Pryce and Pantelides for preprocessing differential algebraic equations for numerical solution.

[1]  Jonathan D. Hauenstein,et al.  Cell decomposition of almost smooth real algebraic surfaces , 2013, Numerical Algorithms.

[2]  Hoon Hong,et al.  Improvements in cad-based quantifier elimination , 1990 .

[3]  Lihong Zhi,et al.  Determining singular solutions of polynomial systems via symbolic-numeric reduction to geometric involutive forms , 2012, J. Symb. Comput..

[4]  A. Seidenberg.  A NEW DECISION METHOD FOR ELEMENTARY ALGEBRA , 1954 .

[5]  Monique Laurent,et al.  A prolongation-projection algorithm for computing the finite real variety of an ideal , 2008, Theor. Comput. Sci..

[6]  Charles W. Wampler,et al.  Finding All Real Points of a Complex Curve , 2006 .

[7]  James H. Davenport,et al.  Real Quantifier Elimination is Doubly Exponential , 1988, J. Symb. Comput..

[8]  George E. Collins,et al.  Hauptvortrag: Quantifier elimination for real closed fields by cylindrical algebraic decomposition , 1975, Automata Theory and Formal Languages.

[9]  Changbo Chen,et al.  Computing cylindrical algebraic decomposition via triangular decomposition , 2009, ISSAC '09.

[10]  Eberhard Freitag,et al.  Analytic Functions of Several Complex Variables , 2011 .

[11]  J. Hauenstein.  Numerically Computing Real Points on Algebraic Sets , 2011, Acta Applicandae Mathematicae.

[12]  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..

[13]  Bruno Buchberger,et al.  Bruno Buchberger's PhD thesis 1965: An algorithm for finding the basis elements of the residue class ring of a zero dimensional polynomial ideal , 2006, J. Symb. Comput..

[14]  Andrew J. Sommese,et al.  The numerical solution of systems of polynomials - arising in engineering and science , 2005 .

[15]  A. Tarski.  A Decision Method for Elementary Algebra and Geometry , 2023 .

[16]  Changbo Chen,et al.  Triangular decomposition of semi-algebraic systems , 2013, J. Symb. Comput..

[17]  Anton Leykin,et al.  Higher-Order Deflation for Polynomial Systems With Isolated Singular Solutions , 2006, math/0602031.

[18]  A. Sard,et al.  The measure of the critical values of differentiable maps , 1942 .

[19]  Nicolai Vorobjov,et al.  Counting connected components of a semialgebraic set in subexponential time , 1992, computational complexity.

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

[21]  一松 信,et al.  R.C. Gunning and H.Rossi: Analytic Functions of Several Complex Variables, Prentice-Hall, Englewood Cliffs, N.J., 1965, 317頁, 15×23cm, $12.50. , 1965 .

[22]  Lihong Zhi,et al.  Computing the multiplicity structure from geometric involutive form , 2008, ISSAC '08.

[23]  Robert M. Corless,et al.  Numerical solutions of index-1 differential algebraic equations can be computed in polynomial time , 2006, Numerical Algorithms.

[24]  Ye Lu.  Finding all real solutions of polynomial systems , 2006 .

[25]  Wenyuan Wu,et al.  Implicit Riquier Bases for PDAE and their semi-discretizations , 2009, J. Symb. Comput..

[26]  C. Pantelides.  The consistent intialization of differential-algebraic systems , 1988 .