Abstract
A new approach to the integration of vertex singularities is described. This approach is based on a non-uniform subdivision of the region of integration and the technique fits well to the subdivision strategy used in many adaptive algorithms. A nice feature with this approach is that it can be used in any dimension and on any region of integration which can be subdivided into subregions of the same form. The strategy can be applied both to vertex singularities and internal point singularities. In the latter case this can be done without an initial subdivision of the region in order to put the singular point in a vertex. It turns out that the technique has excellent numerical stability properties.
Similar content being viewed by others
References
J. Berntsen and T. O. Espelid.Algorithm 706: DCUTRI: An algorithm for adaptive cubature over a collection of triangles. ACM Trans. Math. Softw., 18(3): 329–342, 1992.
R. Cariño, I. Robinson, and E. de Doncker.An algorithm for automatic integration of certain singular functions over a triangle. In T. O. Espelid and A. Genz, editors,Numerical Integration, Recent Developments, Software and Applications, NATO ASI Series C: Math. and Phys. Sciences Vol. 357, pages 295–304, Dordrecht, The Netherlands, 1992. Kluwer Academic Publishers.
E. de Doncker.New Euler-Maclaurin expansions and their application to quadrature over the s-dimensional simplex. Math. Comp., 33: 1003–1018, 1979.
E. de Doncker.Numerical integration and asymptotic expansions., 1980. Ph.D. thesis, Katholieke Universiteit, Leuven, Berlgium.
J. N. Lyness.Applications of extrapolation techniques to multidimensional quadrature of some integrand functions with a singularity. J. Comp. Phys., 20(3):346–364, 1976.
J. N. Lyness.An error functional expansion for n-dimensional quadrature with an integrand function singular at a point. Math. Comp., 30(133):1–23, 1976.
J. N. Lyness.Quadrature over a simplex: Part 1.A representation for the integrand function.SIAM J. Numer. Anal., 15(1):122–133, 1978.
J. N. Lyness,Quadrature over a simplex: Part 2.A representation for the error functional. SIAM J. Numer. Anal., 15(5):870–887, 1978.
J. N. Lyness.On handling singularities in finite elements. In T. O. Espelid and A. Genz, editors,Numerical Integration, Recent Developments, Software and Applications, NATO ASI Series C: Math. and Phys. Sciences Vol. 357, pages 219–233, Dordrecht, The Netherlands, 1992. Kluwer Academic Publishers.
J. N. Lyness and G. Monegato.Quadrature error functional expansions for the simplex when the integrand function has singularities at vertices. Math. Comp., 34(149):213–225, 1980.
A. Sidi.Euler-McLaurin expansions for integrals over triangles of functions having algebraic/logarithmic singularities along an edge. J. Approx. Th., 39(1):39–53, 1983.
W. Squire.Partition-extrapolation methods for numerical quadratures. Int. J. Comp. Math., 5(1):81–91, 1975.
Author information
Authors and Affiliations
Additional information
Dedicated to Carl-Erik Fröberg on the occasion of his 75th birthday.
This work was supported by The Norwegian Research Council for Science and the Humanities.
Rights and permissions
About this article
Cite this article
Espelid, T.O. On integrating vertex singularities using extrapolation. BIT 34, 62–79 (1994). https://doi.org/10.1007/BF01935016
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01935016