Solving parametric systems of polynomial equations over the reals through Hermite matrices - Archive ouverte HAL Accéder directement au contenu
Article Dans Une Revue Journal of Symbolic Computation Année : 2022

Solving parametric systems of polynomial equations over the reals through Hermite matrices

Huu Phuoc Le
Mohab Safey El Din

Résumé

We design a new algorithm for solving parametric systems of equations having finitely many complex solutions for generic values of the parameters. More precisely, let $\mathbf{f} = (f_1, \ldots, f_m)\subset \mathbb{Q}[\mathbf{y}][\mathbf{x}]$ with $\mathbf{y} = (y_1, \ldots, y_{t})$ and $\mathbf{x} = (x_1, \ldots, x_{n})$, $\mathcal{V}\subset \mathbb{C}^t \times \mathbb{C}^n$ be the algebraic set defined by the simultaneous vanishing of the $f_i$'s and $\pi$ be the projection $(\mathbf{y}, \mathbf{x}) \to \mathbf{y}$. Under the assumptions that $\mathbf{f}$ admits finitely many complex solutions when specializing $\mathbf{y}$ to generic values and that the ideal generated by $\mathbf{f}$ is radical, we solve the following algorithmic problem. On input $\mathbf{f}$, we compute {\em semi-algebraic formulas} defining open semi-algebraic sets $\mathcal{S}_1, \ldots, \mathcal{S}_{\ell}$ in the parameters' space $\mathbb{R}^t$ such that $\cup_{i=1}^{\ell} \mathcal{S}_i$ is dense in $\mathbb{R}^t$ and, for $1\leq i \leq \ell$, the number of real points in $\mathcal{V}\cap \pi^{-1}(\eta)$ is invariant when $\eta$ ranges over $\mathcal{S}_i$. This algorithm exploits special properties of some well chosen monomial bases in the quotient algebra $\mathbb{Q}(\mathbf{y})[\mathbf{x}] / I$ where $I\subset \mathbb{Q}(\mathbf{y})[\mathbf{x}]$ is the ideal generated by $\mathbf{f}$ in $\mathbb{Q}(\mathbf{y})[\mathbf{x}]$ as well as the specialization property of the so-called Hermite matrices which represent Hermite's quadratic forms. This allows us to obtain ``compact'' representations of the semi-algebraic sets $\mathcal{S}_i$ by means of semi-algebraic formulas encoding the signature of a given symmetric matrix. When $\mathbf{f}$ satisfies extra genericity assumptions (such as regularity), we use the theory of Gr\"obner bases to derive complexity bounds both on the number of arithmetic operations in $\mathbb{Q}$ and the degree of the output polynomials. More precisely, letting $d$ be the maximal degrees of the $f_i$'s and $\mathfrak{D} = n(d-1)d^n$, we prove that, on a generic input $\mathbf{f}=(f_1,\ldots,f_n)$, one can compute those semi-algebraic formulas with $O\ {\widetilde{~}}\left (\binom{t+\mathfrak{D}}{t}\ 2^{3t}\ n^{2t+1} d^{3nt+2(n+t)+1} \right )$ arithmetic operations in $\mathbb{Q}$ and that the polynomials involved in these formulas have degree bounded by $\mathfrak{D}$. We report on practical experiments which illustrate the efficiency of this algorithm, both on generic parametric systems and parametric systems coming from applications since it allows us to solve systems which are out of reach on the current state-of-the-art.
Fichier principal
Vignette du fichier
mainRRC.pdf (573.56 Ko) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)

Dates et versions

hal-03029441 , version 1 (28-11-2020)
hal-03029441 , version 2 (16-12-2021)

Identifiants

Citer

Huu Phuoc Le, Mohab Safey El Din. Solving parametric systems of polynomial equations over the reals through Hermite matrices. Journal of Symbolic Computation, In press, 112, pp.25-61. ⟨10.1016/j.jsc.2021.12.002⟩. ⟨hal-03029441v2⟩
268 Consultations
341 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More