Published in

Springer Verlag, Numerical Algorithms

DOI: 10.1007/s11075-016-0115-x

Links

Tools

Export citation

Search in Google Scholar

Sharp error bounds for complex floating-point inversion

This paper was not found in any repository, but could be made available legally by the author.
This paper was not found in any repository, but could be made available legally by the author.

Full text: Unavailable

Green circle
Preprint: archiving allowed
Green circle
Postprint: archiving allowed
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

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).