Skip to Main content Skip to Navigation
Preprints, Working Papers, ...

Sharp error bounds for complex floating-point inversion

Abstract : We study the accuracy of the classic algorithm for inverting a complex number given by its real and imaginary parts as floating-point numbers. Our analyses are done in binary floating-point arithmetic with an unbounded exponent range in precision p, and we assume that the elementary arithmetic operations (+, −, ×, /) are rounded to nearest, so that the roundoff unit is u = 2 −p. We prove the componentwise relative error bound 3u for the complex inversion algorithm (assuming p 4), and we show that this bound is asymptotically optimal (as p → ∞) when p is even, and reasonably sharp when using one of the basic IEEE 754 binary formats with an odd precision (p = 53, 113). This componentwise bound obviously leads to the same bound 3u for the normwise relative error. However we prove that the significantly smaller bound 2.707131u holds (assuming p 24) for the normwise relative error, and we illustrate the sharpness of this bound using numerical examples for the basic IEEE 754 binary formats (p = 24, 53, 113).
Document type :
Preprints, Working Papers, ...
Complete list of metadata
Contributor : Jean-Michel Muller <>
Submitted on : Tuesday, September 8, 2015 - 11:28:26 AM
Last modification on : Friday, September 10, 2021 - 2:34:05 PM
Long-term archiving on: : Wednesday, December 9, 2015 - 10:54:20 AM


Files produced by the author(s)


  • HAL Id : ensl-01195625, version 1


Claude-Pierre Jeannerod, Nicolas Louvet, Jean-Michel Muller, Antoine Plet. Sharp error bounds for complex floating-point inversion. 2015. ⟨ensl-01195625v1⟩



Record views


Files downloads