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.

A new polynomial approximation for Jν Bessel functions.

A new polynomial approximation for Jv Bessel functions

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

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

K -Functional, Weighted Moduli of Smoothness, and Best Weighted Polynomial Approximation on a Simplex

Approximate linear relations for Bessel functions

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

On the zeros of derivatives of Bessel functions

A note on trigonometric approximations of Bessel functions of the first kind and trigonometric power sums

Error bounds for a uniform asymptotic approximation of the zeros of the Bessel function $J_ν(x)$

Approximations of the modified Bessel functions of the second kind $K_ν$. Applications in random field generation

Approximation by the modified $ \lambda $-Bernstein-polynomial in terms of basis function

Integral Representations for the Lagrange Polynomials, Shively’s Pseudo-Laguerre Polynomials, and the Generalized Bessel Polynomials

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

Numerical Approximations to the Nonlinear Fractional-Order Logistic Population Model with Fractional-Order Bessel and Legendre Bases

The Solution of a New Kind of Integrals of Bessel Functions Used in the Power Calculation

A new product formula involving Bessel functions

An Introductory Overview of Bessel Polynomials, the Generalized Bessel Polynomials and the Q-Bessel Polynomials

New approximations to the principal real-valued branch of the Lambert W-function

Approximation of Functions of Several Variables by Multidimensional A- and J-fractions with Independent Variables

An Approximate Solution for a Mixed Linear Volterra–Fredholm Integral Equation