A short FE implementation for a 2d homogeneous Dirichlet problem of a Fractional Laplacian
For mathematicians and scientists needing a practical, modifiable code for 2D fractional Laplacian problems.
The paper provides a simple 2D MATLAB finite element code for the homogeneous Dirichlet problem of a fractional Laplacian, filling a gap for researchers needing numerical approximations of nonlocal problems.
In \cite{AcostaBorthagaray}, a complete $n$-dimensional finite element analysis of the homogeneous Dirichlet problem associated to a fractional Laplacian was presented. Here we provide a comprehensive and simple 2D {\it MATLAB}\textsuperscript{\textregistered} finite element code for such a problem. The code is accompanied with a basic discussion of the theory relevant in the context. The main program is written in about 80 lines and can be easily modified to deal with other kernels as well as with time dependent problems. The present work fills a gap by providing an input for a large number of mathematicians and scientists interested in numerical approximations of solutions of a large variety of problems involving nonlocal phenomena in two-dimensional space.