More Instruction Level Parallelism Explains the Actual Efficiency of Compensated Algorithms

Philippe Langlois (1) , Nicolas Louvet (1)

1 LP2A - Laboratoire de Physique Appliquée et d'Automatique

The compensated Horner algorithm and the Horner algorithm with double-double arithmetic improve the accuracy of polynomial evaluation in IEEE-754 floating point arithmetic. Both yield a polynomial evaluation as accurate as if it was computed with the classic Horner algorithm in twice the working precision. Both algorithms also share the same low-level computation of the floating point rounding errors and cost a similar number of floating point operations. We report numerical experiments to exhibit that the compensated algorithm runs at least twice as fast as the double-double one on modern processors. We propose to explain such efficiency by identifying more instruction level parallelism in the compensated implementation. Such property also applies to other compensated algorithms for summation, dot product and triangular linear system solving. More generally this paper illustrates how this kind of performance analysis may be useful to highlight the actual efficiency of numerical algorithms.

Mots clés: Accurate polynomial evaluation; Horner algorithm; compensated Horner algorithm; floating point arithmetic; IEEE-754 standard; instruction level parallelism; performance evaluation.; performance evaluation

Domaines: Logiciel mathématique [cs.MS]

Philippe Langlois, Nicolas Louvet. More Instruction Level Parallelism Explains the Actual Efficiency of Compensated Algorithms. 2007. ⟨hal-00165020⟩