Name
Title | ||
|---|---|---|
On the computation of the reciprocal of floating point expansions using an adapted Newton-Raphson iteration |
Abstract | ||
|---|---|---|
Many numerical problems require a higher computing precision than that offered by common floating point (FP) formats. One common way of extending the precision is to represent numbers in a multiple component format. With so-called floating point expansions, numbers are represented as the unevaluated sum of standard machine precision FP numbers. This format offers the simplicity of using directly available and highly optimized FP operations and is used by multiple-precisions libraries such as Bailey's QD or the analogue Graphics Processing Units tuned version, GQD. In this article we present a new algorithm for computing the reciprocal FP expansion a-1 of a FP expansion a. Our algorithm is based on an adapted Newton-Raphson iteration where we use “truncated” operations (additions, multiplications) involving FP expansions. The thorough error analysis given shows that our algorithm allows for computations of very accurate quotients. Precisely, after q ≤ 0 iterations, the computed FP expansion x = x0 + ... + x2q-1 satisfies the relative error bound |x-a-1/a-1|≤2-2q(p-3)-1, where p > 2 is the precision of the FP representation used (p = 24 for single precision and p = 53 for double precision). |
Year | DOI | Venue |
|---|---|---|
2014 | 10.1109/ASAP.2014.6868632 | Application-specific Systems, Architectures and Processors |
Keywords | Field | DocType |
Newton-Raphson method,floating point arithmetic,FP operations,FP representation precision,adapted Newton-Raphson iteration,double-precision,error analysis,floating point expansion reciprocal,multiple component format,multiple-precision libraries,numerical problems,relative error bound,single-precision,standard machine precision FP numbers,truncated addition operation,truncated multiplication operation,Newton-Raphson iteration,division,floating-point arithmetic,floating-point expansions,high precision arithmetic,multiple-precision arithmetic | Single-precision floating-point format,Discrete mathematics,Arbitrary-precision arithmetic,Computer science,Floating point,Parallel computing,Quotient,Double-precision floating-point format,Algorithm,Machine epsilon,Computation,Extended precision | Conference |
ISSN | Citations | PageRank |
2160-0511 | 2 | 0.43 |
References | Authors | |
4 | 3 | |
Name | Order | Citations | PageRank |
|---|---|---|---|
| Mioara Joldeş | 1 | 110 | 11.53 |
| Jean-Michel Muller | 2 | 466 | 66.61 |
| Valentina Popescu | 3 | 2 | 0.43 |