The Laurent–Horner method for validated evaluation of Chebyshev expansions

Abstract

We develop a simple two-step algorithm for enclosing Chebyshev expansions whose cost is linear in terms of the polynomial degree. The algorithm first transforms the expansion from Chebyshev to the Laurent basis and then applies the interval Horner method. It outperforms the existing eigenvalue-based methods if the degree is high or the evaluation point is close to the boundaries of the domain.