Abstract:This paper deals with the problem of numerically computing the roots of polynomials $p_k(x)$, $k=1,2,\ldots$, of degree $n=2^k-1$ recursively defined by $p_1(x)=x+1$, $p_k(x)=xp_{k-1}(x)^2+1$. An algorithm based on the Ehrlich-Aberth simultaneous iterations complemented by the Fast Multi-pole Method and the fast search of near neighbors of a set of complex numbers is provided. The algorithm, which relies on a specific strategy of selecting initial approximations, costs $O(n\log n)$ arithmetic operations per step. A Fortran 95 implementation is given and numerical experiments are carried out. Experimentally, it turns out that the number of iterations needed to arrive at numerical convergence is $O(\log n)$. This allows us to compute the roots of $p_k(x)$ up to degree $n=2^{24}-1$ in about 16 minutes on a laptop with 16 GB RAM, and up to degree $n=2^{28}-1$ in about 69 minutes on a machine with 256 GB RAM. The case of degree $n=2^{30}-1$ would require higher memory and higher precision to separate the roots. With a suitable adaptation of FMM to the limit of 256 GB RAM and by performing the computation in extended precision (i.e. with 10-byte floating point representation) we were able to compute all the roots in about two weeks of CPU time for $n=2^{30}-1$. From the experimental analysis, explicit asymptotic expressions of the real roots of $p_k(x)$ and an explicit expression of $\min_{i\ne j}|\xi_i^{(k)}-\xi_j^{(k)}|$ for the roots $\xi_i^{(k)}$ of $p_k(x)$ are deduced. The approach is effectively applied to general classes of polynomials defined by a doubling recurrence.

The Newman algorithm for constructing polynomials with restricted coefficients and many real roots

Simultaneous Primitive Root Values Of Polynomials

Polynomials whose coefficients are generalized Leonardo numbers

Numerical computation of the roots of Mandelbrot polynomials: an experimental analysis

The number of real zeros of polynomials with constrained coefficients

Extensions of the Bloch-Pólya Theorem on the number of real zeros of polynomials (II)

Continued Fraction Expansion of Real Roots of Polynomial Systems

Real-rootedness of the type A minuscule polynomials

An Elementary Problem in Galois Theory about the Roots of Irreducible Polynomials

A rectangular additive convolution for polynomials

On random polynomials with an intermediate number of real roots

Counting integer polynomials with several roots of maximal modulus

A polynomial analogue of Jacobsthal function

A strong law of large numbers for real roots of random polynomials

In search of Newton-type inequalities

Primitive Normal Polynomials with Multiple Coefficients Prescribed: an Asymptotic Result

Universality theorems for zeros of random real polynomials with fixed coefficients

Non universality for the variance of the number of real roots of random trigonometric polynomials

The Furstenberg-Sárközy theorem for polynomials in one or more prime variables

Explicit Criterion to Determine the Number of Positive Roots of a Polynomial1)