Results 1  10
of
41
SOLVING SYSTEMS OF POLYNOMIAL EQUATIONS
, 2002
"... These are the lecture notes for ten lectures to be given at the CBMS ..."
Abstract

Cited by 150 (12 self)
 Add to MetaCart
These are the lecture notes for ten lectures to be given at the CBMS
Numerical Decomposition of the Solution Sets of Polynomial Systems into Irreducible Components
, 2001
"... In engineering and applied mathematics, polynomial systems arise whose solution sets contain components of different dimensions and multiplicities. In this article we present algorithms, based on homotopy continuation, that compute much of the geometric information contained in the primary decomposi ..."
Abstract

Cited by 56 (26 self)
 Add to MetaCart
In engineering and applied mathematics, polynomial systems arise whose solution sets contain components of different dimensions and multiplicities. In this article we present algorithms, based on homotopy continuation, that compute much of the geometric information contained in the primary decomposition of the solution set. In particular, ignoring multiplicities, our algorithms lay out the decomposition of the set of solutions into irreducible components, by finding, at each dimension, generic points on each component. As byproducts, the computation also determines the degree of each component and an upper bound on itsmultiplicity. The bound issharp (i.e., equal to one) for reduced components. The algorithms make essential use of generic projection and interpolation, and can, if desired, describe each irreducible component precisely as the common zeroesof a finite number of polynomials.
Numerical Homotopies to compute generic Points on positive dimensional Algebraic Sets
 Journal of Complexity
, 1999
"... Many applications modeled by polynomial systems have positive dimensional solution components (e.g., the path synthesis problems for fourbar mechanisms) that are challenging to compute numerically by homotopy continuation methods. A procedure of A. Sommese and C. Wampler consists in slicing the com ..."
Abstract

Cited by 50 (24 self)
 Add to MetaCart
Many applications modeled by polynomial systems have positive dimensional solution components (e.g., the path synthesis problems for fourbar mechanisms) that are challenging to compute numerically by homotopy continuation methods. A procedure of A. Sommese and C. Wampler consists in slicing the components with linear subspaces in general position to obtain generic points of the components as the isolated solutions of an auxiliary system. Since this requires the solution of a number of larger overdetermined systems, the procedure is computationally expensive and also wasteful because many solution paths diverge. In this article an embedding of the original polynomial system is presented, which leads to a sequence of homotopies, with solution paths leading to generic points of all components as the isolated solutions of an auxiliary system. The new procedure significantly reduces the number of paths to solutions that need to be followed. This approach has been implemented and applied to...
Newton’s method with deflation for isolated singularities of polynomial systems
 Theor. Comp. Sci. 359
"... We present a modification of Newton’s method to restore quadratic convergence for isolated singular solutions of polynomial systems. Our method is symbolicnumeric: we produce a new polynomial system which has the original multiple solution as a regular root. We show that the number of deflation sta ..."
Abstract

Cited by 29 (10 self)
 Add to MetaCart
We present a modification of Newton’s method to restore quadratic convergence for isolated singular solutions of polynomial systems. Our method is symbolicnumeric: we produce a new polynomial system which has the original multiple solution as a regular root. We show that the number of deflation stages is bounded by the multiplicity of the isolated root. Our implementation performs well on a large class of applications. 2000 Mathematics Subject Classification. Primary 65H10. Secondary 14Q99, 68W30. Key words and phrases. Newton’s method, deflation, numerical homotopy algorithms, symbolicnumeric computations. 1
Newton’s method for overdetermined systems of equations
 Mathematics of Computation 69 (2000), 1099–1115. MR 2000j:65133
"... Abstract. Complexity theoretic aspects of continuation methods for the solution of square or underdetermined systems of polynomial equations have been studied by various authors. In this paper we consider overdetermined systems where there are more equations than unknowns. We study Newton’s method f ..."
Abstract

Cited by 25 (3 self)
 Add to MetaCart
Abstract. Complexity theoretic aspects of continuation methods for the solution of square or underdetermined systems of polynomial equations have been studied by various authors. In this paper we consider overdetermined systems where there are more equations than unknowns. We study Newton’s method for such a system. I.
A rankrevealing method with updating, downdating and applications
 SIAM J. Matrix Anal. Appl
"... Abstract. A new rank revealing method is proposed. For a given matrix and a threshold for nearzero singular values, by employing a globally convergent iterative scheme as well as a deflation technique the method calculates approximate singular values below the threshold one by one and returns the a ..."
Abstract

Cited by 24 (8 self)
 Add to MetaCart
Abstract. A new rank revealing method is proposed. For a given matrix and a threshold for nearzero singular values, by employing a globally convergent iterative scheme as well as a deflation technique the method calculates approximate singular values below the threshold one by one and returns the approximate rank of the matrix along with an orthonormal basis for the approximate null space. When a row or column is inserted or deleted, algorithms for updating/downdating the approximate rank and null space are straightforward, stable and efficient. Numerical results exhibiting the advantages of our code over existing packages based on twosided orthogonal rankrevealing decompositions are presented. Also presented are applications of the new algorithm in numerical computation of the polynomial GCD as well as identification of nonisolated zeros of polynomial systems.
Numerical Evidence For A Conjecture In Real Algebraic Geometry
, 1998
"... Homotopies for polynomial systems provide computational evidence for a challenging instance of a conjecture about whether all solutions are real. The implementation of SAGBI homotopies involves polyhedral continuation, flat deformation and cheater's homotopy. The numerical difficulties are overcome ..."
Abstract

Cited by 23 (4 self)
 Add to MetaCart
Homotopies for polynomial systems provide computational evidence for a challenging instance of a conjecture about whether all solutions are real. The implementation of SAGBI homotopies involves polyhedral continuation, flat deformation and cheater's homotopy. The numerical difficulties are overcome if we work in the true synthetic spirit of the Schubert calculus by selecting the numerically most favorable equations to represent the geometric problem. Since a wellconditioned polynomial system allows perturbations on the input data without destroying the reality of the solutions we obtain not just one instance, but a whole manifold of systems that satisfy the conjecture. Also an instance that involves totally positive matrices has been verified. The optimality of the solving procedure is a promising first step towards the development of numerically stable algorithms for the pole placement problem in linear systems theory.
Numerical Irreducible Decomposition using PHCpack
, 2003
"... Homotopy continuation methods have proven to be reliable and efficient to approximate all isolated solutions of polynomial systems. In this paper we show how we can use this capability as a blackbox device to solve systems which have positive dimensional components of solutions. We indicate how the ..."
Abstract

Cited by 21 (14 self)
 Add to MetaCart
Homotopy continuation methods have proven to be reliable and efficient to approximate all isolated solutions of polynomial systems. In this paper we show how we can use this capability as a blackbox device to solve systems which have positive dimensional components of solutions. We indicate how the software package PHCpack can be used in conjunction with Maple and programs written in C. We describe a numerically stable algorithm for decomposing positive dimensional solution sets of polynomial systems into irreducible components.
Polyhedral End Games for Polynomial Continuation
 Numerical Algorithms
, 1998
"... Bernshtein's theorem provides a generically exact upper bound on the number of isolated solutions a sparse polynomial system can have in (C ) n , with C = C n f0g. When a sparse polynomial system has fewer than this number of isolated solutions some face system must have solutions in (C ..."
Abstract

Cited by 21 (8 self)
 Add to MetaCart
Bernshtein's theorem provides a generically exact upper bound on the number of isolated solutions a sparse polynomial system can have in (C ) n , with C = C n f0g. When a sparse polynomial system has fewer than this number of isolated solutions some face system must have solutions in (C ) n . In this paper we address the process of recovering a certificate of deficiency from a diverging solution path. This certificate takes the form of a face system along with approximations of its solutions. We apply extrapolation to estimate the cycle number and the face normal. Applications illustrate the practical usefulness of our approach. keywords : homotopy continuation, polynomial systems, Newton polytopes, Bernshtein bound, cycle number. AMS(MOS) Classification : 14Q99, 52A39, 52B20, 65H10. 1 Introduction All isolated complex solutions to polynomial systems can be approximated numerically by homotopy continuation methods. The strategy is to set up a collection of implicitly d...
Homotopies for intersecting solution components of polynomial systems
 SIAM J. Numer. Anal
, 2004
"... Abstract. We show how to use numerical continuation to compute the intersection C = A∩B of two algebraic sets A and B, where A, B, and C are numerically represented by witness sets. Enroute to this result, we first show how to find the irreducible decomposition of a system of polynomials restricted ..."
Abstract

Cited by 19 (13 self)
 Add to MetaCart
Abstract. We show how to use numerical continuation to compute the intersection C = A∩B of two algebraic sets A and B, where A, B, and C are numerically represented by witness sets. Enroute to this result, we first show how to find the irreducible decomposition of a system of polynomials restricted to an algebraic set. The intersection of components A and B then follows by considering the decomposition of the diagonal system of equations u − v = 0 restricted to {u, v} ∈ A × B. One offshoot of this new approach is that one can solve a large system of equations by finding the solution components of its subsystems and then intersecting these. It also allows one to find the intersection of two components of the two polynomial systems, which is not possible with any previous numerical continuation approach.