Construction of spherical spline quasi-interpolants based on blossoming

A general theory of quasi-interpolants based on quadratic spherical Powell-Sabin splines on spherical triangulations of a sphere-like surface S is developed by using polar forms. As application, various families of discrete and differential quasi-interpolants reproducing quadratic spherical Bézier-Bernstein polynomials or the whole space of the spherical Powell-Sabin quadratic splines of class C¹ are presented.     

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.     

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.     

Bivariate spline interpolation at grid points

Summary. We describe algorithms for constructing point sets at which interpolation by spaces of bivariate splines of arbitrary degree and smoothness is possible. The splines are defined on rectangular partitions adding one or two diagonals to each rectangle. The interpolation sets are selected in such a way that the grid points of the partition are contained in these sets, and no large linear systems have to be solved. Our method is to generate a net of line segments and to choose point sets in these segments which satisfy the Schoenberg-Whitney condition for certain univariate spline spaces such that a principle of degree reduction can be applied. In order to include the grid points in the interpolation sets, we give a sufficient Schoenberg-Whitney type condition for interpolation by bivariate splines supported in certain cones. This approach is completely different from the known interpolation methods for bivariate splines of degree at most three. Our method is illustrated by some numerical examples. Received October 5, 1992 / Revised version received May 13, 1994     

Gregory type quadrature based on quadratic nodal spline interpolation

Summary. Using a method based on quadratic nodal spline interpolation, we define a quadrature rule with respect to arbitrary nodes, and which in the case of uniformly spaced nodes corresponds to the Gregory rule of order two, i.e. the Lacroix rule, which is an important example of a trapezoidal rule with endpoint corrections. The resulting weights are explicitly calculated, and Peano kernel techniques are then employed to establish error bounds in which the associated error constants are shown to grow at most linearly with respect to the mesh ratio parameter. Specializing these error estimates to the case of uniform nodes, we deduce non-optimal order error constants for the Lacroix rule, which are significantly smaller than those calculated by cruder methods in previous work, and which are shown here to compare favourably with the corresponding error constants for the Simpson rule. Received July 27, 1998/ Revised version received February 22, 1999 / Published online January 27, 2000     

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

Barycentric formulae for cardinal (SINC-)interpolants

Summary We present a barycentric representation of cardinal interpolants, as well as a weighted barycentric formula for their efficient evaluation. We also propose a rational cardinal function which in some cases agrees with the corresponding cardinal interpolant and, in other cases, is even more accurate.In numerical examples, we compare the relative accuracy of those various interpolants with one another and with a rational interpolant proposed in former work.Dedicated to the memory of Peter HenriciThis work was done at the University of California at San Diego, La Jolla     

Exact bounds for the uniform approximation by spline interpolants

Summary. This paper completes a result of Reimer (1984) concerning -th-degree cardinal and -periodic interpolation. The method of proof is not restricted to the case of and being odd and seems to be more elementary. Received February 1, 1993 / Revised version received September 14, 1993     

A recursive method for computing interpolants

In this paper, we describe a recursive method for computing interpolants defined in a space spanned by a finite number of continuous functions in R^d${\mathbb{R}}^d$. We apply this method to construct several interpolants such as spline interpolants, tensor product interpolants and multivariate polynomial interpolants. We also give a simple algorithm for solving a multivariate polynomial interpolation problem and constructing the minimal interpolation space for a given finite set of interpolation points.     

Computing near-best fixed pole rational interpolants

We study rational interpolation formulas on the interval [−1,1] for a given set of real or complex conjugate poles outside this interval. Interpolation points which are near-best in a Chebyshev sense were derived in earlier work. The present paper discusses several computation aspects of the interpolation points and the corresponding interpolants. We also study a related set of points (that includes the end points), which is more suitable for applications in rational spectral methods. Some examples are given at the end of this paper.     

An error term and uniqueness for Hermite–Birkhoff interpolation involving only function values and/or first derivative values

This paper discusses various aspects of Hermite–Birkhoff interpolation that involve prescribed values of a function and/or its first derivative. An algorithm is given that finds the unique polynomial satisfying the given conditions if it exists. A mean value type error term is developed which illustrates the ill-conditioning present when trying to find a solution to a problem that is close to a problem that does not have a unique solution. The interpolants we consider and the associated error term may be useful in the development of continuous approximations for ordinary differential equations that allow asymptotically correct defect control. Expressions in the algorithm are also useful in determining whether certain specific types of problems have unique solutions. This is useful, for example, in strategies involving approximations to solutions of boundary value problems by collocation.     

Thinning algorithms for scattered data interpolation

Multistep interpolation of scattered data by compactly supported radial basis functions requires hierarchical subsets of the data. This paper analyzes thinning algorithms for generating evenly distributed subsets of scattered data in a given domain in ℝ ^d .     

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.     

Natural bicubic spline fractal interpolation

Fractal Interpolation functions provide natural deterministic approximation of complex phenomena. Cardinal cubic splines are developed through moments (i.e. second derivative of the original function at mesh points). Using tensor product, bicubic spline fractal interpolants are constructed that successfully generalize classical natural bicubic splines. An upper bound of the difference between the natural cubic spline blended fractal interpolant and the original function is deduced. In addition, the convergence of natural bicubic fractal interpolation functions towards the original function providing the data is studied.     

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.     

A family of Hermite interpolants by bisection algorithms

A two point subdivision scheme with two parameters is proposed to draw curves corresponding to functions that satisfy Hermite conditions on [a, b]. We build two functionsf andf ¹ on dyadic numbers and for some values of the parameters,f is in ¹ withf ¹=f. Examples are provided which show how different the curves can be.     

A note on the optimal stability of bases of univariate functions

This note is concerned with the characterizations and uniqueness of bases of finite dimensional spaces of univariate continuous functions which are optimally stable for evaluation with respect to bases whose elements have no sign changes.     

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.