Forward stable computation of roots of real polynomials with only real distinct roots

Ivan Slapničar, Nevena Jakovčević Stor

Publication date: 21 September 2015

Abstract: As showed in (Fiedler, 1990), any polynomial can be expressed as a characteristic polynomial of a complex symmetric arrowhead matrix. This expression is not unique. If the polynomial is real with only real distinct roots, the matrix can be chosen real. By using accurate forward stable algorithm for computing eigenvalues of real symmetric arrowhead matrices from (Jakovcevic Stor, Slapnicar, Barlow, 2015), we derive a forward stable algorithm for computation of roots of such polynomials in

O (n^{2})

operations. The algorithm computes each root to almost full accuracy. In some cases, the algorithm invokes extended precision routines, but only in the non-iterative part. Our examples include numerically difficult problems, like the well-known Wilkinson's polynomials. Our algorithm compares favourably to other method for polynomial root-finding, like MPSolve or Newton's method.

Has companion code repository: https://github.com/UnofficialJuliaMirror/Arrowhead.jl-1a029416-38b0-5245-b7ed-67249edcb994

Mathematics Subject Classification ID

Numerical computation of eigenvalues and eigenvectors of matrices (65F15) Roundoff error (65G50) Software, source code, etc. for problems pertaining to linear algebra (15-04) Special matrices (15B99)

This page was built for publication: Forward stable computation of roots of real polynomials with only real distinct roots

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