A fast convolution method for the fractional Laplacian in \mathbb{R}
From MaRDI portal
Publication:6420101
arXiv2212.05143MaRDI QIDQ6420101FDOQ6420101
Jorge Cayama, Carlos J. García-Cervera, Francisco de la Hoz, C. M. Cuesta
Publication date: 9 December 2022
Abstract: In this article, we develop a new method to approximate numerically the fractional Laplacian of functions defined on , as well as some more general singular integrals. After mapping into a finite interval, we discretize the integral operator using a modified midpoint rule. The result of this procedure can be cast as a discrete convolution, which can be evaluated efficiently using the Fast-Fourier Transform (FFT). The method provides an efficient, second order accurate, approximation to the fractional Laplacian, without the need to truncate the domain. We first prove that the method gives a second-order approximation for the fractional Laplacian and other related singular integrals; then, we detail the implementation of the method using the fast convolution, and give numerical examples that support its efficacy and efficiency; finally, as an example of its applicability to an evolution problem, we employ the method for the discretization of the nonlocal part of the one-dimensional cubic fractional Schr"odinger equation in the focusing case.
Numerical quadrature and cubature formulas (65D32) Fractional derivatives and integrals (26A33) Fractional partial differential equations (35R11)
This page was built for publication: A fast convolution method for the fractional Laplacian in $\mathbb{R}$
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6420101)
