A Unified Meshfree Pseudospectral Method for Solving Both Classical and Fractional Pdes


In this paper, we propose a meshfree method based on the Gaussian radial basis function (RBF) to solve both classical and fractional PDEs. The proposed method takes advantage of the analytical Laplacian of Gaussian functions so as to accommodate the discretization of the classical and fractional Laplacians in a single framework and avoid the large computational cost for numerical evaluation of the fractional derivatives. These important merits distinguish our method from other existing methods for fractional PDEs. Moreover, our method is simple and easy when handling complex geometries and local refinements, and its computer program implementation re- mains the same for any dimension d ≥ 1. Extensive numerical experiments are provided to study the performance of our method in both approximating the Dirichlet Laplace operators and solving classical and fractional PDE problems. We show that our method has spectral accuracy and can achieve good approximation even with a small number of points. Compared to the recently proposed Wendland RBF method, our method exactly incorporates the Dirichlet boundary conditions into the scheme and is free of the Gibbs phenomenon as observed in the wendland method.


Mathematics and Statistics

Research Center/Lab(s)

Center for High Performance Computing Research

Keywords and Phrases

Classical Laplacian; Fractional Laplacian; Hypergeometric functions; Meshfree method; Pseudospectral method; Radial basis functions

International Standard Serial Number (ISSN)

1064-8275; 1095-7197

Document Type

Article - Journal

Document Version


File Type





© 2021 Society for Industrial and Applied Mathematics (SIAM), All rights reserved.

Publication Date

01 Jan 2021