11,135 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
Spectral Solution with a Subtraction Method to Improve Accuracy
This work addresses the solution to a Dirichlet boundary value problem for the Poisson equation in 1-D, d2u/dx2 = f using a numerical Fourier collocation approach. The order of accuracy of this approach can be increased by modifying f so the periodic extension of the right hand side is suffciently smooth. A proof for the order is given by Sköllermo. This work introduces a subtraction technique to modify the function\u27s right hand side to reduce the discontinuities or improve the smoothness of its periodic extension. This subtraction technique consists of cosine polynomials found by using boundary derivatives. We subtract cosine polynomials to match boundary values and derivatives of f. The derivatives need only be calculated numerically and approximately represent derivatives at the boundaries. Increasing the number of cosine polynomials in the subtraction technique increases the order of accuracy of the solution. The use of cosine polynomials matches well with the Fourier transform approach and is computationally efficient. Implementation of this technique results in a solution with variable accuracy depending on the number of collocation points and approximated boundary derivatives. Results show that the technique can be up to 14th order accurate
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
List decoding of a class of affine variety codes
Consider a polynomial in variables and a finite point ensemble . When given the leading monomial of with respect to
a lexicographic ordering we derive improved information on the possible number
of zeros of of multiplicity at least from . We then use this
information to design a list decoding algorithm for a large class of affine
variety codes.Comment: 11 pages, 5 table
New Acceleration of Nearly Optimal Univariate Polynomial Root-findERS
Univariate polynomial root-finding has been studied for four millennia and is
still the subject of intensive research. Hundreds of efficient algorithms for
this task have been proposed. Two of them are nearly optimal. The first one,
proposed in 1995, relies on recursive factorization of a polynomial, is quite
involved, and has never been implemented. The second one, proposed in 2016,
relies on subdivision iterations, was implemented in 2018, and promises to be
practically competitive, although user's current choice for univariate
polynomial root-finding is the package MPSolve, proposed in 2000, revised in
2014, and based on Ehrlich's functional iterations. By proposing and
incorporating some novel techniques we significantly accelerate both
subdivision and Ehrlich's iterations. Moreover our acceleration of the known
subdivision root-finders is dramatic in the case of sparse input polynomials.
Our techniques can be of some independent interest for the design and analysis
of polynomial root-finders.Comment: 89 pages, 5 figures, 2 table
- …