Polynomial approximation via de la Vallée Poussin means

Size: px
Start display at page:

Download "Polynomial approximation via de la Vallée Poussin means"

Transcription

1 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

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

3 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

4 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

5 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

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 6

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 7

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: { 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

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 9

10 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

11 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

12 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

13 APPROXIMATION IN C u Theorem : Let w = v α,β L [, ] and u = v γ,δ C[, ] satisfy α 2 4 β 2 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

14 Theorem 2: Let w = v α,β (α, β > ) and u = v γ,δ (γ, δ ) satisfy α ν < γ α ν, β ν < δ β ν, 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

15 COMPARISON WITH LAGRANGE INTERPOLATION Let w = v α,β L [, ] and u = v γ,δ C[, ] be such that: α β γ α , δ β 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

16 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 α ν < γ + p α ν β ν < δ + p β ν 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

17 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

18 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

19 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

20 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

21 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

22 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

23 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

24 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

25 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

26 Φ 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

27 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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

110 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 , 998, pp M.R.Capobianco, W.Themistoclakis, On the boundedness of some de la Vallée Poussin operators, East J. Approx., 7, n. 4 (2), Corrigendum in East J. Approx., 3, n. 2 (27), G.Criscuolo, G.Mastroianni, Fourier and Lagrange operators in some Sobolev type spaces, Acta Sci. Math. (Szeged), 6 (995), 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), G. Mastroianni, W. Themistoclakis, De la Vallée Poussin means and Jackson theorem, Acta Sci. Math. (Szeged), 74 (28), W.Themistoclakis, Some interpolating operators of de la Vallée Poussin type, Acta Math. Hungar. 84 No.3 (999), N.Trefethen, Approximation Theory and Approximation Practice (In preparation) Woula Themistoclakis, Chemnitz, September 26-3, 2

Uniform approximation on the sphere by least squares polynomials

Uniform approximation on the sphere by least squares polynomials Uniform approximation on the sphere by least squares polynomials Woula Themistoclakis Marc Van Barel August 7, 8 Abstract The paper concerns the uniform polynomial approximation of a function f, continuous

More information

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

Advanced Computational Fluid Dynamics AA215A Lecture 2 Approximation Theory. Antony Jameson Advanced Computational Fluid Dynamics AA5A Lecture Approximation Theory Antony Jameson Winter Quarter, 6, Stanford, CA Last revised on January 7, 6 Contents Approximation Theory. Least Squares Approximation

More information

Simultaneous Gaussian quadrature for Angelesco systems

Simultaneous Gaussian quadrature for Angelesco systems for Angelesco systems 1 KU Leuven, Belgium SANUM March 22, 2016 1 Joint work with Doron Lubinsky Introduced by C.F. Borges in 1994 Introduced by C.F. Borges in 1994 (goes back to Angelesco 1918). Introduced

More information

References 167. dx n x2 =2

References 167. dx n x2 =2 References 1. G. Akemann, J. Baik, P. Di Francesco (editors), The Oxford Handbook of Random Matrix Theory, Oxford University Press, Oxford, 2011. 2. G. Anderson, A. Guionnet, O. Zeitouni, An Introduction

More information

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.

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. 56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =

More information

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

Some lecture notes for Math 6050E: PDEs, Fall 2016 Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.

More information

On Marcinkiewicz-Zygmund-Type Inequalities

On Marcinkiewicz-Zygmund-Type Inequalities On Marcinkiewicz-Zygmund-Type Inequalities H. N. Mhaskar Department of Mathematics, California State University, Los Angeles, California, 90032, U.S.A. J. Prestin FB Mathematik Universität Rostock, Universitätsplatz

More information

Mean Convergence of Interpolation at Zeros of Airy Functions

Mean Convergence of Interpolation at Zeros of Airy Functions Mean Convergence of Interpolation at Zeros of Airy Functions D. S. Lubinsky Abstract The classical Erdős-Turán theorem established mean convergence of Lagrange interpolants at zeros of orthogonal polynomials.

More information

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

Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of

More information

HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS

HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 7, July 1996 HERMITE MULTIPLIERS AND PSEUDO-MULTIPLIERS JAY EPPERSON Communicated by J. Marshall Ash) Abstract. We prove a multiplier

More information

SHARP L p WEIGHTED SOBOLEV INEQUALITIES

SHARP L p WEIGHTED SOBOLEV INEQUALITIES Annales de l Institut de Fourier (3) 45 (995), 6. SHARP L p WEIGHTED SOBOLEV INEUALITIES Carlos Pérez Departmento de Matemáticas Universidad Autónoma de Madrid 28049 Madrid, Spain e mail: cperezmo@ccuam3.sdi.uam.es

More information

Green s Functions and Distributions

Green s Functions and Distributions CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where

More information

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

arxiv:math/ v2 [math.ap] 3 Oct 2006 THE TAYLOR SERIES OF THE GAUSSIAN KERNEL arxiv:math/0606035v2 [math.ap] 3 Oct 2006 L. ESCAURIAZA From some people one can learn more than mathematics Abstract. We describe a formula for the Taylor series

More information

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

Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform Convergence and Stability of a New Quadrature Rule for Evaluating Hilbert Transform Maria Rosaria Capobianco CNR - National Research Council of Italy Institute for Computational Applications Mauro Picone,

More information

Multiple Orthogonal Polynomials

Multiple Orthogonal Polynomials Summer school on OPSF, University of Kent 26 30 June, 2017 Introduction For this course I assume everybody is familiar with the basic theory of orthogonal polynomials: Introduction For this course I assume

More information

Stability of Kernel Based Interpolation

Stability of Kernel Based Interpolation Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen

More information

Hermite Interpolation and Sobolev Orthogonality

Hermite Interpolation and Sobolev Orthogonality Acta Applicandae Mathematicae 61: 87 99, 2000 2000 Kluwer Academic Publishers Printed in the Netherlands 87 Hermite Interpolation and Sobolev Orthogonality ESTHER M GARCÍA-CABALLERO 1,, TERESA E PÉREZ

More information

Approximation theory

Approximation theory Approximation theory Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 1 1.3 6 8.8 2 3.5 7 10.1 Least 3squares 4.2

More information

High order geometric smoothness for conservation laws

High order geometric smoothness for conservation laws High order geometric smoothness for conservation laws Martin Campos-Pinto, Albert Cohen, and Pencho Petrushev Abstract The smoothness of the solutions of 1D scalar conservation laws is investigated and

More information

Numerical integration formulas of degree two

Numerical integration formulas of degree two Applied Numerical Mathematics 58 (2008) 1515 1520 www.elsevier.com/locate/apnum Numerical integration formulas of degree two ongbin Xiu epartment of Mathematics, Purdue University, West Lafayette, IN 47907,

More information

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

Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Distributional Operators Noname manuscript No. (will be inserted by the editor) Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Distributional Operators Gregory E. Fasshauer Qi Ye Abstract

More information

Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity

Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Approximating Infinity-Dimensional Stochastic Darcy s Equations without Uniform Ellipticity Marcus Sarkis Jointly work with Juan Galvis USC09 Marcus Sarkis (WPI) SPDE without Uniform Ellipticity USC09

More information

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

(0, 0), (1, ), (2, ), (3, ), (4, ), (5, ), (6, ). 1 Interpolation: The method of constructing new data points within the range of a finite set of known data points That is if (x i, y i ), i = 1, N are known, with y i the dependent variable and x i [x

More information

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

BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH BIVARIATE LAGRANGE INTERPOLATION AT THE PADUA POINTS: THE IDEAL THEORY APPROACH LEN BOS, STEFANO DE MARCHI, MARCO VIANELLO, AND YUAN XU Abstract The Padua points are a family of points on the square [,

More information

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

Vector Spaces. Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. Vector Spaces Vector space, ν, over the field of complex numbers, C, is a set of elements a, b,..., satisfying the following axioms. For each two vectors a, b ν there exists a summation procedure: a +

More information

On non negative solutions of some quasilinear elliptic inequalities

On non negative solutions of some quasilinear elliptic inequalities On non negative solutions of some quasilinear elliptic inequalities Lorenzo D Ambrosio and Enzo Mitidieri September 28 2006 Abstract Let f : R R be a continuous function. We prove that under some additional

More information

Stability constants for kernel-based interpolation processes

Stability constants for kernel-based interpolation processes Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento

More information

Orthogonality of hat functions in Sobolev spaces

Orthogonality of hat functions in Sobolev spaces 1 Orthogonality of hat functions in Sobolev spaces Ulrich Reif Technische Universität Darmstadt A Strobl, September 18, 27 2 3 Outline: Recap: quasi interpolation Recap: orthogonality of uniform B-splines

More information

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

L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety

More information

Approximation by Conditionally Positive Definite Functions with Finitely Many Centers

Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal

More information

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

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 Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )

More information

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

A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations Acta Appl Math (21) 19: 861 873 DOI 1.17/s144-8-9351-y A Product Integration Approach Based on New Orthogonal Polynomials for Nonlinear Weakly Singular Integral Equations M. Rasty M. Hadizadeh Received:

More information

Bernstein-Durrmeyer operators with arbitrary weight functions

Bernstein-Durrmeyer operators with arbitrary weight functions Bernstein-Durrmeyer operators with arbitrary weight functions Elena E. Berdysheva German University of Technology in Oman Muscat, Sultanate of Oman Bernstein-Durrmeyer operators with arbitrary weight functions

More information

ENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS. Lecture 3

ENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS. Lecture 3 ENTIRE FUNCTIONS AND COMPLETENESS PROBLEMS A. POLTORATSKI Lecture 3 A version of the Heisenberg Uncertainty Principle formulated in terms of Harmonic Analysis claims that a non-zero measure (distribution)

More information

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

ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist

More information

Markov Sonin Gaussian rule for singular functions

Markov Sonin Gaussian rule for singular functions Journal of Computational and Applied Mathematics 169 (2004) 197 212 www.elsevier.com/locate/cam Markov Sonin Gaussian rule for singular functions G.Mastroianni, D.Occorsio Dipartimento di Matematica, Universita

More information

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

Both these computations follow immediately (and trivially) from the definitions. Finally, observe that if f L (R n ) then we have that. Lecture : One Parameter Maximal Functions and Covering Lemmas In this first lecture we start studying one of the basic and fundamental operators in harmonic analysis, the Hardy-Littlewood maximal function.

More information

Journal of Inequalities in Pure and Applied Mathematics

Journal of Inequalities in Pure and Applied Mathematics Journal of Inequalities in Pure and Applied Mathematics NEWER APPLICATIONS OF GENERALIZED MONOTONE SEQUENCES L. LEINDLER Bolyai Institute, University of Szeged, Aradi vértanúk tere 1 H-6720 Szeged, Hungary

More information

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

Hilbert Spaces. Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Hilbert Spaces Hilbert space is a vector space with some extra structure. We start with formal (axiomatic) definition of a vector space. Vector Space. Vector space, ν, over the field of complex numbers,

More information

Boot camp - Problem set

Boot camp - Problem set Boot camp - Problem set Luis Silvestre September 29, 2017 In the summer of 2017, I led an intensive study group with four undergraduate students at the University of Chicago (Matthew Correia, David Lind,

More information

On Multivariate Newton Interpolation at Discrete Leja Points

On Multivariate Newton Interpolation at Discrete Leja Points On Multivariate Newton Interpolation at Discrete Leja Points L. Bos 1, S. De Marchi 2, A. Sommariva 2, M. Vianello 2 September 25, 2011 Abstract The basic LU factorization with row pivoting, applied to

More information

HOMEOMORPHISMS OF BOUNDED VARIATION

HOMEOMORPHISMS OF BOUNDED VARIATION HOMEOMORPHISMS OF BOUNDED VARIATION STANISLAV HENCL, PEKKA KOSKELA AND JANI ONNINEN Abstract. We show that the inverse of a planar homeomorphism of bounded variation is also of bounded variation. In higher

More information

PART IV Spectral Methods

PART IV Spectral Methods PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,

More information

On Orthogonal Polynomials in Several Variables

On Orthogonal Polynomials in Several Variables Fields Institute Communications Volume 00, 0000 On Orthogonal Polynomials in Several Variables Yuan Xu Department of Mathematics University of Oregon Eugene, Oregon 97403-1222 Abstract. We report on the

More information

Recurrence Relations and Fast Algorithms

Recurrence Relations and Fast Algorithms Recurrence Relations and Fast Algorithms Mark Tygert Research Report YALEU/DCS/RR-343 December 29, 2005 Abstract We construct fast algorithms for decomposing into and reconstructing from linear combinations

More information

Introduction to Orthogonal Polynomials: Definition and basic properties

Introduction to Orthogonal Polynomials: Definition and basic properties Introduction to Orthogonal Polynomials: Definition and basic properties Prof. Dr. Mama Foupouagnigni African Institute for Mathematical Sciences, Limbe, Cameroon and Department of Mathematics, Higher Teachers

More information

Real Variables # 10 : Hilbert Spaces II

Real Variables # 10 : Hilbert Spaces II randon ehring Real Variables # 0 : Hilbert Spaces II Exercise 20 For any sequence {f n } in H with f n = for all n, there exists f H and a subsequence {f nk } such that for all g H, one has lim (f n k,

More information

arxiv: v2 [math.fa] 17 May 2016

arxiv: v2 [math.fa] 17 May 2016 ESTIMATES ON SINGULAR VALUES OF FUNCTIONS OF PERTURBED OPERATORS arxiv:1605.03931v2 [math.fa] 17 May 2016 QINBO LIU DEPARTMENT OF MATHEMATICS, MICHIGAN STATE UNIVERSITY EAST LANSING, MI 48824, USA Abstract.

More information

MATH 590: Meshfree Methods

MATH 590: Meshfree Methods MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu

More information

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

Interpolation. 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter Key References: Interpolation 1. Judd, K. Numerical Methods in Economics, Cambridge: MIT Press. Chapter 6. 2. Press, W. et. al. Numerical Recipes in C, Cambridge: Cambridge University Press. Chapter 3

More information

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

Free energy estimates for the two-dimensional Keller-Segel model Free energy estimates for the two-dimensional Keller-Segel model dolbeaul@ceremade.dauphine.fr CEREMADE CNRS & Université Paris-Dauphine in collaboration with A. Blanchet (CERMICS, ENPC & Ceremade) & B.

More information

Vectors in Function Spaces

Vectors in Function Spaces Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also

More information

Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials

Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Interpolation and Cubature at Geronimus Nodes Generated by Different Geronimus Polynomials Lawrence A. Harris Abstract. We extend the definition of Geronimus nodes to include pairs of real numbers where

More information

Regularity of Weak Solution to Parabolic Fractional p-laplacian

Regularity of Weak Solution to Parabolic Fractional p-laplacian Regularity of Weak Solution to Parabolic Fractional p-laplacian Lan Tang at BCAM Seminar July 18th, 2012 Table of contents 1 1. Introduction 1.1. Background 1.2. Some Classical Results for Local Case 2

More information

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

The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge Vladimir Kozlov (Linköping University, Sweden) 2010 joint work with A.Nazarov Lu t u a ij

More information

Outline of Fourier Series: Math 201B

Outline of Fourier Series: Math 201B Outline of Fourier Series: Math 201B February 24, 2011 1 Functions and convolutions 1.1 Periodic functions Periodic functions. Let = R/(2πZ) denote the circle, or onedimensional torus. A function f : C

More information

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

On the solution of integral equations of the first kind with singular kernels of Cauchy-type International Journal of Mathematics and Computer Science, 7(202), no. 2, 29 40 M CS On the solution of integral equations of the first kind with singular kernels of Cauchy-type G. E. Okecha, C. E. Onwukwe

More information

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

Entrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems September, 27 Solve exactly 6 out of the 8 problems. Prove by denition (in ɛ δ language) that f(x) = + x 2 is uniformly continuous in (, ). Is f(x) uniformly continuous in (, )? Prove your conclusion.

More information

Nonlinear aspects of Calderón-Zygmund theory

Nonlinear aspects of Calderón-Zygmund theory Ancona, June 7 2011 Overture: The standard CZ theory Consider the model case u = f in R n Overture: The standard CZ theory Consider the model case u = f in R n Then f L q implies D 2 u L q 1 < q < with

More information

Strong and Weak Weighted Convergence of Jacobi Series

Strong and Weak Weighted Convergence of Jacobi Series journal of approximation theory 88, 127 (1997) article no. AT9635 Strong and Weak Weighted Convergence of Jacobi Series R. A. Kerman* Department of Mathematics, Brock University, St. Catharines, Ontario,

More information

Class notes: Approximation

Class notes: Approximation Class notes: Approximation Introduction Vector spaces, linear independence, subspace The goal of Numerical Analysis is to compute approximations We want to approximate eg numbers in R or C vectors in R

More information

Interpolation & Polynomial Approximation. Hermite Interpolation I

Interpolation & Polynomial Approximation. Hermite Interpolation I Interpolation & Polynomial Approximation Hermite Interpolation I Numerical Analysis (th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University c 2011

More information

Nonlinear tensor product approximation

Nonlinear tensor product approximation ICERM; October 3, 2014 1 Introduction 2 3 4 Best multilinear approximation We are interested in approximation of a multivariate function f (x 1,..., x d ) by linear combinations of products u 1 (x 1 )

More information

FFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations

FFTs in Graphics and Vision. Homogenous Polynomials and Irreducible Representations FFTs in Graphics and Vision Homogenous Polynomials and Irreducible Representations 1 Outline The 2π Term in Assignment 1 Homogenous Polynomials Representations of Functions on the Unit-Circle Sub-Representations

More information

Hilbert Transforms and Orthonormal Expansions for Exponential Weights

Hilbert Transforms and Orthonormal Expansions for Exponential Weights Hilbert Transforms and Orthonormal Expansions for Exponential Weights S. B. Damelin Abstract. We establish the uniform boundedness of the weighted Hilbert transform for a general class of symmetric and

More information

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

Natural boundary and Zero distribution of random polynomials in smooth domains arxiv: v1 [math.pr] 2 Oct 2017 Natural boundary and Zero distribution of random polynomials in smooth domains arxiv:1710.00937v1 [math.pr] 2 Oct 2017 Igor Pritsker and Koushik Ramachandran Abstract We consider the zero distribution

More information

The heat equation for the Hermite operator on the Heisenberg group

The heat equation for the Hermite operator on the Heisenberg group Hokkaido Mathematical Journal Vol. 34 (2005) p. 393 404 The heat equation for the Hermite operator on the Heisenberg group M. W. Wong (Received August 5, 2003) Abstract. We give a formula for the one-parameter

More information

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

S. Lototsky and B.L. Rozovskii Center for Applied Mathematical Sciences University of Southern California, Los Angeles, CA RECURSIVE MULTIPLE WIENER INTEGRAL EXPANSION FOR NONLINEAR FILTERING OF DIFFUSION PROCESSES Published in: J. A. Goldstein, N. E. Gretsky, and J. J. Uhl (editors), Stochastic Processes and Functional Analysis,

More information

12.0 Properties of orthogonal polynomials

12.0 Properties of orthogonal polynomials 12.0 Properties of orthogonal polynomials In this section we study orthogonal polynomials to use them for the construction of quadrature formulas investigate projections on polynomial spaces and their

More information

Gaussian interval quadrature rule for exponential weights

Gaussian interval quadrature rule for exponential weights Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice

More information

NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS

NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS NEW BOUNDS FOR TRUNCATION-TYPE ERRORS ON REGULAR SAMPLING EXPANSIONS Nikolaos D. Atreas Department of Mathematics, Aristotle University of Thessaloniki, 54006, Greece, e-mail:natreas@auth.gr Abstract We

More information

Wavelets and modular inequalities in variable L p spaces

Wavelets and modular inequalities in variable L p spaces Wavelets and modular inequalities in variable L p spaces Mitsuo Izuki July 14, 2007 Abstract The aim of this paper is to characterize variable L p spaces L p( ) (R n ) using wavelets with proper smoothness

More information

Hardy-Littlewood maximal operator in weighted Lorentz spaces

Hardy-Littlewood maximal operator in weighted Lorentz spaces Hardy-Littlewood maximal operator in weighted Lorentz spaces Elona Agora IAM-CONICET Based on joint works with: J. Antezana, M. J. Carro and J. Soria Function Spaces, Differential Operators and Nonlinear

More information

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

Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation

More information

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

SELF-ADJOINTNESS OF SCHRÖDINGER-TYPE OPERATORS WITH SINGULAR POTENTIALS ON MANIFOLDS OF BOUNDED GEOMETRY Electronic Journal of Differential Equations, Vol. 2003(2003), No.??, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) SELF-ADJOINTNESS

More information

Optimal global rates of convergence for interpolation problems with random design

Optimal global rates of convergence for interpolation problems with random design Optimal global rates of convergence for interpolation problems with random design Michael Kohler 1 and Adam Krzyżak 2, 1 Fachbereich Mathematik, Technische Universität Darmstadt, Schlossgartenstr. 7, 64289

More information

Regularity of the p-poisson equation in the plane

Regularity of the p-poisson equation in the plane Regularity of the p-poisson equation in the plane Erik Lindgren Peter Lindqvist Department of Mathematical Sciences Norwegian University of Science and Technology NO-7491 Trondheim, Norway Abstract We

More information

Nonnegative linearization for little q-laguerre polynomials and Faber basis

Nonnegative linearization for little q-laguerre polynomials and Faber basis Journal of Computational and Applied Mathematics 199 (2007) 89 94 www.elsevier.com/locate/cam Nonnegative linearization for little q-laguerre polynomials and Faber basis Josef Obermaier a,, Ryszard Szwarc

More information

Supplement: Universal Self-Concordant Barrier Functions

Supplement: Universal Self-Concordant Barrier Functions IE 8534 1 Supplement: Universal Self-Concordant Barrier Functions IE 8534 2 Recall that a self-concordant barrier function for K is a barrier function satisfying 3 F (x)[h, h, h] 2( 2 F (x)[h, h]) 3/2,

More information

Spectral Gap and Concentration for Some Spherically Symmetric Probability Measures

Spectral Gap and Concentration for Some Spherically Symmetric Probability Measures Spectral Gap and Concentration for Some Spherically Symmetric Probability Measures S.G. Bobkov School of Mathematics, University of Minnesota, 127 Vincent Hall, 26 Church St. S.E., Minneapolis, MN 55455,

More information

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

DETERMINATION OF THE BLOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION DETERMINATION OF THE LOW-UP RATE FOR THE SEMILINEAR WAVE EQUATION y FRANK MERLE and HATEM ZAAG Abstract. In this paper, we find the optimal blow-up rate for the semilinear wave equation with a power nonlinearity.

More information

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

Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis

More information

Moreover this binary operation satisfies the following properties

Moreover this binary operation satisfies the following properties Contents 1 Algebraic structures 1 1.1 Group........................................... 1 1.1.1 Definitions and examples............................. 1 1.1.2 Subgroup.....................................

More information

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

LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011 LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic

More information

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:

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: 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 2 1 + u 2 2 = u 2. Geometric Meaning: u v = u v cos θ. u θ v 1 Reason: The opposite side is given by u v. u v 2 =

More information

Error formulas for divided difference expansions and numerical differentiation

Error formulas for divided difference expansions and numerical differentiation Error formulas for divided difference expansions and numerical differentiation Michael S. Floater Abstract: We derive an expression for the remainder in divided difference expansions and use it to give

More information

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

Scientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1 Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient

More information

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

LECTURE 16 GAUSS QUADRATURE In general for Newton-Cotes (equispaced interpolation points/ data points/ integration points/ nodes). CE 025 - Lecture 6 LECTURE 6 GAUSS QUADRATURE In general for ewton-cotes (equispaced interpolation points/ data points/ integration points/ nodes). x E x S fx dx hw' o f o + w' f + + w' f + E 84 f 0 f

More information

A trigonometric orthogonality with respect to a nonnegative Borel measure

A trigonometric orthogonality with respect to a nonnegative Borel measure Filomat 6:4 01), 689 696 DOI 10.98/FIL104689M Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A trigonometric orthogonality with

More information

Subelliptic mollifiers and a basic pointwise estimate of Poincaré type

Subelliptic mollifiers and a basic pointwise estimate of Poincaré type Math. Z. 226, 147 154 (1997) c Springer-Verlag 1997 Subelliptic mollifiers and a basic pointwise estimate of Poincaré type Luca Capogna, Donatella Danielli, Nicola Garofalo Department of Mathematics, Purdue

More information

Lecture 7: Kernels for Classification and Regression

Lecture 7: Kernels for Classification and Regression Lecture 7: Kernels for Classification and Regression CS 194-10, Fall 2011 Laurent El Ghaoui EECS Department UC Berkeley September 15, 2011 Outline Outline A linear regression problem Linear auto-regressive

More information

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

On the constant in the reverse Brunn-Minkowski inequality for p-convex balls. On the constant in the reverse Brunn-Minkowski inequality for p-convex balls. A.E. Litvak Abstract This note is devoted to the study of the dependence on p of the constant in the reverse Brunn-Minkowski

More information

A collocation method for solving some integral equations in distributions

A collocation method for solving some integral equations in distributions A collocation method for solving some integral equations in distributions Sapto W. Indratno Department of Mathematics Kansas State University, Manhattan, KS 66506-2602, USA sapto@math.ksu.edu A G Ramm

More information

ON A CLASS OF LINEAR POSITIVE BIVARIATE OPERATORS OF KING TYPE

ON A CLASS OF LINEAR POSITIVE BIVARIATE OPERATORS OF KING TYPE STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume LI, Number 4, December 2006 ON A CLASS OF LINEAR POSITIVE BIVARIATE OPERATORS OF KING TYPE OCTAVIAN AGRATINI Dedicated to Professor Gheorghe Coman at his

More information

Math The Laplacian. 1 Green s Identities, Fundamental Solution

Math The Laplacian. 1 Green s Identities, Fundamental Solution Math. 209 The Laplacian Green s Identities, Fundamental Solution Let be a bounded open set in R n, n 2, with smooth boundary. The fact that the boundary is smooth means that at each point x the external

More information

Second Order Elliptic PDE

Second Order Elliptic PDE Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic

More information

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

POTENTIAL LANDESMAN-LAZER TYPE CONDITIONS AND. 1. Introduction We investigate the existence of solutions for the nonlinear boundary-value problem Electronic Journal of Differential Equations, Vol. 25(25), No. 94, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) POTENTIAL

More information

Recent Results for Moving Least Squares Approximation

Recent Results for Moving Least Squares Approximation Recent Results for Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract. We describe two experiments recently conducted with the approximate moving least squares (MLS) approximation

More information

Probability and Measure

Probability and Measure Part II Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 84 Paper 4, Section II 26J Let (X, A) be a measurable space. Let T : X X be a measurable map, and µ a probability

More information

Singular Integrals. 1 Calderon-Zygmund decomposition

Singular Integrals. 1 Calderon-Zygmund decomposition Singular Integrals Analysis III Calderon-Zygmund decomposition Let f be an integrable function f dx 0, f = g + b with g Cα almost everywhere, with b

More information