Multiple precision computation of interval elementary functions

Fabien Rico

doi:10.1117/12.367648

2 November 1999 Multiple precision computation of interval elementary functions

Fabien Rico

Proceedings Volume 3807, Advanced Signal Processing Algorithms, Architectures, and Implementations IX; (1999) https://doi.org/10.1117/12.367648
Event: SPIE's International Symposium on Optical Science, Engineering, and Instrumentation, 1999, Denver, CO, United States

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.

Citation Download Citation

Fabien Rico "Multiple precision computation of interval elementary functions", Proc. SPIE 3807, Advanced Signal Processing Algorithms, Architectures, and Implementations IX, (2 November 1999); https://doi.org/10.1117/12.367648

ACCESS THE FULL ARTICLE

INSTITUTIONAL
Select your institution to access the SPIE Digital Library.

SELECT YOUR INSTITUTION

PERSONAL
Sign in with your SPIE account to access your personal subscriptions or to use specific features such as save to my library, sign up for alerts, save searches, etc.

PERSONAL SIGN IN

No SPIE Account? Create one

PURCHASE THIS CONTENT

SUBSCRIBE TO DIGITAL LIBRARY

50 downloads per 1-year subscription

Members: $195

Non-members: $335 ADD TO CART

25 downloads per 1 - year subscription

Members: $145

Non-members: $250 ADD TO CART

PURCHASE SINGLE ARTICLE

Includes PDF, HTML & Video, when available

Members: $17.00

Non-members: $21.00 ADD TO CART

PROCEEDINGS
12 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

RIGHTS & PERMISSIONS

Get copyright permission Get copyright permission on Copyright Marketplace

KEYWORDS

Calcium

Error control coding

Evolutionary algorithms

Radon

Aluminum

Computing systems

Lithium

Show All Keywords

Keywords/Phrases

Search In:

Publication Years