Numerical Quadrature for Bessel Transformations with High Oscillations

Shuhuang Xiang

DOI: https://doi.org/10.1007/978-3-642-00464-3_69

2009-01-01

Abstract:We explore higher order numerical quadrature for the integration of systems containing Bessel functions such as $\int_a^b f(x)J_{\nu}(rx)dx$ and $\int_a^b f(x)\cos(r_1x)J_{\nu}(rx)dx$. The decay of the error of the these methods drastically improves as frequency grows.

What problem does this paper attempt to address?

Computational-Complexity on Numerical Integrals

XH WANG,DF HAN

1991-01-01

Abstract:This paper studies the computational complexity on numerical integrals initiated by S. Smale. For a general quadrature functional with k - 1 algebraic precision, the average error of a complexification quadrature formula in the Sobolev space H(k) is obtained. The approach applied in this paper may be also available to other linear approximate problems.
The numerical integration of dual hankel transformation

Jun HUA,Yan-sheng JIANG,Wen-bing WANG

DOI: https://doi.org/10.3969/j.issn.1001-1986.2001.03.019

2001-01-01

Abstract:Due to strong oscillation and slow decay of product term in dual Bessel integration,it is difficult to use ordinary numerical method for the quadrature of dual Bessel.In this paper we divide ［0,∞) into ［0,λ0］ and ［λ0,∞);In ［λ0,∞),the integration can be written as Fourier cosine and sine transformation by asymptotic expression of Hankel Function and it can be evaluated by fast algorthm.In ［0,λ0］,the integration can be computed accurately through direct quadrature.But when we calculate a quantity the related dual Bessel integration with common argument,the algorithm mentioned above seem not to be efficient enough.In this case,according to derivative relation of Bessel function,the constant transformation of the integration in ［0,λ0］ can be done and the changed into the integration of the Bessel function itself which can only be calculated only once.This algorithm has improved efficiency of calculation obviously for dual Bessel integration.
On quadrature of highly oscillatory Bessel function via asymptotic analysis of simplex integrals

Yongxiong Zhou,Ruyun Chen

DOI: https://doi.org/10.1016/j.cam.2024.116239

IF: 2.872

2024-08-31

Journal of Computational and Applied Mathematics

Abstract:In this article, two methods for evaluating highly oscillatory Bessel integrals are explored. Firstly, a polynomial is analyzed as an effective approximation of the simplex integral of a highly oscillatory Bessel function based on Laplace transform, and its error rapidly decreases as the frequency increases. Furthermore, the inner product of f and highly oscillatory Bessel function can be approximated by two other forms of inner product by which one depends on a polynomial and the higher derivatives of f , another depends on Bessel function and the interpolation polynomial of f . In addition, three issues related to highly oscillatory Bessel integrals have also been discussed: inequalities for the convergence rate of Filon-type methods, evaluation of Cauchy principal values, and simplified evaluation on infinite intervals. Through some preliminary numerical experiments, our theoretical analysis has been preliminarily confirmed, and the proposed numerical method is accurate and effective.

mathematics, applied
Quadrature Methods for Highly Oscillatory Singular Integrals

Jing Gao,Marissa Condon,Arieh Iserles,Benjamin Gilvey,Jon Trevelyan

DOI: https://doi.org/10.4208/jcm.1911-m2019-0044

2020-01-01

Journal of Computational Mathematics

Abstract:We address the evaluation of highly oscillatory integrals, with power-law and logarithmic singularities. Such problems arise in numerical methods in engineering. Notably, the evaluation of oscillatory integrals dominates the run-time for wave-enriched boundary integral formulations for wave scattering, and many of these exhibit singularities. We show that the asymptotic behaviour of the integral depends on the integrand and its derivatives at the singular point of the integrand, the stationary points and the endpoints of the integral. A truncated asymptotic expansion achieves an error that decays faster for increasing frequency. Based on the asymptotic analysis, a Filon-type method is constructed to approximate the integral. Unlike an asymptotic expansion, the Filon method achieves high accuracy for both small and large frequency. Complex-valued quadrature involves interpolation at the zeros of polynomials orthogonal to a complex weight function. Numerical results indicate that the complex-valued Gaussian quadrature achieves the highest accuracy when the three methods are compared. However, while it achieves higher accuracy for the same number of function evaluations, it requires significant additional cost of computation of orthogonal polynomials and their zeros.
An Efficient Quadrature Rule for the Oscillatory Infinite Generalized Bessel Transform with a General Oscillator and Its Error Analysis

Hongchao Kang,Hong Wang

DOI: https://doi.org/10.1007/s10915-022-02081-6

2022-12-31

Journal of Scientific Computing

Abstract:Our recent work (Kang and Wang in J Sci Comput 82:1–33, 2020) performed a complete asymptotic analysis and proposed a modified Filon-type method for a class of oscillatory infinite Bessel transform with a general oscillator. In this paper, we present and analyze a different method by converting the integration path to the complex plane for this class of oscillatory infinite Bessel transform. In particular, we establish a series of new quadrature rules for this transform and carry out rigorous analysis, including the cases that the oscillator g ( x ) has either zeros or stationary points. The error analysis shows the advantages that this approach exhibits high asymptotic order, and the accuracy improves significantly as either the frequency or the number of nodes n increases. Furthermore, the constructed method shows higher accuracy and error order by comparing with the existing modified Filon-type method in our recent work (Kang and Wang 2020) at the same computational cost. Some numerical experiments are provided to verify the theoretical results and demonstrate the efficiency of the proposed method.

mathematics, applied
On the Convolution Quadrature Rule for Integral Transforms with Oscillatory Bessel Kernels

Junjie Ma,Huilan Liu

DOI: https://doi.org/10.3390/sym10070239

2018-06-25

Symmetry

Abstract:Lubich’s convolution quadrature rule provides efficient approximations to integrals with special kernels. Particularly, when it is applied to computing highly oscillatory integrals, numerical tests show it does not suffer from fast oscillation. This paper is devoted to studying the convergence property of the convolution quadrature rule for highly oscillatory problems. With the help of operational calculus, the convergence rate of the convolution quadrature rule with respect to the frequency is derived. Furthermore, its application to highly oscillatory integral equations is also investigated. Numerical results are presented to verify the effectiveness of the convolution quadrature rule in solving highly oscillatory problems. It is found from theoretical and numerical results that the convolution quadrature rule for solving highly oscillatory problems is efficient and high-potential.

multidisciplinary sciences
Analysis of Numerical Integration Error for Bessel Integral Identity in Fast Multipole Method for 2D Helmholtz Equation

Hai-jun Wu,Wei-kang Jiang,Yi-jun Liu

DOI: https://doi.org/10.1007/s12204-010-1070-7

2010-01-01

Abstract:In 2D fast multipole method for scattering problems, square quadrature rule is used to discretize the Bessel integral identity for diagonal expansion of 2D Helmholtz kernel, and numerical integration error is introduced. Taking advantage of the relationship between Euler-Maclaurin formula and trapezoidal quadrature rule, and the relationship between trapezoidal and square quadrature rule, sharp computable bound with analytical form on the error of numerical integration of Bessel integral identity by square quadrature rule is derived in this paper. Numerical experiments are presented at the end to demonstrate the accuracy of the sharp computable bound on the numerical integration error.
An Efficient Quadrature Rule for Highly Oscillatory Integrals with Airy Function

Guidong Liu,Zhenhua Xu,Bin Li

DOI: https://doi.org/10.3390/math12030377

IF: 2.4

2024-01-24

Mathematics

Abstract:In this work, our primary focus is on the numerical computation of highly oscillatory integrals involving the Airy function. Specifically, we address integrals of the form ∫0bxαf(x)Ai(−ωx)dx over a finite or semi-infinite interval, where the integrand exhibits rapid oscillations when ω≫1. The inherent high oscillation and algebraic singularity of the integrand make traditional quadrature rules impractical. In view of this, we strategically partition the interval into two segments: [0,1] and [1,b]. For integrals over the interval [0,1], we introduce a Filon-type method based on a two-point Taylor expansion. In contrast, for integrals over [1,b], we transform the Airy function into the first kind of Bessel function. By applying Cauchy's integration theorem, the integral is then reformulated into several non-oscillatory and exponentially decaying integrals over [0,+∞), which can be accurately approximated by the generalized Gaussian quadrature rule. The proposed methods are accompanied by rigorous error analyses to establish their reliability. Finally, we present a series of numerical examples that not only validate the theoretical results but also showcase the accuracy and efficacy of the proposed method.

mathematics
On the Numerical Evaluation of a Class of Oscillatory Integrals in Worldline Variational Calculations

R. Rosenfelder

DOI: https://doi.org/10.48550/arXiv.hep-ph/0603161

2006-03-21

High Energy Physics - Phenomenology

Abstract:Filon-Simpson quadrature rules are derived for integrals of the type \int_a^b dx f(x) sin(xy)/(xy) and \int_a^b dx f(x) 4 sin^2(xy/2)/(xy)^2 which are needed in applications of the worldline variational approach to Quantum Field Theory. These new integration rules reduce to the standard Simpson rule for y = 0 and are exact for y \to \infty when a = 0 and f(0) \ne 0.The subleading term in the asymptotic expansion is also reproduced more and more precisely when the number of integration points is increased. Tests show that the numerical results are indeed stable over a wide range of y-values whereas usual Gauss-Legendre quadrature rules are more precise at low y but fail completely for large values of y. The associated Filon-Simpson weights are given in terms of sine and cosine integrals and have to be evaluated for each value of y. A Fortran program to calculate them in a fast and accurate manner is available. A detailed comparison with the double exponential method of Ooura and Mori is made.
High-Order Numerical Integration on Domains Bounded by Intersecting Level Sets

Lauritz Beck,Florian Kummer

2023-08-21

Abstract:We present a high-order method that provides numerical integration on volumes, surfaces, and lines defined implicitly by two smooth intersecting level sets. To approximate the integrals, the method maps quadrature rules defined on hypercubes to the curved domains of the integrals. This enables the numerical integration of a wide range of integrands since integration on hypercubes is a well known problem. The mappings are constructed by treating the isocontours of the level sets as graphs of height functions. Numerical experiments with smooth integrands indicate a high-order of convergence for transformed Gauss quadrature rules on domains defined by polynomial, rational, and trigonometric level sets. We show that the approach we have used can be combined readily with adaptive quadrature methods. Moreover, we apply the approach to numerically integrate on difficult geometries without requiring a low-order fallback method.

Numerical Analysis
Computation of highly oscillatory integrals using a Fourier extension approximation

Akash Anand,Damini Dhiman

2024-08-30

Abstract:The numerical evaluation of integrals of the form \begin{align*} \int_a^b f(x) e^{ikg(x)}\,dx \end{align*} is an important problem in scientific computing with significant applications in many branches of applied mathematics, science and engineering. The numerical approximation of such integrals using classical quadratures can be prohibitively expensive at high oscillation frequency ($k \gg 1$) as the number of quadrature points needed for achieving a reasonable accuracy must grow proportionally to $k$. To address this significant computational challenge, starting with Filon in 1930, several specialized quadratures have been developed to compute such oscillatory integrals efficiently. A crucial element in such Filon-type quadrature is the accurate evaluation of certain moments which poses a significant challenge when non-linear phase functions $g$ are involved. In this paper, we propose an equispaced-grid Filon-type quadrature for computing such highly oscillatory integrals that utilizes a Fourier extension of the slowly varying envelope $f$. This strategy is primarily aimed at significantly simplifying the moment calculations, even when the phase function has stationary points. Moreover, the proposed approach can also handle certain integrable singularities in the integrand. We analyze the scheme to theoretically establish high-order convergence rates. We also include a wide variety of numerical experiments, including oscillatory integrals with algebraic and logarithmic singularities, to demonstrate the performance of the quadrature.

Numerical Analysis
Numerical integration of ordinary differential equations with rapidly oscillatory factors

J.E. Bunder,A.J. Roberts

DOI: https://doi.org/10.1016/j.cam.2014.12.033

IF: 2.872

2015-07-01

Journal of Computational and Applied Mathematics

Abstract:We present a methodology for numerically integrating ordinary differential equations containing rapidly oscillatory terms. This challenge is distinct from that for differential equations which have rapidly oscillatory solutions: here the differential equation itself has the oscillatory terms. Our method generalises Filon quadrature for integrals, and is analogous to integral techniques designed to solve stochastic differential equations and, as such, is applicable to a wide variety of ordinary differential equations with rapidly oscillating factors. The proposed method flexibly achieves varying levels of accuracy depending upon the truncation of the expansion of certain integrals. Users will choose the level of truncation to suit the parameter regime of interest in their numerical integration.

mathematics, applied
Numerical evaluation of oscillatory integrals via automated steepest descent contour deformation

A. Gibbs,D.P. Hewett,D. Huybrechs

DOI: https://doi.org/10.1016/j.jcp.2024.112787

IF: 4.645

2024-01-25

Journal of Computational Physics

Abstract:Steepest descent methods combining complex contour deformation with numerical quadrature provide an efficient and accurate approach for the evaluation of highly oscillatory integrals. However, unless the phase function governing the oscillation is particularly simple, their application requires a significant amount of a priori analysis and expert user input, to determine the appropriate contour deformation, and to deal with the non-uniformity in the accuracy of standard quadrature techniques associated with the coalescence of stationary points (saddle points) with each other, or with the endpoints of the original integration contour. In this paper we present a novel algorithm for the numerical evaluation of oscillatory integrals with general polynomial phase functions, which automates the contour deformation process and avoids the difficulties typically encountered with coalescing stationary points and endpoints. The inputs to the algorithm are simply the phase and amplitude functions, the endpoints and orientation of the original integration contour, and a small number of numerical parameters. By a series of numerical experiments we demonstrate that the algorithm is accurate and efficient over a large range of frequencies, even for examples with a large number of coalescing stationary points and with endpoints at infinity. As a particular application, we use our algorithm to evaluate cuspoid canonical integrals from scattering theory. A Matlab implementation of the algorithm is made available and is called PathFinder.

physics, mathematical,computer science, interdisciplinary applications
Numerical integrator for highly oscillatory differential equations based on the Neumann series

DOI: https://doi.org/10.1007/s11075-024-01841-9

IF: 2.37

2024-04-26

Numerical Algorithms

Abstract:We propose a third-order numerical integrator based on the Neumann series and the Filon quadrature, designed mainly for highly oscillatory partial differential equations. The method can be applied to equations that exhibit small or moderate oscillations; however, counter-intuitively, large oscillations increase the accuracy of the scheme. With the proposed approach, the convergence order of the method can be easily improved. Error analysis of the method is also performed. We consider linear evolution equations involving first- and second-time derivatives that feature elliptic differential operators, such as the heat equation or the wave equation. Numerical experiments consider the case in which the space dimension is greater than one and confirm the theoretical study.

mathematics, applied
Gaussian quadrature rules for composite highly oscillatory integrals

Menghan Wu,Haiyong Wang

2023-06-08

Abstract:Highly oscillatory integrals of composite type arise in electronic engineering and their calculations is a challenging problem. In this paper, we propose two Gaussian quadrature rules for computing such integrals. The first one is constructed based on the classical theory of orthogonal polynomials and its nodes and weights can be computed efficiently by using tools of numerical linear algebra. We show that the rate of convergence of this rule depends solely on the regularity of the non-oscillatory part of the integrand. The second one is constructed with respect to a sign-changing function and the classical theory of Gaussian quadrature can not be used anymore. We explore theoretical properties of this Gaussian quadrature, including the trajectories of the quadrature nodes and the convergence rate of these nodes to the endpoints of the integration interval, and prove its asymptotic error estimate under suitable hypotheses. Numerical experiments are presented to demonstrate the performance of the proposed methods.

Numerical Analysis
The Numerical Steepest Descent Path Method for Calculating Physical Optics Integrals on Smooth Conducting Quadratic Surfaces

Yu Mao Wu,Li Jun Jiang,Wei E. I. Sha,Weng Cho Chew

DOI: https://doi.org/10.1109/tap.2013.2259788

IF: 5.7

2013-01-01

IEEE Transactions on Antennas and Propagation

Abstract:In this paper, we use the numerical steepest descent path (NSDP) method to analyze the highly oscillatory physical optics (PO) integral on smooth conducting parabolic surfaces, including both monostatic and bistatic cases. Quadratic variations of the amplitude and phase functions are used to approximate the integrand of PO integral. Then the surface PO integral is reduced into several highly oscillatory line integrals. By invoking the NSDP method, these highly oscillatory PO line integrals are defined on the corresponding NSDPs. Furthermore, the critical point contributions for the PO integral are exactly extracted and represented based on the NSDPs. The proposed NSDP method for calculating the PO integral on the smooth conducting surfaces is frequency-independent and error-controllable. Compared with the traditional asymptotic expansion approach, the NSDP method significantly improves the PO integral accuracy by around two digits when the working wave frequencies are not extremely large. Numerical results are given to validate the NSDP method.
On quadrature for singular integral operators with complex symmetric quadratic forms

Jeremy Hoskins,Manas Rachh,Bowei Wu

2023-12-12

Abstract:This paper describes a trapezoidal quadrature method for the discretization of weakly singular, singular and hypersingular boundary integral operators with complex symmetric quadratic forms. Such integral operators naturally arise when complex coordinate methods or complexified contour methods are used for the solution of time-harmonic acoustic and electromagnetic interface problems in three dimensions. The quadrature is an extension of a locally corrected punctured trapezoidal rule in parameter space wherein the correction weights are determined by fitting moments of error in the punctured trapezoidal rule, which is known analytically in terms of the Epstein zeta function. In this work, we analyze the analytic continuation of the Epstein zeta function and the generalized Wigner limits to complex quadratic forms; this analysis is essential to apply the fitting procedure for computing the correction weights. We illustrate the high-order convergence of this approach through several numerical examples.

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

James Bremer

DOI: https://doi.org/10.48550/arXiv.1705.07820

2017-05-22

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.

Numerical Analysis,Classical Analysis and ODEs
High asymptotic order numerical methods for highly oscillatory ODEs with large initial data

Zhongli Liu,Hongjiong Tian

DOI: https://doi.org/10.1016/j.aml.2024.109365

IF: 4.294

2024-11-11

Applied Mathematics Letters

Abstract:In this paper, we propose high asymptotic order numerical methods for solving highly oscillatory second order ODEs with large initial data, where the total energy of the system becomes unbounded as the oscillation frequency grows. The existing asymptotic-numerical solvers are especially designed for the classical energy bounded oscillatory equations, offering no insight into their performance with energy unbounded case. Based on the asymptotic expansion of the solution in the inverse powers of the oscillatory parameter, we propose an asymptotic numerical integrator to solve this class of highly oscillatory ODEs and discuss the computational efficiency for the case of polynomials. One numerical example is given to show the efficiency and accuracy of our proposed asymptotic-numerical solver.

mathematics, applied

Numerical Quadrature for Bessel Transformations with High Oscillations

Computational-Complexity on Numerical Integrals

The numerical integration of dual hankel transformation

On quadrature of highly oscillatory Bessel function via asymptotic analysis of simplex integrals

Quadrature Methods for Highly Oscillatory Singular Integrals

An Efficient Quadrature Rule for the Oscillatory Infinite Generalized Bessel Transform with a General Oscillator and Its Error Analysis

On the Convolution Quadrature Rule for Integral Transforms with Oscillatory Bessel Kernels

Analysis of Numerical Integration Error for Bessel Integral Identity in Fast Multipole Method for 2D Helmholtz Equation

An Efficient Quadrature Rule for Highly Oscillatory Integrals with Airy Function

On the Numerical Evaluation of a Class of Oscillatory Integrals in Worldline Variational Calculations

High-Order Numerical Integration on Domains Bounded by Intersecting Level Sets

Computation of highly oscillatory integrals using a Fourier extension approximation

Numerical integration of ordinary differential equations with rapidly oscillatory factors

Numerical evaluation of oscillatory integrals via automated steepest descent contour deformation

Numerical integrator for highly oscillatory differential equations based on the Neumann series

Gaussian quadrature rules for composite highly oscillatory integrals

The Numerical Steepest Descent Path Method for Calculating Physical Optics Integrals on Smooth Conducting Quadratic Surfaces

On quadrature for singular integral operators with complex symmetric quadratic forms

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

High asymptotic order numerical methods for highly oscillatory ODEs with large initial data