Symmetric Functions Applied to Decomposing Solution Sets of Polynomial Systems

  • Authors:
  • Andrew J. Sommese;Jan Verschelde;Charles W. Wampler

  • Affiliations:
  • -;-;-

  • Venue:
  • SIAM Journal on Numerical Analysis
  • Year:
  • 2002

Quantified Score

Hi-index 0.01

Visualization

Abstract

Many polynomial systems have solution sets comprised of multiple irreducible components, possibly of different dimensions. A fundamental problem of numerical algebraic geometry is to decompose such a solution set, using floating-point numerical processes, into its components. Prior work has shown how to generate sets of generic points guaranteed to include points from every component. Furthermore, we have shown how monodromy can be used to efficiently predict the partition of these points by membership in the components. However, confirmation of this prediction required an expensive procedure of sampling each component to find an interpolating polynomial that vanishes on it. This paper proves theoretically and demonstrates in practice that linear traces suffice for this verification step, which gives great improvement in both computational speed and numerical stability. Moreover, in the case that one may still wish to compute an interpolating polynomial, we show how to do so more efficiently by building a structured grid of samples, using divided differences, and applying symmetric functions. Several test problems illustrate the effectiveness of the new methods.