Stability results for scattered‐data interpolation on Euclidean spheres

Let denote the unit sphere in and the geodesic distance in . A spherical‐basis function approximant is a function of the form , where are real constants, is a fixed function, and is a set of distinct points in . It is known that if is a strictly positive definite function in , then the interpolation matrix is positive definite, hence invertible, for every choice of distinct points and every positive integer M. The paper studies a salient subclass of such functions , and provides stability estimates for the associated interpolation matrices. This revised version was published online in June 2006 with corrections to the Cover Date.     

A new non‐polynomial univariate interpolation formula of Hermite type

A new C ^∞ interpolant is presented for the univariate Hermite interpolation problem. It differs from the classical solution in that the interpolant is of non‐polynomial nature. Its basis functions are a set of simple, compact support, transcendental functions. The interpolant can be regarded as a truncated Multipoint Taylor series. It has essential singularities at the sample points, but is well behaved over the real axis and satisfies the given functional data. The interpolant converges to the underlying real‐analytic function when (i) the number of derivatives at each point tends to infinity and the number of sample points remains finite, and when (ii) the spacing between sample points tends to zero and the number of specified derivatives at each sample point remains finite. A comparison is made between the numerical results achieved with the new method and those obtained with polynomial Hermite interpolation. In contrast with the classical polynomial solution, the new interpolant does not suffer from any ill conditioning, so it is always numerically stable. In addition, it is a much more computationally efficient method than the polynomial approach. This revised version was published online in June 2006 with corrections to the Cover Date.     

Gaussian radial‐basis functions: Cardinal interpolation of ℓp and power‐growth data

Suppose λ is a positive number. Basic theory of cardinal interpolation ensures the existence of the Gaussian cardinal function \(L_\lambda (x) = \sum\nolimits_{k \in \mathbb{Z}} {c_k \exp ( - \lambda (x - k)^2 ),x \in \mathbb{R}} ,\) satisfying the interpolatory conditions \(L_\lambda (j) = \delta _{0j} ,j \in \mathbb{Z}.\) The paper considers the Gaussian cardinal interpolation operator

$(\mathcal{L}_\lambda {\text{y}})(x): = \sum\limits_{k \in \mathbb{Z}} {y_k L_\lambda (x - k),{\text{ y}} = (y_k )_{k \in \mathbb{Z}} ,{\text{ }}x \in \mathbb{R}} ,$

as a linear mapping from ℓ^p(ℤ) into L ^p(ℝ), 1≤ p ∞, and in particular, its behaviour as λ→0⁺. It is shown that \(\left\| {\mathcal{L}_\lambda } \right\|_p \) is uniformly bounded (in λ) for 1 < p < ∞, and that \(\left\| {\mathcal{L}_\lambda } \right\|_1 \asymp \log (1/\lambda )\) as λ→0⁺. The limiting behaviour is seen to be that of the classical Whittaker operator

$\mathcal{W}:{\text{y}} \mapsto \sum\limits_{k \in \mathbb{Z}} {y_k \frac{{\sin \pi (x - k)}}{{\pi (x - k)}}} ,$

in that \(\lim _{\lambda \to 0^ + } \left\| {\mathcal{L}_\lambda {\text{y}} - \mathcal{W}{\text{y}}} \right\|_p = 0,\) for every \({\text{y}} \in \ell ^p (\mathbb{Z}){\text{ and }}1 < p < \infty .\) It is further shown that the Gaussian cardinal interpolants to a function f which is the Fourier transform of a tempered distribution supported in (-π,π) converge locally uniformly to f as λ→0⁺. Multidimensional extensions of these results are also discussed.

     

The C r-fundamental splines of Clough–Tocher and Powell–Sabin types for Lagrange interpolation on a three direction mesh

Let Δ⁽¹⁾ be the uniform three direction mesh of the plane whose vertices are integer points of .Let (respectively of degree d=3r (respectively d=3r+1 ) for r odd (respectively even) on the triangulation , and of degree d=2r (respectively d=2r+1) for r odd (respectively even) on the triangulation . Using linear combinations of translates of these splines we obtain Lagrange interpolants whose corresponding order of approximation is optimal. This revised version was published online in June 2006 with corrections to the Cover Date.     

<Emphasis Type="Italic">C</Emphasis><Superscript>1</Superscript> positivity preserving scattered data interpolation using rational Bernstein-Bézier triangular patch

A local C ¹ positivity preserving scheme is developed using Bernstein-Bézier rational cubic function. The domain is triangulated by Delaunay triangulation method. Simple sufficient conditions are derived on the inner and boundary Bézier ordinates to preserve the shape of positive data. These inner and boundary Bézier ordinates involve weights in their definition. In any triangular patch if the Bézier ordinates do not satisfy the derived conditions of positivity, then these are modified by the weights (free parameters) involved in the construction of Bernstein-Bézier rational cubic function to preserve the shape of positive scattered data.