Related papers: Computation of highly oscillatory Bessel transform…
In a recent article the authors showed that the radiative Transfer equations with multiple frequencies and scattering can be formulated as a nonlinear integral system. In the present article, the formulation is extended to handle reflective…
We present a fast direct solver for boundary integral equations on complex surfaces in three dimensions using an extension of the recently introduced recursive strong skeletonization scheme. For problems that are not highly oscillatory, our…
It is demonstrated is this letter that linear multistep methods for integrating ordinary differential equations can be used to develop a family of fast forward scattering algorithms with higher orders of convergence. Excluding the cost of…
Orthogonal systems in $\mathrm{L}_2(\mathbb{R})$, once implemented in spectral methods, enjoy a number of important advantages if their differentiation matrix is skew-symmetric and highly structured. Such systems, where the differentiation…
We present numerical methods based on the fast Fourier transform (FFT) to solve convolution integral equations on a semi-infinite interval (Wiener-Hopf equation) or on a finite interval (Fredholm equation). We extend and improve a FFT-based…
We implement an efficient method of computation of two dimensional Fourier-type integrals based on approximation of the integrand by Gaussian radial basis functions, which constitute a standard tool in approximation theory. As a result, we…
A formula of Doetsch ({\em Math. Zeitschr.} {\bf 42}, 263 (1937)) is generalized and used to numerically invert the one-sided Laplace transform ${\hat C}(\beta)$. The necessary input is only the values of ${\hat C}(\beta)$ on the positive…
We develop a numerical method for computing with orthogonal polynomials that are orthogonal on multiple, disjoint intervals for which analytical formulae are currently unknown. Our approach exploits the Fokas--Its--Kitaev Riemann--Hilbert…
We present algorithms for computing weakly singular and near-singular integrals arising when solving the 3D Helmholtz equation with curved boundary elements. These are based on the computation of the preimage of the singularity in the…
We introduce a fast, high-precision algorithm for calculating intersections between great circle arcs and lines of constant latitude on the unit sphere. We first propose a simplified intersection point formula with improved speed and…
We introduce a fast Fourier spectral method for the spatially homogeneous Boltzmann equation with non-cutoff collision kernels. Such kernels contain non-integrable singularity in the deviation angle which arise in a wide range of…
Let $X$ be a complete, simply connected harmonic manifold of purely exponential volume growth. This class contains all non-flat harmonic manifolds of non-positive curvature and, in particular all known examples of harmonic manifolds except…
In this article, we prove certain Weber-Schafheitlin type integral formulae for Bessel functions over complex numbers. A special case is a formula for the Fourier transform of regularized Bessel functions on complex numbers. This is applied…
A proof is found for the elliptic integral evaluation of the Bessel moment $$M:=\int_0^\infty t I_0^2(t)K_0^2(t)K_0(2t) {\rm d}t ={1/12} {\bf K}(\sin(\pi/12)){\bf K}(\cos(\pi/12)) =\frac{\Gamma^6(\frac13)}{64\pi^22^{2/3}}$$ resulting from…
Morlet or Gabor wavelet transforms as well as Gaussian smoothing, are widely used in signal processing and image processing. However, the computational complexity of their direct calculations is proportional not only to the number of data…
Approximate $p$-point Leibniz derivation formulas as well as interpolatory Simpson quadrature sums adapted to oscillatory functions are discussed. Both theoretical considerations and numerical evidence concerning the dependence of the…
We present a general diagrammatic approach to the construction of efficient algorithms for computing a Fourier transform on a semisimple algebra. This extends previous work wherein we derive best estimates for the computation of a Fourier…
The paper is concerned with efficient numerical methods for solving a linear system $\phi(A) x= b$, where $\phi(z)$ is a $\phi$-function and $A\in \mathbb R^{N\times N}$. In particular in this work we are interested in the computation of…
The study of turbulent flows calls for measurements with high resolution both in space and in time. We propose a new approach to reconstruct High-Temporal-High-Spatial resolution velocity fields by combining two sources of information that…
The variational discrete element method developed in [28] for dynamic elasto-plastic computations is adapted to compute the deformation of elastic Cosserat materials. In addition to cellwise displacement degrees of freedom (dofs), cellwise…