Preconditioning for radial basis function partition of unity methods

Alfa Heryudono, Elisabeth Larsson, Alison Ramage, Lina von Sydow

Research output: Contribution to journalArticlepeer-review

37 Citations (Scopus)
180 Downloads (Pure)


Meshfree radial basis function (RBF) methods are of interest for solving partial differential equations due to attractive convergence properties, flexibility with respect to geometry, and ease of implementation. For global RBF methods, the computational cost grows rapidly with dimension and problem size, so localised approaches, such as partition of unity or stencil based RBF methods, are currently being developed. An RBF partition of unity method (RBF--PUM) approximates functions through a combination of local RBF approximations. The linear systems that arise are locally unstructured, but with a global structure due to the partitioning of the domain. Due to the sparsity of the matrices, for large scale problems, iterative solution methods are needed both for computational reasons and to reduce memory requirements. In this paper we implement and test different algebraic preconditioning strategies based on the structure of the matrix in combination with incomplete factorisations. We compare their performance for different orderings and problem settings and find that a no-fill incomplete factorisation of the central band of the original discretisation matrix provides a robust and efficient preconditioner.
Original languageEnglish
Pages (from-to)1089-1109
Number of pages21
JournalJournal of Scientific Computing
Issue number3
Publication statusPublished - 19 Oct 2015


  • radial basis function
  • partition of unity
  • iterative method
  • preconditioning
  • algebraic preconditioner


Dive into the research topics of 'Preconditioning for radial basis function partition of unity methods'. Together they form a unique fingerprint.

Cite this