English

High-order Corrected Trapezoidal Rules for Functions with Fractional Singularities

Numerical Analysis 2022-03-22 v2 Numerical Analysis

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 i,j{1,2}i,j\in\{1,2\} and ϕCcN\phi\in C_c^N for N2N\geq 2. 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 2p+4α2p+4-\alpha, where pN0p\in\mathbb{N}_{0} is associated with total number of correction weights. Although we work in 2D setting, we formulate definitions and theorems in nNn\in\mathbb{N} dimensions when appropriate for the sake of generality. We present numerical experiments to validate the order of convergence of the proposed modified quadrature rules.

Keywords

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

R2 v1 2026-06-24T06:43:28.497Z