3,999 research outputs found
Solving polynomial eigenvalue problems by means of the Ehrlich-Aberth method
Given the matrix polynomial , we
consider the associated polynomial eigenvalue problem. This problem, viewed in
terms of computing the roots of the scalar polynomial , is treated
in polynomial form rather than in matrix form by means of the Ehrlich-Aberth
iteration. The main computational issues are discussed, namely, the choice of
the starting approximations needed to start the Ehrlich-Aberth iteration, the
computation of the Newton correction, the halting criterion, and the treatment
of eigenvalues at infinity. We arrive at an effective implementation which
provides more accurate approximations to the eigenvalues with respect to the
methods based on the QZ algorithm. The case of polynomials having special
structures, like palindromic, Hamiltonian, symplectic, etc., where the
eigenvalues have special symmetries in the complex plane, is considered. A
general way to adapt the Ehrlich-Aberth iteration to structured matrix
polynomial is introduced. Numerical experiments which confirm the effectiveness
of this approach are reported.Comment: Submitted to Linear Algebra App
Multiplicity estimates, analytic cycles and Newton polytopes
We consider the problem of estimating the multiplicity of a polynomial when
restricted to the smooth analytic trajectory of a (possibly singular)
polynomial vector field at a given point or points, under an assumption known
as the D-property. Nesterenko has developed an elimination theoretic approach
to this problem which has been widely used in transcendental number theory.
We propose an alternative approach to this problem based on more local
analytic considerations. In particular we obtain simpler proofs to many of the
best known estimates, and give more general formulations in terms of Newton
polytopes, analogous to the Bernstein-Kushnirenko theorem. We also improve the
estimate's dependence on the ambient dimension from doubly-exponential to an
essentially optimal single-exponential.Comment: Some editorial modifications to improve readability; No essential
mathematical change
The complexity and geometry of numerically solving polynomial systems
These pages contain a short overview on the state of the art of efficient
numerical analysis methods that solve systems of multivariate polynomial
equations. We focus on the work of Steve Smale who initiated this research
framework, and on the collaboration between Stephen Smale and Michael Shub,
which set the foundations of this approach to polynomial system--solving,
culminating in the more recent advances of Carlos Beltran, Luis Miguel Pardo,
Peter Buergisser and Felipe Cucker
Counting Solutions of a Polynomial System Locally and Exactly
We propose a symbolic-numeric algorithm to count the number of solutions of a
polynomial system within a local region. More specifically, given a
zero-dimensional system , with
, and a polydisc
, our method aims to certify the existence
of solutions (counted with multiplicity) within the polydisc.
In case of success, it yields the correct result under guarantee. Otherwise,
no information is given. However, we show that our algorithm always succeeds if
is sufficiently small and well-isolating for a -fold
solution of the system.
Our analysis of the algorithm further yields a bound on the size of the
polydisc for which our algorithm succeeds under guarantee. This bound depends
on local parameters such as the size and multiplicity of as well
as the distances between and all other solutions. Efficiency of
our method stems from the fact that we reduce the problem of counting the roots
in of the original system to the problem of solving a
truncated system of degree . In particular, if the multiplicity of
is small compared to the total degrees of the polynomials ,
our method considerably improves upon known complete and certified methods.
For the special case of a bivariate system, we report on an implementation of
our algorithm, and show experimentally that our algorithm leads to a
significant improvement, when integrated as inclusion predicate into an
elimination method
Near Optimal Subdivision Algorithms for Real Root Isolation
We describe a subroutine that improves the running time of any subdivision
algorithm for real root isolation. The subroutine first detects clusters of
roots using a result of Ostrowski, and then uses Newton iteration to converge
to them. Near a cluster, we switch to subdivision, and proceed recursively. The
subroutine has the advantage that it is independent of the predicates used to
terminate the subdivision. This gives us an alternative and simpler approach to
recent developments of Sagraloff (2012) and Sagraloff-Mehlhorn (2013), assuming
exact arithmetic.
The subdivision tree size of our algorithm using predicates based on
Descartes's rule of signs is bounded by , which is better by
compared to known results. Our analysis differs in two key
aspects. First, we use the general technique of continuous amortization from
Burr-Krahmer-Yap (2009), and second, we use the geometry of clusters of roots
instead of the Davenport-Mahler bound. The analysis naturally extends to other
predicates.Comment: 19 pages, 3 figure
Multiplicity Estimates: a Morse-theoretic approach
The problem of estimating the multiplicity of the zero of a polynomial when
restricted to the trajectory of a non-singular polynomial vector field, at one
or several points, has been considered by authors in several different fields.
The two best (incomparable) estimates are due to Gabrielov and Nesterenko.
In this paper we present a refinement of Gabrielov's method which
simultaneously improves these two estimates. Moreover, we give a geometric
description of the multiplicity function in terms certain naturally associated
polar varieties, giving a topological explanation for an asymptotic phenomenon
that was previously obtained by elimination theoretic methods in the works of
Brownawell, Masser and Nesterenko. We also give estimates in terms of Newton
polytopes, strongly generalizing the classical estimates.Comment: Minor revision; To appear in Duke Math. Journa
- …