55 research outputs found
Computing the homology of basic semialgebraic sets in weak exponential time
We describe and analyze an algorithm for computing the homology (Betti
numbers and torsion coefficients) of basic semialgebraic sets which works in
weak exponential time. That is, out of a set of exponentially small measure in
the space of data the cost of the algorithm is exponential in the size of the
data. All algorithms previously proposed for this problem have a complexity
which is doubly exponential (and this is so for almost all data)
Complexity of Computing the Local Dimension of a Semialgebraic Set
AbstractThe paper describes several algorithms related to a problem of computing the local dimension of a semialgebraic set. Let a semialgebraic set V be defined by a system of k inequalities of the formf≥ 0 with f∈R [ X1,⋯ ,Xn ], deg(f) <d , andx∈V . An algorithm is constructed for computing the dimension of the Zariski tangent space to V at x in time (kd)O(n). Let x belong to a stratum of codimension lxin V with respect to a smooth stratification ofV . Another algorithm computes the local dimension dimx(V) with the complexity (k(lx+ 1)d)O(lx2n). Ifl=maxx∈Vlx, and for every connected component the local dimension is the same at each point, then the algorithm computes the dimension of every connected component with complexity (k(l+ 1)d)O(l2n). If V is a real algebraic variety defined by a system of equations, then the complexity of the algorithm is less thankdO(l2n) , and the algorithm also finds the dimension of the tangent space to V at x in time kdO(n). Whenl is fixed, like in the case of a smooth V , the complexity bounds for computing the local dimension are (kd)O(n)andkdO(n) respectively. A third algorithm finds the singular locus ofV in time (kd)O(n2)
Polar Varieties and Efficient Real Elimination
Let be a smooth and compact real variety given by a reduced regular
sequence of polynomials . This paper is devoted to the
algorithmic problem of finding {\em efficiently} a representative point for
each connected component of . For this purpose we exhibit explicit
polynomial equations that describe the generic polar varieties of . This
leads to a procedure which solves our algorithmic problem in time that is
polynomial in the (extrinsic) description length of the input equations and in a suitably introduced, intrinsic geometric parameter, called
the {\em degree} of the real interpretation of the given equation system .Comment: 32 page
CAD Adjacency Computation Using Validated Numerics
We present an algorithm for computation of cell adjacencies for well-based
cylindrical algebraic decomposition. Cell adjacency information can be used to
compute topological operations e.g. closure, boundary, connected components,
and topological properties e.g. homology groups. Other applications include
visualization and path planning. Our algorithm determines cell adjacency
information using validated numerical methods similar to those used in CAD
construction, thus computing CAD with adjacency information in time comparable
to that of computing CAD without adjacency information. We report on
implementation of the algorithm and present empirical data.Comment: 20 page
Cylindrical Algebraic Decomposition Using Local Projections
We present an algorithm which computes a cylindrical algebraic decomposition
of a semialgebraic set using projection sets computed for each cell separately.
Such local projection sets can be significantly smaller than the global
projection set used by the Cylindrical Algebraic Decomposition (CAD) algorithm.
This leads to reduction in the number of cells the algorithm needs to
construct. We give an empirical comparison of our algorithm and the classical
CAD algorithm
Polar Varieties, Real Equation Solving and Data-Structures: The hypersurface case
In this paper we apply for the first time a new method for multivariate
equation solving which was developed in \cite{gh1}, \cite{gh2}, \cite{gh3} for
complex root determination to the {\em real} case. Our main result concerns the
problem of finding at least one representative point for each connected
component of a real compact and smooth hypersurface. The basic algorithm of
\cite{gh1}, \cite{gh2}, \cite{gh3} yields a new method for symbolically solving
zero-dimensional polynomial equation systems over the complex numbers. One
feature of central importance of this algorithm is the use of a
problem--adapted data type represented by the data structures arithmetic
network and straight-line program (arithmetic circuit). The algorithm finds the
complex solutions of any affine zero-dimensional equation system in non-uniform
sequential time that is {\em polynomial} in the length of the input (given in
straight--line program representation) and an adequately defined {\em geometric
degree of the equation system}. Replacing the notion of geometric degree of the
given polynomial equation system by a suitably defined {\em real (or complex)
degree} of certain polar varieties associated to the input equation of the real
hypersurface under consideration, we are able to find for each connected
component of the hypersurface a representative point (this point will be given
in a suitable encoding). The input equation is supposed to be given by a
straight-line program and the (sequential time) complexity of the algorithm is
polynomial in the input length and the degree of the polar varieties mentioned
above.Comment: Late
Constraint Satisfaction Problems over Numeric Domains
We present a survey of complexity results for constraint satisfaction problems (CSPs) over the integers, the rationals, the reals, and the complex numbers. Examples of such problems are feasibility of linear programs, integer linear programming, the max-atoms problem, Hilbert\u27s tenth problem, and many more. Our particular focus is to identify those CSPs that can be solved in polynomial time, and to distinguish them from CSPs that are NP-hard. A very helpful tool for obtaining complexity classifications in this context is the concept of a polymorphism from universal algebra
- …