arrow
Return

Arithmetic Algorithms for Extended Precision Using Floating-Point Expansions

delete2016-04-01
delete18
delete
OA
AI
M
Mioara Joldeş *
O
Olivier Marty *
J
Jean‐Michel Muller *
V
Valentina Popescu *
DOI:10.1109/TC.2015.2441714delete
deleteOriginal
deleteShare
deleteSave
View PDF
Abstract

Abstract

En 中文
Many numerical problems require a higher computing precision than the one offered by standard floating-point (FP) formats. One common way of extending the precision is to represent numbers in a multiple component format. By using the so-called floating-point expansions, real numbers are represented as the unevaluated sum of standard machine precision FP numbers. This representation offers the simplicity of using directly available, hardware implemented and highly optimized, FP operations. It is used by multiple-precision libraries such as Bailey's QD or the analogue Graphics Processing Units (GPU) tuned version, GQD. In this article we briefly revisit algorithms for adding and multiplying FP expansions, then we introduce and prove new algorithms for normalizing, dividing and square rooting of FP expansions. The new method used for computing the reciprocal a(-1) and the square root root a of a FP expansion a is based on an adapted Newton-Raphson iteration where the intermediate calculations are done using truncated operations (additions, multiplications) involving FP expansions. We give here a thorough error analysis showing that it allows very accurate computations. More precisely, after q iterations, the computed FP expansion x = x(0) + ... + x(2q-1) satisfies, for the reciprocal algorithm, the relative error bound: vertical bar(x - a(-1))/a(-1)vertical bar <= 2(-2q(p-3)-1) and, respectively, for the square root one: vertical bar x - 1/root a vertical bar <= 2(-2q(p-3)-1)/root a, where p > 2 is the precision of the FP representation used (p = 24 for single precision and p = 53 for double precision).
Keywords:
Floating-point arithmetic
floating-point expansions
high precision arithmetic
multiple-precision arithmetic
division
reciprocal
square root
Newton-Raphson iteration
AI Summary

AI Summary

Key information extracted from the uploaded paper, including a brief overview, abstract, background, key highlights, visual analysis, and future outlook.

Journal

IEEE Transactions on Computers cover
IEEE Transactions on Computers
IF:
3.8
Papers:
5.3K
Citations:
9.8K

Organization

C
centre national de la recherche scientifique (cnrs)
Scholars:
24.5W
Papers: 18.2W
Citations: 279
E
ecole normale superieure de lyon (ens de lyon)
Scholars:
5.1K
Papers: 3.6K
Citations: 6