Service interruption on Monday 11 July from 12:30 to 13:00: all the sites of the CCSD (HAL, EpiSciences, SciencesConf, AureHAL) will be inaccessible (network hardware connection).
Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

Optimizing polynomials for floating-point implementation

Abstract : The floating-point implementation of a function on an interval often reduces to polynomial approximation, the polynomial being typically provided by Remez algorithm. However, the floating-point evaluation of a Remez polynomial sometimes leads to catastrophic cancellations. This happens when some of the polynomial coefficients are very small in magnitude with respects to others. In this case, it is better to force these coefficients to zero, which also reduces the operation count. This technique, classically used for odd or even functions, may be generalized to a much larger class of functions. An algorithm is presented that forces to zero the smaller coefficients of the initial polynomial thanks to a modified Remez algorithm targeting an incomplete monomial basis. One advantage of this technique is that it is purely numerical, the function being used as a numerical black box. This algorithm is implemented within a larger polynomial implementation tool that is demonstrated on a range of examples, resulting in polynomials with less coefficients than those obtained the usual way.
Complete list of metadata

Cited literature [19 references]  Display  Hide  Download
Contributor : Christoph Lauter Connect in order to contact the contributor
Submitted on : Tuesday, March 4, 2008 - 2:18:24 PM
Last modification on : Saturday, September 11, 2021 - 3:17:02 AM
Long-term archiving on: : Thursday, May 20, 2010 - 11:52:49 PM


Files produced by the author(s)


  • HAL Id : ensl-00260563, version 1
  • ARXIV : 0803.0439



Florent de Dinechin, Christoph Lauter. Optimizing polynomials for floating-point implementation. 2008. ⟨ensl-00260563⟩



Record views


Files downloads