Skip to Main content Skip to Navigation
Journal articles

Algorithms for Computing Cubatures Based on Moment Theory

Mathieu Collowald 1 Evelyne Hubert 2
1 GALAAD2 - Géométrie , Algèbre, Algorithmes
CRISAM - Inria Sophia Antipolis - Méditerranée
2 AROMATH - AlgebRe, geOmetrie, Modelisation et AlgoriTHmes
CRISAM - Inria Sophia Antipolis - Méditerranée , NKUA - National and Kapodistrian University of Athens
Abstract : Quadrature is an approximation of the definite integral of a function by a weighted sum of function values at specified points, or nodes, within the domain of integration. Gaussian quadratures are constructed to yield exact results for any polynomials of degree 2r-1 or less by a suitable choice of r nodes and weights. Cubature is a generalization of quadrature in higher dimension. In this article we elaborate algorithms to compute all minimal cubatures for a given domain and a given degree. We propose first an algorithm in symbolic computation to characterize all cubatures of a given degree with a fixed number of nodes. The determination of the nodes and weights is then left to the computation of the eigenvectors of the matrix identified at the characterization stage and can be performed numerically. The characterisation of cubatures on which our algorithms are based stems from moment theory. We formulate the results there in a basis independent way : Rather than considering the moment matrix, the central object in moment problems, we introduce the underlying linear map from the polynomial ring to its dual, the Hankel operator. This makes natural the use of bases of polynomials other than the monomial basis, and proves to be computationally relevant, either for numerical properties or to exploit symmetry.
Complete list of metadata

Cited literature [77 references]  Display  Hide  Download
Contributor : Evelyne Hubert Connect in order to contact the contributor
Submitted on : Thursday, September 13, 2018 - 10:13:25 AM
Last modification on : Thursday, November 26, 2020 - 3:50:03 PM
Long-term archiving on: : Friday, December 14, 2018 - 1:46:12 PM


Files produced by the author(s)




Mathieu Collowald, Evelyne Hubert. Algorithms for Computing Cubatures Based on Moment Theory. Studies in Applied Mathematics, Wiley-Blackwell, 2018, Orthogonal Polynomials, Special Functions, and Applications, 141 (4), pp.501-546. ⟨10.1111/sapm.12240⟩. ⟨hal-01873077v2⟩