Exit criteria and monotonicity in compound quadrature of Gaussian type

Summary In this short note, for compound quadrature rules of Gaussian type, we prove stopping rules and monotonicity results based on Peano-kernel methods.     

Cubature rule associated with a discrete blending sum of quadratic spline quasi-interpolants

A new cubature rule for a parallelepiped domain is defined by integrating a discrete blending sum of C¹ quadratic spline quasi-interpolants in one and two variables. We give the weights and the nodes of this cubature rule and we study the associated error estimates for smooth functions. We compare our method with cubature rules based on the tensor products of spline quadratures and classical composite Simpson’s rules.     

Superconvergence and ultraconvergence of Newton-Cotes rules for supersingular integrals

In this article, the general (composite) Newton-Cotes rules for evaluating Hadamard finite-part integrals with third-order singularity (which is also called “supersingular integrals”) are investigated and the emphasis is placed on their pointwise superconvergence and ultraconvergence. The main error of the general Newton-Cotes rules is derived, which is shown to be determined by a certain function . Based on the error expansion, the corresponding modified quadrature rules are also proposed. At last, some numerical experiments are carried out to validate the theoretical analysis.     

An evaluation of Clenshaw–Curtis quadrature rule for integration w.r.t. singular measures

This work is devoted to the study of quadrature rules for integration with respect to (w.r.t.) general probability measures with known moments. Automatic calculation of the Clenshaw–Curtis rules is considered and analyzed. It is shown that it is possible to construct these rules in a stable manner for quadrature w.r.t. balanced measures. In order to make a comparison Gauss rules and their stable implementation for integration w.r.t. balanced measures are recalled. Convergence rates are tested in the case of binomial measures.     

Extended product integration rules

In this paper we introduce a class of extended product quadrature rules to associate with the corresponding standard product rules, and present an algorithm for their construction. A general discussion on the convergence of such formulas is then given. Finally some examples and applications are considered.Work sponsored by the Ministero della Pubblica Istruzione of Italy.     

The superconvergence of the composite midpoint rule for the finite-part integral

The composite midpoint rule is probably the simplest one among the Newton-Cotes rules for Riemann integral. However, this rule is divergent in general for Hadamard finite-part integral. In this paper, we turn this rule to a useful one and, apply it to evaluate Hadamard finite-part integral as well as to solve the relevant integral equation. The key point is based on the investigation of its pointwise superconvergence phenomenon, i.e., when the singular point coincides with some a priori known point, the convergence rate of the midpoint rule is higher than what is globally possible. We show that the superconvergence rate of the composite midpoint rule occurs at the midpoint of each subinterval and obtain the corresponding superconvergence error estimate. By applying the midpoint rule to approximate the finite-part integral and by choosing the superconvergence points as the collocation points, we obtain a collocation scheme for solving the finite-part integral equation. More interesting is that the inverse of the coefficient matrix of the resulting linear system has an explicit expression, by which an optimal error estimate is established. Some numerical examples are provided to validate the theoretical analysis.     

Error estimates for Filon quadrature formulae

In this paper, upper bounds for the error of (generalized) Filon quadrature formulae are stated. Furthermore, the main term of this error is derived, yielding simple modified quadrature rules of higher asymptotical precision.     

Romberg quadrature using the Bulirsch sequence

Summary. In this note, we prove a conjecture of Bulirsch concerning the definiteness of the Romberg quadrature rules using the Bulirsch sequence. We compare these rules with the classical Romberg scheme and the Gaussian rules. Received May 16, 2000 / Published online May 30, 2001     

Modified optimal quadrature extensions

Summary Optimal extensions of quadrature rules are of importance in the construction of automatic integrators but many sequences fail to exist in usable form. The paper considers some techniques for overcoming the problem of inextensibility with a minimal effect on integrating efficiency. A modification to the extension procedure proposed recently by Begumisa and Robinson is shown to be just a special case of the standard theory for quadrature extension. Some illustrative examples are included.     

The superconvergence of Newton–Cotes rules for the Hadamard finite-part integral on an interval

We study the general (composite) Newton–Cotes rules for the computation of Hadamard finite-part integral with the second-order singularity and focus on their pointwise superconvergence phenomenon, i.e., when the singular point coincides with some a priori known point, the convergence rate is higher than what is globally possible. We show that the superconvergence rate of the (composite) Newton–Cotes rules occurs at the zeros of a special function and prove the existence of the superconvergence points. Several numerical examples are provided to validate the theoretical analysis. The work of J. Wu was partially supported by the National Natural Science Foundation of China (No. 10671025) and a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (No. CityU 102507). The work of W. Sun was supported in part by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (No. City U 102507) and the National Natural Science Foundation of China (No. 10671077).     

Suboptimal Kronrod extension formulae for numerical quadrature

Summary We consider cases where the Stieltjes polynomial associated with a Gaussian quadrature formula has complex zeros. In such cases a Kronrod extension of the Gaussian rule does not exist. A method is described for modifying the Stieltjes polynomial so that the resulting polynomial has no complex zeros. The modification is performed in such a way that the Kronrod-type extension rule resulting from the addition of then+1 zeros of the modified Stieltjes polynomial to the original knots of the Gaussian rule has only slightly lower degree of precision than normally achieved when the Kronrod extension rule exists. As examples of the use of the method, we present some new formulae extending the classical Gauss-Hermite quadrature rules. We comment on the limited success of the method in extending Gauss-Laguerre rules and show that several modified extensions of the Gauss Gegenbauer formulae exist in cases where the standard Kronrod extension does not.     

The Existence of Good Extensible Polynomial Lattice Rules

Extensible (polynomial) lattice rules have been introduced recently and they are convenient tools for quasi-Monte Carlo integration. It is shown in this paper that for suitable infinite families of polynomial moduli there exist generating parameters for extensible rank-1 polynomial lattice rules such that for all these infinitely many moduli and all dimensions s the quantity R ^(s) and the star discrepancy are small. The case of Korobov-type polynomial lattice rules is also considered.Received April 30, 2002; in revised form August 21, 2002 Published online April 4, 2003     

Error bounds for Romberg quadrature

Denote by the error of a Romberg quadrature rule applied to the function f. We determine approximately the constants in the bounds of the types and for all classical Romberg rules. By a comparison with the corresponding constants of the Gaussian rule we give the statement “The Gaussian quadrature rule is better than the Romberg method” a precise meaning. Received September 10, 1997 / Revised version received February 16, 1998     

From approximating subdivision schemes for exponential splines to high-performance interpolating algorithms

In this work we construct three novel families of approximating subdivision schemes that generate piecewise exponential polynomials and we show how to convert these into interpolating schemes of great interest in curve design for their ability to reproduce important analytical shapes and to provide highly smooth limit curves with a controllable tension.     

Two new formulas for the numerical evaluation of the Hilbert Transform

We develop two algorithms for the numerical evaluation of the semi-infinite Hilbert Transform of functions with a given algebraic behaviour at the origin and at infinity. The first algorithm is connected with a Gauss-Jacobi type quadrature formula for unbounded intervals; the second is based on a rational Bernstein-type operator. Error estimates for different classes of functions are shown. Finally numerical examples are given, comparing the rules among themselves.     

Shifted lattice rules based on a general weighted discrepancy for integrals over Euclidean space

We approximate weighted integrals over Euclidean space by using shifted rank-1 lattice rules with good bounds on the “generalised weighted star discrepancy”. This version of the discrepancy corresponds to the classic L_∞ weighted star discrepancy via a mapping to the unit cube. The weights here are general weights rather than the product weights considered in earlier works on integrals over R^d. Known methods based on an averaging argument are used to show the existence of these lattice rules, while the component-by-component technique is used to construct the generating vector of these shifted lattice rules. We prove that the bound on the weighted star discrepancy considered here is of order O(n^−1+δ) for any δ>0 and with the constant involved independent of the dimension. This convergence rate is better than the O(n^−1/2) achieved so far for both Monte Carlo and quasi-Monte Carlo methods.     

Quadrature rules based on partial fraction expansions

Quadrature rules are typically derived by requiring that all polynomials of a certain degree be integrated exactly. The nonstandard issue discussed here is the requirement that, in addition to the polynomials, the rule also integrates a set of prescribed rational functions exactly. Recurrence formulas for computing such quadrature rules are derived. In addition, Fejér's first rule, which is based on polynomial interpolation at Chebyshev nodes, is extended to integrate also rational functions with pre-assigned poles exactly. Numerical results, showing a favorable comparison with similar rules that have been proposed in the literature, are presented. An error analysis of a representative test problem is given. This revised version was published online in June 2006 with corrections to the Cover Date.     

Onp-generator fully symmetric quadrature rules

Summary We consider fully symmetric quadrature rules for fully symmetricn-dimensional integration regions. When the region is a product region it is well known that product Gaussian rules exist. These obtain an approximation of polynomial degree 4p+1 based on (2p+1) ⁿ function values arranged on a rectangular grid. We term rules using such a grid,p-generator rules. In this paper we determine the necessary conditions on the region of integration forp-generator rules of degree 4p+1 to exist. Regions with this property are termed PropertyQ regions and besides product spaces, this class includes the hypersphere and other related regions.Work performed under the auspices of the U.S. Energy Research and Development Administration     

Approximation by interpolating variational splines

In this paper we present an approximation problem of parametric curves and surfaces from a Lagrange or Hermite data set. In particular, we study an interpolation problem by minimizing some functional on a Sobolev space that produces the new notion of interpolating variational spline. We carefully establish a convergence result. Some specific cases illustrate the generality of this work.     

Optimization of parameters for curve interpolation by cubic splines

In this paper, we present an interpolation method for curves from a data set by means of the optimization of the parameters of a quadratic functional in a space of parametric cubic spline functions. The existence and the uniqueness of this problem are shown. Moreover, a convergence result of the method is established in order to justify the method presented. The aforementioned functional involves some real non-negative parameters; the optimal parametric curve is obtained by the suitable optimization of these parameters. Finally, we analyze some numerical and graphical examples in order to show the efficiency of our method.