Algorithms and Codes for the Macdonald Function: Recent Progress and Comparisons

Published

February 1, 2003

Author(s)

B R. Fabijonas, Daniel W. Lozier, J M. Rappoport

Abstract

The modified Bessel function K_iŅχ, also known as the Macdonald function, finds application in the Kontorovich-Lebedev integral transform when χ and Ņ are real and positive. In this paper, a comparison of three codes for computing this function is made. These codes differ in algorithmic approach, timing, and regions of validity. One of them can be tested independently of the other two through Wronskian checks, and therefore is used as a standard against which the others are compared.

Citation

NIST Interagency/Internal Report (NISTIR) - 6596

Report Number

6596

Volume

161

NIST Pub Series

NIST Interagency/Internal Report (NISTIR)

Pub Type

NIST Pubs

Keywords

Bessel functions, computation of special functions, Macdonald function

Citation

Fabijonas, B. , Lozier, D. and Rappoport, J. (2003), Algorithms and Codes for the Macdonald Function: Recent Progress and Comparisons, NIST Interagency/Internal Report (NISTIR), National Institute of Standards and Technology, Gaithersburg, MD (Accessed October 16, 2025)