Skip to Main content Skip to Navigation
New interface
Conference papers

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 metadata

Cited literature [16 references]  Display  Hide  Download
Contributor : Marc Mezzarobba Connect in order to contact the contributor
Submitted on : Sunday, April 28, 2013 - 12:16:22 PM
Last modification on : Tuesday, October 25, 2022 - 4:24:42 PM
Long-term archiving on: : Tuesday, April 4, 2017 - 1:30:19 AM



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⟩



Record views


Files downloads