On Markov operators preserving polynomials

The paper is concerned with a special class of positive linear operators acting on the space C(K)$C(K)$ of all continuous functions defined on a convex compact subset K   of R^d${\mathbf{R}}^d$, d?1$d?1$, having non-empty interior. Actually, this class consists of all positive linear operators T   on C(K)$C(K)$ which leave invariant the polynomials of degree at most 1 and which, in addition, map polynomials into polynomials of the same degree. Among other things, we discuss the existence of such operators in the special case where K is strictly convex by also characterizing them within the class of positive projections. In particular we show that such operators exist if and only if ∂K   is an ellipsoid. Furthermore, a characterization of balls of R^d${\mathbf{R}}^d$ in terms of a special class of them is furnished. Additional results and illustrative examples are presented as well.     

Minimizing the quasi-interpolation error for bivariate discrete quasi-interpolants

We define a class of discrete quasi-interpolants based on bivariate box splines by imposing the exactness on a space of polynomials of total degree, depending on the box spline and minimizing a constant appearing in the leading term of an appropriate quasi-interpolation error estimate. We give some C¹$C^1$ quadratic and C²$C^2$ quartic examples and compare them with other well-known quasi-interpolants.     

Geometric construction of quintic parametric B-splines

The aim of this paper is to present a new class of B-spline-like functions with tension properties. The main feature of these basis functions consists in possessing C³$C^3$ or even C⁴$C^4$ continuity and, at the same time, being endowed by shape parameters that can be easily handled. Therefore they constitute a useful tool for the construction of curves satisfying some prescribed shape constraints. The construction is based on a geometric approach which uses parametric curves with piecewise quintic components.     

KAM theory and the 3D Euler equation

We prove that the dynamical system defined by the hydrodynamical Euler equation on any closed Riemannian 3-manifold M   is not mixing in the C^k$C^k$ topology (k>4$k>4$ and non-integer) for any prescribed value of helicity and sufficiently large values of energy. This can be regarded as a 3D version of Nadirashvili's and Shnirelman's theorems showing the existence of wandering solutions for the 2D Euler equation. Moreover, we obtain an obstruction for the mixing under the Euler flow of C^k$C^k$-neighborhoods of divergence-free vectorfields on M  . On the way we construct a family of functionals on the space of divergence-free C¹$C^1$ vectorfields on the manifold, which are integrals of motion of the 3D Euler equation. Given a vectorfield these functionals measure the part of the manifold foliated by ergodic invariant tori of fixed isotopy types. We use the KAM theory to establish some continuity properties of these functionals in the C^k$C^k$-topology. This allows one to get a lower bound for the C^k$C^k$-distance between a divergence-free vectorfield (in particular, a steady solution) and a trajectory of the Euler flow.     

A note on the zeros of Freud–Sobolev orthogonal polynomials

We prove that the zeros of a certain family of Sobolev orthogonal polynomials involving the Freud weight function e^-x4${\mathrm{e}}^{-x^4}$ on R$\mathbb{R}$ are real, simple, and interlace with the zeros of the Freud polynomials, i.e., those polynomials orthogonal with respect to the weight function e^-x4${\mathrm{e}}^{-x^4}$. Some numerical examples are shown.     

Lifespan of classical discontinuous solutions to general quasilinear hyperbolic systems of conservation laws with small BV initial data: Rarefaction waves

In the present paper the author investigates the global structure stability of Riemann solutions for general quasilinear hyperbolic systems of conservation laws under small BV perturbations of the initial data, where the Riemann solution contains rarefaction waves, while the perturbations are in BV but they are assumed to be C¹$C^1$-smooth, with bounded and possibly large C¹$C^1$-norms. Combining the techniques employed by Li–Kong with the modified Glimm’s functional, the author obtains a lower bound of the lifespan of the piecewise C¹$C^1$ solution to a class of generalized Riemann problems, which can be regarded as a small BV perturbation of the corresponding Riemann problem. This result is also applied to the system of traffic flow on a road network using the Aw–Rascle model.     

Existence of global smooth solution to the initial boundary value problem for p-system with damping

This paper is concerned with the initial boundary value problem for the p$p$-system with damping. We prove the existence of the global smooth solution under the assumption that only the C⁰$C^0$-norm of the derivative of the initial data is sufficiently small, while the C⁰$C^0$-norm of the initial data is not necessarily small. The proof is based on several key a priori estimates, the maximum principle and the characteristic method.     

Semilocal convergence of a continuation method with Hölder continuous second derivative in Banach spaces

In this paper, the semilocal convergence of a continuation method combining the Chebyshev method and the convex acceleration of Newton’s method used for solving nonlinear equations in Banach spaces is established by using recurrence relations under the assumption that the second Frëchet derivative satisfies the Hölder continuity condition. This condition is mild and works for problems in which the second Frëchet derivative fails to satisfy Lipschitz continuity condition. A new family of recurrence relations are defined based on two constants which depend on the operator. The existence and uniqueness regions along with a closed form of the error bounds in terms of a real parameter α∈[0,1]$\alpha \in [0,1]$ for the solution x^∗$x^{\ast }$ is given. Two numerical examples are worked out to demonstrate the efficacy of our approach. On comparing the existence and uniqueness regions for the solution obtained by our analysis with those obtained by using majorizing sequences under Hölder continuity condition on F^″$F^{″}$, it is found that our analysis gives improved results. Further, we have observed that for particular values of the α$\alpha$, our analysis reduces to those for the Chebyshev method (α=0$\alpha =0$) and the convex acceleration of Newton’s method (α=1)$(\alpha =1)$ respectively with improved results.     

On the inverse limit stability of endomorphisms

We present several results suggesting that the concept of C¹$C^1$-inverse (limit structural) stability is free of singularity theory. An example of a robustly transitive, C¹$C^1$-inverse stable endomorphism with a persistent critical set is given. We show that every C¹$C^1$-inverse stable, axiom A endomorphism satisfies a certain strong transversality condition (T). We prove that every attractor–repeller endomorphism satisfying axiom A and condition (T  ) is C¹$C^1$-inverse stable. The latter is applied to Hénon maps, rational functions and others. This leads us to conjecture that C¹$C^1$-inverse stable endomorphisms are exactly those which satisfy axiom A and condition (T).     

High-order compact solution of the one-dimensional heat and advection–diffusion equations

In this work, we propose a high-order accurate method for solving the one-dimensional heat and advection–diffusion equations. We apply a compact finite difference approximation of fourth-order for discretizing spatial derivatives of these equations and the cubic C¹$C^1$-spline collocation method for the resulting linear system of ordinary differential equations. The cubic C¹$C^1$-spline collocation method is an A-stable method for time integration of parabolic equations. The proposed method has fourth-order accuracy in both space and time variables, i.e. this method is of order O(h⁴,k⁴)$O(h^4\text{,}{k}^4)$. Additional to high-order of accuracy, the proposed method is unconditionally stable which will be proved in this paper. Numerical results show that the compact finite difference approximation of fourth-order and the cubic C¹$C^1$-spline collocation method give an efficient method for solving the one-dimensional heat and advection–diffusion equations.     

Convergence analysis for least-squares finite element approximations of second-order two-point boundary value problems

In this paper, a C⁰$C^0$ least-squares finite element method for second-order two-point boundary value problems is considered. The problem is recast as a first-order system. Standard and improved optimal error estimates in maximum-norms are established. Superconvergence estimates at interelement, Lobatto, and Gauss points are developed. Numerical experiments are given to illustrate theoretical results.     

To the theory of viscosity solutions for uniformly elliptic Isaacs equations

We show how a theorem about the solvability in C^1,1$C^{1,1}$ of special Isaacs equations can be used to obtain existence and uniqueness of viscosity solutions of general uniformly nondegenerate Isaacs equations. We apply it also to establish the C^1+χ$C^{1+\chi }$ regularity of viscosity solutions and show that finite-difference approximations have an algebraic rate of convergence. The main coefficients of the Isaacs equations are supposed to be in C^γ$C^{\gamma }$ with γ slightly less than 1/2.     

Weak Dirichlet processes with a stochastic control perspective

The motivation of this paper is to prove verification theorems for stochastic optimal control of finite dimensional diffusion processes without control in the diffusion term, in the case where the value function is assumed to be continuous in time and once differentiable in the space variable (C^0,1$C^{0,1}$) instead of once differentiable in time and twice in space (C^1,2$C^{1,2}$), like in the classical results. For this purpose, the replacement tool of the Itô formula will be the Fukushima–Dirichlet decomposition for weak Dirichlet processes. Given a fixed filtration, a weak Dirichlet process is the sum of a local martingale M$M$ plus an adapted process A$A$ which is orthogonal, in the sense of covariation, to any continuous local martingale. The decomposition mentioned states that a C^0,1$C^{0,1}$ function of a weak Dirichlet process with finite quadratic variation is again a weak Dirichlet process. That result is established in this paper and it is applied to the strong solution of a Cauchy problem with final condition.     

New representation and factorizations of the higher-order ultraspherical-type differential equations

The paper deals with the class of linear differential equations of any even order 2α+4$2\alpha +4$, α∈_N0$\alpha \in {\mathbb{N}}_0$, which are associated with the so-called ultraspherical-type polynomials. These polynomials form an orthogonal system on the interval [−1,1]$[-1,1]$ with respect to the ultraspherical weight function (1−x²)^α${(1-x^2)}^{\alpha }$ and additional point masses of equal size at the two endpoints. The differential equations of “ultraspherical-type” were developed by R. Koekoek in 1994 by utilizing special function methods. In the present paper, a new and completely elementary representation of these higher-order differential equations is presented. This result is used to deduce the orthogonality relation of the ultraspherical-type polynomials directly from the differential equation property. Moreover, we introduce two types of factorizations of the corresponding differential operators of order 2α+4$2\alpha +4$ into a product of α+2$\alpha +2$ linear second-order operators.     

Nonuniform continuity of the solution map to the two component Camassa–Holm system

We prove that the solution map of the two-component Camassa–Holm system is not uniformly continuous as a map from a bounded subset of the Sobolev space H^s(T)×H^r(T)$H^s(\mathbb{T})\times H^r(\mathbb{T})$ to C([0,1],H^s(T)×H^r(T))$C([0,1],H^s(\mathbb{T})\times H^r(\mathbb{T}))$ when s?1$s?1$ and r?0$r?0$. We also demonstrate the nonuniform continuous property in the continuous function space C¹(T)×C¹(T)$C^1(\mathbb{T})\times C^1(\mathbb{T})$.