High-order Corrected Trapezoidal Rules for Functions with Fractional Singularities
Abstract
In this paper, we introduce and analyze arbitrarily high-order quadrature rules for evaluating the two-dimensional singular integrals of the forms \begin{align} I_{i,j} = \int_{\mathbb{R}^2}\phi(x)\frac{x_ix_j}{|x|^{2+\alpha}} \d x, \quad 0< \alpha < 2 \end{align} where and for . This type of singular integrals and its quadrature rule appear in the numerical discretization of fractional Laplacian in non-local Fokker-Planck Equations in 2D. The quadrature rules are trapezoidal rules equipped with correction weights for points around singularity. We prove the order of convergence is , where is associated with total number of correction weights. Although we work in 2D setting, we formulate definitions and theorems in dimensions when appropriate for the sake of generality. We present numerical experiments to validate the order of convergence of the proposed modified quadrature rules.
Cite
@article{arxiv.2110.03838,
title = {High-order Corrected Trapezoidal Rules for Functions with Fractional Singularities},
author = {Senbao Jiang and Xiaofan Li},
journal= {arXiv preprint arXiv:2110.03838},
year = {2022}
}
Comments
Updated version