On the componentwise accuracy of complex floating-point division with an FMA - Archive ouverte HAL Access content directly
Conference Papers Year : 2013

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

(1, 2) , (1, 2) , (1, 2)
1
2

Abstract

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.
Fichier principal
Vignette du fichier
Jeannerod_Louvet_Muller_final.pdf (160.33 Ko) Télécharger le fichier
Origin : Files produced by the author(s)
Loading...

Dates and versions

ensl-00734339 , version 1 (21-09-2012)
ensl-00734339 , version 2 (07-02-2013)

Identifiers

  • HAL Id : ensl-00734339 , version 2

Cite

Claude-Pierre Jeannerod, Nicolas Louvet, Jean-Michel Muller. On the componentwise accuracy of complex floating-point division with an FMA. 21st IEEE Symposium on Computer Arithmetic, Apr 2013, Austin, TX, United States. pp.83--90. ⟨ensl-00734339v2⟩
532 View
750 Download

Share

Gmail Facebook Twitter LinkedIn More