Transfinite interpolation on the medians of a triangle and best L 1-approximation

Let be a triangle in and let be the set of its three medians. We construct interpolants to smooth functions using transfinite (or blending) interpolation on The interpolants are of type f(₁)+g(₂)+h(₃), where (₁,₂,₃) are the barycentric coordinates with respect to the vertices of . Based on an error representation formula, we prove that the interpolant is the unique best L¹-approximant by functions of this type subject the function to be approximated is from a certain convexity cone in C³().Received: 17 December 2003     

Interpolation lattices in several variables

Principal lattices are classical simplicial configurations of nodes suitable for multivariate polynomial interpolation in n dimensions. A principal lattice can be described as the set of intersection points of n + 1 pencils of parallel hyperplanes. Using a projective point of view, Lee and Phillips extended this situation to n + 1 linear pencils of hyperplanes. In two recent papers, two of us have introduced generalized principal lattices in the plane using cubic pencils. In this paper we analyze the problem in n dimensions, considering polynomial, exponential and trigonometric pencils, which can be combined in different ways to obtain generalized principal lattices.We also consider the case of coincident pencils. An error formula for generalized principal lattices is discussed. Partially supported by the Spanish Research Grant BFM2003-03510, by Gobierno de Aragón and Fondo Social Europeo.     

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.     

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.     

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     

On Lagrange and Hermite interpolation in R k

Summary A method for the construction of a set of data of interpolation in several variables is given. The resulting data, which are either function values or directional derivatives values, give rise to a space of polynomials, in such a way that unisolvence is guaranteed. The interpolating polynomial is calculated using a procedure which generalizes the Newton divided differences formula for a single variable.     

Lattices on simplicial partitions

In this paper, (d+1)-pencil lattices on simplicial partitions in R^d are studied. The barycentric approach naturally extends the lattice from a simplex to a simplicial partition, providing a continuous piecewise polynomial interpolant over the extended lattice. The number of degrees of freedom is equal to the number of vertices of the simplicial partition. The constructive proof of this fact leads to an efficient computer algorithm for the design of a lattice.     

Lattices on simplicial partitions which are not simply connected

In this paper, (d+1)-pencil lattices on simplicial partitions in R^d, which are not simply connected, are studied. It is shown, how the fact that a partition is not simply connected can be used to increase the flexibility of a lattice. A local modification algorithm is developed also to deal with slight partition topology changes that may appear afterwards a lattice has already been constructed.     

Multivariate polynomial interpolation under projectivities II: Neville-Aitken formulas

This is the second part of a note on interpolation by real polynomials of several real variables. For certain regular knot systems (geometric or regular meshes, tensor product grids), Neville-Aitken algorithms are derived explicitly. By application of a projectivity they can be extended in a simple way to arbitrary (k+1)-pencil lattices as recently introduced by Lee and Phillips. A numerical example is given.Partially supported by CICYT Res. Grant PS87-0060.Travel Grant Programa Europa 1991, CAI Zaragoza.     

On the bivariate Shepard-Lidstone operators

We propose a new combination of the bivariate Shepard operators (Coman and Trîmbi?a?, 2001 [2]) by the three point Lidstone polynomials introduced in Costabile and Dell’Accio (2005) [7]. The new combination inherits both degree of exactness and Lidstone interpolation conditions at each node, which characterize the interpolation polynomial. These new operators find application to the scattered data interpolation problem when supplementary second order derivative data are given (Kraaijpoel and van Leeuwen, 2010 [13]). Numerical comparison with other well known combinations is presented.     

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.     

Bivariate Hermite-Birkhoff interpolation and Vandermonde determinants

The concepts of Vandermonde determinant and confluent Vandermonde determinant are extended to the multidimensional setting by relating them to multivariate interpolation problems. With an approach different from that of other recent papers on this subject, the values of these determinants are computed, recovering and extending the results of those papers.Partially supported by Research Grant PS900121 DGICYT.     

On the least squares fit by radial functions to multidimensional scattered data

Summary This paper investigates some aspects of discrete least squares approximation by translates of certain classes of radial functions. Its specific aims are (i) to provide conditions under which the associated least squares matrix is invertible and (ii) to give upper bounds for the Euclidean norms of the inverses of these matrices (when they exist).The second named author was supported by the National Science Foundation under grant number DMS-8901345     

On a class of weak Tchebycheff systems

In this paper we study the approximation power, the existence of a normalized B-basis and the structure of a degree-raising process for spaces of the formrequiring suitable assumptions on the functions u and v. The results about degree raising are detailed for special spaces of this form which have been recently introduced in the area of CAGD.     

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.     

Transfinite mean value interpolation in general dimension

Mean value interpolation is a simple, fast, linearly precise method of smoothly interpolating a function given on the boundary of a domain. For planar domains, several properties of the interpolant were established in a recent paper by Dyken and the second author, including: sufficient conditions on the boundary to guarantee interpolation for continuous data; a formula for the normal derivative at the boundary; and the construction of a Hermite interpolant when normal derivative data is also available. In this paper we generalize these results to domains in arbitrary dimension.     

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 Power of Standard Information for Weighted Approximation

We study weighted approximation of multivariate functions for classes of standard and linear information in the worst case and average case settings. Under natural assumptions, we show a relation between n th minimal errors for these two classes of information. This relation enables us to infer convergence and error bounds for standard information, as well as the equivalence of tractability and strong tractability for the two classes. April 11, 2001. Final version received: May 29, 2001.