Share Email Print

Proceedings Paper

Multiple precision computation of interval elementary functions
Author(s): Fabien Rico
Format Member Price Non-Member Price
PDF $17.00 $21.00

Paper Abstract

The control of error during numerical evaluation of an expression is a crucial problem. A natural solution is to use interval arithmetic during calculations. For this reason, we need algorithms able to efficiently determine the image of an interval by the elementary functions. Adapting classical algorithms to interval arithmetic leads to a large overestimation of the result. In this paper, we propose a method based on Taylor approximation which evaluates functions on the bounds of the interval and deduces the resulting image. The computation of each bound in done with the Taylor approximation evaluated by the Smith scheme. Indeed, polynomial approximation seems to be the most adapted to our range of precision (several hundreds of digits). This approach is better than the classical algorithm using the Horner scheme and than the other polynomial approximations (e.g. minimax and Chebyshev polynomials) which are more precise but also more complicated. Furthermore, we present a rigorous evaluation of the error. Thus, the desired accuracy of the result is given, even on the numerically instable points (big number, (pi) /2 ...). For the computation of a point, our algorithm implemented with GMP is up to 2 times faster than the corresponding computation with softwares like MAPLE or MUPAD which don't even guarantee the relative precision. We show that the control of the error is not costly. And since the Taylor approximation is just a bit less precise but far more simple than the other approximations, it can be rapidly determined and, all in all, is the most efficient.

Paper Details

Date Published: 2 November 1999
PDF: 12 pages
Proc. SPIE 3807, Advanced Signal Processing Algorithms, Architectures, and Implementations IX, (2 November 1999); doi: 10.1117/12.367648
Show Author Affiliations
Fabien Rico, Univ. de Provence (France)

Published in SPIE Proceedings Vol. 3807:
Advanced Signal Processing Algorithms, Architectures, and Implementations IX
Franklin T. Luk, Editor(s)

© SPIE. Terms of Use
Back to Top
Sign in to read the full article
Create a free SPIE account to get access to
premium articles and original research
Forgot your username?