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.

Approximate linear relations for Bessel functions

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

A new polynomial approximation for Jv Bessel functions

Almost Exact Boundary Condition for One-Dimensional Schrödinger Equations

A new polynomial approximation for Jν Bessel functions.

ALmost EXact Boundary Conditions for Transient Schrödinger–Poisson System

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

The Approximation of the Exact Boundary Conditions at an Artificial Boundary for Linear Elastic Equations and Its Application

Certain Classes Of Finite-Series Relationships And Generating Functions Involving The Generalized Bessel Polynomials

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

Generating Generalized Bessel Equations by Virtue of Bose Operator Algebra and Entangled State Representations

Asymptotic boundary conditions and numerical methods for nonlinear elliptic problems on unbounded domains

Stable and Accurate Least Squares Radial Basis Function Approximations on Bounded Domains

An analytical approximation scheme to two point boundary value problems of ordinary differential equations

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

Sufficient Conditions for Linear Operators Related to Confluent Hypergeometric Function and Generalized Bessel Function of the First Kind to Belong to a Certain Class of Analytic Functions

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

The variation of Barnes and Bessel zeta functions

An Approximate Solution of Nonlinear Wave Equation

Asymptotic Matching the Self-Consistent Expansion to Approximate the Modified Bessel Functions of the Second Kind

On the Inverse of the Linearization Coefficients of Bessel Polynomials