Multiple precision evaluation of the Airy Ai function with reduced cancellation

Sylvain Chevillard 1 Marc Mezzarobba 2, *
* Corresponding author
2 ARIC - Arithmetic and Computing
Inria Grenoble - Rhône-Alpes, LIP - Laboratoire de l'Informatique du Parallélisme
Abstract : The series expansion at the origin of the Airy function Ai(x) is alternating and hence problematic to evaluate for x > 0 due to cancellation. Based on a method recently proposed by Gawronski, Müller, and Reinhard, we exhibit two functions F and G, both with nonnegative Taylor expansions at the origin, such that Ai(x) = G(x)/F(x). The sums are now well-conditioned, but the Taylor coefficients of G turn out to obey an ill-conditioned three-term recurrence. We use the classical Miller algorithm to overcome this issue. We bound all errors and our implementation allows an arbitrary and certified accuracy, that can be used, e.g., for providing correct rounding in arbitrary precision.
Complete list of metadatas

Cited literature [16 references]  Display  Hide  Download

https://hal.inria.fr/hal-00767085
Contributor : Marc Mezzarobba <>
Submitted on : Sunday, April 28, 2013 - 12:16:22 PM
Last modification on : Thursday, February 7, 2019 - 4:33:34 PM
Long-term archiving on : Tuesday, April 4, 2017 - 1:30:19 AM

Identifiers

Collections

Citation

Sylvain Chevillard, Marc Mezzarobba. Multiple precision evaluation of the Airy Ai function with reduced cancellation. 21st IEEE Symposium on Computer Arithmetic, 2013, Austin, TX, United States. pp.175-182, ⟨10.1109/ARITH.2013.33⟩. ⟨hal-00767085v2⟩

Share

Metrics

Record views

705

Files downloads

851