Curvature variation minimizing cubic Hermite interpolants

In this paper, planar parametric Hermite cubic interpolants with small curvature variation are studied. By minimization of an appropriate approximate functional, it is shown that a unique solution of the interpolation problem exists, and has a nice geometric interpretation. The best solution of such a problem is a quadratic geometric interpolant. The optimal approximation order 4 of the solution is confirmed. The approach is combined with strain energy minimization in order to obtain G¹ cubic interpolatory spline.     

Hermite pseudospectral approximations. An error estimate

We find an error bound for the pseudospectral approximation of a function in terms of Hermite functions, hn(x)=e−x²Hn(x), in certain weighted Sobolev spaces. We use properties of Hermite polynomials, as well as an asymptotic expression for their largest zeroes to achieve the mentioned estimate.     

A Priori Error Estimates for Finite Element Methods for H (2,1)-Elliptic Equations

The convergence of finite element methods for linear elliptic boundary value problems of second and forth order is well understood. In this article, we introduce finite element approximations of some linear semi-elliptic boundary value problem of mixed order on a two-dimensional rectangular domain Q. The equation is of second order in one direction and forth order in the other and appears in the optimal control of parabolic partial differential equations if one eliminates the control and the state (or the adjoint state) in the first order optimality conditions. We establish a regularity result and estimate for the finite element error of conforming approximations of this equation. The finite elements in use have a tensor product structure, in one dimension we use linear, quadratic or cubic Lagrange elements in the other dimension cubic Hermite elements. For these elements, we prove the error bound O(h ² + τ ^k ) in the energy norm and O((h ² + τ ^k )(h ² + τ)) in the L ²(Q)-norm.     

Planar G2 Hermite interpolation with some fair,C-shaped curves

G² Hermite data consists of two points, two unit tangent vectors at those points, and two signed curvatures at those points. The planar G² Hermite interpolation problem is to find a planar curve matching planar G² Hermite data. In this paper, a C-shaped interpolating curve made of one or two spirals is sought. Such a curve is considered fair because it comprises a small number of spirals. The C-shaped curve used here is made by joining a circular arc and a conic in a G² manner. A curve of this type that matches given G² Hermite data can be found by solving a quadratic equation. The new curve is compared to the cubic Bézier curve and to a curve made from a G² join of a pair of quadratics. The new curve covers a much larger range of the G² Hermite data that can be matched by a C-shaped curve of one or two spirals than those curves cover.     

Comportement semi-classique pour l'opérateur de Schrödinger à potentiel périodique

In this paper, the Schrödinger operator with a periodic potential is considered. Let V a smooth periodic function, we study the semi-classical behavior for a continuum spectrum of −h² Δ + V (h → 0). We are interested by localization and width of bands. We give the interaction matrix up to an exponentially small error, measured by Agmon's distance between the wells. A detailed investigation of the spectrum is made for the case where V has one nondegenerate minima per unit cell. We also investigate the spectral properties of −h² Δ + V + ΔV, where ΔV is a smooth positive perturbation with compact support.     

Monotone and Convex C 1 Hermite Interpolants Generated by a Subdivision Scheme

Abstract. We propose C ¹ Hermite interpolants generated by the general subdivision scheme introduced by Merrien [17] and satisfying monotonicity or convexity constraints. For arbitrary values and slopes of a given function f at the end-points of a bounded interval, which are compatible with the contraints, the given algorithms construct shape-preserving interpolants. Moreover, these algorithms are quite simple and fast as well as adapted to CAGD. We also give error estimates in the case of interpolation of smooth functions.     

On using cubic spline for the solution of problems in calculus of variations

Two different approaches based on cubic B-spline are developed to approximate the solution of problems in calculus of variations. Both direct and indirect methods will be described. It is known that, when using cubic spline for interpolating a function g∈C⁴[a,b] on a uniform partition with the step size h, the optimal order of convergence derived is O(h⁴). In Zarebnia and Birjandi (J. Appl. Math. 1–10, 2012) a non-optimal O(h²) method based on cubic B-spline has been used to solve the problems in calculus of variations. In this paper at first we will obtain an optimal O(h⁴) indirect method using cubic B-spline to approximate the solution. The convergence analysis will be discussed in details. Also a locally superconvergent O(h⁶) indirect approximation will be describe. Finally the direct method based on cubic spline will be developed. Some test problems are given to demonstrate the efficiency and applicability of the numerical methods.     

Summability of Hermite Series

Let f∈L ^p (R), 1≤p≤t8, and c _j be the inner product of f and the Hermite function h _j . Assume that c _j 's satisfy $\left| {c_j } \right| \cdot f = o\left( 1 \right)\;\quad as\;j \to \infty $ If r=5/4, then the Hermite series Σc _j h _j conerges to f almost everywhere. If r=9/4-1/p, the Σ c _j h _j converges to f in L ^p (R).     

New error bounds for the penalty method and extrapolation

New error bounds are obtained for the Babu?ka penalty method which justify the use of extrapolation. For the problemΔu=f in Ω,u=g on ?Ω we show that, for a particular choice of boundary weight, repeated extrapolation yields a quasioptimal approximate solution. For example, the error in the second extrapolate (using cubic spline approximants) isO (h ³) when measured in the energy norm. Nearly optimalL ₂ error estimates are also obtained.     

Stability and convergence of fully discrete finite element schemes for the acoustic wave equation

In this paper, we investigate the stability and convergence of some fully discrete finite element schemes for solving the acoustic wave equation where a discontinuous Galerkin discretization in space is used. We first review and compare conventional time-stepping methods for solving the acoustic wave equation. We identify their main properties and investigate their relationship. The study includes the Newmark algorithm which has been used extensively in applications. We present a rigorous stability analysis based on the energy method and derive sharp stability results covering some well-known CFL conditions. A convergence analysis is carried out and optimal a priori error estimates are obtained. For sufficiently smooth solutions, we demonstrate that the maximal error in the L ²-norm error over a finite time interval converges optimally as O(h ^p+1+??t ^s ), where p denotes the polynomial degree, s=1 or 2, h the mesh size, and ??t the time step.     

On Hermite–Hadamard inequality for $${\varvec{}}$$-convex stocastic processes

In this paper, we investigate the Hermite–Hadamard type inequality for the class of some h-convex stochastic processes, which is an extension of the Hermite–Hadamard inequality given by Barráez et al. (Math. Æterna 5:571–581, 2015). We also provide the estimates of both sides of the Hermite–Hadamard type inequality for h-convex stochastic processes, where h is any non-negative function with \(h(t)+h(1-t)\le 1\) for \(0\le t\le 1\).     

Cubic Splines with Minimal Norm

Natural cubic interpolatory splines are known to have a minimal L ₂-norm of its second derivative on the C ² (or W ₂ ² ) class of interpolants. We consider cubic splines which minimize some other norms (or functionals) on the class of interpolatory cubic splines only. The cases of classical cubic splines with defect one (interpolation of function values) and of Hermite C ¹ splines (interpolation of function values and first derivatives) with spline knots different from the points of interpolation are discussed.     

Superconvergence analysis and extrapolation of quasi-Wilson nonconforming finite element method for nonlinear Sobolev equations

Quasi-Wilson nonconforming finite element approximation for a class of nonlinear Sobolev equa-tions is discussed on rectangular meshes. We first prove that this element has two special characters by novel approaches. One is that (▽h ( u-Ihu )1, ▽hvh) h may be estimated as order O ( h2 ) when u ∈ H3 (Ω), where Ihu denotes the bilinear interpolation of u , vh is a polynomial belongs to quasi-Wilson finite element space and ▽h denotes the piecewise defined gradient operator, h is the mesh size tending to zero. The other is that the consistency error of this element is of order O ( h2 ) /O ( h3 ) in broken H 1-norm, which is one/two order higher than its interpolation error when u ∈ H3 (Ω) /H4 (Ω). Then we derive the optimal order error estimate and su- perclose property via mean-value method and the known high accuracy result of bilinear element. Furthermore, we deduce the global superconvergence through interpolation post processing technique. At last, an extrapola- tion result of order O ( h3 ), two order higher than traditional error estimate, is obtained by constructing a new suitable extrapolation scheme.     

Curve design with rational Pythagorean-hodograph curves

The dual Bézier representation offers a simple and efficient constructive approach to rational curves with rational offsets (rational PH curves). Based on the dual form, we develop geometric algorithms for approximating a given curve with aG ² piecewise rational PH curve. The basic components of the algorithms are an appropriate geometric segmentation andG ² Hermite interpolation. The solution involves rational PH curves of algebraic class 4; these curves and important special cases are studied in detail.     

A constructive approach to cubic Hermite Fractal Interpolation Function and its constrained aspects

The theory of splines is a well studied topic, but the kinship of splines with fractals is novel. We introduce a simple explicit construction for a -cubic Hermite Fractal Interpolation Function (FIF). Under some suitable hypotheses on the original function, we establish a priori estimates (with respect to the L _p -norm, 1≤p≤∞) for the interpolation error of the -cubic Hermite FIF and its first derivative. Treating the first derivatives at the knots as free parameters, we derive suitable values for these parameters so that the resulting cubic FIF enjoys global smoothness. Consequently, our method offers an alternative to the standard moment construction of -cubic spline FIFs. Furthermore, we identify appropriate values for the scaling factors in each subinterval and the derivatives at the knots so that the graph of the resulting -cubic FIF lies within a prescribed rectangle. These parameters include, in particular, conditions for the positivity of the cubic FIF. Thus, in the current article, we initiate the study of the shape preserving aspects of fractal interpolation polynomials. We also provide numerical examples to corroborate our results.     

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.     

A class of difference scheme for solving telegraph equation by new non-polynomial spline methods

In this paper, by using a new non-polynomial parameters cubic spline in space direction and compact finite difference in time direction, we get a class of new high accuracy scheme of O(τ⁴ + h²) and O(τ⁴ + h⁴) for solving telegraph equation if we suitably choose the cubic spline parameters. Meanwhile, stability condition of the difference scheme has been carried out. Finally, numerical examples are used to illustrate the efficiency of the new difference scheme.     

Error Estimates for the Solution of the Radial Schr?dinger Equation by the Rayleigh--Ritz Finite Element Method

Approximate eigenvalues and eigenfunctions are obtained forthe radial Schr?dinger equation by applying the Rayleigh—Ritzmethod to a function space consisting of polynomial splinesof odd degree. Computable a posteriori error estimates for theeigenfunction error estimates are obtained. The sharpness ofthese estimates is illustrated for the harmonic oscillator andWoods—Saxon potentials, using both cubic splines and piecewisecubic Hermite polynomials.     

On Quasi-Interpolation with Non-uniformly Distributed Centers on Domains and Manifolds

The paper studies quasi-interpolation by scaled shifts of a smooth and rapidly decaying function. The centers are images of a smooth mapping of the hZⁿ-lattice in R^s, sn, and the scaling parameters are proportional to h. We show that for a large class of generating functions the quasi-interpolants provide high order approximations up to some prescribed accuracy. Although in general the approximants do not converge as h tends to zero, the remaining saturation error is negligible in numerical computations if a scalar parameter is suitably chosen. The lack of convergence is compensated for by a greater flexibility in the choice of generating functions used in numerical methods for solving operator equations.     

An approach to geometric interpolation by Pythagorean-hodograph curves

The problem of geometric interpolation by Pythagorean-hodograph (PH) curves of general degree n is studied independently of the dimension d????2. In contrast to classical approaches, where special structures that depend on the dimension are considered (complex numbers, quaternions, etc.), the basic algebraic definition of a PH property together with geometric interpolation conditions is used. The analysis of the resulting system of nonlinear equations exploits techniques such as the cylindrical algebraic decomposition and relies heavily on a computer algebra system. The nonlinear equations are written entirely in terms of geometric data parameters and are independent of the dimension. The analysis of the boundary regions, construction of solutions for particular data and homotopy theory are used to establish the existence and (in some cases) the number of admissible solutions. The general approach is applied to the cubic Hermite and Lagrange type of interpolation. Some known results are extended and numerical examples provided.