A Generalized Fast Multipole Method for Nonoscillatory Kernels

Zydrunas Gimbutas, Vladimir Rokhlin
2003 SIAM Journal on Scientific Computing  
We present a modification of the Fast Multipole Method (FMM) in two dimensions. While previous implementations of the FMM have been designed for harmonic kernels, our algorithm works for a large class of kernels that satisfy fairly general conditions, amounting to the kernel being sufficiently smooth away from the diagonal. Our algorithm approximates appropriately chosen parts of the kernel with "tensor products" of Legendre expansions and uses the Singular Value Decomposition (SVD) to compress
more » ... the resulting representations. The obtained singular function expansions replace the Taylor and Laurent expansions used in the original FMM. The algorithm requires O(N) operations, and is stable and robust. The performance of the algorithm is illustrated with numerical examples.
doi:10.1137/s1064827500381148 fatcat:q4h344mubne37nexjyyq64htre