Related papers: Computation of highly oscillatory Bessel transform…
We present a fast algorithm for computing the diffracted field from arbitrary binary (sharp-edged) planar apertures and occulters in the scalar Fresnel approximation, for up to moderately high Fresnel numbers ($\lesssim 10^3$). It uses a…
The usual fluid equations describing the large-scale evolution of mass density in the universe can be written as local in the density, velocity divergence, and velocity potential fields. As a result, the perturbative expansion in small…
We propose a novel method, called the dimension-changing transformation (DCT), to compute one-loop Feynman integrals and recently introduced fixed-branch integrals to arbitrary orders in $\epsilon$. The DCT relates one-loop Feynman…
We present a new efficient method to compute the angular power spectra of large-scale structure observables that circumvents the numerical integration over Bessel functions, expanding on a recently proposed algorithm based on FFTlog. This…
In this paper, we introduce an algorithm that provides approximate solutions to semi-linear ordinary differential equations with highly oscillatory solutions, which, after an appropriate change of variables, can be rewritten as…
Cosmological perturbation theory is a powerful tool to predict the statistics of large-scale structure in the weakly non-linear regime, but even at 1-loop order it results in computationally expensive mode-coupling integrals. Here we…
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…
We propose two easy-to-implement fast algorithms based on moment-matching to compute the nonlocal potential $\varphi(\textbf{x})=(U\ast \rho)(\textbf{x})$ on bounded domain, where the kernel $U$ is singular at the origin and the density…
We present a new method for the numerical solution of singular integral equations on the real axis. The method's value stems from an explicit formula for the Cauchy integral of a complex exponential multiplied by a rational function.…
Computation of the spherical harmonic rotation coefficients or elements of Wigner's d-matrix is important in a number of quantum mechanics and mathematical physics applications. Particularly, this is important for the Fast Multipole Methods…
The moments of Bessel functions and Bessel-trigonometric functions play a basic role in many practical problems and numerical analysis. This paper presents a complete analysis for these moments based on the recursive relations of Bessel…
Efficient methods are proposed, for computing integrals appeaing in electronic structure calculations. The methods consist of two parts: the first part is to represent the integrals as contour integrals and the second one is to evaluate the…
Higher precision efficient computation of period 1 relaxation oscillations of strongly nonlinear and singularly perturbed Rayleigh equations with external periodic forcing is presented. The computations are performed in the context of…
Fortran 77 programs for the computation of modified Bessel functions of purely imaginary order are presented. The codes compute the functions $K_{ia}(x)$, $L_{ia}(x)$ and their derivatives for real $a$ and positive $x$; these functions are…
A method of numerically evaluating slowly convergent monotone series is described. First, we apply a condensation transformation due to Van Wijngaarden to the original series. This transforms the original monotone series into an alternating…
We introduce two efficient algorithms for computing the partial Fourier transforms in one and two dimensions. Our study is motivated by the wave extrapolation procedure in reflection seismology. In both algorithms, the main idea is to…
Computing resonance and anti-resonance backbone curves in complex nonlinear mechanical systems is of high engineering relevance but remains computationally challenging, especially for large finite-element (FE) models. Existing…
The modified Bessel function of the second kind K$\nu$ appears in a wide variety of applied scientific fields. While its use is greatly facilitated by an implementation in most numerical libraries, overflow issues can be encountered…
Scalar diffraction calculations such as the angular spectrum method (ASM) and Fresnel diffraction, are widely used in the research fields of optics, X-rays, electron beams, and ultrasonics. It is possible to accelerate the calculation using…
A new highly accurate numerical approximation scheme based on a Gauss type Clenshaw-Curtis Quadrature for Fredholm integral equations of the second kind, whose kernel is either discontinuous or not smooth along the main diagonal, is…