A stable algorithm for divergence-free radial basis functions in the flat limit

Kathryn P. Drake, Grady B. Wright

Research output: Contribution to journalArticlepeer-review

5 Scopus citations

Abstract

The direct method used for calculating smooth radial basis function (RBF) interpolants in the flat limit becomes numerically unstable. The RBF-QR algorithm bypasses this ill-conditioning using a clever change of basis technique. We extend this method for computing interpolants involving matrix-valued kernels, specifically surface divergence-free RBFs on the sphere, in the flat limit. Results illustrating the effectiveness of this algorithm are presented for a divergence-free vector field on the sphere from samples at scattered points.

Original languageEnglish
Article number109595
JournalJournal of Computational Physics
Volume417
DOIs
StatePublished - 15 Sep 2020

Keywords

  • Divergence-free
  • Ill-conditioning
  • RBFQR
  • Radial basis functions
  • Shape parameter
  • Sphere

Fingerprint

Dive into the research topics of 'A stable algorithm for divergence-free radial basis functions in the flat limit'. Together they form a unique fingerprint.

Cite this