Construction of iteration functions for the simultaneous computation of the solutions of equations and algebraic systems

We construct iteration functions for the simultaneous computation of the solutions of a system of equations, with local quadratic convergence: they generalize to the multivariate case the well-known Weierstrass function for polynomials, which is expected to be globally convergent except on a zero-measured set of starting points. We clarify these functions using univariate interpolation. Both for polynomials and algebraic systems with real coefficients, we extend the conjecture of global convergence to the research of real roots or solutions.     

On the shape of solution sets of systems of (functional) equations

Solution sets of systems of linear equations over fields are characterized as being affine subspaces. But what can we say about the “shape” of the set of all solutions of other systems of equations? We study solution sets over arbitrary algebraic structures, and we give a necessary condition for a set of n-tuples to be the set of solutions of a system of equations in n unknowns over a given algebra. In the case of Boolean equations we obtain a complete characterization, and we also characterize solution sets of systems of Boolean functional equations.     

Degree formulas and signature formulas for the Euler characteristic of real algebraic sets

We describe formulas, due to several authors, for the computation of the Euler characteristic of real algebraic sets. We also give a new result for real algebraic sets with a compact set of singular points.     

Geometric control of hybrid systems

In this paper, we present a geometric approach for computing controlled invariant sets for hybrid control systems. While the problem is well studied in the ellipsoidal case, this family is quite conservative for constrained or switched linear systems. We reformulate the invariance of a set as an inequality for its support function that is valid for any convex set. This produces novel algebraic conditions for the invariance of sets with polynomial or piecewise quadratic support functions.     

Complexity of solving parametric polynomial systems

In this paper, we present three algorithms: the first one solves zero-dimensional parametric homogeneous polynomial systems within single exponential time in the number n of unknowns; it decomposes the parameter space into a finite number of constructible sets and computes the finite number of solutions by parametric rational representations uniformly in each constructible set. The second algorithm factirizes absolutely multivariate parametic polynomials within single exponential time in n and in the upper bound d on the degree of the factorized polynomials. The third algorithm decomposes algebraic varieties defined by parametric polynomial systems of positive dimension into absolutely irreducible components uniformly in the values of the parameters. The complexity bound for this algorithm is double exponential in n. On the other hand, the lower bound on the complexity of the problem of resolution of parametric polynomial systems is double exponential in n. Bibliography: 72 titles.     

The Hardness of Polynomial Equation Solving

Elimination theory was at the origin of algebraic geometry in the nineteenth century and now deals with the algorithmic solving of multivariate polynomial equation systems over the complex numbers or, more generally, over an arbitrary algebraically closed field. In this paper we investigate the intrinsic sequential time complexity of universal elimination procedures for arbitrary continuous data structures encoding input and output objects of elimination theory (i.e., polynomial equation systems) and admitting the representation of certain limit objects. Our main result is the following: let there be given such a data structure and together with this data structure a universal elimination algorithm, say P, solving arbitrary parametric polynomial equation systems. Suppose that the algorithm P avoids unnecessary branchings and that P admits the efficient computation of certain natural limit objects (as, e.g., the Zariski closure of a given constructible algebraic set or the parametric greatest common divisor of two given algebraic families of univariate polynomials). Then P$ cannot be a polynomial time algorithm. The paper contains different variants of this result and discusses their practical implications.     

A multivariate QD-like algorithm

The problem of constructing a univariate rational interpolant or Padé approximant for given data can be solved in various equivalent ways: one can compute the explicit solution of the system of interpolation or approximation conditions, or one can start a recursive algorithm, or one can obtain the rational function as the convergent of an interpolating or corresponding continued fraction.In case of multivariate functions general order systems of interpolation conditions for a multivariate rational interpolant and general order systems of approximation conditions for a multivariate Padé approximant were respectively solved in [6] and [9]. Equivalent recursive computation schemes were given in [3] for the rational interpolation case and in [5] for the Padé approximation case. At that moment we stated that the next step was to write the general order rational interpolants and Padé approximants as the convergent of a multivariate continued fraction so that the univariate equivalence of the three main defining techniques was also established for the multivariate case: algebraic relations, recurrence relations, continued fractions. In this paper a multivariate qd-like algorithm is developed that serves this purpose.     

Some classes of abstract simplicial complexes motivated by module theory

In this paper we analyze some classes of abstract simplicial complexes relying on algebraic models arising from module theory. To this regard, we consider a left-module on a unitary ring and find models of abstract complexes and related set operators having specific regularity properties, which are strictly interrelated to the algebraic properties of both the module and the ring.Next, taking inspiration from the aforementioned models, we carry out our analysis from modules to arbitrary sets. In such a more general perspective, we start with an abstract simplicial complex and an associated set operator. Endowing such a set operator with the corresponding properties obtained in our module instances, we investigate in detail and prove several properties of three subclasses of abstract complexes.More specifically, we provide uniformity conditions in relation to the cardinality of the maximal members of such classes. By means of the notion of OSS-bijection, we prove a correspondence theorem between a subclass of closure operators and one of the aforementioned families of abstract complexes, which is similar to the classic correspondence theorem between closure operators and Moore systems. Next, we show an extension property of a binary relation induced by set systems when they belong to one of the above families.Finally, we provide a representation result in terms of pairings between sets for one of the three classes of abstract simplicial complexes studied in this work.     

分片代数集合的不可约性和同构

本文利用多元样条函数来定义分片代数集合,讨论了分片代数集合的不可约性和同构问题,给出了分片代数集合不可约的两个等价条件,并把分片代数集合的同构分类问题转化为交换代数的同构分类问题。     

Topology of injective endomorphisms of real algebraic sets

Using only basic topological properties of real algebraic sets and regular morphisms we show that any injective regular self-mapping of a real algebraic set is surjective. Then we show that injective morphisms between germs of real algebraic sets define a partial order on the equivalence classes of these germs divided by continuous semi-algebraic homeomorphisms. We use this observation to deduce that any injective regular self-mapping of a real algebraic set is a homeomorphism. We show also a similar local property. All our results can be extended to arc-symmetric semi-algebraic sets and injective continuous arc-symmetric morphisms, and some results to Euler semi-algebraic sets and injective continuous semi-algebraic morphisms.Mathematics Subject Classification (2000):14Pxx, 14A10, 32B10     

零误差计算

研究采用有误差的数值计算来获得无误差的准确值具有重要的理论价值和应用价值.这种通过近似的数值方法获得准确结果的计算被称为零误差计算.本文首先指出,只有一致离散集合中的数才能够开展零误差计算,即有非零隔离界的数集,这也是"数"可以进行零误差计算的一个充要条件.以此为基本出发点,本文分析代数数零误差计算的最低理论精度,该精度对应于恢复近似代数数的准确值时必要的误差控制条件,但由于所采用恢复算法的局限性,这一理论精度往往不能保证成功恢复出代数数的准确值.为此,本文给出采用PSLQ (partial-sum-LQ-decomposition)算法进行代数数零误差计算所需的精度控制条件,与基于LLL (Lenstra-Lenstra-Lovász)算法相比,该精度控制条件关于代数数次数的依赖程度由二次降为拟线性,从而可降低相应算法的复杂度.最后探讨零误差计算未来的发展趋势.     

The Nother and Riemann-Roch type theorems for piecewise algebraic curve

A piecewise algebraic curve is a curve determined by the zero set of a bivariate spline function. In this paper, the Nother type theorems for Cμpiecewise algebraic curves are obtained. The theory of the linear series of sets of places on the piecewise algebraic curve is also established. In this theory, singular cycles are put into the linear series, and a complete series of the piecewise algebraic curves consists of all effective ordinary cycles in an equivalence class and all effective singular cycles which are equivalent specifically to any effective ordinary cycle in the equivalence class. This theory is a generalization of that of linear series of the algebraic curve. With this theory and the fundamental theory of multivariate splines on smoothing cofactors and global conformality conditions, and the results on the general expression of multivariate splines, we get a formula on the index, the order and the dimension of a complete series of the irreducible Cμpiecewise algebraic curves and the degree, the genus and the smoothness of the curves, hence the Riemann-Roch type theorem of the Cμpiecewise algebraic curve is established.     

Structure indices for multidimensional systems

The structure indices of a one-dimensional system are an importantset of invariants. In this paper we examine a generalizationof this concept to multidimensional linear systems, which correspondsto the algebraic concept of a Hilbert series. We use the standardtheory of the Hilbert series to explain some of the previousID system-theoretic results. We discuss the computation of nDstructure indices from an initial condition set, and the invariantswhich can be derived from these indices.     

Computations in differential and difference modules

Constructive methods based on the Gröbner bases theory have been used many times in commutative algebra over the past 20 years, in particular, they allow the computation of such important invariants of manifolds given by systems of algebraic equations as their Hilbert polynomials. In differential and difference algebra, the analogous roles play characteristic sets.In this paper, algorithms for computations in differential and difference modules, which allow for the computation of characteristic sets (Gröbner bases) in differential, difference, and polynomial modules and differential (difference) dimension polynomials, are described. The algorithms are implemented in the algorithmic language REFAL.     

Means on dyadic symmetrie sets and polar decompositions

In this paper we develop the theory of the geometric mean and the spectral mean on dyadic symmetric sets, an algebraic generalization of symmetric spaces of noncompact type, and apply them to obtain decomposition theorems of involutive systems. In particular we show for involutive dyadic symmetric sets: every involutive dyadic symmetric set admits a canonical polar decomposition with factors the geometric and spectral means.     

Stochastic Arithmetic: s-spaces and Some Applications

It has been recently shown that computation with stochastic numbers as regard to addition and multiplication by scalars can be reduced to computation in familiar vector spaces. In this work we show how this can be used for the algebraic solution of linear systems of equations with stochastic right-hand sides. On several examples we compare the algebraic solution with the simulated solution using the CADNA package.     

The Nöther and Riemann-Roch type theorems for piecewise algebraic curve

A piecewise algebraic curve is a curve determined by the zero set of a bivariate spline function. In this paper, the Nöther type theorems for C ^µ piecewise algebraic curves are obtained. The theory of the linear series of sets of places on the piecewise algebraic curve is also established. In this theory, singular cycles are put into the linear series, and a complete series of the piecewise algebraic curves consists of all effective ordinary cycles in an equivalence class and all effective singular cycles which are equivalent specifically to any effective ordinary cycle in the equivalence class. This theory is a generalization of that of linear series of the algebraic curve. With this theory and the fundamental theory of multivariate splines on smoothing cofactors and global conformality conditions, and the results on the general expression of multivariate splines, we get a formula on the index, the order and the dimension of a complete series of the irreducible C ^µ piecewise algebraic curves and the degree, the genus and the smoothness of the curves, hence the Riemann-Roch type theorem of the C ^µ piecewise algebraic curve is established.     

Isosingular Sets and Deflation

This article introduces the concept of isosingular sets, which are irreducible algebraic subsets of the set of solutions to a system of polynomial equations constructed by taking the closure of points with a common singularity structure. The definition of these sets depends on deflation, a procedure that uses differentiation to regularize solutions. A weak form of deflation has proven useful in regularizing algebraic sets, making them amenable to treatment by the algorithms of numerical algebraic geometry. We introduce a strong form of deflation and define deflation sequences, which, in a different context, are the sequences arising in Thom–Boardman singularity theory. We then define isosingular sets in terms of deflation sequences. We also define the isosingular local dimension and examine the properties of isosingular sets. While isosingular sets are of theoretical interest as constructs for describing singularity structures of algebraic sets, they also expand the kinds of algebraic set that can be investigated with methods from numerical algebraic geometry.     

Geometry and asymptotics of wavefronts

Methods of convex analysis and differential geometry are applied to the study of properties of nonconvex sets in the plane. Constructions of the theory of α-sets are used as a tool for investigation of problems of the control theory and the theory of differential games. The notions of the bisector and of a pseudovertex of a set introduced in the paper, which allow ones to study the geometry of sets and compute their measure of nonconvexity, are of independent interest. These notions are also useful in studies of evolution of sets of attainability of controllable systems and in constructing of wavefronts. In this paper, we develop a numerically-analytical approach to finding pseudovertices of a curve, computation of the measure of nonconvexity of a plane set, and constructing front sets on the basis these data.In the paper, we give the results of numerical constructing of bisectors and wavefronts for plane sets. We demonstrate the relation between nonsmoothness of wavefronts and singularity of the geometry of the initial set. We also single out a class of sets whose bisectors have an asymptote.