We propose a new algorithm fur the classical and still practically important problem of approximating zeros zi of an nth degree polynomial p(x) within error bound 2(-b) max(j) [z(j)]. The algorithm uses O((n(2) log n)...
详细信息
We propose a new algorithm fur the classical and still practically important problem of approximating zeros zi of an nth degree polynomial p(x) within error bound 2(-b) max(j) [z(j)]. The algorithm uses O((n(2) log n) log(bn)) arithmetic operations and comparisons for approximating all the n zeros and O((kn log n) log(bn)) for approximating the it zeros lying in a fixed domain (disc or square) and isolated from the other zeros. Unlike the previous fast algorithms of this kind, the new algorithm has its simple elementary description, is convenient for practical implementation, and allows the users to adapt the computational precision to the current level of approximation achieved in the process of computing and ultimately to the requirements to the output precision for each zero of p(x). The algorithm relies on our novel versions of Weyl's qundtree construction and Newton's iteration. (C) 2000 Academic Press.
We substantially improve the known algorithms for approximating all the complexzeros of an n(th) degree polynomial p(x). Our new algorithms save both Boolean and arithmetic sequential time, versus the previous best a...
详细信息
We substantially improve the known algorithms for approximating all the complexzeros of an n(th) degree polynomial p(x). Our new algorithms save both Boolean and arithmetic sequential time, versus the previous best algorithms of Schonhage [1], Pan [2], and Neff and Pelf [3]. In parallel (NC) implementation, we dramatically decrease the number of processors, versus the parallel algorithm of Neff [4], which was the only NC algorithm known for this problem so far. Specifically, under the simple normalization assumption that the variable a: has been scaled so as to confine the zeros of p(x) to the unit disc {x : \x\ less than or equal to 1}, our algorithms (which promise to be practically effective) approximate all the zeros of p(a) within the absolute error bound 2(-b), by using order of n arithmetic operations and order of (b + n)n(2) Boolean (bitwise) operations (in both cases up to within polylogarithmic factors). The algorithms allow their optimal (work preserving) NC parallelization, so that they can be implemented by using polylogarithmic time and the orders of n arithmetic processors or (b + n)n(2) Boolean processors. All the cited bounds on the computational complexity are within polylogarithmic factors from the optimum (in terms of n and b) under both arithmetic and Boolean models of computation (in the Boolean case, under the additional (realistic) assumption that n = O(b)).
The classical problem of solving an nth degree polynomial equation has substantially influenced the development of mathematics throughout the centuries and still has several important applications to the theory and pr...
详细信息
The classical problem of solving an nth degree polynomial equation has substantially influenced the development of mathematics throughout the centuries and still has several important applications to the theory and practice of present-day computing. We briefly recall the history of the algorithmic approach to this problem and then review some successful solution algorithms. We end by outlining some algorithms of 1995 that solve this problem at a surprisingly low computational cost.
暂无评论