The quadrature error associated with a regular quadrature rule for evaluation of a layer potential increases rapidly when the evaluation point approaches the surface and the integral becomes nearly singular. Error estimates are needed to determine when the accuracy is insufficient and a more costly special quadrature method should be utilized. The final result of this paper are such quadrature error estimates for the composite Gauss-Legendre rule and the global trapezoidal rule, when applied to evaluate layer potentials defined over smooth curved surfaces in R^3. The estimates have no unknown coefficients and can be efficiently evaluated given the discretization of the surface, invoking a local one-dimensional root-finding procedure. They are derived starting with integrals over curves, using complex analysis involving contour integrals, residue calculus and branch cuts. By complexifying the parameter plane, the theory can be used to derive estimates also for curves in in R^3. These results are then used in the derivation of the estimates for integrals over surfaces. In this procedure, we also obtain error estimates for layer potentials evaluated over curves in R^2. Such estimates combined with a local root-finding procedure for their evaluation were earlier derived for the composite Gauss-Legendre rule for layer potentials written on complex form [4]. This is here extended to provide quadrature error estimates for both complex and real formulations of layer potentials, both for the Gauss-Legendre and the trapezoidal rule. Numerical examples are given to illustrate the performance of the quadrature error estimates. The estimates for integration over curves are in many cases remarkably precise, and the estimates for curved surfaces in R^3 are also sufficiently precise, with sufficiently low computational cost, to be practically useful.
翻译:当评价点接近表层和整体体时,与评估层潜力的定期二次曲线规则有关的二次误差会迅速增加,当评价点接近表面和整体体时,估计数会迅速增加。需要错误估计,以确定准确性不足和应使用更昂贵的特殊二次曲线方法时,本文件的最终结果是复合高斯-Legendre规则和全球捕猎规则的二次误差估计数,用于评价平滑曲线表面的平滑曲线表面所定义的层潜力时。 估计数没有未知系数,而且鉴于表层的离散,可以有效评估这些估计数,并援引当地一维根根根调查程序。这些估计数是从曲线的内含整体值开始,使用涉及等离子内、残余的计算和分支切割的复杂分析方法。如果将参数平面规则规则的复合图层的估算也用这些理论来得出估计数,那么,为了对地表表面的精度和层的精确性能的精确性能的估算,这种精确性能的估算是用来对地表层进行精确性推的。