Related papers: Interpolatory quadrature rules for oscillatory int…
We consider the computation of quadrature rules that are exact for a Chebyshev set of linearly independent functions on an interval $[a,b]$. A general theory of Chebyshev sets guarantees the existence of rules with a Gaussian property, in…
We apply a Gaussian state formalism to track fluctuating perturbations that act on the position and momentum quadrature variables of a harmonic oscillator. Following a seminal proposal by Tsang and Caves [Phys. Rev. Lett. 105, 123601…
We investigate estimating scalar oscillatory integrals by integrating by parts in directions based on $(x_1 \partial_{x_1} f(x) ,..., x_n \partial_{x_n}f(x))$, where $f(x)$ is the phase function. We prove a theorem which provides estimates…
In this paper we develop a theory for oscillatory integrals with complex phases. When $f:{\mathbb C}^n \to {\mathbb C}$, we evaluate this phase function on the basic character ${\rm e}(z) := e^{2\pi i x} e^{2\pi i y}$ of ${\mathbb C} \simeq…
A simple quantum mechanical model consisting of a discrete level resonantly coupled to a continuum of finite width, where the coupling can be varied from perturbative to strong (Fano-Anderson model), is considered. The particle is initially…
We present an algorithm for the efficient numerical evaluation of integrals of the form \[ I(\omega) = \int_0^1 F( x,\mathrm e^{\mathrm i \omega x}; \omega) \, \mathrm d x \] for sufficiently smooth but otherwise arbitrary $F$ and $\omega…
We provide tools to help automate the error analysis of algorithms that evaluate simple functions over the floating-point numbers. The aim is to obtain tight relative error bounds for these algorithms, expressed as a function of the unit…
An effective solution to the problem of Hermite $G^1$ interpolation with a clothoid curve is provided. At the beginning the problem is naturally formulated as a system of nonlinear equations with multiple solutions that is generally…
In this work, we develop a fast and accurate method for the scattering of flexural-gravity waves by a thin plate of varying thickness overlying a fluid of infinite depth. This problem commonly arises in the study of sea ice and ice shelves,…
In the framework of a Particle-In-Cell scheme for some 1D Vlasov-Poisson system depending on a small parameter, we propose a time-stepping method which is numerically uniformly accurate when the parameter goes to zero. Based on an…
We propose a Gauss-Legendre quadrature based sampling on the rotation group for the representation of a band-limited signal such that the Fourier transform (FT) of a signal can be exactly computed from its samples. Our figure of merit is…
In this paper, we present a meshless method belonging to the family of element-free Galerkin (EFG) methods. The distinguishing feature of the presented meshless method is that it allows accurate enforcement of essential boundary conditions.…
We propose an exponential integrator Fourier pseudospectral method DEI-FP for solving the "Good" Boussinesq (GB) equation. The numerical scheme is based on a Deuflhard-type exponential integrator and a Fourier pseudospectral method for…
The paper improves the accuracy of the one-dimensional fractional Fourier transform (FRFT) by leveraging closed Newton-Cotes quadrature rules. Using the weights derived from the Composite Newton-Cotes rules of order QN, we demonstrate that…
In this work we develop the Gaussian quadrature rule for weight functions involving fractional powers, exponentials and Bessel functions of the first kind. Besides the computation based on the use of the standard and the modified Chebyshev…
The implementation of the finite element method for linear elliptic equations requires to assemble the stiffness matrix and the load vector. In general, the entries of this matrix-vector system are not known explicitly but need to be…
The integrating factor and exponential time differencing methods are implemented and tested for solving the time-dependent Kohn--Sham equations. Popular time propagation methods used in physics, as well as other robust numerical approaches,…
A general approach was proposed in this article to develop high-order exponentially fitted basis functions for finite element approximations of multi-dimensional drift-diffusion equations for modeling biomolecular electrodiffusion…
The so-called non-uniform FFT (NFFT) is a family of algorithms for efficiently computing the Fourier transform of finite-length signals, whenever the time or frequency grid is nonuniformly spaced. Following the usual classification, there…
Unevenly spaced samples from a periodic function are common in signal processing and can often be viewed as a perturbed equally spaced grid. In this paper, we analyze how the uneven distribution of the samples impacts the quality of…