410 research outputs found
TR-2008003: Unified Nearly Optimal Algorithms for Structured Integer Matrices and Polynomials
We seek the solution of banded, Toeplitz, Hankel, Vandermonde, Cauchy and other structured linear systems of equations with integer coefficients. By combining Hensel’s symbolic lifting with either divide-and-conquer algorithms or numerical iterative refinement, we unify the solution for all these structures. We yield the solution in nearly optimal randomized Boolean time, which covers both solution and its correctness verification. Our algorithms and nearly optimal time bounds are extended to the computation of the determinant of a structured integer matrix, its rank and a basis for its null space as well as to some fundamental computations with univariate polynomials that have integer coefficients. Furthermore, we allow to perform lifting modulo a properly bounded power of two t
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
- …