... In this article, we develop a theory of @-finite sequences and functions which provides a unified framework to express algorithms proving and discovering multivariate identities. This approach is vindicated by an implementation.

Given a system of polynomial equations and inequations with coefficients in the field of rational numbers, we show how to compute a geometric resolution of the set of common roots of the system over the field of complex numbers. A geometric resolution consists of a primitive element of the algebraic extension defined by the set of roots, its minimal polynomial and the parametrizations of the coordinates. Such a representation of the solutions has a long history which goes back to Leopold Kronecker and has been revisited many times in computer algebra. We introduce a new generation of probabilistic algorithms where all the computations use only univariate or bivariate polynomials. We give a new codification of the set of solutions of a positive dimensional algebraic variety relying on a new global version of Newton’s iterator. Roughly speaking the complexity of our algorithm is polynomial in some kind of degree of the system, in its height, and linear in the complexity of evaluation

We discuss the use of a single generic linear combination of multiplication matrices, and its reordered Schur factorization, to find the roots of a system of multivariate polynomial equations. The principal contribution of the paper is to show how to reduce the multivariate problem to a univariate problem, even in the case of multiple roots, in a numerically stable way. 1 Introduction The technique of solving systems of multivariate polynomial equations via eigenproblems has become a topic of active research (with applications in computer-aided design and control theory, for example) at least since the papers [2, 6, 9]. One may approach the problem via various resultant formulations or by Grobner bases. As more understanding is gained, it is becoming clearer that eigenvalue problems are the "weakly nonlinear nucleus to which the original, strongly nonlinear task may be reduced"[13]. Early works concentrated on the case of simple roots. An example of such was the paper [5], which use...

We present an algorithm which converts a given Gröbner basis of a polynomial ideal I to a Gröbner basis of I with respect to another term order. The conversion is done in several steps following a path in the Gröbner fan of I. Each conversion step is based on the computation of a Gröbner basis of a toric degeneration of I. c ○ 1997 Academic Press Limited 1.

There are very close connections between the arithmetic of integer lattices, algebraic properties of the associated ideals, and the geometry and the combinatorics of corresponding polyhedra. In this paper we investigate the generating sets ("Gröbner bases") of integer lattices that correspond to the Gröbner bases of the associated binomial ideals. Extending results by Sturmfels & Thomas, we obtain a geometric characterization of the universal Gröbner basis in terms of the vertices and edges of the associated corner polyhedra. In the special case where the lattice has finite index, the corner polyhedra were studied by Gomory, and there is a close connection to the "group problem in integer programming." We present exponential lower and upper bounds for the maximal size of a reduced Grobner basis. The initial complex of (the ideal of) a lattice is shown to be dual to the boundary of a certain simple polyhedron.

In this paper we describe some experimental findings on selection strategies for Gröbner basis computation with the Buchberger algorithm. In particular, the results suggest that the "sugar flavor" of the "normal selection", implemented first in CoCoA, then in AlPI, and now in SCRATCHPAD-II, is the best choice for a selection strategy. It has to be combined with the "straight-forward" simplification strategy and with a special form of the Gebauer-Möller criteria to obtain the best results. The idea of the "sugar flavor" is the following: the Buchberger algorithm for homogeneous ideals, with degree-compatible term ordering and normal selection strategy, usually works fine. Homogenizing the basis of the ideal is good for the strategy, but bad for the basis to be computed. The sugar flavor computes, for every polynomial in the course of the algorithm, "the degree that it would have if computed with the homogeneous algorithm", and uses this phantom degree (the sugar) only for the selection strategy. We have tested several examples with different selection strategies, and the sugar flavor has proved to be always the best choice or very near to it. The comparison between the different variants of the sugar flavor has been made, but the results are up to now inconclusive. We include a complete deterministic description of the Buchberger algorithm as it was used in our experiments.

Abstract. We present two algorithms for the arithmetic of cubic curves with a totally ramified prime at infinity. The first algorithm, inspired by Cantor’s reduction for hyperelliptic curves, is easily implemented with a few lines of code, making use of a polynomial arithmetic package. We prove explicit reducedness criteria for superelliptic curves of genus 3 and 4, which show the correctness of the algorithm. The second approach, quite general in nature and applicable to further classes of curves, uses the FGLM algorithm for switching between Gröbner bases for different orderings. Carrying out the computations symbolically, we obtain explicit reduction formulae in terms of the input data. 1.

Abstract not found

Given a polynomial system of n equations in n unknowns that depends on some parameters, we de ne the notion of parametric geometric resolution as a means to represent some generic solutions in terms of the parameters. The coefficients of this resolution are rational functions of the parameters; we first show that their degree is bounded by the Bézout number d n , where d is a bound on the degrees of the input system. We then present a probabilistic algorithm to compute such a resolution; in short, its complexity is polynomial in the size of the output and the probability of success is controlled by a quantity polynomial in the Bézout number. We present several applications of this process, to computations in the Jacobian of hyperelliptic curves and to questions of real geometry.

We discuss changing the variable ordering for a regular chain in positive dimension. This quite general question has applications going from implicitization problems to the symbolic resolution of some systems of differential algebraic equations. We propose a modular method, reducing the problem to dimension zero and using Newton-Hensel lifting techniques. The problems raised by the choice of the specialization points, the lack of the (crucial) information of what are the free and algebraic variables for the new ordering, and the efficiency regarding the other methods are discussed. Strong hypotheses (but not unusual) for the initial regular chain are required. Change of ordering in dimension zero is taken as a subroutine.