The eigenvalues of Hermite and rational spectral differentiation matrices

Summary We derive expressions for the eigenvalues of spectral differentiation matrices for unbounded domains. In particular, we consider Galerkin and collocation methods based on Hermite functions as well as rational functions (a Fourier series combined with a cotangent mapping). We show that (i) first derivative matrices have purely imaginary eigenvalues and second derivative matrices have real and negative eigenvalues, (ii) for the Hermite method the eigenvalues are determined by the roots of the Hermite polynomials and for the rational method they are determined by the Laguerre polynomials, and (iii) the Hermite method has attractive stability properties in the sense of small condition numbers and spectral radii.     

Pythagorean-hodograph preserving mappings

We study the scaled Pythagorean-hodograph (PH) preserving mappings. These mappings make offset-rational isothermal surfaces and map PH curves to PH curves. We present a method to produce a great number of the scaled PH preserving mappings. For an application of the PH preserving mappings, we solve the Hermite interpolation problem for PH curves in the space.     

Computation of the differentiation matrix for the Hermite interpolating polynomials

We obtain formulas for computing the elements of the differentiation matrix for special cases of the Hermite interpolating polynomials. They are expressed in terms of the elements of the differentiation matrices of the Lagrange interpolating polynomials in various systems of interpolation nodes, which can easily be calculated on a computer. These formulas find application in numerical realization of collocation finiteelement methods for solving differential problems.Translated fromVychislitel'naya i Prikladnaya Matematika, Issue 71, 1990, pp. 43–49.     

Approximation of trigonometric functions by trigonometric polynomials with interpolation

The paper studies the approximation order of periodic functions by trigonometric polynomials with interpolation in arbitrary set of nodes. A method of construction of Hermite interpolation polynomials is pointed out.     

An efficient algorithm for periodic Hermite spline interpolation with shifted nodes

Generalized Hermite spline interpolation with periodic splines of defect 2 on an equidistant lattice is considered. Then the classic periodic Hermite spline interpolation with shifted interpolation nodes is obtained as a special case.By means of a new generalization of Euler-Frobenius polynomials the symbol of the considered interpolation problem is defined. Using this symbol, a simple representation of the fundamental splines can be given. Furthermore, an efficient algorithm for the computation of the Hermite spline interpolant is obtained, which is mainly based on the fast Fourier transform.     

Numerical factorization of a polynomial by rational Hermite interpolation

We derive a class of iterative formulae to find numerically a factor of arbitrary degree of a polynomialf(x) based on the rational Hermite interpolation. The iterative formula generates the sequence of polynomials which converge to a factor off(x). It has a high convergence order even for a factor which includes multiple zeros. Some numerical examples are also included.     

Using FFT-based techniques in polynomial and matrix computations: recent advances and applicatons

Computations with univariate polynomials, like the evaluation of product, quotient, remainder, greatest common divisor, etc, are closely related to linear algebra computations performed with structured matrices having the Toeplitz-like or the Hankel-like structures.

The discrete Fourier transform, and the FFT algorithms for its computation, constitute a powerful tool for the design and analysis of fast algorithms for solving problems involving polynomials and structured matrices.

We recall the main correlations between polynomial and matrix computations and present two recent results in this field: in particular, we show how Fourier methods can speed up the solution of a wide class of problems arising in queueing theory where certain Markov chains, defined by infinite block Toeplitz matrices in generalized Hessenberg form, must be solved. Moreover, we present a new method for the numerical factorization of polynomials based on a matrix generalization of Koenig's theorem. This method, that relies on the evaluation/interpolation technique at the Fourier points, reduces the problem of polynomial factorization to the computation of the LU decomposition of a banded Toeplitz matrix with its rows and columns suitably permuted. Numerical experiments that show the effectiveness of our algorithms are presented     

Approximation by interpolating variational splines

In this paper we present an approximation problem of parametric curves and surfaces from a Lagrange or Hermite data set. In particular, we study an interpolation problem by minimizing some functional on a Sobolev space that produces the new notion of interpolating variational spline. We carefully establish a convergence result. Some specific cases illustrate the generality of this work.     

Quadrature rules based on partial fraction expansions

Quadrature rules are typically derived by requiring that all polynomials of a certain degree be integrated exactly. The nonstandard issue discussed here is the requirement that, in addition to the polynomials, the rule also integrates a set of prescribed rational functions exactly. Recurrence formulas for computing such quadrature rules are derived. In addition, Fejér's first rule, which is based on polynomial interpolation at Chebyshev nodes, is extended to integrate also rational functions with pre-assigned poles exactly. Numerical results, showing a favorable comparison with similar rules that have been proposed in the literature, are presented. An error analysis of a representative test problem is given. This revised version was published online in June 2006 with corrections to the Cover Date.     

Harmonic interpolation of Hermite type based on Radon projections with constant distances

We study harmonic interpolation of Hermite type of harmonic functions based on Radon projections with constant distances of chords. We show that the interpolation polynomials are continuous with respect to the angles and the distances. When the chords coalesce to some points on the unit circle, we prove that the interpolation polynomials tend to a Hermite interpolation polynomial at the coalescing points.     

Some researches on multivariate Lagrange interpolation along the sufficiently intersected algebraic manifold

In this article, Lagrange interpolation by polynomials in several variables is studied. Particularly on the sufficiently intersected algebraic manifolds, we discuss the dimension about the interpolation space of polynomials. After defining properly posed set of nodes (or PPSN for short) along the sufficiently intersected algebraic manifolds, we prove the existence of PPSN and give the number of points in PPSN of any degree. Moreover, in order to compute the number of points in PPSN concretely, we propose the operator ? _k with reciprocal difference.     

Newton basis for multivariate Birkhoff interpolation

Multivariate Birkhoff interpolation is the most complex polynomial interpolation problem and people know little about it so far. In this paper, we introduce a special new type of multivariate Birkhoff interpolation and present a Newton paradigm for it. Using the algorithms proposed in this paper, we can construct a Hermite system for any interpolation problem of this type and then obtain a Newton basis for the problem w.r.t. the Hermite system.     

Fast computation of determinants of Bézout matrices and application to curve implicitization

When using bivariate polynomial interpolation for computing the implicit equation of a rational plane algebraic curve given by its parametric equations, the generation of the interpolation data is the most costly of the two stages of the process. In this work a new way of generating those interpolation data with less computational cost is presented. The method is based on an efficient computation of the determinants of certain constant Bézout matrices.     

About measures and nodal systems for which the Hermite interpolants uniformly converge to continuous functions on the circle and interval

We obtain the Laurent polynomial of Hermite interpolation on the unit circle for nodal systems more general than those formed by the n-roots of complex numbers with modulus one. Under suitable assumptions for the nodal system, that is, when it is constituted by the zeros of para-orthogonal polynomials with respect to appropriate measures or when it satisfies certain properties, we prove the convergence of the polynomial of Hermite-Fejér interpolation for continuous functions. Moreover, we also study the general Hermite interpolation problem on the unit circle and we obtain a sufficient condition on the interpolation conditions for the derivatives, in order to have uniform convergence for continuous functions.Finally, we obtain some improvements on the Hermite interpolation problems on the interval and for the Hermite trigonometric interpolation.     

Estimation de l'erreur d'interpolation d'Hermite dans ℝ n

Summary This paper is devoted to study the Hermite interpolation error in an open subset of ⁿ .It follows a previous work of Arcangeli and Gout [1]. Like this one, it is based principally on the paper of Ciarlet and Raviart [7].We obtain two kinds of the Hermite interpolation error, the first from the Hermite interpolation polynomial, the other from approximation method using the Taylor polynomial.Finally in the last part we study some numerical examples concerning straight finite element methods: in the first and second examples, we use finite elements which are included in the affine theory, but it is not the case in the last example. However, in this case, it is possible to refer to the affine theory by the way of particular study (cf. Argyris et al. [2]; Ciarlet [6]; ciarlet and Raviart [7]; Raviart [11]).

     

A single‐node characteristic collocation method for unsteady‐state convection‐diffusion equations in three‐dimensional spaces

We develop a nonconventional single‐node characteristic collocation method with piecewise‐cubic Hermite polynomials for the numerical simulation to unsteady‐state advection‐diffusion transport partial differential equations. This method greatly reduces the number of unknowns in the conventional collocation method, and generates accurate numerical solutions even if very large time steps are taken. The reduction of number of nodes has great potential for problems defined on high space dimensions, which appears in such problems as quantification of uncertainties in subsurface porous media. The method developed here is easy to formulate. Numerical experiments are presented to show the strong potential of the method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 786–802, 2011     

Newton interpolation at Leja points

The Newton form is a convenient representation for interpolation polynomials. Its sensitivity to perturbations depends on the distribution and ordering of the interpolation points. The present paper bounds the growth of the condition number of the Newton form when the interpolation points are Leja points for fairly general compact sets K in the complex plane. Because the Leja points are defined recursively, they are attractive to use with the Newton form. If K is an interval, then the Leja points are distributed roughly like Chebyshev points. Our investigation of the Newton form defined by interpolation at Leja points suggests an ordering scheme for arbitrary interpolation points.Research supported in part by NSF under Grant DMS-8704196 and by U.S. Air Force Grant AFSOR-87-0102.On leave from University of Kentucky, Department of Mathematics, Lexington, KY 40506, U.S.A.     

Ein lokal überlinear konvergentes Verfahren zur Bestimmung von Rückkehrpunkten implizit definierter Raumkurven

Summary A numerical technique is described for computing turning points of a space curveL implicitly defined by a nonlinear system ofn equations inn+1 variables. The basic idea is a local parametrization ofL where the parameter that gives the next iterate is determined by applying one step of the well-known method for minimizing a real function using cubic Hermite interpolation with two nodes. The method is shown to convergeQ-super-linearly and withR-order of at least two. A numerical example concerning the analysis of nonlinear resistive circuits shows the algorithm to work effectively on real life problems.

     

Computation of the Riesz-Herglotz transform and its application to quadrature formulas over the unit circle

Summary A method is proposed for the computation of the Riesz-Herglotz transform. Numerical experiments show the effectiveness of this method. We study its application to the computation of integrals over the unit circle in the complex plane of analytic functions. This approach leads us to the integration by Taylor polynomials. On the other hand, with the goal of minimizing the quadrature error bound for analytic functions, in the set of quadrature formulas of Hermite interpolatory type, we found that this minimum is attained by the quadrature formula based on the integration of the Taylor polynomial. These two different approaches suggest the effectiveness of this formula. Numerical experiments comparing with other quadrature methods with the same domain of validity, or even greater such as Szeg? formulas, (traditionally considered as the counterpart of the Gauss formulas for integrals on the unit circle) confirm the superiority of the numerical estimations. This work was supported by the ministry of education and culture of Spain under contract PB96-1029.     

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.