On the computation of the reciprocal of floating point expansions using an adapted Newton-Raphson iteration

Mioara Joldes 1, * Jean-Michel Muller 2, * Valentina Popescu 1, *
* Auteur correspondant
1 LAAS-MAC - Équipe Méthodes et Algorithmes en Commande
LAAS - Laboratoire d'analyse et d'architecture des systèmes [Toulouse]
2 ARIC - Arithmetic and Computing
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
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 \emph{multiple component} format. With so-called \emph{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 $i\geq0$ iterations, the computed FP expansion $x=x_0+\ldots+x_{2^i-1}$ satisfies the relative error bound: $\abs{\frac{x-a^{-1}}{a^{-1}}}\leq 2^{-2^i(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).
Type de document :
Communication dans un congrès
25th IEEE International Conference on Application-specific Systems, Architectures and Processors, ASAP, Jun 2014, Zurich, Switzerland. pp.8, 2014
Liste complète des métadonnées


https://hal.archives-ouvertes.fr/hal-00957379
Contributeur : Mioara Joldes <>
Soumis le : lundi 10 mars 2014 - 12:07:09
Dernière modification le : mardi 10 janvier 2017 - 15:11:31
Document(s) archivé(s) le : mardi 10 juin 2014 - 11:55:13

Fichier

invNewton.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

  • HAL Id : hal-00957379, version 1

Citation

Mioara Joldes, Jean-Michel Muller, Valentina Popescu. On the computation of the reciprocal of floating point expansions using an adapted Newton-Raphson iteration. 25th IEEE International Conference on Application-specific Systems, Architectures and Processors, ASAP, Jun 2014, Zurich, Switzerland. pp.8, 2014. <hal-00957379>

Partager

Métriques

Consultations de
la notice

753

Téléchargements du document

579