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 align Ii,j = ∫R2φ(x)xixj|x|2+α x, 0< α < 2 align where i,j∈\1,2\ and φ∈ CcN for N≥ 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-α, where p∈N0 is associated with total number of correction weights. Although we work in 2D setting, we formulate definitions and theorems in n∈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.

0

Turn this paper into a full lesson

ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.

Discussion (0)

Sign in to join the discussion.

Loading comments…