The functions erf and erfc computed with arbitrary precision and explicit error bounds

Sylvain Chevillard 1, 2, 3, *
* Auteur correspondant
2 CARAMEL - Cryptology, Arithmetic: Hardware and Software
Inria Nancy - Grand Est, LORIA - ALGO - Department of Algorithms, Computation, Image and Geometry
Abstract : The error function erf is a special function. It is widely used in statistical computations for instance, where it is also known as the standard normal cumulative probability. The complementary error function is defined as erfc(x)=erf(x)-1. In this paper, the computation of erf(x) and erfc(x) in arbitrary precision is detailed: our algorithms take as input a target precision t' and deliver approximate values of erf(x) or erfc(x) with a relative error bounded by 2^(-t'). We study three different algorithms for evaluating erf and erfc. These algorithms are completely detailed. In particular, the determination of the order of truncation, the analysis of roundoff errors and the way of choosing the working precision are presented. The scheme used for implementing erf and erfc and the proofs are expressed in a general setting, so they can directly be reused for the implementation of other functions. We implemented the three algorithms and studied experimentally what is the best algorithm to use in function of the point x and the target precision t'.
Type de document :
Article dans une revue
Information and Computation, Elsevier, 2012, 216, pp.72 -- 95
Liste complète des métadonnées


https://hal-ens-lyon.archives-ouvertes.fr/ensl-00356709
Contributeur : Sylvain Chevillard <>
Soumis le : jeudi 27 mai 2010 - 14:46:00
Dernière modification le : mardi 13 décembre 2016 - 15:40:28
Document(s) archivé(s) le : lundi 22 octobre 2012 - 14:30:10

Identifiants

  • HAL Id : ensl-00356709, version 3

Citation

Sylvain Chevillard. The functions erf and erfc computed with arbitrary precision and explicit error bounds. Information and Computation, Elsevier, 2012, 216, pp.72 -- 95. <ensl-00356709v3>

Partager

Métriques

Consultations de
la notice

1470

Téléchargements du document

2371