Least-squares polynomial approximation on weakly admissible meshes: Disk and triangle

We construct symmetric polar WAMs (weakly admissible meshes) with low cardinality for least-squares polynomial approximation on the disk. These are then mapped to an arbitrary triangle. Numerical tests show that the growth of the least-squares projection uniform norm is much slower than the theoretical bound, and even slower than that of the Lebesgue constant of the best known interpolation points for the triangle. As opposed to good interpolation points, such meshes are straightforward to compute for any degree. The construction can be extended to polygons by triangulation.     

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).     

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.     

A note on the discrete approximation of discontinuous curves and surfaces

This is a note on the paper [A. Kouibia, A.J. López-Linares, M. Pasadas, Approximation of discontinuous curves and surfaces with tangent conditions, J. Comput. Appl. Math. 193 (2006) 51–64]. We consider the constructing problem of a discontinuous parametric curve or surface from a finite set of points and tangent conditions. We develop a method based on the theory of discrete smoothing variational splines conveniently adapted to introduce the tangent conditions and the discontinuity set. We give a convergence result and we analyze some numerical and graphical examples in order to illustrate the effectiveness of the presented method.     

The Superconvergence of the Composite Trapezoidal Rule for Hadamard Finite Part Integrals

The composite trapezoidal rule has been well studied and widely applied for numerical integrations and numerical solution of integral equations with smooth or weakly singular kernels. However, this quadrature rule has been less employed for Hadamard finite part integrals due to the fact that its global convergence rate for Hadamard finite part integrals with (p+1)-order singularity is p-order lower than that for the Riemann integrals in general. In this paper, we study the superconvergence of the composite trapezoidal rule for Hadamard finite part integrals with the second-order and the third-order singularity, respectively. We obtain superconvergence estimates at some special points and prove the uniqueness of the superconvergence points. Numerical experiments confirm our theoretical analysis and show that the composite trapezoidal rule is efficient for Hadamard finite part integrals by noting the superconvergence phenomenon. The work of this author was partially supported by the National Natural Science Foundation of China(No.10271019), a grant from the Research Grants Council of the Hong Kong Special Administractive Region, China (Project No. City 102204) and a grant from the Laboratory of Computational Physics The work of this author was supported in part by a grant from the Research Grants Council of the Hong Kong Special Administrative Region, China (Project No. CityU 102204).     

The superconvergence of the Newton-Cotes rule for Cauchy principal value integrals

We consider the general (composite) Newton-Cotes method for the computation of Cauchy principal value integrals and focus on its pointwise superconvergence phenomenon, which means that the rate of convergence of the Newton-Cotes quadrature rule is higher than what is globally possible when the singular point coincides with some a priori known point. The necessary and sufficient conditions satisfied by the superconvergence point are given. Moreover, the superconvergence estimate is obtained and the properties of the superconvergence points are investigated. Finally, some numerical examples are provided to validate the theoretical results.     

Eigenfrequencies of fractal drums

A method for the computation of eigenfrequencies and eigenmodes of fractal drums is presented. The approach involves first conformally mapping the unit disk to a polygon approximating the fractal and then solving a weighted eigenvalue problem on the unit disk by a spectral collocation method. The numerical computation of the complicated conformal mapping was made feasible by the use of the fast multipole method as described in [L. Banjai, L.N. Trefethen, A multipole method for Schwarz–Christoffel mapping of polygons with thousands of sides, SIAM J. Sci. Comput. 25(3) (2003) 1042–1065]. The linear system arising from the spectral discretization is large and dense. To circumvent this problem we devise a fast method for the inversion of such a system. Consequently, the eigenvalue problem is solved iteratively. We obtain eight digits for the first eigenvalue of the Koch snowflake and at least five digits for eigenvalues up to the 20th. Numerical results for two more fractals are shown.     

Algorithms for spline wavelet packets on an interval

The aim of this paper is to describe decomposition and reconstruction algorithms for spline wavelet packets on a closed interval. In order to generate packet spaces of dyadic dimensions, it is necessary to modify the approach for spline wavelets on an interval as studied by Chui, Quak and Weyrich in [3, 11]. The first author was supported by the Department of the Air Force, contract F33600-94-M-2603, and the second author by the Department of Defense, contract H98230-R5-93-9187.     

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.     

Quadrature formulas descending from BS Hermite spline quasi-interpolation

Two new classes of quadrature formulas associated to the BS Boundary Value Methods are discussed. The first is of Lagrange type and is obtained by directly applying the BS methods to the integration problem formulated as a (special) Cauchy problem. The second descends from the related BS Hermite quasi-interpolation approach which produces a spline approximant from Hermite data assigned on meshes with general distributions. The second class formulas is also combined with suitable finite difference approximations of the necessary derivative values in order to define corresponding Lagrange type formulas with the same accuracy.     

Jacobi matrices for sums of weight functions

Orthogonal polynomials are conveniently represented by the tridiagonal Jacobi matrix of coefficients of the recurrence relation which they satisfy. LetJ ₁ andJ ₂ be finite Jacobi matrices for the weight functionsw ₁ andw ₂, resp. Is it possible to determine a Jacobi matrix \(\tilde J\) , corresponding to the weight functions \(\tilde w\) =w ₁+w ₂ using onlyJ ₁ andJ ₂ and if so, what can be said about its dimension? Thus, it is important to clarify the connection between a finite Jacobi matrix and its corresponding weight function(s). This leads to the need for stable numerical processes that evaluate such matrices. Three newO(n ²) methods are derived that “merge” Jacobi matrices directly without using any information about the corresponding weight functions. The first can be implemented using any of the updating techniques developed earlier by the authors. The second new method, based on rotations, is the most stable. The third new method is closely related to the modified Chebyshev algorithm and, although it is the most economical of the three, suffers from instability for certain kinds of data. The concepts and the methods are illustrated by small numerical examples, the algorithms are outlined and the results of numerical tests are reported.     

Fast and accurate interpolation of large scattered data sets on the sphere

In this paper a new efficient algorithm for spherical interpolation of large scattered data sets is presented. The solution method is local and involves a modified spherical Shepard’s interpolant, which uses zonal basis functions as local approximants. The associated algorithm is implemented and optimized by applying a nearest neighbour searching procedure on the sphere. Specifically, this technique is mainly based on the partition of the sphere in a suitable number of spherical zones, the construction of spherical caps as local neighbourhoods for each node, and finally the employment of a spherical zone searching procedure. Computational cost and storage requirements of the spherical algorithm are analyzed. Moreover, several numerical results show the good accuracy of the method and the high efficiency of the proposed algorithm.     

Two shape preserving lagrangeC 2-interpolants

Summary We present a LagrangeC ²-interpolant to scattered convex data which preserves convexity. We also present a LagrangeC ²-interpolant to uniformly spaced monotone data sites which preserves monotonicity. In both cases no further conditions are required on the data values. These interpolants are explicitely described and local. Error isO(h ³) when the function to be interpolated isC ³.     

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.     

The error norm of Gauss-Lobatto quadrature formulae for weight functions of Bernstein-Szegö type

Summary In certain spaces of analytic functions the error term of the Gauss-Lobatto quadrature formula relative to a (nonnegative) weight function is a continuous linear functional. Here we compute the norm of the error functional for the Bernstein-Szegö weight functions consisting of any of the four Chebyshev weights divided by an arbitrary quadratic polynomial that remains positive on [–1, 1]. The norm can subsequently be used to derive bounds for the error functional. The efficiency of these bounds is illustrated with some numerical examples.Work supported in part by a grant from the Research Council of the Graduate School, University of Missouri-Columbia.     

Radial basis functions under tension

We discuss multivariate interpolation with some radial basis function, called radial basis function under tension (RBFT). The RBFT depends on a positive parameter which provides a convenient way of controlling the behavior of the interpolating surface. We show that our RBFT is conditionally positive definite of order at least one and give a construction of the native space, namely a semi-Hilbert space with a semi-norm, minimized by such an interpolant. Error estimates are given in terms of this semi-norm and numerical examples illustrate the behavior of interpolating surfaces.     

Fully symmetric integration rules for the 4-cube

In this paper we discuss fully symmetric integration rules of degree 7 and 9 for the 4-cube. In particular we are interested in good rules. (i.e. rules with all the evaluation points inside the cube and all the weights positive).This work was supported by the Norwegian Research Council for Sciences and Humanities.     

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.     

An IMT-type quadrature formula with the same asymptotic performance as the DE formula

We propose an IMT-type quadrature formula which achieves the same asymptotic error estimate as the DE formula. The point of the idea is to optimize the parameters of the IMT-type transformation depending on the number of sampling points. We also show the performance of our IMT-type quadrature formula by numerical examples.