Abstract
One considers the problem of finding hard to round cases of a periodic function for large floating-point inputs, more precisely when the function cannot be efficiently approximated by a polynomial. This is one of the last few issues that prevents from guaranteeing an efficient computation of correctly rounded transcendentals for the whole IEEE-754 double precision format. The first non-naive algorithm for that problem is presented, with a heuristic complexity of O(2°.⁶⁷⁶p) for a precision of p bits. The efficiency of the algorithm is shown on the largest IEEE-754 double precision binade for the sine function, and some corresponding bad cases are given. We can hope that all the worst cases of the trigonometric functions in their whole domain will be found within a few years, a task that was considered out of reach until now.
Original language | English |
---|---|
Title of host publication | Proceedings |
Subtitle of host publication | 18th IEEE Symposium on Computer Arithmetic : ARITH 18, Montpellier, France, June 25-27, 2007 |
Editors | Peter Kornerup, Jean-Michel Muller |
Place of Publication | Los Alamitos, Calif. |
Publisher | Institute of Electrical and Electronics Engineers (IEEE) |
Pages | 133-140 |
Number of pages | 8 |
ISBN (Print) | 0769528546 |
DOIs | |
Publication status | Published - 2007 |
Event | Symposium on Computer Arithmetic (18th : 2007) - Montpellier, France Duration: 25 Jun 2007 → 27 Jun 2007 |
Conference
Conference | Symposium on Computer Arithmetic (18th : 2007) |
---|---|
City | Montpellier, France |
Period | 25/06/07 → 27/06/07 |