English

Quadrature for Singular Integrals over convex Polytopes

Numerical Analysis 2025-11-19 v1 Numerical Analysis

Abstract

A new algorithm for the efficient numerical approximation of weakly singular integrals over convex polytopes is introduced. Such integrals appear in the Galerkin discretizations of integral equations and nonlocal partial differential equations. The polytope is decomposed into a number of convex hulls of a singular and regular face. This expresses the singularity in a single variable which is effectively handled by Gauss-Jacobi quadrature. The decomposition algorithm is applicable to general finite polytopes. The Cartesian product of two simplices and two cubes will be discussed as special cases and numerical examples will be presented to illustrate the convergence of the resulting quadrature scheme.

Keywords

Cite

@article{arxiv.2511.13974,
  title  = {Quadrature for Singular Integrals over convex Polytopes},
  author = {Johannes Tausch},
  journal= {arXiv preprint arXiv:2511.13974},
  year   = {2025}
}

Comments

20 pages

R2 v1 2026-07-01T07:42:21.997Z