Paper Info
Title
On the Cohen--Olivier algorithm for computing ζ (s): Error analysis in the real case for an arbitrary precision
Abstract
The only published error analysis for an approximation algorithm computing the Riemann zeta-function @z(s), due to Henri Cohen and Michel Olivier, does evaluate the error of the approximation, but is not concerned by the fact that the computations required to calculate this approximation will be carried on with a finite precision arithmetic (by a computer), and thus produce other (rounding) errors. As a first step towards clearing this matter we provide a complete error analysis of the Cohen-Olivier algorithm when s is real with s=1/2, s1. We prove that, if s can be written with D"s bits in base 2, then in order to compute @z(s) in any relative precision P=11, that is, in order to compute a P-bit number @z"P(s) such that |@z"P(s)-@z(s)| is certified to be smaller than the number represented by a ''1'' at the Pth and last significant bit-place of |@z"P(s)|, it is sufficient to perform all the computations (i.e. additions, subtractions, multiplications, divisions, and computation of k^-^s for integers k=2) with an internal precisionD=max(D"s,P+max(14,@?3logP2log2+2.71@?)), and then to round to the nearest P-bits number. For instance if the wanted precision is P=1000 (and if s has no more than 1018 significant bits), then an internal precision D=1018 is sufficient.
Year
DOI
Venue
2007
10.1016/j.aam.2005.10.003
Advances in Applied Mathematics
Keywords
DocType
Volume
relative precision,complete error analysis,P-bit number,real case,error analysis,approximation algorithm,arbitrary precision,finite precision arithmetic,wanted precision,internal precision,Cohen-Olivier algorithm,nearest P-bits number
Journal
38
Issue
ISSN
Citations 
1
0196-8858
0
PageRank 
References 
Authors
0.34
0
2
Name
Order
Citations
PageRank
Y. -F. S. Pétermann100.68
Jean-Luc Rémy241.16