Tangent Graeffe iteration

From MaRDI portal



Abstract: Graeffe iteration was the choice algorithm for solving univariate polynomials in the XIX-th and early XX-th century. In this paper, a new variation of Graeffe iteration is given, suitable to IEEE floating-point arithmetics of modern digital computers. We prove that under a certain generic assumption the proposed algorithm converges. We also estimate the error after N iterations and the running cost. The main ideas from which this algorithm is built are: classical Graeffe iteration and Newton Diagrams, changes of scale (renormalization), and replacement of a difference technique by a differentiation one. The algorithm was implemented successfully and a number of numerical experiments are displayed.


The authors present a version of renormalized Graeffe's algorithm, also considered in their previous paper [J. Complexity 17, No. 3, 541-573 (2001; reviewed above)], for solving real and complex univariate polynomials. Recall that Graeffe iteration maps a degree \(d\) polynomial \(f(x)\) into the degree \(d\) polynomial \[ Gf(x)=(-1)^d f(\sqrt{x}) f(-\sqrt{x}). \] and that, if \(\zeta_1,\zeta_2,\ldots,\zeta_d\) are the roots of \(f\), then the roots of \(Gf\) are \(\zeta_1^2,\zeta_2^2,\ldots,\zeta_d^2.\)











This page was built for publication: Tangent Graeffe iteration

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5952132)