Approximating and preconditioning the stiffness matrix in the GoFD approximation of the fractional Laplacian
arXiv:2407.07199 · doi:10.4208/cicp.OA-2024-0079
Abstract
In the finite difference approximation of the fractional Laplacian the stiffness matrix is typically dense and needs to be approximated numerically. The effect of the accuracy in approximating the stiffness matrix on the accuracy in the whole computation is analyzed and shown to be significant. Four such approximations are discussed. While they are shown to work well with the recently developed grid-over finite difference method (GoFD) for the numerical solution of boundary value problems of the fractional Laplacian, they differ in accuracy, economics to compute, performance of preconditioning, and asymptotic decay away from the diagonal line. In addition, two preconditioners based on sparse and circulant matrices are discussed for the iterative solution of linear systems associated with the stiffness matrix. Numerical results in two and three dimensions are presented.
arXiv comment28 pages
References in corpus (10)
- A fractional Laplace equation: regularity of solutions and Finite Element approximations
- A short FE implementation for a 2d homogeneous Dirichlet problem of a Fractional Laplacian
- 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
- Computing the ground and first excited states of the fractional Schrodinger equation in an infinite potential well
- A unified meshfree pseudospectral method for solving both classical and fractional PDEs
- Local convergence of the FEM for the integral fractional Laplacian
- A grid-overlay finite difference method for the fractional Laplacian on arbitrary bounded domains
- A novel and simple spectral method for nonlocal PDEs with the fractional Laplacian
- Meshfree finite difference solution of homogeneous Dirichlet problems of the fractional Laplacian