Vector cell-average multiresolution based on Hermite interpolation

Harten’s interpolatory multiresolution representation of data has been extended in the case of point-value discretization to include Hermite interpolation by Warming and Beam in [17]. In this work we extend Harten’s framework for multiresolution analysis to the vector case for cell-averaged data, focusing on Hermite interpolatory techniques. *Supported by European Community IHP projects HPRN-CT-2002-00282 and HPRN-CT-2005-00286. **Supported by European Community IHP projects HPRN-CT-2002-00282 and HPRN-CT-2005-00286, and by FPU grant from M.E.C.D. AP2000-1386. ^†Supported by European Community IHP projects HPRN-CT-2002-00282 and HPRN-CT-2005-00286.     

Error estimates of Hermite interpolation

This paper presents a procedure for obtaining error estimates for Hermite interpolation at the Chebyshev nodes {cos ((2j+1)/2n)} _j ^=0n–1 –1x1, for functionsf(x) of various orders of continuity. The procedure is applicable in many cases when the usual Lagrangian error bound is not, and is a better bound, in general, when both are applicable.     

On a Hermite interpolation on the sphere

The purpose of this paper is to put forward a kind of Hermite interpolation scheme on the unit sphere. We prove the superposition interpolation process for Hermite interpolation on the sphere and give some examples of interpolation schemes. The numerical examples shows that this method for Hermite interpolation on the sphere is feasible. And this paper can be regarded as an extension and a development of Lagrange interpolation on the sphere since it includes Lagrange interpolation as a particular case.     

Error estimates for Hermite interpolation on spheres

In this paper, we prove convergence rates for spherical spline Hermite interpolation on the sphere S^d−1 via an error estimate given in a technical report by Luo and Levesley. The functionals in the Hermite interpolation are either point evaluations of pseudodifferential operators or rotational differential operators, the desirable feature of these operators being that they map polynomials to polynomials. Convergence rates for certain derivatives are given in terms of maximum point separation.     

On the Hermite interpolation

Explicit representations for the Hermite interpolation and their derivatives of any order are provided.Furthermore,suppose that the interpolated function f has continuous derivatives of sufficiently high order on some sufficiently small neighborhood of a given point x and any group of nodes are also given on the neighborhood.If the derivatives of any order of the Hermite interpolation polynomial of f at the point x are applied to approximating the corresponding derivatives of the function f(x),the asymptotic representations for the remainder are presented.     

Error inequalities for quintic and biquintic discrete Hermite interpolation

In this paper we shall develop a class of discrete Hermite interpolates in one and two independent variables. Further, we offer explicit error bounds in ?_∞ norm for the quintic and biquintic discrete Hermite interpolates. Some numerical examples are included to illustrate the results obtained.     

Quasi-Newton method by Hermite interpolation

This paper describes a new attempt to solve the problem of computing a local minimizer of a sufficiently often differentiable unconstrained objective function. In every step of the iteration, a special Hermite interpolant is constructed. Old iteration points serve as points of support with the function value and gradient information. This yields a quasi-Newton algorithm with quadratic convergence order.     

A fast algorithm for multivariate Hermite interpolation

Multivariate Hermite interpolation is widely applied in many fields, such as finite element construction, inverse engineering, CAD etc.. For arbitrarily given Hermite interpolation conditions, the typical method is to compute the vanishing ideal I （the set of polynomials satisfying all the homogeneous interpolation conditions are zero） and then use a complete residue system modulo I as the interpolation basis. Thus the interpolation problem can be converted into solving a linear equation system. A generic algorithm was presented in [18], which is a generalization of BM algorithm [22] and the complexity is O（τ^3） where r represents the number of the interpolation conditions. In this paper we derive a method to obtain the residue system directly from the relative position of the points and the corresponding derivative conditions （presented by lower sets） and then use fast GEPP to solve the linear system with O（（τ ＋ 3）τ^2） operations, where τ is the displacement-rank of the coefficient matrix. In the best case τ = 1 and in the worst case τ = [τ/n], where n is the number of variables.     

六次PH曲线C~1 Hermite插值

本文讨论六次PH(pythagorean hodograph)曲线的Hermite插值问题.六次PH曲线可以分为两种类型,本文使用参数曲线的复数表示形式,分别给出这两类曲线的构造方法.在给定C1连续的Hermite条件下,需要指定一个自由参数以确定插值曲线,本文进一步阐述这个自由参数的几何意义.由于六次PH曲线是非正则曲线,对于第一类曲线,不易控制奇异点在曲线中的位置;而对于第二类曲线,奇异点可以在构造过程中显式地被指定,因此可以有效地避免其在特定曲线段上的出现.     

Hermite interpolation of nonsmooth functions preserving boundary conditions

This article is devoted to the construction of a Hermite-type regularization operator transforming functions that are not necessarily into globally finite-element functions that are piecewise polynomials. This regularization operator is a projection, it preserves appropriate first and second order polynomial traces, and it has approximation properties of optimal order. As an illustration, it is used to discretize a nonhomogeneous Navier-Stokes problem, with tangential boundary condition.

     

A wavelets method for plane elasticity problem

A Neumann boundary value problem of plane elasticity problem in the exterior circular domain is reduced into an equivalent natural boundary integral equation and a Poisson integral formula with the DtN method. Using the trigonometric wavelets and Galerkin method, we obtain a fast numerical method for the natural boundary integral equation which has an unique solution in the quotient space. We decompose the stiffness matrix in our numerical method into four circulant and symmetrical or antisymmetrical submatrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT) instead of the inverse matrix. Examples are given for demonstrating our method has good accuracy of our method even though the exact solution is almost singular.     

Cardinal Hermite interpolation using positive definite functions

In this paper, we study cardinal Hermite interpolation by using positive definite functions. Among other things, we establish a procedure that employs the multiquadrics for cardinal Hermite interpolation.     

Hermite interpolation by Pythagorean hodograph curves of degree seven

Polynomial Pythagorean hodograph (PH) curves form a remarkable subclass of polynomial parametric curves; they are distinguished by having a polynomial arc length function and rational offsets (parallel curves). Many related references can be found in the article by Farouki and Neff on Hermite interpolation with PH quintics. We extend the Hermite interpolation scheme by taking additional curvature information at the segment boundaries into account. As a result we obtain a new construction of curvature continuous polynomial PH spline curves. We discuss Hermite interpolation of boundary data (points, first derivatives, and curvatures) with PH curves of degree 7. It is shown that up to eight possible solutions can be found by computing the roots of two quartic polynomials. With the help of the canonical Taylor expansion of planar curves, we analyze the existence and shape of the solutions. More precisely, for Hermite data which are taken from an analytical curve, we study the behaviour of the solutions for decreasing stepsize . It is shown that a regular solution is guaranteed to exist for sufficiently small stepsize , provided that certain technical assumptions are satisfied. Moreover, this solution matches the shape of the original curve; the approximation order is 6. As a consequence, any given curve, which is assumed to be (curvature continuous) and to consist of analytical segments can approximately be converted into polynomial PH form. The latter assumption is automatically satisfied by the standard curve representations of Computer Aided Geometric Design, such as Bézier or B-spline curves. The conversion procedure acts locally, without any need for solving a global system of equations. It produces polynomial PH spline curves of degree 7.

     

Hermite interpolation with radial basis functions on spheres

We show how conditionally negative definite functions on spheres coupled with strictly completely monotone functions (or functions whose derivative is strictly completely monotone) can be used for Hermite interpolation. The classes of functions thus obtained have the advantage over the strictly positive definite functions studied in [17] that closed form representations (as opposed to series expansions) are readily available. Furthermore, our functions include the historically significant spherical multiquadrics. Numerical results are also presented. This revised version was published online in June 2006 with corrections to the Cover Date.     

Strictly positive definiteness of Hermite interpolation on spheres

In this paper, we obtain some sufficient conditions for positive definite kernels to be strictly positive definite and hence well‐posed for Hermite scattered data interpolation on Euclidean unit spheres. This revised version was published online in June 2006 with corrections to the Cover Date.     

Fourier and barycentric formulae for equidistant Hermite trigonometric interpolation

We consider the Hermite trigonometric interpolation problem of order 1 for equidistant nodes, i.e., the problem of finding a trigonometric polynomial t that interpolates the values of a function and of its derivative at equidistant points. We give a formula for the Fourier coefficients of t in terms of those of the two classical trigonometric polynomials interpolating the values and those of the derivative separately. This formula yields the coefficients with a single FFT. It also gives an aliasing formula for the error in the coefficients which, on its turn, yields error bounds and convergence results for differentiable as well as analytic functions. We then consider the Lagrangian formula and eliminate the unstable factor by switching to the barycentric formula. We also give simplified formulae for even and odd functions, as well as consequent formulae for Hermite interpolation between Chebyshev points.     

Periodic interpolation and wavelets on sparse grids

Nested spaces of multivariate periodic functions forming a non-stationary multiresolution analysis are investigated. The scaling functions of these spaces are fundamental polynomials of Lagrange interpolation on a sparse grid. The approach based on Boolean sums leads to sample and wavelet spaces of significantly lower dimension and good approximation order. The algorithms for complete decomposition and reconstruction are of simple structure and low complexity. This revised version was published online in June 2006 with corrections to the Cover Date.     

Local mass conservative Hermite interpolation for the solution of flow problems by a multi-domain boundary element approach

This paper presents a local Hermite radial basis function interpolation scheme for the velocity and pressure fields. The interpolation for velocity satisfies the continuity equation (mass conservative interpolation) while the pressure interpolation obeys the pressure equation. Additionally, the Dual Reciprocity Boundary Element method (DRBEM) is applied to obtain an integral representation of the Navier-Stokes equations. Then, the proposed local interpolation is used to obtain the values of the field variables and their partial derivatives at the boundary of the sub-domains. This interpolation allows one to obtain the boundary values needed for the integral formulas for velocity and pressure at some nodes within the sub-domains. In the proposed approach the boundary elements are merely used to parameterize the geometry, but not for the evaluation of the integrals as it is usually done. The presented multi-domain approach is different from the traditional ones in boundary elements because the resulting integral equations are non singular and the boundary data needed for the boundary integrals are approximated using a local interpolation. Some accurate results for simple Stokes problems and for the Navier-Stokes equations at low Reynolds numbers up to Re = 400 were obtained.     

Hermite interpolation and its numerical differentiation formula involving symmetric functions

By utilizing symmetric functions,this paper presents explicit representations for Hermite interpolation and its numerical differentiation formula.And the corresponding error estimates are also provided.