Matrix decomposition RBF algorithm for solving 3D elliptic problems

In this study, we propose an efficient algorithm for the evaluation of the particular solutions of three-dimensional inhomogeneous elliptic partial differential equations using radial basis functions. The collocation points are placed on concentric spheres and thus the resulting global matrix possesses a block circulant structure. This structure is exploited to develop an efficient matrix decomposition algorithm for the solution of the resulting system. Further savings in the matrix decomposition algorithm are obtained by the use of fast Fourier transforms. The proposed algorithm is used, in conjunction with the method of fundamental solutions for the solution of three-dimensional inhomogeneous elliptic boundary value problems.