Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

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

Abstract : 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.
Complete list of metadata
Contributor : Huu Phuoc Le <>
Submitted on : Saturday, November 28, 2020 - 2:21:16 PM
Last modification on : Tuesday, March 23, 2021 - 9:28:03 AM


Files produced by the author(s)


  • HAL Id : hal-03029441, version 1


Huu Phuoc Le, Mohab Safey El Din. Solving parametric systems of polynomial equations over the reals through Hermite matrices. 2020. ⟨hal-03029441⟩



Record views


Files downloads