A high-order fast method for computing convolution integral with smooth kernel

作者：Qiang Ji^*

来源：Computer Physics Communications, 2010, 181(2): 313-316.

DOI：10.1016/j.cpc.2009.10.005

• Convolution integral
• Green's function
• Simpson rule
• FFT

In this paper we report on a high-order fast method to numerically calculate convolution integral with smooth non-periodic kernel. This method is based on the Newton-Cotes quadrature rule for the integral approximation and an FFT method for discrete summation. The method can have an arbitrarily high-order accuracy in principle depending on the number of points used in the integral approximation and a computational cost of O (N log(N)), where N is the number of grid points. For a three-point Simpson rule approximation, the method has an accuracy of O(h(4)), where h is the size of the computational grid. Applications of the Simpson rule based algorithm to the calculation of a one-dimensional continuous Gauss transform and to the calculation of a two-dimensional electric field from a charged beam are also presented.

• 出版日期2010-2