Efficient Monte Carlo Method for Integral Fractional Laplacian in Multiple Dimensions

by   Changtao Sheng, et al.

In this paper, we develop a Monte Carlo method for solving PDEs involving an integral fractional Laplacian (IFL) in multiple dimensions. We first construct a new Feynman-Kac representation based on the Green function for the fractional Laplacian operator on the unit ball in arbitrary dimensions. Inspired by the "walk-on-spheres" algorithm proposed in [24], we extend our algorithm for solving fractional PDEs in the complex domain. Then, we can compute the expectation of a multi-dimensional random variable with a known density function to obtain the numerical solution efficiently. The proposed algorithm finds it remarkably efficient in solving fractional PDEs: it only needs to evaluate the integrals of expectation form over a series of inside ball tangent boundaries with the known Green function. Moreover, we carry out the error estimates of the proposed method for the n-dimensional unit ball. Finally, ample numerical results are presented to demonstrate the robustness and effectiveness of this approach for fractional PDEs in unit disk and complex domains, and even in ten-dimensional unit balls.


page 19

page 20


Fast Fourier-like Mapped Chebyshev Spectral-Galerkin Methods for PDEs with Integral Fractional Laplacian in Unbounded Domains

In this paper, we propose a fast spectral-Galerkin method for solving PD...

Numerical Approximation of the Fractional Laplacian on R Using Orthogonal Families

In this paper, using well-known complex variable techniques, we compute ...

Extension of δ_-ziti method in the unit ball: Numerical integration, resolution of Poisson's problem and Heat transfer

Inspired by the Galerkin and particular method, a new approximation appr...

A universal solution scheme for fractional and classical PDEs

We propose a unified meshless method to solve classical and fractional P...

Fractional Laplacian - Quadrature rules for singular double integrals in 3D

In this article, quadrature rules for the efficient computation of the s...

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...