921 research outputs found
Irredundant Triangular Decomposition
Triangular decomposition is a classic, widely used and well-developed way to
represent algebraic varieties with many applications. In particular, there
exist sharp degree bounds for a single triangular set in terms of intrinsic
data of the variety it represents, and powerful randomized algorithms for
computing triangular decompositions using Hensel lifting in the
zero-dimensional case and for irreducible varieties. However, in the general
case, most of the algorithms computing triangular decompositions produce
embedded components, which makes it impossible to directly apply the intrinsic
degree bounds. This, in turn, is an obstacle for efficiently applying Hensel
lifting due to the higher degrees of the output polynomials and the lower
probability of success. In this paper, we give an algorithm to compute an
irredundant triangular decomposition of an arbitrary algebraic set defined
by a set of polynomials in C[x_1, x_2, ..., x_n]. Using this irredundant
triangular decomposition, we were able to give intrinsic degree bounds for the
polynomials appearing in the triangular sets and apply Hensel lifting
techniques. Our decomposition algorithm is randomized, and we analyze the
probability of success
On the complexity of computing with zero-dimensional triangular sets
We study the complexity of some fundamental operations for triangular sets in
dimension zero. Using Las-Vegas algorithms, we prove that one can perform such
operations as change of order, equiprojectable decomposition, or quasi-inverse
computation with a cost that is essentially that of modular composition. Over
an abstract field, this leads to a subquadratic cost (with respect to the
degree of the underlying algebraic set). Over a finite field, in a boolean RAM
model, we obtain a quasi-linear running time using Kedlaya and Umans' algorithm
for modular composition. Conversely, we also show how to reduce the problem of
modular composition to change of order for triangular sets, so that all these
problems are essentially equivalent. Our algorithms are implemented in Maple;
we present some experimental results
Reverse engineering of CAD models via clustering and approximate implicitization
In applications like computer aided design, geometric models are often
represented numerically as polynomial splines or NURBS, even when they
originate from primitive geometry. For purposes such as redesign and
isogeometric analysis, it is of interest to extract information about the
underlying geometry through reverse engineering. In this work we develop a
novel method to determine these primitive shapes by combining clustering
analysis with approximate implicitization. The proposed method is automatic and
can recover algebraic hypersurfaces of any degree in any dimension. In exact
arithmetic, the algorithm returns exact results. All the required parameters,
such as the implicit degree of the patches and the number of clusters of the
model, are inferred using numerical approaches in order to obtain an algorithm
that requires as little manual input as possible. The effectiveness, efficiency
and robustness of the method are shown both in a theoretical analysis and in
numerical examples implemented in Python
Computing Equilibria of Semi-algebraic Economies Using Triangular Decomposition and Real Solution Classification
In this paper, we are concerned with the problem of determining the existence
of multiple equilibria in economic models. We propose a general and complete
approach for identifying multiplicities of equilibria in semi-algebraic
economies, which may be expressed as semi-algebraic systems. The approach is
based on triangular decomposition and real solution classification, two
powerful tools of algebraic computation. Its effectiveness is illustrated by
two examples of application.Comment: 24 pages, 5 figure
- …