Local MFS Matrix Decomposition Algorithms for Elliptic BVPs in Annuli

Abstract

We apply the local method of fundamental solutions (LMFS) to boundary value problems (BVPs) for the Laplace and homogeneous biharmonic equations in annuli. By appropriately choosing the collocation points, the LMFS discretization yields sparse block circulant system matrices. As a result, matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs) can be used for the solution of the systems resulting in considerable savings in both computational time and storage requirements. The accuracy of the method and its ability to solve large scale problems are demonstrated by applying it to several numerical experiments.