EN FR
Homepage Inria website


Section: New Results

Solving systems over the reals and applications

It is well-known that every non-negative univariate real polynomial can be written as the sum of two polynomial squares with real coefficients. When one allows a weighted sum of finitely many squares instead of a sum of two squares, then one can choose all coefficients in the representation to lie in the field generated by the coefficients of the polynomial. In particular, this allows an effective treatment of polynomials with rational coefficients. In [11], we describe, analyze and compare both from the theoretical and practical points of view, two algorithms computing such a weighted sums of squares decomposition for univariate polynomials with rational coefficients. The first algorithm, due to the third author relies on real root isolation, quadratic approximations of positive polynomials and square-free decomposition but its complexity was not analyzed. We provide bit complexity estimates, both on the runtime and the output size of this algorithm. They are exponential in the degree of the input univariate polynomial and linear in the maximum bitsize of its complexity. This analysis is obtained using quantifier elimination and root isolation bounds. The second algorithm, due to Chevillard, Harrison, Joldes and Lauter, relies on complex root isolation and square-free decomposition and has been introduced for certifying positiveness of poly-nomials in the context of computer arithmetics. Again, its complexity was not analyzed. We provide bit complexity estimates, both on the runtime and the output size of this algorithm, which are polynomial in the degree of the input polynomial and linear in the maximum bitsize of its complexity. This analysis is obtained using Vieta's formula and root isolation bounds. Finally, we report on our implementations of both algorithms and compare them in practice on several application benchmarks. While the second algorithm is, as expected from the complexity result, more efficient on most of examples, we exhibit families of non-negative polynomials for which the first algorithm is better.

[9] describes our freely distributed Maple library spectra , for Semidefinite Programming solved Exactly with Computational Tools of Real Algebra. It solves linear matrix inequalities with symbolic computation in exact arithmetic and it is targeted to small-size, possibly degenerate problems for which symbolic infeasibility or feasibility certificates are required.

Let Sn be a compact basic semi-algebraic set defined as the real solution set of multivariate polynomial inequalities with rational coefficients. In [19], we design an algorithm which takes as input a polynomial system defining S and an integer p0 and returns the n-dimensional volume of S at absolute precision 2-p. Our algorithm relies on the relationship between volumes of semi-algebraic sets and periods of rational integrals. It makes use of algorithms computing the Picard-Fuchs differential equation of appropriate periods, properties of critical points, and high-precision numerical integration of differential equations. The algorithm runs in essentially linear time with respect to p. This improves upon the previous exponential bounds obtained by Monte-Carlo or moment-based methods. Assuming a conjecture of Dimca, the arithmetic cost of the algebraic subroutines for computing Picard-Fuchs equations and critical points is singly exponential in n and polynomial in the maximum degree of the input.

Let 𝐟=(f1,...,fs) be a sequence of polynomials in [X1,...,Xn] of maximal degree D and Vn be the algebraic set defined by 𝐟 and r be its dimension. The real radical 𝐟 associated to 𝐟 is the largest ideal which defines the real trace of V. When V is smooth, we show in [13], that 𝐟re, has a finite set of generators with degrees bounded by degV. Moreover, we present a probabilistic algorithm of complexity (snDn)O(1) to compute the minimal primes of 𝐟re. When V is not smooth, we give a probabilistic algorithm of complexity sO(1)(nD)O(nr2r) to compute rational parametrizations for all irreducible components of the real algebraic set Vn.

Let (g1,...,gp) in [X1,...,Xn] and S be the basic closed semi-algebraic set defined by g10,...,gp0. The S-radical of 𝐟, which is denoted by 𝐟S, is the ideal associated to the Zariski closure of VS. We give a probabilistic algorithm to compute rational parametrizations of all irreducible components of that Zariski closure, hence encoding 𝐟S. Assuming now that D is the maximum of the degrees of the fi's and the gi's, this algorithm runs in time 2p(s+p)O(1)(nD)O(rn2r).

Experiments are performed to illustrate and show the efficiency of our approaches on computing real radicals.

In [14], we consider the second-order discontinuous differential equation y''+η sgn (y)=θy+αsin(βt) where the parameters η,θ,α,β are real. The main goal is to discuss the existence of periodic solutions. Under explicit conditions, the number of such solutions is given. Furthermore, for each of these periodic solutions, an explicit formula is provided.