A Fast Treecode for Multiquadric Interpolation with Varying Shape Parameters

  • Authors:
  • Quan Deng;Tobin A. Driscoll

  • Affiliations:
  • dengquan@udel.edu and driscoll@udel.edu;-

  • Venue:
  • SIAM Journal on Scientific Computing
  • Year:
  • 2012

Quantified Score

Hi-index 0.00

Visualization

Abstract

A treecode algorithm is presented for the fast evaluation of multiquadric radial basis function (RBF) approximations. The method is a dual approach to one presented by Krasny and Wang, which applies far-field expansions to clusters of RBF centers (source points). The new approach clusters evaluation points instead and is therefore easily able to cope with basis functions that have different multiquadric shape parameters. The new treecode is able to evaluate an approximation on $N$ centers at $M$ points in $O((N+M) \log M)$ time in the ideal case when evaluation points are uniformly distributed. When coupled with a two-level restricted additive Schwarz preconditioner for GMRES iterations, the treecode is well suited for use within an adaptive RBF iteration, previously described by Driscoll and Heryudono, as is demonstrated by experiments on test functions.