Error bounds on complex floating-point multiplication

Source abstract

Given floating-point arithmetic with t t -digit base- β \beta significands in which all arithmetic operations are performed as if calculated to infinite precision and rounded to a nearest representable value, we prove that the product of complex values z 0 z_0 and z 1 z_1 can be computed with maximum absolute error | z 0 ‖ z 1 | 1 2 β 1 − t 5 |z_0\|z_1| \frac {1}{2} \beta ^{1 - t} \sqrt {5} . In particular, this provides relative error bounds of 2 − 24 5 2^{-24} \sqrt {5} and 2 − 53 5 2^{-53} \sqrt {5} for IEEE 754 single and double precision arithmetic respectively, provided that overflow, underflow, and denormals do not occur. We also provide the numerical worst cases for IEEE 754 single and double precision arithmetic.