Computing points in connected components defined by a real inequation: algorithms, complexity and implementations, Part I

TL;DR

This paper introduces a probabilistic algorithm for computing sample points in each connected component of semi-algebraic sets defined by polynomial inequalities, with complexity estimates and practical benchmarks.

Contribution

It presents a new probabilistic approach based on critical points and reductions to zero-dimensional systems, with refined complexity analysis and real-world experiments.

Findings

Algorithm has cubic complexity in the B{é}zout bound for exact points.

Cost is quartic in the B{é}zout bound for rational approximations.

Practical experiments outperform existing methods on complex polynomials.

Abstract

We consider the problem of computing sample points in each connected component of a semi-algebraic set defined by the non-vanishing or the positivity of an n-variate polynomial of degree d, with rational coefficients of bit size bounded by $\tau$. Such a problem is a basic routine in effective real algebraic geometry, used in higher-level algorithms for solving polynomial systems over the reals and finds many applications in sciences. We design a probabilistic algorithm for solving this problem, which is based on reductions to different routines for solving zero-dimensional polynomial systems. It assumes that the input polynomial satisfies sufficiently generic properties (namely, smoothness of its defining hypersurface). This is done through the computations of critical points of well-chosen maps to capture the connected components of the semi-algebraic set under study. We derive a bit complexity estimate for the cost of this algorithm, which is, in terms of the B{\'e}zout bound d(d -1)^{n-1}, essentially cubic for obtaining parametrisations of the sought-for real points. Moreover, we also consider the case of obtaining rational approximations of those points, which are precise enough to lie in the same connected components as their exact counterparts, which yields a cost that is essentially quartic in the B{\'e}zout bound. In these complexity estimates, we take into account the degree structure of the input polynomial and its partial derivatives, allowing for a more refined bit complexity when the partial derivative of the input polynomial have degree lower than expected. We also analyse the probability of success of those algorithms. We report on practical experiments, benchmarking with random dense input polynomials as well as polynomials coming from applications, which were out of reach of the state-of-the-art implementations, and hence illustrate the practical efficiency of these new algorithms.