A simple solver for the fractional Laplacian in multiple dimensions

Abstract

We present a simple discretization scheme for the hypersingular integral representation of the fractional Laplace operator and solver for the corresponding fractional Laplacian problem. Through singularity subtraction, we obtain a regularized integrand that is amenable to the trapezoidal rule with equispaced nodes, assuming a high degree of regularity in the underlying function (i.e., $u\in C^6(R^d)$). The resulting quadrature scheme gives a discrete operator on a regular grid that is translation-invariant and thus can be applied quickly with the fast Fourier transform. For discretizations of problems related to space-fractional diffusion on bounded domains, we observe that the underlying linear system can be efficiently solved via preconditioned Krylov methods with a preconditioner based on the finite-difference (non-fractional) Laplacian. We show numerical results illustrating the error of our simple scheme as well the efficiency of our preconditioning approach, both for the elliptic (steady-state) fractional diffusion problem and the time-dependent problem.

27 pages, 3 figures

References in corpus (7)

• Aspects of an adaptive finite element method for the fractional Laplacian: a priori and a posteriori error estimates, efficient implementation and multigrid solver
• Towards an Efficient Finite Element Method for the Integral Fractional Laplacian on Polygonal Domains
• What Is the Fractional Laplacian?
• Regularity theory and high order numerical methods for the (1d)-Fractional Laplacian
• Hybrid Finite Element - Spectral Method for the Fractional Laplacian: Approximation Theory and Efficient Solver
• Finite difference methods for fractional Laplacians
• Spectral Method for the Fractional Laplacian in 2D and 3D

Cited by in corpus (2)

• Numerical methods for nonlocal and fractional models
• Efficient Numerical Method for Models Driven by Lévy Process via Hierarchical Matrices