A new class of quadrature rules for estimating the error in Gauss quadrature
Aleksandar V. Pejčev,Lothar Reichel,Miodrag M. Spalević,Stefan M. Spalević
DOI: https://doi.org/10.1016/j.apnum.2024.06.011
IF: 2.994
2024-06-21
Applied Numerical Mathematics
Abstract:The need to evaluate Gauss quadrature rules arises in many applications in science and engineering. It often is important to be able to estimate the quadrature error when applying an l -point Gauss rule, Gl(f) , where f is an integrand of interest. Such an estimate often is furnished by applying another quadrature rule, Qk(f) , with k>l nodes, and using the difference Qk(f)−Gl(f) or its magnitude as an estimate for the quadrature error in Gl(f) or its magnitude. The classical approach to estimate the error in Gl(f) is to let Qk(f) , with k=2l+1 , be the Gauss-Kronrod quadrature rule associated with Gl(f) . However, it is well known that the Gauss-Kronrod rule associated with a Gauss rule Gl(f) might not exist for certain measures that determine the Gauss rule and for certain numbers of nodes. This prompted M. M. Spalević [1] to develop generalized averaged Gauss rules, Gˆ2l+1 , with 2l+1 nodes for estimating the error in Gl(f) . Similarly as for (2l+1) -node Gauss-Kronrod rules, l nodes of the rule Gˆ2l+1 agree with the nodes of Gl . However, generalized averaged Gauss rules are not internal for some measures. They therefore may not be applicable when the integrand only is defined on the convex hull of the support of the measure. This paper describes a new kind of quadrature rules that may be internal also when generalized averaged quadrature rules are not. The construction of the new quadrature rules is based on theory developed by Peherstorfer [2] . Their application is particularly attractive when the rule Gˆ2l+1 is not internal, the integrand cannot be evaluated at all its nodes, and the integrand is inexpensive to evaluate at the quadrature points. Computed examples that illustrate the performance of the new quadrature rules introduced in this paper are presented.
mathematics, applied