Year: 2021
Author: Yong Zhang, Chijie Zhuang, Shidong Jiang
Communications in Computational Physics, Vol. 29 (2021), Iss. 5 : pp. 1570–1582
Abstract
Based on the recently-developed sum-of-exponential (SOE) approximation, in this article, we propose a fast algorithm to evaluate the one-dimensional convolution potential $φ(x)=K∗ρ=∫^1_{0}K(x−y)ρ(y)dy$ at (non)uniformly distributed target grid points {$x_i$}$^M_{i=1}$, where the kernel $K(x)$ might be singular at the origin and the source density function $ρ(x)$ is given on a source grid ${{{y_i}}}^N_{j=1}$ which can be different from the target grid. It achieves an optimal accuracy, inherited from the interpolation of the density $ρ(x)$, within $\mathcal{O}(M+N)$ operations. Using the kernel's SOE approximation $K_{ES}$, the potential is split into two integrals: the exponential convolution $φ_{ES}$=$K_{ES}∗ρ$ and the local correction integral $φ_{cor}=(K−K_{ES})∗ρ$. The exponential convolution is evaluated via the recurrence formula that is typical of the exponential function. The local correction integral is restricted to a small neighborhood of the target point where the kernel singularity is considered. Rigorous estimates of the optimal accuracy are provided. The algorithm is ideal for parallelization and favors easy extensions to complicated kernels. Extensive numerical results for different kernels are presented.
You do not have full access to this article.
Already a Subscriber? Sign in as an individual or via your institution
Journal Article Details
Publisher Name: Global Science Press
Language: English
DOI: https://doi.org/10.4208/cicp.OA-2020-0116
Communications in Computational Physics, Vol. 29 (2021), Iss. 5 : pp. 1570–1582
Published online: 2021-01
AMS Subject Headings: Global Science Press
Copyright: COPYRIGHT: © Global Science Press
Pages: 13
Keywords: One dimensional convolution sum of exponentials singular kernel discrete density.
Author Details
-
Fast $$\theta $$-Maruyama scheme for stochastic Volterra integral equations of convolution type: mean-square stability and strong convergence analysis
Wang, Mengjie | Dai, Xinjie | Yu, Yanyan | Xiao, AiguoComputational and Applied Mathematics, Vol. 42 (2023), Iss. 3
https://doi.org/10.1007/s40314-023-02248-3 [Citations: 2] -
Discrete Lehmann representation of imaginary time Green's functions
Kaye, Jason | Chen, Kun | Parcollet, OlivierPhysical Review B, Vol. 105 (2022), Iss. 23
https://doi.org/10.1103/PhysRevB.105.235115 [Citations: 29] -
Decomposing Imaginary-Time Feynman Diagrams Using Separable Basis Functions: Anderson Impurity Model Strong-Coupling Expansion
Kaye, Jason | Huang, Zhen | Strand, Hugo U. R. | Golež, DenisPhysical Review X, Vol. 14 (2024), Iss. 3
https://doi.org/10.1103/PhysRevX.14.031034 [Citations: 1] -
A Kernel-Independent Sum-of-Exponentials Method
Gao, Zixuan | Liang, Jiuyang | Xu, ZhenliJournal of Scientific Computing, Vol. 93 (2022), Iss. 2
https://doi.org/10.1007/s10915-022-01999-1 [Citations: 2] -
A Unified Structure Preserving Scheme for a Multispecies Model with a Gradient Flow Structure and Nonlocal Interactions via Singular Kernels
Zhang, Yong | Zhao, Yu | Zhou, ZhennanSIAM Journal on Scientific Computing, Vol. 43 (2021), Iss. 3 P.B539
https://doi.org/10.1137/20M1348911 [Citations: 1]