• DocumentCode
    610868
  • Title

    Multiple-Precision Evaluation of the Airy Ai Function with Reduced Cancellation

  • Author

    Chevillard, S. ; Mezzarobba, M.

  • Author_Institution
    Apics Project-Team, Inria, Sophia Antipolis, France
  • fYear
    2013
  • fDate
    7-10 April 2013
  • Firstpage
    175
  • Lastpage
    182
  • 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 Rein hard, 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.
  • Keywords
    differential equations; series (mathematics); Airy Ai function; Taylor coefficients; cancellation reduction; classical Miller algorithm; ill-conditioned three-term recurrence; linear ordinary differential equation; multiple-precision evaluation; nonnegative Taylor expansions; series expansion; Accuracy; Algorithm design and analysis; Approximation algorithms; Approximation methods; Equations; Shape; Taylor series; Miller method; Special functions; algorithm; arbitrary precision; asymptotics; correct rounding; error bounds; numerical evaluation;
  • fLanguage
    English
  • Publisher
    ieee
  • Conference_Titel
    Computer Arithmetic (ARITH), 2013 21st IEEE Symposium on
  • Conference_Location
    Austin, TX
  • ISSN
    1063-6889
  • Print_ISBN
    978-1-4673-5644-2
  • Type

    conf

  • DOI
    10.1109/ARITH.2013.33
  • Filename
    6545905