Related papers: The Spectral Ewald method for singly periodic doma…
We propose an efficient and accurate solver for the nonlocal potential in the Davey-Stewartson equation using nonuniform FFT (NUFFT). A discontinuity in the Fourier transform of the nonlocal potential causes accuracy locking if the…
We apply a new numerical method, the singular Fourier-Pad\'e (SFP) method invented by Driscoll and Fornberg (2001, 2011), to price European-type options in L\'evy and affine processes. The motivation behind this application is to reduce the…
Pseudospectral time domain (PSTD) methods are widely used in many branches of acoustics for the numerical solution of the wave equation, including biomedical ultrasound and seismology. The use of the Fourier collocation spectral method in…
Spectral embedding uses eigenfunctions of the discrete Laplacian on a weighted graph to obtain coordinates for an embedding of an abstract data set into Euclidean space. We propose a new pre-processing step of first using the eigenfunctions…
In this paper, we propose a computationally efficient iterative algorithm for proper orthogonal decomposition (POD) using random sampling based techniques. In this algorithm, additional rows and columns are sampled and a merging technique…
In this paper, we propose a computationally efficient iterative algorithm for proper orthogonal decomposition (POD) using random sampling based techniques. In this algorithm, additional rows and columns are sampled and a merging technique…
The standard particle-in-cell algorithm suffers from grid heating. There exists a gridless alternative which bypasses the deposition step and calculates each Fourier mode of the charge density directly from the particle positions. We show…
In solid-state physics, energies of crystals are usually computed with a plane-wave discretization of Kohn-Sham equations. However the presence of Coulomb singularities requires the use of large plane-wave cut-offs to produce accurate…
The accurate calculation of guided electromagnetic modes in optical nanofibres is critical for applications in nanophotonics, from quantum interfaces to vectorial light sensing. Standard textbook methods rely on solving a $4\times4$ matrix…
abridged] A method to rapidly estimate the Fourier power spectrum of a point distribution is presented. This method relies on a Taylor expansion of the trigonometric functions. It yields the Fourier modes from a number of FFTs, which is…
Fast Multipole Methods (FMMs) based on the oscillatory Helmholtz kernel can reduce the cost of solving N-body problems arising from Boundary Integral Equations (BIEs) in acoustic or electromagnetics. However, their cost strongly increases…
This paper presents an efficient parallel direct algorithm with near-optimal complexity for the compact fourth and sixth-order approximation of the three-dimensional Helmholtz equations [1] with the problem coefficient depending on only one…
In this article, we present an $O(N \log N)$ rapidly convergent algorithm for the numerical approximation of the convolution integral with radially symmetric weakly singular kernels and compactly supported densities. To achieve the reduced…
We extend the publicly available quantumfdtd code. It was originally intended for solving the time-independent three-dimensional Schr\"odinger equation via the finite-difference time-domain (FDTD) method and for extracting the ground,…
The plane wave method is most widely used for solving the Kohn-Sham equations in first-principles materials science computations. In this procedure, the three-dimensional (3-dim) trial wave functions' fast Fourier transform (FFT) is a…
If the phase retrieval problem can be solved by a method similar to that of solving a system of linear equations under the context of FFT, the time complexity of computer based phase retrieval algorithm would be reduced. Here I present such…
We discuss a method based on a segmentary approximation of solutions of the Schr\"odinger by quadratic splines, for which the coefficients are determined by a variational method that does not require the resolution of complicated algebraic…
We present here a Finite Element Method devoted to the simulation of 3D periodic structures of arbitrary geometry. The numerical method based on ARPACK and PARDISO libraries, is discussed with the aim of extracting the eigenmodes of…
In this article, we first introduce a singular fractional Sturm-Liouville eigen-problems (SFSLP) on unbounded domain. The associated fractional differential operators in these problems are both Weyl and Caputo type . The properties of…
In this paper, we investigate the performance of pseudo-spectral methods in computing nearly singular solutions of fluid dynamics equations. We consider two different ways of removing the aliasing errors in a pseudo-spectral method. The…