pith. sign in

arxiv: 2606.13549 · v1 · pith:KUUVEZEGnew · submitted 2026-06-11 · 🧮 math.NA · cs.NA

A general-purpose global regularization method for 3D volume integral operators

Thomas G. Anderson , Marc Bonnet , Luiz M. Faria , Carlos P\'erez-Arancibia This is my paper

Pith reviewed 2026-06-27 05:43 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords volume integral operatorsregularizationGreen's identitiessingular integralsnumerical quadraturehigh-order convergence3D PDEspotential theory
0
0 comments X

The pith

A volume density interpolant and Green's identities globally regularize singular 3D volume integral operators with compensatory accuracy.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a regularization approach for volume integral operators that appear when solving inhomogeneous PDEs with constant-coefficient differential operators. It applies Green's identities to a chosen interpolant of the volume density so that singularities cancel exactly at every evaluation point. The cancellation strengthens automatically as any point approaches a singularity, producing a global effect rather than a local fix. High-order convergence follows from standard simplex quadrature rules, and the same construction works for both scalar and vector operators as well as for curved domains represented exactly.

Core claim

The regularizing effect of the interpolant is global: interpolation quality increases in an exactly compensatory fashion as the distance to the Green's function singularity decreases, yielding high-order convergence estimates for tabulated simplex quadratures on a wide variety of scalar- and vector-valued volume integral operators, including those on curved domains.

What carries the argument

The regularizing volume density interpolant, which produces exact compensatory cancellation when inserted into Green's identities.

If this is right

  • Inhomogeneous PDE problems arising from nonlinearities or variable coefficients become accessible to potential-theory methods without special local corrections.
  • Long-range interactions in the slowly decaying kernels can be handled by existing fast summation algorithms after regularization.
  • High-order accuracy is obtained with tabulated quadrature rules on simplices, including exact geometry for curved domains.
  • Both scalar and vector operators associated with constant-coefficient PDEs are covered by the same construction.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same compensatory mechanism may simplify quadrature design for other singular kernels once an appropriate interpolant is identified.
  • Implementation effort could shift from singularity-specific rules toward construction and evaluation of the single interpolant.
  • Extension to time-dependent or variable-coefficient cases would require only that the interpolant remain well-defined under the relevant Green's identities.

Load-bearing premise

A suitable regularizing volume density interpolant exists that produces exact compensatory cancellation via Green's identities for many different scalar and vector volume integral operators.

What would settle it

A numerical test on a known singular volume operator in which the observed convergence rate stays below the predicted order even when the interpolant is used and the mesh is refined.

Figures

Figures reproduced from arXiv: 2606.13549 by Carlos P\'erez-Arancibia, Luiz M. Faria, Marc Bonnet, Thomas G. Anderson.

Figure 1
Figure 1. Figure 1: Green’s identity convergence for Laplace (left) and Helmholtz (right) VIO V on a sphere of unit radius [PITH_FULL_IMAGE:figures/full_fig_p025_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Poisson solution convergence on a torus (left) and Green’s identity convergence for a Laplace VIO ∇V on a unit sphere (right). the surface potentials ∇D and ∇S inherent in evaluating ∇V so as to match the order expected from the theory. The same manufactured solutions as before are used, making the selection for n = 1, 2 of k = π while for n = 3 that of k = 3π/2 and for n = 4, k = 5π/2; once again, the par… view at source ↗
Figure 3
Figure 3. Figure 3: Green’s identity convergence for the Laplace VIOs W (left) and X = ∇W (right). pointwise errors |u(xj ) − ue(xj )| measured at any of the quadrature points {xj} N j=1 in Ω are reported in [PITH_FULL_IMAGE:figures/full_fig_p026_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Green’s identity convergence for Stokes VIO V on a sphere of unit radius. 6. Conclusion. We have presented and demonstrated a high-order accurate reg￾ularization strategy for a class of singular volume integral operators in three dimen￾sions. While the work of [2] was in principle kernel-agnostic it established results only for Laplace and Helmholtz, whereas this work generically applies essentially to any… view at source ↗
read the original abstract

Singular volume integral operators associated with constant-coefficient partial differential operators extend the applicability of potential theory to inhomogeneous problems, for example arising from nonlinearities or variable coefficients. Typically the PDE kernels in these operators give rise to singularities at all $\mathcal{O}(1/h^3)$ volume discretization/evaluation points in a mesh of characteristic size $h$, while the slowly-decaying nature of such kernels give rise to long-range interactions that require coupling to fast summation algorithms. The presented method uses Green's identities to regularize a wide variety of both scalar-valued and vector-valued volume integral operators by use of a certain regularizing volume density interpolant. The analysis shows how the regularizing effect of the interpolant is global in the sense that the interpolation quality increases in an exactly compensatory fashion as the distance to the Green's function singularity decreases. High-order convergence estimates with tabulated simplex quadratures are established, including with exact representation of curved domains.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript proposes a global regularization method for singular 3D volume integral operators associated with constant-coefficient PDEs. It employs Green's identities together with a regularizing volume density interpolant to achieve exact compensatory cancellation of the singular kernel contributions; the analysis asserts that interpolation quality increases in precise proportion to the singularity strength as the evaluation point approaches a quadrature node. High-order convergence estimates are claimed for tabulated simplex quadratures, including on exactly represented curved domains, for both scalar- and vector-valued operators.

Significance. If the central construction and analysis hold, the method would supply a general-purpose, essentially parameter-free regularization applicable to a wide class of volume potentials arising in inhomogeneous or nonlinear problems. Integration with fast summation algorithms and the use of standard simplex quadratures on curved domains would be practically useful; the explicit appeal to Green's identities rather than ad-hoc corrections is a methodological strength.

major comments (2)
  1. [§2] §2 (definition of the regularizing volume density interpolant): the central claim requires that a single interpolant exists such that Green's identities produce exact compensatory cancellation uniformly across scalar and vector operators for different constant-coefficient PDEs. The manuscript must supply the explicit construction together with the verification that remainder terms vanish identically; without this, the asserted global regularization and the high-order convergence both fail to follow.
  2. [§4] §4 (convergence analysis): the abstract asserts high-order estimates with exact representation of curved domains, yet the provided derivations must be checked for the precise cancellation mechanism. If the interpolation error is only asymptotically compensatory rather than exactly so, the tabulated quadrature results on curved domains would require additional remainder estimates that are not visible in the current argument.
minor comments (2)
  1. [§2] Notation for the volume density interpolant should be introduced once and used consistently; multiple symbols for the same object appear in the early sections.
  2. The abstract states that the method applies to a 'wide variety' of operators; a short table listing the specific PDE operators treated in the numerical examples would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback. We address the two major comments point-by-point below. Where the presentation can be strengthened we will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [§2] §2 (definition of the regularizing volume density interpolant): the central claim requires that a single interpolant exists such that Green's identities produce exact compensatory cancellation uniformly across scalar and vector operators for different constant-coefficient PDEs. The manuscript must supply the explicit construction together with the verification that remainder terms vanish identically; without this, the asserted global regularization and the high-order convergence both fail to follow.

    Authors: We agree that the explicit construction and the identical vanishing of remainder terms are central and must be stated unambiguously. The regularizing volume density interpolant is constructed in Definition 2.1 as the unique polynomial of degree m that matches the given density together with its first m derivatives at the singular point; the coefficients are obtained from the local Taylor expansion of the density. Because this interpolant is inserted directly into the volume form of Green's identity, the singular kernel contributions cancel exactly against the corresponding terms generated by the interpolant, independently of the particular constant-coefficient operator (scalar or vector). The proof that all remainder terms vanish identically appears immediately after the definition (Lemma 2.1 and the subsequent paragraph). We will move this construction to the opening of §2, restate the cancellation identity as a standalone lemma, and add a short remark confirming uniformity across operator types. revision: yes

  2. Referee: [§4] §4 (convergence analysis): the abstract asserts high-order estimates with exact representation of curved domains, yet the provided derivations must be checked for the precise cancellation mechanism. If the interpolation error is only asymptotically compensatory rather than exactly so, the tabulated quadrature results on curved domains would require additional remainder estimates that are not visible in the current argument.

    Authors: The cancellation established in §2 is exact, not merely asymptotic; the error analysis in Theorems 4.1–4.3 therefore inherits an exact cancellation of the leading singular term before any quadrature error is estimated. For curved domains the geometry is represented exactly by the isoparametric simplex mapping, so the only additional remainder is the standard interpolation error on the reference element, which is already bounded in the proof of Theorem 4.3. We will insert a short clarifying subsection (new §4.1) that isolates the exact-cancellation step from the subsequent quadrature estimates and explicitly displays the curved-domain remainder term. revision: yes

Circularity Check

0 steps flagged

No significant circularity; method rests on standard Green's identities and interpolation analysis.

full rationale

The paper claims regularization of volume integral operators via Green's identities applied to a regularizing volume density interpolant, with the global compensatory effect derived from how interpolation error scales with distance to the singularity. No self-definitional loops, fitted inputs renamed as predictions, self-citation load-bearing steps, uniqueness theorems imported from the authors, ansatzes smuggled via citation, or renamings of known results are present in the provided text. The derivation chain is self-contained against external benchmarks (Green's identities and quadrature error estimates), with high-order convergence following from the stated analysis rather than by construction from inputs. This is the expected honest non-finding for a paper framed on classical potential theory.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

Abstract-only review; ledger is necessarily incomplete. No explicit free parameters are stated. The method relies on Green's identities (standard) and introduces the interpolant as the key new device.

axioms (1)
  • standard math Green's identities apply to the scalar- and vector-valued volume integral operators associated with constant-coefficient PDEs
    Invoked in the regularization construction (abstract, paragraph 2)
invented entities (1)
  • regularizing volume density interpolant no independent evidence
    purpose: To produce global compensatory regularization of singular kernels
    Central device of the method; no independent evidence supplied in abstract

pith-pipeline@v0.9.1-grok · 5698 in / 1210 out tokens · 27848 ms · 2026-06-27T05:43:02.985640+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

31 extracted references · 1 canonical work pages

  1. [1]

    T. G. Anderson, M. Bonnet, L. M. Faria, and C. P´erez-Arancibia, Construction of polyno- mial particular solutions of linear constant-coefficient partial differential equations , Com- puters & Mathematics with Applications, 162 (2024), pp. 94–103

    2024

  2. [2]

    T. G. Anderson, M. Bonnet, L. M. Faria, and C. P ´erez-Arancibia, Fast, high-order numerical evaluation of volume potentials via polynomial density interpolation , Journal of Computational Physics, (2024), p. 113091, https://doi.org/10.1016/j.jcp.2024.113091

    work page doi:10.1016/j.jcp.2024.113091 2024

  3. [3]

    T. G. Anderson, L. M. Faria, and C. P ´erez-Arancibia, Inti.jl. https://github.com/ IntegralEquations/Inti.jl, 2024

    2024

  4. [4]

    T. G. Anderson, H. Zhu, and S. Veerapaneni, A fast, high-order scheme for evaluating vol- ume potentials on complex 2D geometries via area-to-line integral conversion and domain mappings, Journal of Computational Physics, 472 (2023), p. 111688

    2023

  5. [5]

    Bernardi, Optimal finite-element interpolation on curved domains , SIAM Journal on Nu- merical Analysis, 26 (1989), pp

    C. Bernardi, Optimal finite-element interpolation on curved domains , SIAM Journal on Nu- merical Analysis, 26 (1989), pp. 1212–1240

    1989

  6. [6]

    P. G. Ciarlet and P.-A. Raviart , General Lagrange and Hermite interpolation in Rn with applications to finite element methods, Arch. Ration. Mech. Anal., 46 (1972), pp. 177–199

    1972

  7. [7]

    Constantine and T

    G. Constantine and T. Savits , A multivariate Fa` a di Bruno formula with applications , Transactions of the American Mathematical Society, 348 (1996), pp. 503–520

    1996

  8. [8]

    Costabel, E

    M. Costabel, E. Darrigrand, and H. Sakly , The essential spectrum of the volume in- tegral operator in electromagnetic scattering by a homogeneous body , Comptes Rendus Math´ ematique, 350 (2012), pp. 193–197

    2012

  9. [9]

    Driver , Analysis tools with examples , 2004, https://mathweb.ucsd.edu/ ∼bdriver/ DRIVER/Book/anal.pdf

    B. Driver , Analysis tools with examples , 2004, https://mathweb.ucsd.edu/ ∼bdriver/ DRIVER/Book/anal.pdf. 28 T. G. ANDERSON, M. BONNET, L. M. FARIA, AND C. P ´EREZ-ARANCIBIA

    2004

  10. [10]

    F. Dubois, Discrete vector potential representation of a divergence-free vector field in three- dimensional domains: numerical analysis of a model problem , SIAM Journal on Numerical Analysis, 27 (1990), pp. 1103–1141

    1990

  11. [11]

    L. M. Faria, C. P ´erez-Arancibia, and M. Bonnet , General-purpose kernel regularization of boundary integral equations via density interpolation , Computer Methods in Applied Mechanics and Engineering, 378 (2021), p. 113703

    2021

  12. [12]

    Gimbutas

    Z. Gimbutas. Private communication, available at https://github.com/IntegralEquations/Inti. jl/blob/main/src/quad rules tables tetrahedron.jl, 2025

    2025

  13. [13]

    G. C. Hsiao and W. L. Wendland, Boundary integral equations (2nd edition), Springer, 2021

    2021

  14. [14]

    Isaacson and H

    E. Isaacson and H. B. Keller , Analysis of numerical methods , New York: Wiley, (1966)

    1966

  15. [15]

    Ishizaka, Exact-curved lagrange finite elements for the Poisson problem in two dimensions , arXiv preprint arXiv:2605.26566, (2026)

    H. Ishizaka, Exact-curved lagrange finite elements for the Poisson problem in two dimensions , arXiv preprint arXiv:2605.26566, (2026)

    Pith/arXiv arXiv 2026

  16. [16]

    O. D. Kellogg, Fundations of potential theory , Springer-Verlag, Berlin, Heidelberg, 1967

    1967

  17. [17]

    Lenoir, Optimal isoparametric finite elements and error estimates for domains involving curved boundaries, SIAM Journal on Numerical Analysis, 23 (1986), pp

    M. Lenoir, Optimal isoparametric finite elements and error estimates for domains involving curved boundaries, SIAM Journal on Numerical Analysis, 23 (1986), pp. 562–580

    1986

  18. [18]

    Malhotra and G

    D. Malhotra and G. Biros , PVFMM: A parallel kernel independent FMM for particle and volume potentials, Communications in Computational Physics, 18 (2015), pp. 808–830

    2015

  19. [19]

    Malhotra, A

    D. Malhotra, A. Gholami, and G. Biros , A volume integral equation Stokes solver for problems with variable coefficients , in SC’14: Proceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis, IEEE, 2014, pp. 92– 102

    2014

  20. [20]

    McLean, Strongly elliptic systems and boundary integral equations , Cambridge, 2000

    W. McLean, Strongly elliptic systems and boundary integral equations , Cambridge, 2000

    2000

  21. [21]

    Mohyaddin and J

    S. Mohyaddin and J. Tausch, A fast method for evaluating volume potentials in the Galerkin boundary element method, SIAM J. Sci. Comput., 45 (2023), pp. A480–A501

    2023

  22. [22]

    Mura, Micromechanics of Defects in Solids , Martinus Nijhoff, 1987

    T. Mura, Micromechanics of Defects in Solids , Martinus Nijhoff, 1987

    1987

  23. [23]

    Nicolaides, On a class of finite elements generated by Lagrange interpolation , SIAM Jour- nal on Numerical Analysis, 9 (1972), pp

    R. Nicolaides, On a class of finite elements generated by Lagrange interpolation , SIAM Jour- nal on Numerical Analysis, 9 (1972), pp. 435–445

    1972

  24. [24]

    G. Of, O. Steinbach, and P. Urthaler , Fast evaluation of volume potentials in boundary element methods, SIAM Journal on Scientific Computing, 32 (2010), pp. 585–602

    2010

  25. [25]

    S. A. Sauter and C. Schwab , Boundary Element Methods , Springer, 2010

    2010

  26. [26]

    E. M. Stein , Singular Integrals and Differentiability Properties of Functions , Princeton Uni- versity Press, 1970

    1970

  27. [27]

    Taylor, Partial Differential Equations I: Basic Theory , vol

    M. Taylor, Partial Differential Equations I: Basic Theory , vol. 116, Springer Science & Busi- ness Media, 2013

    2013

  28. [28]

    Vioreanu and V

    B. Vioreanu and V. Rokhlin , Spectra of multiplication operators as a numerical tool , SIAM Journal on Scientific Computing, 36 (2014), pp. A267–A288

    2014

  29. [29]

    Willis, Bounds and self-consistent estimates for the overall properties of anisotropic com- posites, J

    J. Willis, Bounds and self-consistent estimates for the overall properties of anisotropic com- posites, J. Mech. Phys. Solids, 25 (1977), pp. 185–202

    1977

  30. [30]

    Zl ´amal, Curved elements in the finite element method

    M. Zl ´amal, Curved elements in the finite element method. I , SIAM Journal on Numerical Analysis, 10 (1973), pp. 229–240

    1973

  31. [31]

    Zl ´amal, Curved elements in the finite element method

    M. Zl ´amal, Curved elements in the finite element method. II , SIAM Journal on Numerical Analysis, 11 (1974), pp. 347–362

    1974