Computation of the error functions erf & erfc in arbitrary precision with correct rounding

Abstract : In this paper, the computation of erf(x) in arbitrary precision is detailed. A feature of our implementation is correct rounding: the returned result is the exact result (as if it were computed with infinite precision) rounded according to the specified rounding mode. The algorithm that computes the correctly rounded value of erf(x) for any argument x is detailed in this paper. In particular, the choice of the approximation formula, the determination of the truncation rank and of the computing precision are presented. When the current truncation rank and computing precision do not suffice to determine the correctly rounded value, they are increased and the computation is done again. The optimal strategy to adapt them is given in this paper. Finally, timings on some experiments are given, and the implementation of the complementary error function erfc is then sketched.
Type de document :
[Research Report] RR-6465, INRIA. 2008
Liste complète des métadonnées

Littérature citée [16 références]  Voir  Masquer  Télécharger
Contributeur : Nathalie Revol <>
Soumis le : mercredi 3 août 2011 - 16:50:28
Dernière modification le : mardi 16 janvier 2018 - 15:34:30
Document(s) archivé(s) le : dimanche 4 décembre 2016 - 04:59:04


Fichiers produits par l'(les) auteur(s)


  • HAL Id : inria-00261360, version 3



Sylvain Chevillard, Nathalie Revol. Computation of the error functions erf & erfc in arbitrary precision with correct rounding. [Research Report] RR-6465, INRIA. 2008. 〈inria-00261360v3〉



Consultations de la notice


Téléchargements de fichiers