On the relative error of computing complex square roots in floating-point arithmetic

Claude-Pierre Jeannerod (1) , Jean-Michel Muller (2, 1)

1 ARIC - Arithmetic and Computing
2 CNRS - Centre National de la Recherche Scientifique

We study the accuracy of a classical approach to computing complex square-roots in floating-point arithmetic. Our analyses are done in binary floating-point arithmetic in precision p, and we assume that the (real) arithmetic operations +, −, ×, ÷, √ are rounded to nearest, so the unit roundoff is u = 2^−p. We show that in the absence of underflow and overflow, the componentwise and normwise relative errors of this approach are at most 7 / 2 u and u √ 37/2, respectively, and this without having to neglect terms of higher order in u. We then provide some input examples showing that these bounds are reasonably sharp for the three basic binary interchange formats (binary32, binary64, and binary128) of the IEEE 754 standard for floating-point arithmetic.

Domaines: Autre [cs.OH]; Arithmétique des ordinateurs

Claude-Pierre Jeannerod, Jean-Michel Muller. On the relative error of computing complex square roots in floating-point arithmetic. ACSSC 2017 - 51st Asilomar Conference on Signals, Systems, and Computers, Oct 2017, Pacific Grove, United States. pp.737-740, ⟨10.1109/ACSSC.2017.8335442⟩. ⟨ensl-01780265⟩