Asymptotic expansions and fast computation of oscillatory Hilbert transforms

In this paper, we study the asymptotics and fast computation of the one-sided oscillatory Hilbert transforms of the form $$H^{+}(f(t)e^{i\omega t})(x)=-int_{0}^{\infty}e^{i\omega t}\frac{f(t)}{t-x}dt,\qquad \omega>0,\qquad x\geq 0,$$ where the bar indicates the Cauchy principal value and $f$ is a real-valued function with analytic continuation in the first quadrant, except possibly a branch point of algebraic type at the origin. When $x=0$, the integral is interpreted as a Hadamard finite-part integral, provided it is divergent. Asymptotic expansions in inverse powers of $\omega$ are derived for each fixed $x\geq 0$, which clarify the large $\omega$ behavior of this transform. We then present efficient and affordable approaches for numerical evaluation of such oscillatory transforms. Depending on the position of $x$, we classify our discussion into three regimes, namely, $x=\mathcal{O}(1)$ or $x\gg1$, $0<x\ll 1$ and $x=0$. Numerical experiments show that the convergence of the proposed methods greatly improve when the frequency $\omega$ increases. Some extensions to oscillatory Hilbert transforms with Bessel oscillators are briefly discussed as well.Comment: 32 pages, 6 figures, 4 table

Semi-analytical description of the modulator section of the coherent electron cooling

In the coherent electron cooling, the modern hadron beam cooling technique, each hadron receives an individual kick from the electric field of the amplified electron density perturbation created in the modulator by this hadron in a co-propagating electron beam. We developed a method for computing the dynamics of these density perturbations in an infinite electron plasma with any equilibrium velocity distribution -- a possible model for the modulator. We derived analytical expressions for the dynamics of the density perturbations in the Fourier-Laplace domain for a variety of 1D, 2D, and 3D equilibrium distributions of the electron beam. To obtain the space-time dynamics, we employed the fast Fourier transform (FFT) algorithm. We also found an analytical solution in the space-time domain for the 1D Cauchy equilibrium distribution, which serves as a benchmark for our general approach based on numerical evaluation of the integral transforms and as a fast alternative to the numerical computations. We tested the method for various distributions and initial conditions.Comment: 11 pages, 5 figure

A fast and well-conditioned spectral method for singular integral equations

We develop a spectral method for solving univariate singular integral equations over unions of intervals by utilizing Chebyshev and ultraspherical polynomials to reformulate the equations as almost-banded infinite-dimensional systems. This is accomplished by utilizing low rank approximations for sparse representations of the bivariate kernels. The resulting system can be solved in ${\cal O}(m^2n)$ operations using an adaptive QR factorization, where $m$ is the bandwidth and $n$ is the optimal number of unknowns needed to resolve the true solution. The complexity is reduced to ${\cal O}(m n)$ operations by pre-caching the QR factorization when the same operator is used for multiple right-hand sides. Stability is proved by showing that the resulting linear operator can be diagonally preconditioned to be a compact perturbation of the identity. Applications considered include the Faraday cage, and acoustic scattering for the Helmholtz and gravity Helmholtz equations, including spectrally accurate numerical evaluation of the far- and near-field solution. The Julia software package SingularIntegralEquations.jl implements our method with a convenient, user-friendly interface

Propagation Speed of the Maximum of the Fundamental Solution to the Fractional Diffusion-Wave Equation

In this paper, the one-dimensional time-fractional diffusion-wave equation with the fractional derivative of order $1 \le \alpha \le 2$ is revisited. This equation interpolates between the diffusion and the wave equations that behave quite differently regarding their response to a localized disturbance: whereas the diffusion equation describes a process, where a disturbance spreads infinitely fast, the propagation speed of the disturbance is a constant for the wave equation. For the time fractional diffusion-wave equation, the propagation speed of a disturbance is infinite, but its fundamental solution possesses a maximum that disperses with a finite speed. In this paper, the fundamental solution of the Cauchy problem for the time-fractional diffusion-wave equation, its maximum location, maximum value, and other important characteristics are investigated in detail. To illustrate analytical formulas, results of numerical calculations and plots are presented. Numerical algorithms and programs used to produce plots are discussed.Comment: 22 pages 6 figures. This paper has been presented by F. Mainardi at the International Workshop: Fractional Differentiation and its Applications (FDA12) Hohai University, Nanjing, China, 14-17 May 201

The exponentially convergent trapezoidal rule

It is well known that the trapezoidal rule converges geometrically when applied to analytic functions on periodic intervals or the real line. The mathematics and history of this phenomenon are reviewed and it is shown that far from being a curiosity, it is linked with computational methods all across scientific computing, including algorithms related to inverse Laplace transforms, special functions, complex analysis, rational approximation, integral equations, and the computation of functions and eigenvalues of matrices and operators

Conditioning moments of singular measures for entropy optimization. I

In order to process a potential moment sequence by the entropy optimization method one has to be assured that the original measure is absolutely continuous with respect to Lebesgue measure. We propose a non-linear exponential transform of the moment sequence of any measure, including singular ones, so that the entropy optimization method can still be used in the reconstruction or approximation of the original. The Cauchy transform in one variable, used for this very purpose in a classical context by A.\ A.\ Markov and followers, is replaced in higher dimensions by the Fantappi\`{e} transform. Several algorithms for reconstruction from moments are sketched, while we intend to provide the numerical experiments and computational aspects in a subsequent article. The essentials of complex analysis, harmonic analysis, and entropy optimization are recalled in some detail, with the goal of making the main results more accessible to non-expert readers. Keywords: Fantappi\`e transform; entropy optimization; moment problem; tube domain; exponential transformComment: Submitted to Indagnationes Mathematicae, I. Gohberg Memorial issu