skip to main content
research-article

Algorithm 877: A Subroutine Package for Cylindrical Functions of Complex Order and Nonnegative Argument

Published: 01 July 2008 Publication History

Abstract

The algorithm presented provides a package of subroutines for calculating the cylindrical functions Jν(x), Nν(x), Hν(1)(x), Hν(2)(x) where the order ν is complex and the real argument x is nonnegative. The algorithm is written in Fortran 95 and calculates the functions using single, double, or quadruple precision according to the value of a parameter defined in the algorithm. The methods of calculating the functions are based on a series expansion, Debye's asymptotic expansions, Olver's asymptotic expansions, and recurrence methods (Miller's algorithms). The relative errors of the functional values computed by this algorithm using double precision are less than 2.4×10 − 13 in the region 0 ≤ Re ν ≤ 64, 0 ≤ Im ν ≤ 63, 0.024 ≤ x ≤ 97.

Supplementary Material

ZIP File (877.zip)
Software for A Subroutine Package for Cylindrical Functions of Complex Order and Nonnegative Argument

References

[1]
Amos, D. E. 1983a. Computation of Bessel functions of complex argument. Tech. rep. SAND83--0086, Sandia National Laboratories.
[2]
Amos, D. E. 1983b. Computation of Bessel functions of complex argument and large order. Tech. rep. SAND83--0643, Sandia National Laboratories.
[3]
Char, B. W., Geddes, K. O., Gonnet, G. H., Leong, B. L., Monagan, M. B., and Watt, S. M. 1991. Maple V Library Reference Manual. Springer-Verlag, Berlin, Germany.
[4]
Fabijonas, B. R. 2004. Algorithm 838: Airy functions. ACM Trans. Math. Soft. 30, 4, 491--501.
[5]
Felsen, L. B. and Marcuvitz, N. 1994. Radiation and Scattering of Waves. IEEE Press, Los Alamitos, CA, 323--328.
[6]
Gil, A., Segura, J., and Temme, N. M. 2002a. Computing complex Airy functions by numerical quadrature. Numer. Algor. 30, 11--23.
[7]
Gil, A., Segura, J., and Temme, N. M. 2002b. Evaluation of the modified Bessel function of the third kind of imaginary orders. J. Comput. Phys. 175, 398--411.
[8]
Gil, A., Segura, J., and Temme, N. M. 2003. Computation of the modified Bessel function of the third kind of imaginary orders: uniform Airy-type asymptotic expansion. J. Comput. Appl. Math. 153, 225--234.
[9]
Gil, A., Segura, J., and Temme, N. M. 2004. Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments. ACM Trans. Math. Soft. 30, 2, 145--158.
[10]
Goldstein, M. and Thaler, R. M. 1959. Recurrence techniques for the calculation of Bessel functions. Math. Tables Aids Comp. 13, 102--108.
[11]
Kodama, M. 1995. Numerical calculation of the Bessel function of complex order using the recurrence method. IECE Trans. Fundamentals E78-A, 4, 506--516.
[12]
Kodama, M. 1996. New algorithm on the recurrence method for numerical calculation of Bessel function of complex order. IECE Trans. Fundamentals E79-A, 4, 621--623.
[13]
Kodama, M., Yamasato, M., and Yamashioro, S. 1995. Numerical calculation of the Neumann function N ν(x) of complex order ν. IECE Trans. Fundamentals E78-A, 6, 727--736.
[14]
Kuki, H. 1972. Complex gamma function with error control. Comm. ACM 15, 4, 262--267.
[15]
Lozier, D. W. and Olver, F. W. J. 2000. Numerical Evaluation of Special Functions. homepage: http://math.nist.gov/mcsd/Reports/2001/nesf/.
[16]
Metcalf, M. and Reid, J. 1999. Fortran 90/95 Explained 2nd Ed. Oxford University Press, Oxford, England, UK.
[17]
Olver, F. W. 1970. Bessel functions of integer order. In Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, M. Abramowitz and I. A. Stegun, Eds. Dover Publications, Inc., New York, NY, Chapter 9, 355--433.
[18]
Olver, F. W. J. 1975. Asymptotics and Special Functions. Academic Press, New York, NY.
[19]
Olver, F. W. J. and Sookne, D. J. 1972. Note on backward recurrence algorithms. Math. Comput. 26, 120, 941--947.
[20]
Rappoport, J. M. 2001. Canonical vector polynomials for the computation of complex order Bessel functions with the tau method. Comput. Math. Appl. 41, 399--406.
[21]
Rashid, M. A. and Kodama, M. 2000. Numerical calculation of cylindrical functions of complex order using Debye's asymptotic series. IECE Trans. Fundamentals E83-A, 12, 2664--2671.
[22]
Sensiper, S. 1957. Cylindrical radio waves. IRE Trans. Antennas Propag. AP-5, 1, 56--70.
[23]
Stoer, J. and Bulirsch, R. 1991. Introduction to Numerical Analysis 2nd Ed. Springer-Verlag, Berlin, Germany.
[24]
Temme, N. M. 1975. On the numerical evaluation of the modified Bessel function of the third kind. J. Comput. Phys. 19, 324--337.
[25]
Temme, N. M. 1976. On the numerical evaluation of the ordinary Bessel function of the second kind. J. Comput. Phys. 21, 343--350.
[26]
Thompson, I. J. and Barnett, A. R. 1985. COULCC: A continued-fraction algorithm for coulomb functions of complex order with complex arguments. Comput. Phys. Comm. 36, 363--372.
[27]
Thompson, I. J. and Barnett, A. R. 1986. Coulomb and Bessel functions of complex arguments and order. J. Comput. Phys. 64, 490--509.
[28]
Watson, G. N. 1944. A treatise on the theory of Bessel functions, Second ed. Cambridge University Press, Cambridge, UK.
[29]
Wolfram, S. 2003. The Mathematica Book 5th Ed. Wolfram Media, Champaign, IL.

Cited By

View all
  • (2023)An Accurate and Numerically Stable Formulation for Computing the Electromagnetic Fields in Uniform Bend Rectangular WaveguidesIEEE Transactions on Microwave Theory and Techniques10.1109/TMTT.2022.321983071:3(988-996)Online publication date: Mar-2023
  • (2017)Tartarus: A relativistic Green’s function quantum average atom codeHigh Energy Density Physics10.1016/j.hedp.2017.06.00224(33-38)Online publication date: Sep-2017
  • (2013)Coupling Analysis in Concentric Ring WaveguidesJournal of Lightwave Technology10.1109/JLT.2013.226363331:13(2140-2145)Online publication date: Jul-2013
  • Show More Cited By

Index Terms

  1. Algorithm 877: A Subroutine Package for Cylindrical Functions of Complex Order and Nonnegative Argument

    Recommendations

    Comments

    Information & Contributors

    Information

    Published In

    cover image ACM Transactions on Mathematical Software
    ACM Transactions on Mathematical Software  Volume 34, Issue 4
    July 2008
    152 pages
    ISSN:0098-3500
    EISSN:1557-7295
    DOI:10.1145/1377596
    Issue’s Table of Contents
    Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected]

    Publisher

    Association for Computing Machinery

    New York, NY, United States

    Publication History

    Published: 01 July 2008
    Accepted: 01 June 2007
    Revised: 01 June 2007
    Received: 01 March 2006
    Published in TOMS Volume 34, Issue 4

    Permissions

    Request permissions for this article.

    Check for updates

    Badges

    Author Tags

    1. Cylindrical functions
    2. complex order
    3. nonnegative argument
    4. numerical calculation

    Qualifiers

    • Research-article
    • Research
    • Refereed

    Contributors

    Other Metrics

    Bibliometrics & Citations

    Bibliometrics

    Article Metrics

    • Downloads (Last 12 months)4
    • Downloads (Last 6 weeks)0
    Reflects downloads up to 06 Oct 2024

    Other Metrics

    Citations

    Cited By

    View all
    • (2023)An Accurate and Numerically Stable Formulation for Computing the Electromagnetic Fields in Uniform Bend Rectangular WaveguidesIEEE Transactions on Microwave Theory and Techniques10.1109/TMTT.2022.321983071:3(988-996)Online publication date: Mar-2023
    • (2017)Tartarus: A relativistic Green’s function quantum average atom codeHigh Energy Density Physics10.1016/j.hedp.2017.06.00224(33-38)Online publication date: Sep-2017
    • (2013)Coupling Analysis in Concentric Ring WaveguidesJournal of Lightwave Technology10.1109/JLT.2013.226363331:13(2140-2145)Online publication date: Jul-2013
    • (2011)Algorithm 912ACM Transactions on Mathematical Software (TOMS)10.1145/1916461.191647137:4(1-25)Online publication date: 1-Feb-2011
    • (2011)Basic Methods for Computing Special FunctionsRecent Advances in Computational and Applied Mathematics10.1007/978-90-481-9981-5_4(67-121)Online publication date: 2011

    View Options

    Get Access

    Login options

    Full Access

    View options

    PDF

    View or Download as a PDF file.

    PDF

    eReader

    View online with eReader.

    eReader

    Media

    Figures

    Other

    Tables

    Share

    Share

    Share this Publication link

    Share on social media