Fast convolution with radial kernels at nonequispaced knots

  • Authors:
  • Daniel Potts;Gabriele Steidl;Arthur Nieslony

  • Affiliations:
  • Institute of Mathematics, University of Lübeck, 23560, Lübeck, Germany;Faculty of Mathematics and Computer Science, University of Mannheim, 68131, Mannheim, Germany;Faculty of Mathematics and Computer Science, University of Mannheim, 68131, Mannheim, Germany

  • Venue:
  • Numerische Mathematik
  • Year:
  • 2004

Quantified Score

Hi-index 0.01

Visualization

Abstract

We develop a new algorithm for the fast evaluation of linear combinations of radial functions * based on the recently developed fast Fourier transform at nonequispaced knots. For smooth kernels, e.g. the Gaussian, our algorithm requires * arithmetic operations. In case of singular kernels an additional regularization procedure must be incorporated and the algorithm has the arithmetic complexity * if either the points yj or the points xk are “reasonably uniformly distributed”. We prove error estimates to obtain clues about the choice of the involved parameters and present numerical examples for various singular and smooth kernels in two dimensions.