English

SinCoTrap: A High-Order Locally Corrected Trapezoidal Rule for Periodic Singular Integrals in Arbitrary Dimensions

Numerical Analysis 2026-07-14 v1

Abstract

We present SinCoTrap (Singularity-Corrected Trapezoidal Rule), a high-order locally corrected trapezoidal method for periodic singular integrals in arbitrary dimension dd with kernel xs|\boldsymbol{x}|^{-s}, 0<s<d0<s<d. The scheme preserves the uniform tensor grid and modifies only a fixed, small stencil of weights near the singularity. For a correction order pp, the resulting quadrature attains the error rate O(h2p+2+ds)O(h^{2p+2+d-s}). We derive explicit, mesh-independent limiting correction weights via analytic continuation of a special generalization of the Riemann zeta function, yielding rapidly computable formulas that can be pretabulated for each (d,s,p)(d,s,p). This makes SinCoTrap both efficient in application and robust for high-order accuracy across a broad class of periodic singular integrals.

Keywords

Cite

@article{arxiv.2607.12390,
  title  = {SinCoTrap: A High-Order Locally Corrected Trapezoidal Rule for Periodic Singular Integrals in Arbitrary Dimensions},
  author = {Hengzhun Chen and Yingzhou Li},
  journal= {arXiv preprint arXiv:2607.12390},
  year   = {2026}
}