On a generalization of the Bessel function Neumann expansion
This is an incremental theoretical contribution for researchers working on series expansions and iterative methods.
The authors propose a generalization of the Bessel-Neumann expansion using basis functions that satisfy an infinite-dimensional linear ODE associated with a Hessenberg matrix, and provide a procedure to compute the basis functions and coefficients. They show that non-standard basis functions can achieve faster convergence than standard Bessel functions.
The Bessel-Neumann expansion (of integer order) of a function $g:\mathbb{C}\rightarrow\mathbb{C}$ corresponds to representing $g$ as a linear combination of basis functions $ϕ_0,ϕ_1,\ldots$, i.e., $g(z)=\sum_{\ell = 0}^\infty w_\ell ϕ_\ell(s)$, where $ϕ_i(z)=J_i(z)$, $i=0,\ldots$, are the Bessel functions. In this work, we study an expansion for a more general class of basis functions. More precisely, we assume that the basis functions satisfy an infinite dimensional linear ordinary differential equation associated with a Hessenberg matrix, motivated by the fact that these basis functions occur in certain iterative methods. A procedure to compute the basis functions as well as the coefficients is proposed. Theoretical properties of the expansion are studied. We illustrate that non-standard basis functions can give faster convergence than the Bessel functions.