Affiliation:
1. Computer Centre, Australian National University, Box 4, Canberra, ACT 2600, Australia
Abstract
Let ƒ(
x
) be one of the usual elementary functions (exp, log, artan, sin, cosh, etc.), and let
M
(
n
) be the number of single-precision operations required to multiply
n
-bit integers. It is shown that ƒ(
x
) can be evaluated, with relative error
Ο
(2
-
n
), in
Ο
(
M
(
n
)log (
n
)) operations as
n
→ ∞, for any floating-point number
x
(with an
n
-bit fraction) in a suitable finite interval. From the Schönhage-Strassen bound on
M
(
n
), it follows that an
n
-bit approximation to ƒ(
x
) may be evaluated in
Ο
(
n
log
2
(
n
) log log(
n
)) operations. Special cases include the evaluation of constants such as π,
e
, and
e
π
. The algorithms depend on the theory of elliptic integrals, using the arithmetic-geometric mean iteration and ascending Landen transformations.
Publisher
Association for Computing Machinery (ACM)
Subject
Artificial Intelligence,Hardware and Architecture,Information Systems,Control and Systems Engineering,Software
Cited by
253 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献
1. Fuchsian holonomic sequences;Applicable Algebra in Engineering, Communication and Computing;2023-09-24
2. Computing error bounds for asymptotic expansions of regular P-recursive sequences;Mathematics of Computation;2023-08-16
3. Order bounds for C2-finite sequences;Proceedings of the 2023 International Symposium on Symbolic and Algebraic Computation;2023-07-24
4. On the Order of Power Series and the Sum of Square Roots Problem;Proceedings of the 2023 International Symposium on Symbolic and Algebraic Computation;2023-07-24
5. Unscented Kalman filter with stable embedding for simple, accurate, and computationally efficient state estimation of systems on manifolds in Euclidean space;International Journal of Robust and Nonlinear Control;2022-11-29