A bound for the error term in the Brent-McMillan algorithm (original) (raw)
Abstract:The Brent-McMillan algorithm B3 (1980), when implemented with binary splitting, is the fastest known algorithm for high-precision computation of Euler's constant. However, no rigorous error bound for the algorithm has ever been published. We provide such a bound and justify the empirical observations of Brent and McMillan. We also give bounds on the error in the asymptotic expansions of functions related to modified Bessel functions.
Submission history
From: Richard Brent [view email]
[v1] Fri, 29 Nov 2013 23:47:32 UTC (8 KB)