Results 1  10
of
238
Recursive Markov chains, stochastic grammars, and monotone systems of nonlinear equations
 IN STACS
, 2005
"... We define Recursive Markov Chains (RMCs), a class of finitely presented denumerable Markov chains, and we study algorithms for their analysis. Informally, an RMC consists of a collection of finitestate Markov chains with the ability to invoke each other in a potentially recursive manner. RMCs offer ..."
Abstract

Cited by 67 (11 self)
 Add to MetaCart
We define Recursive Markov Chains (RMCs), a class of finitely presented denumerable Markov chains, and we study algorithms for their analysis. Informally, an RMC consists of a collection of finitestate Markov chains with the ability to invoke each other in a potentially recursive manner. RMCs offer a natural abstract model for probabilistic programs with procedures. They generalize, in a precise sense, a number of well studied stochastic models, including Stochastic ContextFree Grammars (SCFG) and MultiType Branching Processes (MTBP). We focus on algorithms for reachability and termination analysis for RMCs: what is the probability that an RMC started from a given state reaches another target state, or that it terminates? These probabilities are in general irrational, and they arise as (least) fixed point solutions to certain (monotone) systems of nonlinear equations associated with RMCs. We address both the qualitative problem of determining whether the probabilities are 0, 1 or inbetween, and
Semidefinite Representations for Finite Varieties
 MATHEMATICAL PROGRAMMING
, 2002
"... We consider the problem of minimizing a polynomial over a semialgebraic set defined by polynomial equalities and inequalities. When the polynomial equalities have a finite number of complex solutions and define a radical ideal we can reformulate this problem as a semidefinite programming prob ..."
Abstract

Cited by 37 (7 self)
 Add to MetaCart
We consider the problem of minimizing a polynomial over a semialgebraic set defined by polynomial equalities and inequalities. When the polynomial equalities have a finite number of complex solutions and define a radical ideal we can reformulate this problem as a semidefinite programming problem. This semidefinite program involves combinatorial moment matrices, which are matrices indexed by a basis of the quotient vector space R[x 1 , . . . , x n ]/I. Our arguments are elementary and extend known facts for the grid case including 0/1 and polynomial programming. They also relate to known algebraic tools for solving polynomial systems of equations with finitely many complex solutions. Semidefinite approximations can be constructed by considering truncated combinatorial moment matrices; rank conditions are given (in a grid case) that ensure that the approximation solves the original problem at optimality.
Towards an open curved kernel
 In Proc. Annual ACM Symp. on Computational Geometry
, 2004
"... Our work goes towards answering the growing need for the robust and efficient manipulation of curved objects in numerous applications. The kernel of the cgal library provides several functionalities which are, however, mostly restricted to linear objects. We focus here on the arrangement of conic ar ..."
Abstract

Cited by 31 (14 self)
 Add to MetaCart
Our work goes towards answering the growing need for the robust and efficient manipulation of curved objects in numerous applications. The kernel of the cgal library provides several functionalities which are, however, mostly restricted to linear objects. We focus here on the arrangement of conic arcs in the plane. Our first contribution is the design, implementation and testing of a kernel for computing arrangements of circular arcs. A preliminary C++ implementation exists also for arbitrary conic curves. We discuss the representation and predicates of the geometric objects. Our implementation is targeted for inclusion in the cgal library. Our second contribution concerns exact and efficient algebraic algorithms for the case of conics. They treat all inputs, including degeneracies, and they are implemented as part of the library synaps 2.1. Our tools include Sturm sequences, resultants, Descartes ’ rule, and isolating points. Thirdly, our experiments on circular arcs show that our ∗ Work partially supported by the IST Programme of the EU as a
Almost tight recursion tree bounds for the Descartes method
 In Proc. Int. Symp. on Symbolic and Algebraic Computation
, 2006
"... We give a unified (“basis free”) framework for the Descartes method for real root isolation of squarefree real polynomials. This framework encompasses the usual Descartes ’ rule of sign method for polynomials in the power basis as well as its analog in the Bernstein basis. We then give a new bound ..."
Abstract

Cited by 30 (3 self)
 Add to MetaCart
We give a unified (“basis free”) framework for the Descartes method for real root isolation of squarefree real polynomials. This framework encompasses the usual Descartes ’ rule of sign method for polynomials in the power basis as well as its analog in the Bernstein basis. We then give a new bound on the size of the recursion tree in the Descartes method for polynomials with real coefficients. Applied to polynomials A(X) = P n i=0 aiXi with integer coefficients ai  < 2 L, this yields a bound of O(n(L + log n)) on the size of recursion trees. We show that this bound is tight for L = Ω(log n), and we use it to derive the best known bit complexity bound for the integer case.
Real Algebraic Numbers: Complexity Analysis and Experimentation
 RELIABLE IMPLEMENTATIONS OF REAL NUMBER ALGORITHMS: THEORY AND PRACTICE, LNCS (TO APPEAR
, 2006
"... We present algorithmic, complexity and implementation results concerning real root isolation of a polynomial of degree d, with integer coefficients of bit size ≤ τ, using Sturm (Habicht) sequences and the Bernstein subdivision solver. In particular, we unify and simplify the analysis of both metho ..."
Abstract

Cited by 30 (17 self)
 Add to MetaCart
We present algorithmic, complexity and implementation results concerning real root isolation of a polynomial of degree d, with integer coefficients of bit size ≤ τ, using Sturm (Habicht) sequences and the Bernstein subdivision solver. In particular, we unify and simplify the analysis of both methods and we give an asymptotic complexity bound of eOB(d 4 τ 2). This matches the best known bounds for binary subdivision solvers. Moreover, we generalize this to cover the non squarefree polynomials and show that within the same complexity we can also compute the multiplicities of the roots. We also consider algorithms for sign evaluation, comparison of real algebraic numbers and simultaneous inequalities, and we improve the known bounds at least by a factor of d. Finally, we present our C++ implementation in synaps and some preliminary experiments on various data sets.
M.F.: Bernstein’s basis and real root isolation
, 2004
"... Abstract. In this mostly expository paper we explain how the Bernstein basis, widely used in computeraided geometric design, provides an efficient method for real root isolation, using de Casteljau’s algorithm. We discuss the link between this approach and more classical methods for real root isola ..."
Abstract

Cited by 27 (8 self)
 Add to MetaCart
Abstract. In this mostly expository paper we explain how the Bernstein basis, widely used in computeraided geometric design, provides an efficient method for real root isolation, using de Casteljau’s algorithm. We discuss the link between this approach and more classical methods for real root isolation. We also present a new improved method for isolating real roots in the Bernstein basis inspired by Roullier and Zimmerman.
Fast and Exact Geometric Analysis of Real Algebraic Plane Curves
 Proceedings of the International Symposium on Symbolic and Algebraic Computation (ISSAC’07), 151–158
, 2007
"... An algorithm is presented for the geometric analysis of an algebraic curve f(x, y) = 0 in the real affine plane. It computes a cylindrical algebraic decomposition (CAD) of the plane, augmented with adjacency information. The adjacency information describes the curve’s topology by a topologically eq ..."
Abstract

Cited by 27 (12 self)
 Add to MetaCart
An algorithm is presented for the geometric analysis of an algebraic curve f(x, y) = 0 in the real affine plane. It computes a cylindrical algebraic decomposition (CAD) of the plane, augmented with adjacency information. The adjacency information describes the curve’s topology by a topologically equivalent planar graph. The numerical data in the CAD gives an embedding of the graph. The algorithm is designed to provide the exact result for all inputs but to perform only few symbolic operations for the sake of efficiency. In particular, the roots of f(α, y) at a critical xcoordinate α are found with adaptiveprecision arithmetic in all cases, using a variant of the Bitstream Descartes method (Eigenwillig et al., 2005). The algorithm may choose a generic coordinate system for parts of the analysis but provides its result in the original system. The algorithm is implemented as C++ library AlciX in the EXACUS project. Running time comparisons with top by GonzalezVega and Necula (2002), and with cad2d by Brown demonstrate its efficiency. Categories and Subject Descriptors: I.1.4 [Symbolic and Algebraic Manipulation]: Applications;
Current Issues in SamplingBased Motion Planning
, 2003
"... In this paper, we discuss the field of samplingbased motion planning. In contrast to methods that construct boundary representations of configuration space obstacles, samplingbased methods use only information from a collision detector as they search the configuration space. The simplicity of this ..."
Abstract

Cited by 25 (1 self)
 Add to MetaCart
In this paper, we discuss the field of samplingbased motion planning. In contrast to methods that construct boundary representations of configuration space obstacles, samplingbased methods use only information from a collision detector as they search the configuration space. The simplicity of this approach, along with increases in computation power and the development of efficient collision detection algorithms, has resulted in the introduction of a number of powerful motion planning algorithms, capable of solving challenging problems with many degrees of freedom. First, we trace how samplingbased motion planning has developed. We then discuss a variety of important issues for samplingbased motion planning, including uniform and regular sampling, topological issues, and search philosophies. Finally, we address important issues regarding the role of randomization in samplingbased motion planning.
Counting complexity classes for numeric computations II: Algebraic and semialgebraic sets (Extended Abstract)
 J. COMPL
, 2004
"... We define counting #P classes #P ¡ and in the BlumShubSmale setting of computations over the real or complex numbers, respectively. The problems of counting the number of solutions of systems of polynomial inequalities over ¢ , or of systems of polynomial equalities over £ , respectively, turn ou ..."
Abstract

Cited by 21 (13 self)
 Add to MetaCart
We define counting #P classes #P ¡ and in the BlumShubSmale setting of computations over the real or complex numbers, respectively. The problems of counting the number of solutions of systems of polynomial inequalities over ¢ , or of systems of polynomial equalities over £ , respectively, turn out to be natural complete problems in these classes. We investigate to what extent the new counting classes capture the complexity of computing basic topological invariants of semialgebraic sets (over ¢ ) and algebraic sets (over £). We prove that the problem to compute the (modified) Euler characteristic of semialgebraic sets is FP #P¤complete, and that the problem to compute the geometric degree of complex algebraic sets is FP #P¥complete. We also define new counting complexity classes GCR and GCC in the classical Turing model via taking Boolean parts of the classes above, and show that the problems to compute the Euler characteristic and the geometric degree of (semi)algebraic sets given by integer polynomials are complete in these classes. We complement the results in the Turing model by proving, for all k ¦ ∈ , the FPSPACEhardness of the problem of computing the kth Betti number of the set of real zeros of a given integer polynomial. This holds with respect to the singular homology as well as for the BorelMoore homology.
Revisiting Two Theorems of Curto and Fialkow on Moment Matrices
, 2004
"... We revisit two results of Curto and Fialkow on moment matrices. The first result asserts that every sequence... ..."
Abstract

Cited by 19 (4 self)
 Add to MetaCart
We revisit two results of Curto and Fialkow on moment matrices. The first result asserts that every sequence...