Polynomial approximation via de la Vallée Poussin means

Similar documents
Uniform approximation on the sphere by least squares polynomials

Advanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson

Simultaneous Gaussian quadrature for Angelesco systems

References 167. dx n x2 =2

256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.

Some lecture notes for Math 6050E: PDEs, Fall 2016

On Marcinkiewicz-Zygmund-Type Inequalities

Mean Convergence of Interpolation at Zeros of Airy Functions

Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION

HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS

SHARP L p WEIGHTED SOBOLEV INEQUALITIES

Green s Functions and Distributions

arxiv:math/ v2 [math.ap] 3 Oct 2006

Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform

Multiple Orthogonal Polynomials

Stability of Kernel Based Interpolation

Hermite Interpolation and Sobolev Orthogonality

Approximation theory

High order geometric smoothness for conservation laws

Numerical integration formulas of degree two

Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Distributional Operators

Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity

(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ).

BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH

Vector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms.

On non negative solutions of some quasilinear elliptic inequalities

Stability constants for kernel-based interpolation processes

Orthogonality of hat functions in Sobolev spaces

L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS

Approximation by Conditionally Positive Definite Functions with Finitely Many Centers

Finite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product

A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations

Bernstein-Durrmeyer operators with arbitrary weight functions

ENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS. Lecture 3

ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction

Markov Sonin Gaussian rule for singular functions

Both these computations follow immediately (and trivially) from the definitions. Finally, observe that if f L (R n ) then we have that.

Journal of Inequalities in Pure and Applied Mathematics

Hilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space.

Boot camp - Problem set

On Multivariate Newton Interpolation at Discrete Leja Points

HOMEOMORPHISMS OF BOUNDED VARIATION

PART IV Spectral Methods

On Orthogonal Polynomials in Several Variables

Recurrence Relations and Fast Algorithms

Introduction to Orthogonal Polynomials: Definition and basic properties

Real Variables # 10 : Hilbert Spaces II

arxiv: v2 [math.fa] 17 May 2016

MATH 590: Meshfree Methods

Interpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter

Free energy estimates for the two-dimensional Keller-Segel model

Vectors in Function Spaces

Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials

Regularity of Weak Solution to Parabolic Fractional p-laplacian

The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge

Outline of Fourier Series: Math 201B

On the solution of integral equations of the first kind with singular kernels of Cauchy-type

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems

Nonlinear aspects of Calderón-Zygmund theory

Strong and Weak Weighted Convergence of Jacobi Series

Class notes: Approximation

Interpolation & Polynomial Approximation. Hermite Interpolation I

Nonlinear tensor product approximation

FFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations

Hilbert Transforms and Orthonormal Expansions for Exponential Weights

Natural boundary and Zero distribution of random polynomials in smooth domains arxiv: v1 [math.pr] 2 Oct 2017

The heat equation for the Hermite operator on the Heisenberg group

S. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA

12.0 Properties of orthogonal polynomials

Gaussian interval quadrature rule for exponential weights

NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS

Wavelets and modular inequalities in variable L p spaces

Hardy-Littlewood maximal operator in weighted Lorentz spaces

Error estimates for the Raviart-Thomas interpolation under the maximum angle condition

SELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY

Optimal global rates of convergence for interpolation problems with random design

Regularity of the p-poisson equation in the plane

Nonnegative linearization for little q-laguerre polynomials and Faber basis

Supplement: Universal Self-Concordant Barrier Functions

Spectral Gap and Concentration for Some Spherically Symmetric Probability Measures

DETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains

Moreover this binary operation satisfies the following properties

LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011

Recall: Dot product on R 2 : u v = (u 1, u 2 ) (v 1, v 2 ) = u 1 v 1 + u 2 v 2, u u = u u 2 2 = u 2. Geometric Meaning:

Error formulas for divided difference expansions and numerical differentiation

Scientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1

LECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes).

A trigonometric orthogonality with respect to a nonnegative Borel measure

Subelliptic mollifiers and a basic pointwise estimate of Poincaré type

Lecture 7: Kernels for Classification and Regression

On the constant in the reverse Brunn-Minkowski inequality for p-convex balls.

A collocation method for solving some integral equations in distributions

ON A CLASS OF LINEAR POSITIVE BIVARIATE OPERATORS OF KING TYPE

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Second Order Elliptic PDE

POTENTIAL LANDESMAN-LAZER TYPE CONDITIONS AND. 1. Introduction We investigate the existence of solutions for the nonlinear boundary-value problem

Recent Results for Moving Least Squares Approximation

Probability and Measure

Singular Integrals. 1 Calderon-Zygmund decomposition

Transcription:

Summer School on Applied Analysis 2 TU Chemnitz, September 26-3, 2 Polynomial approximation via de la Vallée Poussin means Lecture 2: Discrete operators Woula Themistoclakis CNR - National Research Council of Italy Institute for Computational Applications Mauro Picone, Naples, Italy. Woula Themistoclakis, Chemnitz, September 26-3, 2

De la Vallée Poussin means: V n,m (w, f, x) := 2m where S n (w, f, x) := c k (w, f) := n c k (w, f)p k (w, x) k= p k (w, y)f(y)w(y)dy n+m k=n m S k (w, f, x) Fourier sum Fourier coefficients w(x) := v α,β (x) = ( x) α ( + x) β, α, β >, being a Jacobi weight and {p j (w, x)} j the corresponding system of orthonormal Jacobi polynomials Quasi-projection: V n,m (w, P ) = P, P P n m Near best polynomial: (f V n,m (w, f))u p CE n m (f) u,p holds for any p, m = θn (θ ], [ fixed) and suitable u, w. Woula Themistoclakis, Chemnitz, September 26-3, 2 2

Other forms: V n,m (w, f, x) = V n,m (w, f, x) = n+m k= H n,m (w, x, y)f(y)w(y)dy () µ m n,k c k (w, f)p k (w, x) (2) n+m H n,m (w, x, y) := K r (w, x, y) 2m r=n m r K r (w, x, y) := p j (w, x)p j (w, y) j= de la Vallée Poussin kernel Darboux kernel c k (w, f) := µ m n,k := p k (w, y)f(y)w(y)dy if k n m, n + m k 2m if n m < k < n + m. Fourier coefficients Woula Themistoclakis, Chemnitz, September 26-3, 2 3

Lagrange interpolation at the zeros of orthogonal polynomials Take the Fourier sum: S n (w, f, x) := K n (w, x, y)f(y)w(y)dy Apply the Gaussian rule: g(x)w(x)dx = n λ n,j g(x n,j ), j= g P 2n In this way we obtain the Lagrange polynomial of degree n interpolating f at the n zeros {x n,j } j of p n (w), i.e. L n (w, f, x) := n λ n,j K n (w, x, x n,j )f(x n,j ) j= with λ n,j K n (w, x n,i, x n,j ) = δ i,j, j =,..., n. Woula Themistoclakis, Chemnitz, September 26-3, 2 4

DISCRETE DE LA VALLEE POUSSIN MEANS By applying the previous Gaussian rule to the integrals in () V n,m (w, f, x) = (2) V n,m (w, f, x) = n+m k= H n,m (w, x, y)f(y)w(y)dy µ m n,k c k (w, f)p k (w, x) we obtain the following discrete operator: ( ) Ṽ n,m (w, f, x) = (2 ) Ṽ n,m (w, f, x) = n λ n,j H n,m (w, x, x n,j )f(x n,j ) j= n+m k= µ m n,k c n,k (w, f) p k (w, x) Discrete Fourier coefficients: c n,k (w, f) := n j= λ n,jf(x n,j )p k (w, x n,j ) Woula Themistoclakis, Chemnitz, September 26-3, 2 5

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) = n j= λ n,j K n (w, x, x n,j )f(x n,j ), Ṽ n,m (w, f, x) = n j= λ n,j H n,m (w, x, x n,j )f(x n,j ), m < n Computation: Both computable from the data f(x n,j ), j =,.., n Invariance: { Ln (w, P ) = S n (w, P ) = P, P P n Ṽ n,m (w, P ) = V n,m (w, P ) = P, P P n m Approximation: { L n (w) C u C u C log n Is sup n,m Ṽn,m(w) C u C u for any u < true for suitable u?? Interpolation: { Ln (w, f, x n,j ) = f(x n,j ), j =,..., n Is it true that Ṽn,m(w, f, x n,j ) = f(x n,j )?? Woula Themistoclakis, Chemnitz, September 26-3, 2 6

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) = n j= λ n,j K n (w, x, x n,j )f(x n,j ), Ṽ n,m (w, f, x) = n j= λ n,j H n,m (w, x, x n,j )f(x n,j ), m < n Computation: Both computable from the data f(x n,j ), j =,.., n Invariance: { Ln (w, P ) = S n (w, P ) = P, P P n Ṽ n,m (w, P ) = V n,m (w, P ) = P, P P n m Approximation: { L n (w) C u C u C log n Is sup n,m Ṽn,m(w) C u C u for any u < true for suitable u?? Interpolation: { Ln (w, f, x n,j ) = f(x n,j ), j =,..., n Is it true that Ṽn,m(w, f, x n,j ) = f(x n,j )?? Woula Themistoclakis, Chemnitz, September 26-3, 2 7

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) = n j= λ n,j K n (w, x, x n,j )f(x n,j ), Ṽ n,m (w, f, x) = n j= λ n,j H n,m (w, x, x n,j )f(x n,j ), m < n Computation: Both computable from the data f(x n,j ), j =,.., n Invariance: { Ln (w, P ) = S n (w, P ) = P, P P n Ṽ n,m (w, P ) = V n,m (w, P ) = P, P P n m Approximation: { L n (w) C u C u C log n Is sup n,m Ṽn,m(w) C u C u for any u < true for suitable u?? Interpolation: { Ln (w, f, x n,j ) = f(x n,j ), j =,..., n Is it true that Ṽn,m(w, f, x n,j ) = f(x n,j )?? Woula Themistoclakis, Chemnitz, September 26-3, 2 8

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) = n j= λ n,j K n (w, x, x n,j )f(x n,j ), Ṽ n,m (w, f, x) = n j= λ n,j H n,m (w, x, x n,j )f(x n,j ), m < n Computation: Both computable from the data f(x n,j ), j =,.., n Invariance: { Ln (w, P ) = S n (w, P ) = P, P P n Ṽ n,m (w, P ) = V n,m (w, P ) = P, P P n m Approximation: { Ln (w) C u Cu C log n for any u Is sup n,m Ṽn,m(w) C u Cu < true for suitable u?? Interpolation: { Ln (w, f, x n,j ) = f(x n,j ), j =,..., n Is it true that Ṽn,m(w, f, x n,j ) = f(x n,j )?? Woula Themistoclakis, Chemnitz, September 26-3, 2 9

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) = n j= λ n,j K n (w, x, x n,j )f(x n,j ), Ṽ n,m (w, f, x) = n j= λ n,j H n,m (w, x, x n,j )f(x n,j ), m < n Computation: Both computable from the data f(x n,j ), j =,.., n Invariance: { Ln (w, P ) = S n (w, P ) = P, P P n Ṽ n,m (w, P ) = V n,m (w, P ) = P, P P n m Approximation: { Ln (w) C u Cu C log n for any u Is sup n,m Ṽn,m(w) C u Cu < true for suitable u?? Interpolation: { Ln (w, f, x n,j ) = f(x n,j ), j =,..., n Is it true that Ṽn,m(w, f, x n,j ) = f(x n,j )?? Woula Themistoclakis, Chemnitz, September 26-3, 2

Connection with the continuous de la Vallée Poussin operator Case p < : We have Ṽn,m(w, f)u p C V n,m (w) L p w u L p w u where p + p = and λ n(v, x n,k ) := [ n k= ( n k= p 2 k(v, x n,k ) λ n (u p, x n,k ) f(x n,k ) p ) p ] v(x n,k ) x 2 n,k n Case p = : We have Ṽn,m(w, f)u C V n,m (w) C u C u ( ) max f(x n,k) u(x n,k ) k n C > being independent of n, m, f in both the cases. Woula Themistoclakis, Chemnitz, September 26-3, 2

Proof for p =. By λ n,k = λ n (w, x n,k ) w(x n,k ) Marcinkiewicz inequality, we get x 2 n,k n, and [ ] n Ṽn,m(w, f)u = max u(x) λ n (w, x n,k )H n,m (w, x, x n,k )f(x n,k ) x k= [ ] n ( w ) C max u(x) λ n x u, x n,k H n,m (w, x, x n,k )(fu)(x n,k ) k= ( ) [ ] n ( w ) C max (fu)(x n,k) u(x) λ n k n u, x n,k H n,m (w, x, x n,k ) C = C ( max k n (fu)(x n,k) ) ( ) max (fu)(x n,k) k n max x max x [ u(x) k= V n,m (w) C u C u. H n,m (w, x, y) w(y) ] u(y) dy Woula Themistoclakis, Chemnitz, September 26-3, 2 2

APPROXIMATION IN C u Theorem : Let w = v α,β L [, ] and u = v γ,δ C[, ] satisfy α 2 4 β 2 4 < γ α 2 + 5 4, < δ β 2 + 5 4 and γ δ α β 2 Then for all integers n and any m = θn with < θ < arbitrarily fixed, the map Ṽn,m(w) : C u C u is uniformly bounded w.r.t. n, m and E n+m (f) u, [f Ṽn,m(w, f)]u CE n m (f) u, holds for every f C u, C > being independent of f, n, m. Woula Themistoclakis, Chemnitz, September 26-3, 2 3

Theorem 2: Let w = v α,β (α, β > ) and u = v γ,δ (γ, δ ) satisfy α 2 + 4 ν < γ α 2 + 5 4 ν, β 2 + 4 ν < δ β 2 + 5 4 ν, for some ν 2 Then for all integers n and m = θn with < θ < arbitrarily fixed, the map Ṽn,m(w) : C u C u is uniformly bounded w.r.t. n, m and E n+m (f) u, [f Ṽn,m(w, f)]u CE n m (f) u, holds for every f C u, C > being independent of f, n, m. Woula Themistoclakis, Chemnitz, September 26-3, 2 4

COMPARISON WITH LAGRANGE INTERPOLATION Let w = v α,β L [, ] and u = v γ,δ C[, ] be such that: α 2 + 4 β 2 + 4 γ α 2 + 5 4, δ β 2 + 5 4. Then for all sufficiently large pair of integers n and m = θn ( < θ < fixed) and for each f C u, we have Lagrange error: [f L n (w, f)]u C log n E n (f) u, De la V.P. error: [f Ṽ m n (w, f)]u CE n m (f) u, where C > is independent of n, f in both the cases. Woula Themistoclakis, Chemnitz, September 26-3, 2 5

APPROXIMATION IN L p u Theorem 3: Let p < and assume that w = v α,β L [, ] and u = v γ,δ L p [, ], with w u Lp [, ], satisfy the bounds α 2 + 4 ν < γ + p α 2 + 5 4 ν β 2 + 4 ν < δ + p β 2 + 5 4 ν for some ν 2 Then for all n, m N with m = θn ( < θ < fixed) and each f L p u (everywhere defined on ], [), we have Ṽn,m(w, f)u p C ( n k= where we recall that λ n (u p, x n,k ) u p (x n,k ) λ n (u p, x n,k ) f(x n,k ) p ) p x 2 n,k n. Lagrange case L n (w, f): The same estimate holds with p and ν =. Woula Themistoclakis, Chemnitz, September 26-3, 2 6

Error estimates in Sobolev type spaces Wr p (u) := { } f L p u : f (r ) AC loc, and f (r) ϕ r L p u f p Wr (u) := fu p + f (r) ϕ r u p, ϕ(x) := x 2 Note: For any f W p r (u), we have E n (f) u,p C n r f (r) ϕ r u p Theorem 4: Under the assumptions of Theorem 3, for all f W p r (u), we have [f Ṽn,m(w, f)]u p C n r f (r) ϕ r u p, f Ṽn,m(w, f) W p s (u) C n r s f W p r (u), < s r where C > is independent of f, n, m and p (setting L u := C u). Lagrange case: L n (w, f) verifies the same estimates, but for p / {, } and ν =. Woula Themistoclakis, Chemnitz, September 26-3, 2 7

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) := Ṽ n,m (w, f, x) := n λ n,k K n (w, x, x n,k )f(x n,k ), k= n λ n,k H n,m (w, x, x n,k )f(x n,k ), k= n > m Invariance: Ṽn,m(w) : f Ṽn,m(w, f) P n+m is a quasi projection Approximation: Ṽn,m(w, f) solves the critical cases p =,. Interpolation: Is it true that Ṽn,m(w, f, x n,k ) = f(x n,k ), k =,.., n?? Woula Themistoclakis, Chemnitz, September 26-3, 2 8

Theorem 5 Let w be such that, for all x ], [, we have (3) p n+s (w, x)+p n s (w, x) = p n (w, x)q(x), deg(q) s < n, Then Ṽn,m(w, f, x n,i ) = f(x n,i ), i =,.., n, holds for all n m >. Examples: Bernstein Szego weights defined by w(x) := ( x) α ( + x) β, α = β =, Chebyshev weights 2 w(x) := w(x) :=, w(x) := p(x) x 2 p(x) x2, deg(p) 2 p(x) x + x, deg(p) provide polynomials satisfying (3). Woula Themistoclakis, Chemnitz, September 26-3, 2 9

Proof. Note that we can write H n,m (w, x n,i, x n,j ) = 2m where: m r= K n+r (w, x n,i, x n,j ) = K n (w, x n,i, x n,j ) + [ Kn+r (w, x n,i, x n,j ) + K n (r+) (w, x n,i, x n,j ) ] r p n+s (w, x n,i )p n+s (w, x n,j ), + = + + r K n (r+) (w, x n,i, x n,j ) = K n (w, x n,i, x n,j ) p n s (w, x n,i )p n s (w, x n,j ). Hence by (3) we get p n+s (w, x n,i ) = p n s (w, x n,i ), i =,.., n, and the kernel in n Ṽ n,m (w, f, x n,i ) = λ n,j H n,m (w, x n,i, x n,j )f(x n,j ) j= s= s= reduces to H n,m (w, x n,i, x n,j ) = K n (w, x n,i, x n,j ) = δ i,j [λ n,j ]. Woula Themistoclakis, Chemnitz, September 26-3, 2 2

COMPARISON WITH LAGRANGE INTERPOLATION L n (w, f, x) := Ṽ n,m (w, f, x) := n λ n,k K n (w, x, x n,k )f(x n,k ), k= n λ n,k H n,m (w, x, x n,k )f(x n,k ), k= n > m Invariance: { Ln (w) : f L n (w, f) P n projection Ṽ n,m (w) : f Ṽn,m(w, f) P n+m quasi projection Approximation: { [f Ln (w, f)]u C log n E n (f) u, [f Ṽn,m(w, f)]u CE n m (f) u, Interpolation: { Ln (w, f, x n,k ) = f(x n,k ), for all w = v α,β Ṽ n,m (w, f, x n,k ) = f(x n,k ), α = β = 2, n m Woula Themistoclakis, Chemnitz, September 26-3, 2 2

De la Vallée Poussin type polynomial spaces DEF: S n,m (w) := span{ λ n,k H n,m (w, x, x n,k ) : k =,..., n} Interpolation property: λ n,k H n,m (w, x n,h, x n,k ) = δ h,k dim S n,m (w) = n Invariance property: Ṽn,m(w, P ) = P, deg(p ) n m P n m S n,m (w) P n+m Theorem 6: In the interpolating case, w = v α,β with α = β = 2, the operator Ṽn,m(w) : f Ṽn,m(w, f) is a projection on S n,m (w), i.e. we have f S n,m (w) f = Ṽn,m(w, f) Woula Themistoclakis, Chemnitz, September 26-3, 2 22

INTERPOLATING BASIS OF S n,m (w): { S n,m (w) := span Φ m n,k (w, x) := λ n,k H n,m (w, x, x n,k ), } k =,..., n De la V. P. interpolating polynomial: Ṽ n,m (w, f, x) = n f(x n,k )Φ m n,k(w, x) k= Under the assumptions of Theorem 3, for all a k R, k =,.., n, we have ( ) ( n p n u λ a k Φn,k(w)) m k (u p, x n,k ) a k p if p < k= k= p max a k u(x n,k ) if p = k n i.e. {Φ m n,k (w)} k is a Marcinkiewicz basis in L p u, for all p. Woula Themistoclakis, Chemnitz, September 26-3, 2 23

q k (w) := ORTHOGONAL BASIS OF S n,m (w) p k (w) if k n m m + n k 2m p k(w) m n + k 2m p 2n k(w) if n m < k < n Theorem 7: The set {q k (w)} k is an orthogonal basis of S n,m (w), i.e. we have S n,m (w) := span{q k (w) : k =,,..., n } with q h (w, x)q k (w, x)w(x)dx = δ h,k if k n m m 2 +(n k) 2 De la Vallée Poussin interpolating polynomial: [ n n ] Ṽ n,m (w, f, x) = λ n,i p k (w, x n,i )f(x n,i ) q k (w, x) k= i= 2m 2 if n m < k < n Woula Themistoclakis, Chemnitz, September 26-3, 2 24

Proof. We are going to state the basis transformation n Φ m n,k(w, x) = λ n,k j= p j (w, x n,k )q j (w, x), k =,..., n. Recall that n+m Φ m n,k(w, s x) := λ n,k H n,m (w, x, x n,k ) = λ n,k µ m n,j p j (w, x n,k )p j (w, x) = λ n,k + n m j= n+m j=n+ n j= n + m j p j (w, x n,k )p j (w, x) + 2m p j(w, x n,k )p j (w, x) j=n m+ n + m j 2m p j(w, x n,k )p j (w, x) i.e., by changing the summation variables, we have Woula Themistoclakis, Chemnitz, September 26-3, 2 25

Φ m n,k = λ n,k + m s= n m j= p j (w, x n,k )p j (w) + m s= ] m s 2m p n+s(w, x n,k )p n+s (w) and using p n+s (w, x n,k ) = p n s (w, x n,k ), we get n m Φ m n,k s = λ n,k j= m + λ n,k n = λ n,k j= s= p j (w, x n,k )p j (w) + p n s (w, x n,k ) p j (w, x n,k )q j (w). m + s 2m p n s(w, x n,k )p n s (w) [ m + s 2m p n s(w) m s ] 2m p n+s(w) Woula Themistoclakis, Chemnitz, September 26-3, 2 26

COMPARISON WITH LAGRANGE INTERPOLATION Chebyshev case: L n (w, f, x) := Ṽ n,m (w, f, x) := n k= n k= c n,k (w, f) p k (w, x), c n,k (w, f) q k (w, x), n > m Interpolation: { Ln (w, f, x n,k ) = f(x n,k ), k =,..., n Ṽ n,m (w, f, x n,k ) = f(x n,k ), k =,..., n Invariance: { Ln (w) : f L n (w, f) P n projection Ṽ n,m (w) : f Ṽn,m(w, f) S n,m (w) projection Approximation: { [f Ln (w, f)]u C log n E n (f) u, [f Ṽn,m(w, f)]u CE n m (f) u, Woula Themistoclakis, Chemnitz, September 26-3, 2 27

.5 Lagrange interpolation for n = 5.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 28

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =..5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 29

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 3

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.3.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 3

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.4.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 32

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.5.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 33

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.6.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 34

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.7.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 35

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.8.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 36

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.9.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 37

.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 38

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =..8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 39

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.2.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 4

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.3.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 4

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.4.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 42

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.5.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 43

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.6.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 44

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.7.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 45

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.8.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 46

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.9.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 47

Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.8.6.4.2.2.4.6.8 Woula Themistoclakis, Chemnitz, September 26-3, 2 48

.2 Lagrange interpolation for n = 2.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 49

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =..8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 5

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.2.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 5

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.3.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 52

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.4.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 53

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.5.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 54

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.6.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 55

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.7.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 56

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.8.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 57

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.9.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 58

.2 De la V.P. interpolation for n = 2 and m = θ n, with θ =.8.6.4.2.2 Woula Themistoclakis, Chemnitz, September 26-3, 2 59

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =..3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 6

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.2.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 6

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.3.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 62

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.4.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 63

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.5.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 64

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.6.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 65

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.7.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 66

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.8.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 67

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.9.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 68

.4 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.3.2...2.3.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 69

.2 Fundamental Lagrange polynomial for n = 5.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 7

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =..8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 7

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.2.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 72

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.3.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 73

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.4.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 74

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.5.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 75

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.6.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 76

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.7.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 77

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.8.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 78

.2 De la VP kernel polynomials for n = 5 and m = θ n, with θ =.9.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 79

.2 Fejer kernel polynomials for n = 5 and m = θ n, with θ =.8.6.4.2.2.4 Woula Themistoclakis, Chemnitz, September 26-3, 2 8

2.5 Lagrange interpolation for n = 5 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 8

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =. 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 82

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.2 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 83

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.3 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 84

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.4 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 85

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.5 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 86

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.6 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 87

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.7 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 88

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.8 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 89

2.5 De la V.P. interpolation for n = 5 and m = θ n, with θ =.9 2.5.5.5.5 2 2.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 9

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =..5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 9

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 92

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.3.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 93

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.4.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 94

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.5.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 95

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.6.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 96

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.7.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 97

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.8.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 98

.5 Error curves by Lagrange (in red) and de la V.P. interpolation for θ =.9.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 99

3.5 Lagrange interpolation for n = 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2

3.5 De la V.P. interpolation for n = and m = θ n, with θ =. 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.2 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 2

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.3 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 3

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.4 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 4

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.5 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 5

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.6 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 6

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.7 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 7

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.8 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 8

3.5 De la V.P. interpolation for n = and m = θ n, with θ =.9 3 2.5 2.5.5.5 Woula Themistoclakis, Chemnitz, September 26-3, 2 9

SOME REFERENCES: M.R.Capobianco, G.Criscuolo, G.Mastroianni, Special Lagrange and Hermite interpolation processes in: Approximation Theory and Applications, T.M.Rassias Editor, Hadronic Press, Palm Harbor, USA, ISBN -57485-4-5, 998, pp.37 62 M.R.Capobianco, W.Themistoclakis, On the boundedness of some de la Vallée Poussin operators, East J. Approx., 7, n. 4 (2), 47 444. Corrigendum in East J. Approx., 3, n. 2 (27), 223 226. G.Criscuolo, G.Mastroianni, Fourier and Lagrange operators in some Sobolev type spaces, Acta Sci. Math. (Szeged), 6 (995), 3 48. G.Mastroianni, G.Milovanovic, Interpolation processes. Basic theory and applications. Springer Monographs in Mathematics. Springer-Verlag, Berlin, 28. G. Mastroianni, M.G.Russo, Lagrange interpolation in weighted Besov spaces, Constr. Approx., 5 (999), 257 289. G. Mastroianni, W. Themistoclakis, De la Vallée Poussin means and Jackson theorem, Acta Sci. Math. (Szeged), 74 (28), 47 7. W.Themistoclakis, Some interpolating operators of de la Vallée Poussin type, Acta Math. Hungar. 84 No.3 (999), 22-235 N.Trefethen, Approximation Theory and Approximation Practice (In preparation) http://www2.maths.ox.ac.uk/chebfun/atap/ Woula Themistoclakis, Chemnitz, September 26-3, 2