On the computation of Gaussian quadrature rules for Chebyshev sets of linearly independent functions

DOI10.1137/21M1456935zbMATH Open1507.65065arXiv1710.11244WikidataQ114073949 ScholiaQ114073949MaRDI QIDQ6293265FDOQ6293265

Daan Huybrechs

Publication date: 30 October 2017

Abstract: We consider the computation of quadrature rules that are exact for a Chebyshev set of linearly independent functions on an interval $[a,b]$. A general theory of Chebyshev sets guarantees the existence of rules with a Gaussian property, in the sense that $2l$ basis functions can be integrated exactly with just $l$ points and weights. Moreover, all weights are positive and the points lie inside the interval $[a,b]$. However, the points are not the roots of an orthogonal polynomial or any other known special function as in the case of regular Gaussian quadrature. The rules are characterized by a nonlinear system of equations, and earlier numerical methods have mostly focused on finding suitable starting values for a Newton iteration to solve this system. In this paper we describe an alternative scheme that is robust and generally applicable for so-called complete Chebyshev sets. These are ordered Chebyshev sets where the first $k$ elements also form a Chebyshev set for each $k$. The points of the quadrature rule are computed one by one, increasing exactness of the rule in each step. Each step reduces to finding the unique root of a univariate and monotonic function. As such, the scheme of this paper is guaranteed to succeed. The quadrature rules are of interest for integrals with non-smooth integrands that are not well approximated by polynomials.

Mathematics Subject Classification ID

Numerical quadrature and cubature formulas (65D32) Approximate quadratures (41A55)

This page was built for publication: On the computation of Gaussian quadrature rules for Chebyshev sets of linearly independent functions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6293265)