34328 articles – 26891 Notices  [english version]

ensl-00734339, version 1

On the componentwise accuracy of complex floating-point division with an FMA

Claude-Pierre Jeannerod () 12, Nicolas Louvet () 12, Jean-Michel Muller (, http://perso.ens-lyon.fr/jean-michel.muller/) 12

21st IEEE Symposium on Computer Arithmetic (2013) 8

Résumé : This paper deals with the accuracy of complex division in radix-two floating-point arithmetic. Assuming that a fused multiply-add (FMA) instruction is available and that no underflow/overflow occurs, we study how to ensure high relative accuracy in the componentwise sense. Since this essentially reduces to evaluating accurately three expressions of the form ac+bd, an obvious approach would be to perform three calls to Kahan's compensated algorithm for 2 by 2 determinants. However, in the context of complex division, two of those expressions are such that ac and bd have the same sign, suggesting that cheaper schemes should be used here (since cancellation cannot occur). We first give a detailed accuracy analysis of such schemes for the sum of two nonnegative products, providing not only sharp bounds on both their absolute and relative errors, but also sufficient conditions for the output of one of them to coincide with the output of Kahan's algorithm. By combining Kahan's algorithm with this particular scheme, we then deduce two new division algorithms. Our first algorithm is a straight-line program whose componentwise relative error is always at most 5u+13u^2 with u the unit roundoff; we also provide examples of inputs for which the error of this algorithm approaches 5u, thus showing that our upper bound is essentially the best possible. When tests are allowed we show with a second algorithm that the bound above can be further reduced to (9/2) u+9u^2, and that this improved bound is reasonably sharp.

  • 1 :  ARIC (Inria Grenoble Rhône-Alpes / LIP Laboratoire de l'Informatique du Parallélisme)
  • INRIA – CNRS : UMR5668 – Université Claude Bernard - Lyon I (UCBL) – École Normale Supérieure (ENS) - Lyon
  • 2 :  Laboratoire de l'Informatique du Parallélisme (LIP)
  • PRES Université de Lyon – CNRS : UMR5668 – INRIA – École Normale Supérieure (ENS) - Lyon – Université Claude Bernard - Lyon I (UCBL)
  • Domaine : Informatique/Autre
  • Mots-clés : Floating-point arithmetic – complex division – fused multiply-add (FMA) – rounding error analysis
  • Versions disponibles :  v1 (21-09-2012) v2 (07-02-2013)
  • ensl-00734339, version 1
  • oai:hal-ens-lyon.archives-ouvertes.fr:ensl-00734339
  • Contributeur : 
  • Soumis le : Vendredi 21 Septembre 2012, 14:43:51
  • Dernière modification le : Lundi 10 Décembre 2012, 11:02:53