skip to main content
10.1145/1073884.1073898acmconferencesArticle/Chapter ViewAbstractPublication PagesissacConference Proceedingsconference-collections
Article

Hybrid symbolic-numeric integration in multiple dimensions via tensor-product series

Published: 24 July 2005 Publication History

Abstract

We present a new hybrid symbolic-numeric method for the fast and accurate evaluation of definite integrals in multiple dimensions. This method is well-suited for two classes of problems: (1) analytic integrands over general regions in two dimensions, and (2) families of analytic integrands with special algebraic structure over hyperrectangular regions in higher dimensions.The algebraic theory of multivariate interpolation via natural tensor product series was developed in the doctoral thesis by Chapman, who named this broad new scheme of bilinear series expansions "Geddes series" in honour of his thesis supervisor. This paper describes an efficient adaptive algorithm for generating bilinear series of Geddes-Newton type and explores applications of this algorithm to multiple integration. We will present test results demonstrating that our new adaptive integration algorithm is effective both in high dimensions and with high accuracy. For example, our Maple implementation of the algorithm has successfully computed nontrivial integrals with hundreds of dimensions to 10-digit accuracy, each in under 3 minutes on a desktop computer.Current numerical multiple integration methods either become very slow or yield only low accuracy in high dimensions, due to the necessity to sample the integrand at a very large number of points. Our approach overcomes this difficulty by using a Geddes-Newton series with a modest number of terms to construct an accurate tensor-product approximation of the integrand. The partial separation of variables achieved in this way reduces the original integral to a manageable bilinear combination of integrals of essentially half the original dimension. We continue halving the dimensions recursively until obtaining one-dimensional integrals, which are then computed by standard numeric or symbolic techniques.

References

[1]
J. Bernsten, T. O. Espelid, and A. C. Genz. Algorithm 698: DCUHRE - An Adaptive Multidimensional Integration Routine for a Vector of Integrals. ACM Transactions on Mathematical Software, 17:452--456, 1991.
[2]
O. A. Carvajal. A New Hybrid Symbolic-Numeric Method for Multiple Integration Based on Tensor-Product Series Approximations. Master's thesis, Univ of Waterloo, Waterloo, ON, Canada, 2004.
[3]
F. W. Chapman. Generalized Orthogonal Series for Natural Tensor Product Interpolation. PhD thesis, Univ of Waterloo, Waterloo, ON, Canada, 2003.
[4]
W. Cheney and W. Light. A Course in Approximation Theory. The Brooks/Cole Series in Advanced Mathematics. Brooks/Cole Publishing Co., Pacific Grove, California, 2000.
[5]
K. O. Geddes. Numerical Integration in a Symbolic Context. In B. W. Char, editor, Proc of SYMSAC'86, pages 185--191, New York, 1986. ACM Press.
[6]
K. O. Geddes and G. J. Fee. Hybrid Symbolic-Numeric Integration in Maple. In P. Wang, editor, Proc of ISAAC'92, pages 36--41, New York, 1992. ACM Press.
[7]
A. C. Genz and A. A. Malik. An Adaptive Algorithm for Numerical Integration over an N-Dimensional Rectangular Region. Journal of Computational and Applied Mathematics, 6:295--302, 1980.
[8]
J. M. Hammersley and D. C. Handscomb. Monte Carlo Methods. Methuen, 1964.
[9]
T. X. He. Dimensionality Reducing Expansion of Multivariate Integration. Birkhaüser, 2001.
[10]
F. J. Hickernell. What Affects Accuracy of Quasi-Monte Carlo Quadrature? In H. Niederreiter and J. Spanier, editors, Monte Carlo and Quasi-Monte Carlo Methods, pages 16--55. Springer-Verlag, Berlin, 2000.
[11]
C. Lemieux and P. L'Ecuyer. On Selection Criteria for Lattice Rules and Other Quasi-Monte Carlo Point Sets. Mathematics and Computers in Simulation, 55(1-3):139--148, 2001.
[12]
I. H. Sloan and S. Joe. Lattice Methods for Multiple Integration. Oxford University Press, 1994.
[13]
A. H. Stroud. Approximate Calculation of Multiple Integrals. Prentice-Hall, 1971.

Cited By

View all
  • (2020)Optimal Truncations for Multivariate Fourier and Chebyshev Series: Mysteries of the Hyperbolic Cross: Part I: Bivariate CaseJournal of Scientific Computing10.1007/s10915-020-01131-182:2Online publication date: 31-Jan-2020
  • (2014)Continuous analogues of matrix factorizationsProceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences10.1098/rspa.2014.0585471:2173(20140585-20140585)Online publication date: 12-Nov-2014
  • (2014)Comparison of Some Reduced Representation ApproximationsReduced Order Methods for Modeling and Computational Reduction10.1007/978-3-319-02090-7_3(67-100)Online publication date: 2014
  • Show More Cited By

Index Terms

  1. Hybrid symbolic-numeric integration in multiple dimensions via tensor-product series

    Recommendations

    Comments

    Information & Contributors

    Information

    Published In

    cover image ACM Conferences
    ISSAC '05: Proceedings of the 2005 international symposium on Symbolic and algebraic computation
    July 2005
    388 pages
    ISBN:1595930957
    DOI:10.1145/1073884
    Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected]

    Sponsors

    Publisher

    Association for Computing Machinery

    New York, NY, United States

    Publication History

    Published: 24 July 2005

    Permissions

    Request permissions for this article.

    Check for updates

    Author Tags

    1. Geddes series scheme
    2. Geddes-Newton series expansions
    3. approximation of functions
    4. bilinear series
    5. deconstruction/approximation/reconstruction technique (DART)
    6. multiple integration
    7. splitting operator
    8. symbolic-numeric algorithms
    9. tensor products

    Qualifiers

    • Article

    Conference

    ISSAC05
    Sponsor:

    Acceptance Rates

    Overall Acceptance Rate 395 of 838 submissions, 47%

    Contributors

    Other Metrics

    Bibliometrics & Citations

    Bibliometrics

    Article Metrics

    • Downloads (Last 12 months)4
    • Downloads (Last 6 weeks)0
    Reflects downloads up to 03 Mar 2025

    Other Metrics

    Citations

    Cited By

    View all
    • (2020)Optimal Truncations for Multivariate Fourier and Chebyshev Series: Mysteries of the Hyperbolic Cross: Part I: Bivariate CaseJournal of Scientific Computing10.1007/s10915-020-01131-182:2Online publication date: 31-Jan-2020
    • (2014)Continuous analogues of matrix factorizationsProceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences10.1098/rspa.2014.0585471:2173(20140585-20140585)Online publication date: 12-Nov-2014
    • (2014)Comparison of Some Reduced Representation ApproximationsReduced Order Methods for Modeling and Computational Reduction10.1007/978-3-319-02090-7_3(67-100)Online publication date: 2014
    • (2011)High-precision numerical integrationJournal of Symbolic Computation10.1016/j.jsc.2010.08.01046:7(741-754)Online publication date: 1-Jul-2011
    • (2011)ForewordJournal of Symbolic Computation10.1016/j.jsc.2010.08.00946:7(735-740)Online publication date: 1-Jul-2011

    View Options

    Login options

    View options

    PDF

    View or Download as a PDF file.

    PDF

    eReader

    View online with eReader.

    eReader

    Figures

    Tables

    Media

    Share

    Share

    Share this Publication link

    Share on social media