Abstract
Asymptotic expansions related to the integration of well-behaved functions on simplices and cubes have been known for several decades. Extensions of these results to classes of vertex and line singularities are also known. The nature of these expansions justifies extrapolation using the ε-algorithm of Wynn. In principle this requires a uniform subdivision of the region. This was implemented in QUADPACK for finite intervals and in TRIEX for triangles about 15 years ago. In this paper its incorporation in CUBPACK, a software package for automatic integration over a collection of cubes and simplices, is described and some results are reported. We also report on a special subdivision strategy that offers an alternative approach for higher-dimensional problems.
Similar content being viewed by others
References
J. Berntsen and T.O. Espelid, Error estimation in automatic quadrature routines, ACM Trans. Math. Software 17 (1991) 233–252.
J. Berntsen, T.O. Espelid and A. Genz, An adaptive algorithm for the approximate calculation of multiple integrals, ACM Trans. Math. Software 17 (1991) 437–451.
J. Berntsen, T.O. Espelid and A. Genz, Algorithm 698: DCUHRE – an adaptive multidimensional integration routine for a vector of integrals, ACM Trans. Math. Software 17 (1991) 452–456.
R. Cools, Monomial cubature rules since “Stroud”: A compilation – part 2, J. Comput. Appl. Math. 112(1/2) (1999) 21–27.
R. Cools and A. Haegemans, CUBPACK: Progress report, in: Numerical Integration – Recent Developments, Software and Applications, eds. T.O. Espelid and A. Genz, NATO Advanced Science Institutes Series C: Mathematical and Physical Sciences, Vol. 357 (Kluwer Academic, Dordrecht, 1992) pp. 305–315.
R. Cools and A. Haegemans, CUBPACK: A package for automatic cubature; framework description, ACM Trans. Math. Software 29(3) (2003) 287–296; source code available from http://www.cs. kuleuven.ac.be/~nines/research/CUBPACK/.
R. Cools and K.J. Kim, Rotation invariant cubature formulas over the n-dimensional unit cube, J. Comput. Appl. Math. 132 (2001) 15–32.
R. Cools and B. Maerten, A hybrid subdivision strategy for adaptive integration routines, J. Universal Comput. Sci. 4(5) (1998) 485–499.
R. Cools and P. Rabinowitz, Monomial cubature rules since 'stroud': A compilation, J. Comput. Appl. Math. 48 (1993) 309–326.
E. de Doncker and I. Robinson, An algorithm for automatic integration over a triangle using nonlinear extrapolation, ACM Trans. Math. Software 10 (1984) 1–16.
E. de Doncker and I. Robinson, TRIEX: Integration over a TRIangle using nonlinear EXtrapolation, ACM Trans. Math. Software 10 (1984) 17–22.
T.O. Espelid, On the construction of good fully symmetric integration rules, SIAM J. Numer. Anal. 24 (1987) 855–881.
T.O. Espelid and A. Genz, DECUHR: An algorithm for automatic integration of singular functions over a hyperrectangular region, Numer. Algorithms 8 (1994) 201–220.
A. Genz, A package for testing multiple integration subroutines, in: Numerical Integration – Recent Developments, Software and Applications, eds. P. Keast and G. Fairweather, NATO Advanced Science Institutes Series C: Mathematical and Physical Sciences, Vol. 203 (D. Reidel, Dordrecht, 1987) pp. 337–340.
A. Genz, An adaptive numerical integration algorithm for simplices, in: Computing in the 90s – Proc. of the 1st Great Lakes Computer Science Conf., eds. N.A. Sherwani, E. de Doncker and J.A. Kapenga, Lecture Notes in Computer Science, Vol. 507 (Springer, New York, 1991) pp. 279–292.
A.C. Genz and A.A. Malik, An adaptive algorithm for numerical integration over an N-dimensional rectangular region, J. Comput. Appl. Math. 6 (1980) 295–302.
A. Haegemans, Algorithm 34: An algorithm for the automatic integration over a triangle, Computing 19 (1977) 179–187.
J.N. Lyness, An error functional expansion for N-dimensional quadrature with an integrand function singular at a point, Math. Comp. 30 (1976) 1–23.
J.N. Lyness and E. de Doncker-Kapenga, On quadrature error expansions, part I, J. Comput. Appl. Math. 17 (1987) 131–149.
J.N. Lyness and E. de Doncker, Quadrature error expansions. II. The full corner singularity, Numer. Math. 64 (1993) 355–370.
R. Piessens, E. de Doncker-Kapenga, C.W. Ñberhuber and D.K. Kahaner, QUADPACK (Springer, Berlin, 1983).
A. Sidi, Euler–Maclaurin expansions for integrals over triangles and squares of functions having algebraic/logarithmic singularities along an edge, J. Approx. Theory 39 (1983) 39–53.
P. Van Dooren and L. De Ridder, An adaptive algorithm for numerical integration over an n-dimensional cube, J. Comput. Appl. Math. 2 (1976) 207–217.
P. Verlinden and A. Haegemans, An error expansion for cubature with an integrand with homogeneous boundary singularities, Numer. Math. 65 (1993) 383–406.
P. Wynn, On a device for computing the em(sn) transformation, Math. Comp. 10 (1956) 91–96.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Cools, R. Extrapolation and Adaptivity in Software for Automatic Numerical Integration on a Cube. Numerical Algorithms 34, 259–269 (2003). https://doi.org/10.1023/B:NUMA.0000005365.15663.94
Issue Date:
DOI: https://doi.org/10.1023/B:NUMA.0000005365.15663.94