Abstract:We describe a method for the rapid numerical evaluation of the Bessel functions of the first and second kinds of nonnegative real orders and positive arguments. Our algorithm makes use of the well-known observation that although the Bessel functions themselves are expensive to represent via piecewise polynomial expansions, the logarithms of certain solutions of Bessel's equation are not. We exploit this observation by numerically precomputing the logarithms of carefully chosen Bessel functions and representing them with piecewise bivariate Chebyshev expansions. Our scheme is able to evaluate Bessel functions of orders between $0$ and $1\sep,000\sep,000\sep,000$ at essentially any positive real argument. In that regime, it is competitive with existing methods for the rapid evaluation of Bessel functions and has several advantages over them. First, our approach is quite general and can be readily applied to many other special functions which satisfy second order ordinary differential equations. Second, by calculating the logarithms of the Bessel functions rather than the Bessel functions themselves, we avoid many issues which arise from numerical overflow and underflow. Third, in the oscillatory regime, our algorithm calculates the values of a nonoscillatory phase function for Bessel's differential equation and its derivative. These quantities are useful for computing the zeros of Bessel functions, as well as for rapidly applying the Fourier-Bessel transform. The results of extensive numerical experiments demonstrating the efficacy of our algorithm are presented. A Fortran package which includes our code for evaluating the Bessel functions as well as our code for all of the numerical experiments described here is publically available.

[O III]/[N II] as an Abundance Indicator at High Redshift

An algorithm for the rapid numerical evaluation of Bessel functions of real orders and arguments

An Effective Approximation Algorithm for Second-Order Singular Functional Differential Equations

Numerical solutions of Korteweg de Vries and Korteweg de Vries-Burger's equations using computer programming

A New Polynomial Approximation for J&gt;inf&lt;&#957;&gt;/inf&lt;Bessel Functions

Numerical calculation of Bessel, Hankel and Airy functions

Bidimensional Gegenbauer Polynomials for Variable‐Order Time‐Fractional Integro‐Partial Differential Equation With a Weakly Singular Kernel

Non-Standard Finite Difference and Vieta-Lucas Orthogonal Polynomials for the Multi-Space Fractional-Order Coupled Korteweg-de Vries Equation

A Riemann–Hilbert Method to Algebro-Geometric Solutions of the Korteweg–de Vries Equation

Positive Solutions and Their Existence of a Nonlinear Hadamard Fractional-Order Differential Equation with a Singular Source Item Using Spectral Analysis

On a novel iterative method to compute polynomial approximations to Bessel functions of the first kind and its connection to the solution of fractional diffusion/diffusion-wave problems

Derivatives with respect to the degree and order of associated Legendre functions for $|z|>1$ using modified Bessel functions

A new polynomial approximation for Jv Bessel functions

On the computation of fundamental functions and Abelian differentials of the third kind

Some formulas for Bessel functions related to diag(1, −1, −1)-matrices and an intertwining operator

Numerical Solution of Higher-Order Linear and Nonlinear Ordinary Differential Equations with Orthogonal Rational Legendre Functions

Modified Bessel Functions in Analytic Number Theory

Infinite series representations for Bessel functions of the first kind of integer order

Generalized Bessel Quasilinearization Technique Applied to Bratu and Lane-Emden-Type Equations of Arbitrary Order

The Lamé functions and elliptic soliton solutions: Bilinear approach

Numerical Solutions of the Multi-Space Fractional-Order Coupled Korteweg–De Vries Equation with Several Different Kernels