Computational aspects of multivariate polynomial interpolation: Indexing the coefficients

An algorithm is derived for generating the information needed to pass efficiently between multi-indices of neighboring degrees, of use in the construction and evaluation of interpolating polynomials and in the construction of good bases for polynomial ideals. This revised version was published online in June 2006 with corrections to the Cover Date.     

On the optimal stability of bases of univariate functions

Summary. This paper is concerned with bases of finite dimensional spaces of univariate continuous functions which are optimally stable for evaluation. The only bases considered are those whose elements have no sign changes. Among these, an optimally stable basis is characterized under the assumption that the set of points where each basis function is nonzero is an interval. A uniqueness result and many examples of such optimally stable bases are also provided. Received May 26, 2000 / Published online August 17, 2001     

Error estimates in W 2,∞-semi-norms for discrete interpolating D 2-splines

We derive error estimates in W^2,∞-semi-norms for multivariate discrete D²-splines that interpolate an unknown function at the vertices of given triangulations. These results are widely based on the construction of approximation operators and linear projectors onto piecewise polynomial spaces having weakly stable local bases.     

Bases for kernel-based spaces

Since it is well-known (De Marchi and Schaback (2001) [4]) that standard bases of kernel translates are badly conditioned while the interpolation itself is not unstable in function space, this paper surveys the choices of other bases. All data-dependent bases turn out to be defined via a factorization of the kernel matrix defined by these data, and a discussion of various matrix factorizations (e.g. Cholesky, QR, SVD) provides a variety of different bases with different properties. Special attention is given to duality, stability, orthogonality, adaptivity, and computational efficiency. The “Newton” basis arising from a pivoted Cholesky factorization turns out to be stable and computationally cheap while being orthonormal in the “native” Hilbert space of the kernel. Efficient adaptive algorithms for calculating the Newton basis along the lines of orthogonal matching pursuit conclude the paper.     

Riesz bounds of Wilson bases generated byB-splines

In this paper, we are concerned with biorthogonal Wilson bases having B-splines as well as powers of sinc functions as window functions. We prove properties of B-splines and exponential Euler splines and use these properties to estimate the Riesz bounds of the Wilson bases.     

Cubic Spline Prewavelets on the Four-Directional Mesh

Dedicated to Professor M. J. D. Powell on the occasion of his sixty-fifth birthday and his retirement. In this paper, we design differentiable, two-dimensional, piecewise polynomial cubic prewavelets of particularly small compact support. They are given in closed form, and provide stable, orthogonal decompositions of L ² (R ² ) . In particular, the splines we use in our prewavelet constructions give rise to stable bases of spline spaces that contain all cubic polynomials, whereas the more familiar box spline constructions cannot reproduce all cubic polynomials, unless resorting to a box spline of higher polynomial degree.     

A tension approach to controlling the shape of cubic spline surfaces on FVS triangulations

We propose a parametric tensioned version of the FVS macro-element to control the shape of the composite surface and remove artificial oscillations, bumps and other undesired behaviour. In particular, this approach is applied to C¹ cubic spline surfaces over a four-directional mesh produced by two-stage scattered data fitting methods.     

On the construction of local quadratic spline quasi-interpolants on bounded rectangular domains

In this paper local bivariate C¹$C^1$ spline quasi-interpolants on a criss-cross triangulation of bounded rectangular domains are considered and a computational procedure for their construction is proposed. Numerical and graphical tests are provided.     

Least supported bases and local linear independence

Summary. We introduce the concept of least supported basis, which is very useful for numerical purposes. We prove that this concept is equivalent to the local linear independence of the basis. For any given locally linearly independent basis we characterize all the bases of the space sharing the same property. Several examples for spline spaces are given. Received December 4, 1992 / Revised version received March 2, 1993     

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.     

On near-best discrete quasi-interpolation on a four-directional mesh

Spline quasi-interpolants are practical and effective approximation operators. In this paper, we construct QIs with optimal approximation orders and small infinity norms called near-best discrete quasi-interpolants which are based on Ω-splines, i.e. B-splines with octagonal supports on the uniform four-directional mesh of the plane. These quasi-interpolants are exact on some space of polynomials and they minimize an upper bound of their infinity norms depending on a finite number of free parameters. We show that this problem has always a solution, in general nonunique. Concrete examples of such quasi-interpolants are given in the last section.     

A shape preserving area true approximation of histogram by rational splines

For a given histogram, we consider an application of a simple rational spline to a shape preserving area true approximation of the histogram. An algorithm for determination of the spline is as easy as one with a quadratic polynomial spline, while the latter does not always preserve the shape of the histogram. Some numerical examples are given at the end of the paper.     

Fourier Series with Linear Splines on an Interval

We construct orthonormal bases of linear splines on a finite interval [a, b] and then we study the Fourier series associated to these orthonormal bases. For continuous functions defined on [a, b], we prove that the associated Fourier series converges pointwisely on (a, b) and also uniformly on [a, b], if it convergences pointwisely at a and b.     

Error estimates in Sobolev spaces for interpolating thin plate splines under tension

This paper discusses _Lp$L_p$-error estimates for interpolation by thin plate spline under tension of a function in the classical Sobolev space on an open bounded set with a Lipschitz-continuous boundary. A property of convergence is also given when the set of interpolating points becomes more and more dense.     

Spline smoothing under constraints on derivatives

An efficient algorithm for computing a smoothing polynomial splines under inequality constraints on derivatives is introduced where both order and breakpoints ofs can be prescribed arbitrarily. By using the B-spline representation ofs, the original semi-infinite constraints are replaced by stronger finite ones, leading to a least squares problem with linear inequality constraints. Then these constraints are transformed into simple box constraints by an appropriate substitution of variables so that efficient standard techniques for solving such problems can be applied. Moreover, the smoothing term commonly used is replaced by a cheaply computable approximation. All matrix transformations are realized by numerically stable Givens rotations, and the band structure of the problem is exploited as far as possible.     

Interpolation by periodic splines with Birkhoff knots

Summary We give a complete characterization of the Hermite interpolation problem by periodic splines with Birkhoff knots. As a dual result we derive the characterization of the Birkhoff interpolation by periodic splines with multiple knots.Sponsored by the Bulgarian Ministry of Education and Science under Contract No. MM-15     

A nodal spline interpolant for the Gregory rule of even order

Summary The Gregory rule is a well-known example in numerical quadrature of a trapezoidal rule with endpoint corrections of a given order. In the literature, the methods of constructing the Gregory rule have, in contrast to Newton-Cotes quadrature,not been based on the integration of an interpolant. In this paper, after first characterizing an even-order Gregory interpolant by means of a generalized Lagrange interpolation operator, we proceed to explicitly construct such an interpolant by employing results from nodal spline interpolation, as established in recent work by the author and C.H. Rohwer. Nonoptimal order error estimates for the Gregory rule of even order are then easily obtained.     

Fair upper bounds for the curvature in univariate convex interpolation

In convex interpolation the curvature of the interpolants should be as small as possible. We attack this problem by treating interpolation subject to bounds on the curvature. In view of the concexity the lower bound is equal to zero while the upper bound is assumed to be piecewise constant. The upper bounds are called fair with respect to a function class if the interpolation problem becomes solvable for all data sets in strictly convex position. We derive fair a priori bounds for classes of quadraticC ¹, cubicC ², and quarticC ³ splines on refined grids.     

Convex interval interpolation with cubic splines

A necessary and sufficient criterion is presented under which the problem of the convex interval interpolation with cubicC ¹-splines has at least one solution. The criterion is given as an algorithm which turns out to be effective.Dedicated to Professor Julius Albrecht on the occasion of his 60th birthday.