In honour of Professor William Gear
|
|
- Ashlee Flowers
- 5 years ago
- Views:
Transcription
1 Functional Calculus and Numerical Analysis Michel Crouzeix Université de Rennes 1 ICNAAM 2011 In honour of Professor William Gear Halkidiki, September 2011
2 The context Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. We assume D(A) H, with dense imbedding, and the spectrum satisfies σ(a) C. Example 1. A square matrix. H = D(A) = C d, A C d,d.
3 The context Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. We assume D(A) H, with dense imbedding, and the spectrum satisfies σ(a) C. Example 2. Dirichlet Laplace operator. H = L 2 (Ω), V = H0 1 (Ω), a(u, v) = Ω u. v dx, D(A) = {u V ; Au H, Au, v = a(u, v), v V } = H 2 (Ω) H0 1 (Ω), Au = u.
4 The context Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. We assume D(A) H, with dense imbedding, and the spectrum satisfies σ(a) C. Example 3. Robin Laplace operator. H = L 2 (Ω), V = H 1 (Ω), a C( Ω; R). a(u, v) = Ω u. v dx + a u v dσ, D(A) = {u V ; Au H, Au, v = a(u, v), v V } = {u H 2 (Ω); n u + au = 0, on Ω}, Au = u. Ω
5 The context Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. We assume D(A) H, with dense imbedding, and the spectrum satisfies σ(a) C. Example 4. Complex Robin-Laplace operator. H = L 2 (Ω), V = H 1 (Ω), a, b C( Ω; R). a(u, v) = Ω u. v dx + (a+ib) u v dσ, D(A) = {u V ; Au H, Au, v = a(u, v), v V } = {u H 2 (Ω); n u + (a+ib)u = 0, on Ω}, Au = u. Ω
6 The self-adjoint case Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. If A is a self-adjoint operator, spectral theory is very efficient. We know that the spectrum σ(a) is real and we have r(a) = sup r(x), x σ(a) if the rational function r is bounded on σ(a).
7 The self-adjoint case Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. If A is a self-adjoint operator, spectral theory is very efficient. We know that the spectrum σ(a) is real and we have r(a) = sup r(x), x σ(a) if the rational function r is bounded on σ(a). By a density argument, this (in)equality allows to define f(a) L(H), for any function f that is continuous and bounded on σ(a). The previous inequality also holds with f instead of r.
8 The self-adjoint case Let us consider a closed linear operator A L(D(A), H) on a Hilbert space H. If A is a self-adjoint operator, spectral theory is very efficient. We know that the spectrum σ(a) is real and we have f(a) = sup f(x), x σ(a) for all continuous functions f that are bounded on σ(a). Furthermore the map f f(a), is a homomorphism from the algebra C b (σ(a)) into the algebra L(H).
9 Example of application Assume that we have a self-adjoint positive definite operator A L(D(A), H). We can define C(t) = cos(t A), and we are able to show that u(t) = C(t) u 0, with u 0 D(A), is the unique solution in C 2 (R, H) C 0 (R, D(A)) of u (t) + A u(t) = 0, t R, u(0) = u 0, u (0) = 0.
10 Continuity proof. Hyp : u 0 D(A), u(t) = C(t)u 0. With f 1 (z) = cos(t z), f 2 (z) = (cos((t+h) z) cos((t h) z))/z, we have, for z > 0, f 1 (z) 1 and f 2 (z) 2 ht. This allows to set C(t) = f 1 (A), and B(t, h) = f 2 (A), and we have C(t) 1, B(t, h) 2 ht.
11 Continuity proof. Hyp : u 0 D(A), u(t) = C(t)u 0. With f 1 (z) = cos(t z), f 2 (z) = (cos((t+h) z) cos((t h) z))/z, we have, for z > 0, f 1 (z) 1 and f 2 (z) 2 ht. This allows to set C(t) = f 1 (A), and B(t, h) = f 2 (A), and we have C(t) 1, B(t, h) 2 ht. Furthermore B(t, h) = (C(t + h) C(t h))a 1,
12 Continuity proof. Hyp : u 0 D(A), u(t) = C(t)u 0. With f 1 (z) = cos(t z), f 2 (z) = (cos((t+h) z) cos((t h) z))/z, we have, for z > 0, f 1 (z) 1 and f 2 (z) 2 ht. This allows to set C(t) = f 1 (A), and B(t, h) = f 2 (A), and we have C(t) 1, B(t, h) 2 ht. Furthermore B(t, h) = (C(t + h) C(t h))a 1, This yields u(t+h) u(t h) = B(t, h)au 0, and u(t+h) u(t h) 2 ht Au 0.
13 Continuity proof. Hyp : u 0 D(A), u(t) = C(t)u 0. With f 1 (z) = cos(t z), f 2 (z) = (cos((t+h) z) cos((t h) z))/z, we have, for z > 0, f 1 (z) 1 and f 2 (z) 2 ht. This allows to set C(t) = f 1 (A), and B(t, h) = f 2 (A), and we have C(t) 1, B(t, h) 2 ht. Furthermore B(t, h) = (C(t + h) C(t h))a 1, This yields u(t+h) u(t h) = B(t, h)au 0, and u(t+h) u(t h) 2 ht Au 0. This shows that u C(R; H), if u 0 D(A), but also, by density, if u 0 H.
14 Continuity proof. Hyp : u 0 D(A), u(t) = C(t)u 0. With f 1 (z) = cos(t z), f 2 (z) = (cos((t+h) z) cos((t h) z))/z, we have, for z > 0, f 1 (z) 1 and f 2 (z) 2 ht. This allows to set C(t) = f 1 (A), and B(t, h) = f 2 (A), and we have C(t) 1, B(t, h) 2 ht. Furthermore B(t, h) = (C(t + h) C(t h))a 1, This yields u(t+h) u(t h) = B(t, h)au 0, and u(t+h) u(t h) 2 ht Au 0. A similar proof gives that u C(R; D(A)), if u 0 D(A).
15 Differentiability proof. Hyp : u 0 D(A) and u(t) = C(t)u 0. We have cos(t t z)/z = 1/z (t s)cos(s z) ds. 0 Hence (homomorphism of algebra) C(t) A 1 = A 1 and u(t) = u 0 t 0 t 0 (t s)c(s) ds. (t s)au(s) ds.
16 Differentiability proof. Hyp : u 0 D(A) and u(t) = C(t)u 0. We have u(t) = u 0 t 0 (t s)au(s) ds. It follows u(0) = u 0, u (t) = t 0 Au(s) ds, u (0) = 0, u (t) = Au(t), and u C 2 (R; H).
17 Application 2. Heat equation Let A L(D(A), H) a self-adjoint, positive definite operator. We consider the heat equation u (t) + Au(t) = 0, u(0) = u 0 given in H, that we discretize in space by a finite elements method u h (t) + A hu h (t) = 0, u h (0) = u 0 h given in V h H, and then, in time, by the implicit scheme, t n = n t, u n+1 h u n h t + A h u n+1 h = 0.
18 Application 2. Heat equation Let A L(D(A), H) a self-adjoint, positive definite operator. We consider the heat equation u (t) + Au(t) = 0, u(0) = u 0 given in H, that we discretize in space by a finite elements method u h (t) + A hu h (t) = 0, u h (0) = u 0 h given in V h H, and then, in time, by the implicit scheme, t n = n t, u n+1 h u n h t + A h u n+1 h = 0, whence u n h = (r( ta h)) n u 0 h, with r(z) = (1 + z) 1.
19 Application 2. Heat equation Error estimate u(t n ) u n h u(t n) u h (t n ) + u h (t n ) u n h The first term u(t n ) u h (t n ) comes from the finite element method. We only consider the second u h (t n ) u n h = ( exp( n ta h ) r( ta h ) n) u 0 h. If A h is positive definite, (recall that r(x) = (1 + x) 1 ) exp( n ta h ) r( ta h ) n sup (e x ) n r(x) n C/n. x>0
20 Application 2. Heat equation Error estimate u(t n ) u n h u(t n) u h (t n ) + u h (t n ) u n h The first term u(t n ) u h (t n ) comes from the finite element method. We only consider the second u h (t n ) u n h = ( exp( n ta h ) r( ta h ) n) u 0 h. If A h is positive definite, (recall that r(x) = (1 + x) 1 ) exp( n ta h ) r( ta h ) n sup (e x ) n r(x) n C/n. x>0 Hence u(t n ) u n h u(t n) u h (t n ) + C t t n u 0 h.
21 Application 3. Convergence of GMRES. Let us consider a large linear system Ax = b, with b = 1. At the step n, GMRES algorithm builts x (n) Span{b, Ab,..., A n 1 b} element of the Krylov subspace spanned by b that minimizes the residual r (n) = b Ax (n).
22 Application 3. Convergence of GMRES. Let us consider a large linear system Ax = b, with b = 1. At the step n, GMRES algorithm builts x (n) Span{b, Ab,..., A n 1 b} element of the Krylov subspace spanned by b that minimizes the residual r (n) = b Ax (n). Therefore, we have r (n) min{ p(a) ;p(0) = 1, p P n }.
23 Application 3. Convergence of GMRES. Let us consider a large linear system Ax = b, with b = 1. At the step n, GMRES algorithm builts x (n) Span{b, Ab,..., A n 1 b} element of the Krylov subspace spanned by b that minimizes the residual r (n) = b Ax (n). If furthermore, A is self-adjoint with σ(a) [α, β], 0 < α < β, r (n) min{ p(a) ;p(0) = 1, p P n } min{ sup p(x) ;p(0) = 1, p P n }. x [α,β]
24 Application 3. Convergence of GMRES. Let us consider a large linear system Ax = b, with b = 1. At the step n, GMRES algorithm builts x (n) Span{b, Ab,..., A n 1 b} element of the Krylov subspace spanned by b that minimizes the residual r (n) = b Ax (n). If furthermore, A is self-adjoint with σ(a) [α, β], 0 < α < β, r (n) min{ p(a) ;p(0) = 1, p P n } min{ sup p(x) ;p(0) = 1, p P n } x [α,β] 2 β + α γ n + γ n, with γ =. β α
25 Spectral sets, J. von Neumann, 1951 A set X C, satisfying σ(a) X, is called a K-spectral set for the operator A if we have r(a) K sup z X r(z) for all rational functions that are bounded on X. If K = 1 we say that X is a spectral set.
26 Spectral sets A set X C, satisfying σ(a) X, is called a K-spectral set for the operator A if we have r(a) K sup z X r(z) for all rational functions that are bounded on X. Examples 1) If A is self-adjoint, or normal, then the spectrum σ(a) is a spectral set for A. 2) If A 1, then the unit disk { z 1} is a spectral set for A. 3) If Re Av, v 0, v H, then the right half-plane {Re z 0} is a spectral set for A.
27 Spectral sets A set X C, satisfying σ(a) X, is called a K-spectral set for the operator A if we have r(a) K sup z X r(z) Example 4. Let Ω be a domain that contains σ(a). If we may write the Cauchy formula r(a) = 2πi 1 Ω r(σ)(σ A) 1 dσ, then Ω is a K-spectral set for A with K = 1 2π Ω (σ A) 1 dσ.
28 There are many other ways to represent r(a) on an integral form : for instance by using Laplace transform (Hille-Phillips calculus), we can also use the Mellin transform, or the Poisson kernel, or the double layer potential... The following two lemmata may be useful in order to bound such integrals.
29 Lemma 1. Assume that dm(t) is a complex-valued measure, bounded on E, M(t) L(H), M(t) = M (t) 0, in E, r is a rational function bounded by 1 on E. Then E r(t)m(t) dm(t) E M(t) dm(t).
30 Lemma 1. Assume that dm(t) is a complex-valued measure, bounded on E, M(t) L(H), M(t) = M (t) 0, in E, r is a rational function bounded by 1 on E. Then E r(t)m(t) dm(t) Key of the proof. We have E M(t) dm(t) M(t)u, v M(t)u, u 1/2 M(t)v, v 1/2.
31 Lemma 2. Assume that dm(t) is a complex-valued measure, bounded on E, M(t), N(t) L(H), N(t) = N (t), in E, α > 0 Re M(t) = 1 2 (M(t) + M(t) ) N(t) α, in E, r is a rational function bounded by 1 on E. Then E r(t)(m(t)) 1 dm(t) E (N(t)) 1 dm(t).
32 Lemma 2. Assume that r is bounded by 1, α > 0, Re M(t) = 1 2 (M(t) + M(t) ) N(t) α, in E. Then E r(t)(m(t)) 1 dm(t) E (N(t)) 1 dm(t). Key of the proof. M(t) 1 u, v N(t) 1 u, u 1/2 N(t) 1 v, v 1/2.
33 The numerical range. Recall that the numerical range of a matrix A C d,d is defined by W(A) = { Av, v = v Av ; v C d, v = 1}. It is a closed convex subset of C, (Toeplitz & Hausdorff) that contains the spectrum σ(a).
34 The numerical range. Recall that the numerical range of a matrix A C d,d is defined by W(A) = { Av, v = v Av ; v C d, v = 1}. It is a closed convex subset of C, (Toeplitz & Hausdorff) that contains the spectrum σ(a). Let us consider a bounded convex domain Ω. We assume that W(A) Ω. If ν denotes the unit outward normal at the point σ Ω, then ν z Ω is equivalent to Re > 0, σ Ω. σ z
35 ν unit normal at σ Ω, s arclength on Ω, Π σ Ω ν σ dσ = i ν ds z Π σ is equivalent to Re ν σ z > 0.
36 The numerical range. Recall that the numerical range of a matrix A C d,d is defined by W(A) = { Av, v = v Av ; v C d, v = 1}. It is a closed convex subset of C, (Toeplitz & Hausdorff) that contains the spectrum σ(a). Let us consider a bounded convex domain Ω. We assume that W(A) Ω. If ν denotes the unit outward normal in the point σ Ω, then ν z Ω is equivalent to Re > 0, σ Ω. σ z Therefore W(A) Ω is equivalent to ν(σ A) 1 + ν( σ A ) 1 > 0, σ Ω.
37 The numerical range. The numerical range of an operator A L(D(A), H) is defined by W(A) = { Av, v = a(v, v);v D(A), v H = 1}. It is a convex subset of C, its closure contains the spectrum σ(a). If Ω is a convex domain, then W(A) Ω implies ν(σ A) 1 + ν( σ A ) 1 0, σ Ω.
38 Numerical Range and functional calculus An application of the two lemmata is the following theorem Theorem. If Ω is a closed convex domain of the complex plane with a conic boundary and if A is a closed operator with W(A) Ω, then Ω is a 3.2-spectral set for the operator A.
39 Numerical Range and functional calculus Theorem. If Ω is a closed convex domain of the complex plane with a conic boundary and if A is a closed operator with W(A) Ω, then Ω is a 3.2-spectral set for the operator A. Remark. Then f(a) is well defined for all functions that are holomorphic in the interior of Ω, bounded and continuous in Ω. Furthermore f(a) 3.2 sup z Ω f(z).
40 Numerical Range and functional calculus Application Corollary. If there exist a Hilbert space V H, V dense in H, and constants α > 0, λ, M, N R, such that then the problem α v 2 V Re Av, v + λ v 2 H M v 2 V Im Av, v N v H v V u (t) + A u(t) = 0, t R, u(0) = u 0 D(A), u (0) = v 0 D(A 1/2 ). has a unique solution in C 2 (R, H) C 0 (R, D(A)). Proof.
41 Numerical Range and functional calculus Application Corollary. If there exist a Hilbert space V H, V dense in H, and constants α > 0, λ, M, N R, such that then the problem α v 2 V Re Av, v + λ v 2 H M v 2 V Im Av, v N v H v V u (t) + A u(t) = 0, t R, u(0) = u 0 D(A), u (0) = v 0 D(A 1/2 ). has a unique solution in C 2 (R, H) C 0 (R, D(A)). Proof. If v H = 1 and Av, v = x+iy, then we have x+λ α v 2 V α N 2y2,.
42 Numerical Range and functional calculus Application Corollary. If there exist a Hilbert space V H, V dense in H, and constants α > 0, λ, M, N R, such that then the problem α v 2 V Re Av, v + λ v 2 H M v 2 V Im Av, v N v H v V u (t) + A u(t) = 0, t R, u(0) = u 0 D(A), u (0) = v 0 D(A 1/2 ). has a unique solution in C 2 (R, H) C 0 (R, D(A)). Proof. If v H = 1 and Av, v = x+iy, then we have x+λ α v 2 V α N 2y2, i.e. W(A) Ω := {x+iy ; x+λ α N 2y2 }.
43 Numerical Range and functional calculus Application Proof. If v H = 1 and Av, v = x+iy, then we have x+λ α v 2 V α N 2y2, i.e. W(A) Ω := {x+iy ; x+λ α N 2y2 }. It can be seen that z Ω implies Im z 1/2 ω, with ω = max( λ +, N/2 α). Thus cos(t z) is a holomorphic function in Ω with bound cos(t z) e ω t... The same proof as in the self-adjoint context works.
44 Numerical Range and functional calculus Theorem. If Ω is a closed convex domain of the complex plane with a conic boundary and if A is a closed operator with W(A) Ω, then Ω is a 3.2-spectral set for the operator A. Proof in the strip case. We assume Ω = {x+iy; y 2 1 0}. We write A = B+iC, with B and C self-adjoint. The condition W(A) Ω reads C 1 and we have to show r(a) 3.2, if the rational function r is bounded by 1 in the strip Ω.
45 We assume r(z) 1 for z Ω = {x+iy; y 2 1 0}, A = B+iC, with B and C self-adjoint, C 1. We want to show r(a) 3.2.
46 We assume r(z) 1 for z Ω = {x+iy; y 2 1 0}, A = B+iC, with C 1. We want to show r(a) 3.2. Proof. We set R 1 = r(σ) 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ, Ω 2πi R 2 = 1 2πi Ω r(σ)( σ A ) 1 d σ,
47 We assume r(z) 1 for z Ω = {x+iy; y 2 1 0}, A = B+iC, with C 1. We want to show r(a) 3.2. Proof. We set R 1 = r(σ) 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ, Ω 2πi R 2 = 1 2πi Ω r(σ)( σ A ) 1 d σ, and note that, from the Cauchy formula, r(a) = 1 2πi Ω r(σ)(σ A) 1 dσ = R 1 + R 2.
48 We assume r(z) 1 for z Ω = {x+iy; y 2 1 0}, A = B+iC, with C 1. We want to show r(a) 3.2. Proof. We set R 1 = r(σ) 1 ( (σ A) 1 dσ ( σ A ) 1 d σ ), Ω 2πi R 2 = 1 2πi Ω r(σ)( σ A ) 1 d σ, and note that, from the Cauchy formula, r(a) = 1 2πi Ω r(σ)(σ A) 1 dσ = R 1 + R 2. Therefore, it suffices to show that R 1 2 and R
49 Proof of R 1 2. We have R 1 = r(σ) 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ, Ω 2πi and it is easily seen that 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ 0. 2πi
50 Proof of R 1 2. We have R 1 = r(σ) 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ, Ω 2πi and it is easily seen that 1 ( (σ A) 1 dσ ( σ A ) 1 d σ ) 0. 2πi Thus it follows from Lemma 1 that 1 ( R 1 (σ A) 1 dσ ( σ A ) 1 d σ ) = 2. Ω 2πi
51 Proof of R We have R 2 = 1 2πi Ω r(σ)( σ A ) 1 d σ. We note that Ω = Ω + Ω, with Ω + = R+i and Ω = R i. Ω + 0 Ω
52 Proof of R We have R 2 = 1 2πi Ω r(σ)( σ A ) 1 d σ. We note that Ω = Ω + Ω, with Ω + = R+i and Ω = R i. Ω + 0 Ω On Ω +, we have σ = σ 2i and on Ω we have σ = σ+2i.
53 Proof of R We obtain R 2 = 1 2πi = 1 2πi = 1 2πi ( ( Ω r(σ)( σ A ) 1 d σ Ω + r(σ)( σ A ) 1 d σ + Ω + r(σ)(σ 2i A ) 1 dσ + Ω r(σ)( σ A ) 1 d σ ) Ω r(σ)(σ+2i A ) 1 dσ We note that r(σ)(σ±2i A ) 1 is analytic for σ Ω. ).
54 Ω + Γ 0 Γ + We can replace and Ω + f(σ) dσ by Ω f(σ) dσ by + + f(x) dx, f(x) dx.
55 Proof of R This yields R 2 = 1 2πi = 1 2πi = 2 π ( ( Ω + r(σ)(σ 2i A ) 1 dσ + r(x)(x 2i A ) 1 dx + r(x)((x A ) 2 + 4) 1 dx. Ω r(σ)(σ+2i A ) 1 dσ r(x)(x+2i A ) 1 dx ) )
56 Proof of R This yields R 2 = 2 π r(x)((x A ) 2 + 4) 1 dx. Recall that A = B+iC and C 1, thus Re ( (x A ) ) = (x B) 2 C (x B)
57 Proof of R This yields R 2 = 2 π r(x)((x A ) 2 + 4) 1 dx. Recall that A = B+iC and C 1, thus Re ( (x A ) ) = (x B) 2 C (x B) We may apply Lemma 2 and we get R 2 2 π 2 π sup λ R ((x B)2 + 3) 1 dx dy (y λ) =
58 Numerical Range and functional calculus Theorem. If Ω is a convex domain with a conic boundary and if A is a closed operator with W(A) Ω, then Ω is a 3.2-spectral set for the operator A. We have only considered here the case of a strip, but the proof in the general case is very similar...
59 A more general result is : The numerical range of an operator A W(A) := { Av, v ;v D(A), v = 1} is a K-spectral set for A with K
60 Application : Convergence of GMRES. We assume that W(A) Ω and 0 / Ω, Ω bounded convex domain. Let Φ be the Riemann mapping from C \ Ω onto the exterior of the unit disk, satisfying Φ(z) = γz + δ + O(z 1 ), as z, γ > 0. The Faber polynomial of degree n, F n, is defined by (Φ(z)) n = F n (z) + O(z 1 ), as z. Ω Φ
61 Application : Convergence of GMRES. We assume that W(A) Ω and 0 / Ω, Ω bounded convex domain. Let Φ be the Riemann mapping from C \ Ω onto the exterior of the unit disk, satisfying Φ(z) = γz + δ + O(z 1 ), as z, γ > 0. The Faber polynomial of degree n, F n, is defined by It is known that Similarly we have F n (A) = (Φ(z)) n = F n (z) + O(z 1 ), as z. F n (z) = Ω (Φ(σ))n 2πi 1 ( dσ σ z d σ σ z Ω Φ(σ)n 2πi 1 ( (σ A) 1 dσ ( σ A ) 1 d σ ). ).
62 Application : Convergence of GMRES. We have F n (A) = Ω Φ(σ)n 2πi 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ. Recall that ν = 1 i dσ ds, thus ( ) ( 1 2πi (σ A) 1 dσ ( σ A ) 1 d σ = 1 2π ν(σ A) 1 + ν( σ A ) 1 ) ds 0.
63 Application : Convergence of GMRES. We have F n (A) = Ω Φ(σ)n 2πi 1 ( ) (σ A) 1 dσ ( σ A ) 1 d σ. Recall that ν = 1 i dσ ds, thus ( ) ( 1 2πi (σ A) 1 dσ ( σ A ) 1 d σ = 1 2π ν(σ A) 1 + ν( σ A ) 1 ) ds 0. We can apply Lemma 1 F n (A) Ω ( ) 1 2πi (σ A) 1 dσ ( σ A ) 1 d σ = = 2.
64 Application : Convergence of GMRES. We consider a large linear system Ax = b, with b = 1. GMRES algorithm built at step n, an element x (n) that minimizes the residual r (n) = b Ax (n). Therefore, r (n) min{ p(a) ;p(0) = 1, p P n }.
65 Application : Convergence of GMRES. We consider a large linear system Ax = b, with b = 1. GMRES algorithm built at step n, an element x (n) that minimizes the residual r (n) = b Ax (n). Therefore, r (n) min{ p(a) ;p(0) = 1, p P n }. With the choice p(z) = F n (z)/f n (0), we get, r (n) F n (A) / F n (0) 2/ F n (0).
66 Application : Convergence of GMRES. We consider a large linear system Ax = b, with b = 1. GMRES algorithm built at step n, an element x (n) that minimizes the residual r (n) = b Ax (n). Therefore, r (n) min{ p(a) ;p(0) = 1, p P n }. With the choice p(z) = F n (z)/f n (0), we get, r (n) 2 F n (0) 2 γ n γ 1, with γ = Φ(0) > 1, by using an estimate of Kövari-Pommerenke F n (0) Φ(0) n 1/ Φ(0).
67 Application : Convergence of GMRES. Setting γ = Φ(0) > 1, we have shown that, if the matrix satisfies W(A) Ω, then we have the estimate min{ p(a) ;p(0) = 1, p P n } 2 F n (0) 2 γ n γ 1. This generalizes to the matrices an estimate due to Kövari and Pommerenke in the scalar case min{sup p(z) ;p(0) = 1, p P n } z Ω 2 F n (0) 2 γ n γ 1.
68 Application : Convergence of GMRES. Setting γ = Φ(0) > 1, we have shown that, if the matrix satisfies W(A) Ω, then we have the estimate min{ p(a) ;p(0) = 1, p P n } 2 F n (0) 2 γ n γ 1. This generalizes to the matrices an estimate due to Kövari and Pommerenke in the scalar case min{sup p(z) ;p(0) = 1, p P n } z Ω 2 F n (0) 2 γ n γ 1. The previous proof is due to Bernhard Beckermann.
69 B. Beckermann, Image numérique, GMRES et polynômes de Faber, CRAS, série 1, 340 (2006), M. Crouzeix, Numerical range and functional calculus in Hilbert space, J. Funct. Anal., 244 (2007), M. Crouzeix, A functional calculus based on the numerical range. Applications, Linear and Multil. Algebra, 56, 1, (2008), B.&F. Delyon, Generalization of Von Neumann s spectral sets and integral representation of operators, Bull. Soc. Math. France, 1, (1999), J. von Neumann, Eine Spektraltheorie für allgemeine Operatoren eines unitären Raumes, Math. Nachr. 4 (1951),
Banach Algebras where the Singular Elements are Removable Singularities
Banach Algebras where the Singular Elements are Removable Singularities Lawrence A. Harris Department of Mathematics, University of Kentucky, Lexington, Kentucky 40506-0027 E-mail: larry@ms.uky.edu Let
More information1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0
4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for
More informationInvariant subspaces for operators whose spectra are Carathéodory regions
Invariant subspaces for operators whose spectra are Carathéodory regions Jaewoong Kim and Woo Young Lee Abstract. In this paper it is shown that if an operator T satisfies p(t ) p σ(t ) for every polynomial
More informationHolomorphic functions which preserve holomorphic semigroups
Holomorphic functions which preserve holomorphic semigroups University of Oxford London Mathematical Society Regional Meeting Birmingham, 15 September 2016 Heat equation u t = xu (x Ω R d, t 0), u(t, x)
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationIn honor to Professor Francisco J. Lisbona. ZARAGOZA, September 2012
A MATRIX FACTORIZATION BASED ON NUMERICAL RADII Michel Crouzeix Université de Rennes 1 Workshop on Numerical Methods for Ordinary and Partial Di erential Equations and Applications In honor to Professor
More informationMath 123 Homework Assignment #2 Due Monday, April 21, 2008
Math 123 Homework Assignment #2 Due Monday, April 21, 2008 Part I: 1. Suppose that A is a C -algebra. (a) Suppose that e A satisfies xe = x for all x A. Show that e = e and that e = 1. Conclude that e
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationCriteria for Determining If A Subset is a Subspace
These notes closely follow the presentation of the material given in David C. Lay s textbook Linear Algebra and its Applications (3rd edition). These notes are intended primarily for in-class presentation
More informationCharacterization of half-radial matrices
Characterization of half-radial matrices Iveta Hnětynková, Petr Tichý Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague 8, Czech Republic Abstract Numerical radius r(a) is the
More informationComplex Analysis Problems
Complex Analysis Problems transcribed from the originals by William J. DeMeo October 2, 2008 Contents 99 November 2 2 2 200 November 26 4 3 2006 November 3 6 4 2007 April 6 7 5 2007 November 6 8 99 NOVEMBER
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Spring 2018 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Spring 208 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationSecond 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 informationMULTICENTRIC HOLOMORPHIC CALCULUS FOR n TUPLES OF COMMUTING OPERATORS
Adv. Oper. Theory https://doi.org/10.15352/aot.1804-1346 ISSN: 2538-225X (electronic) https://projecteuclid.org/aot MULTICENTRIC HOLOMORPHIC CALCULUS FOR n TUPLES OF COMMUTING OPERATORS DIANA ANDREI Communicated
More informationIntroduction to Index Theory. Elmar Schrohe Institut für Analysis
Introduction to Index Theory Elmar Schrohe Institut für Analysis Basics Background In analysis and pde, you want to solve equations. In good cases: Linearize, end up with Au = f, where A L(E, F ) is a
More informationClassical stuff - title to be changed later
CHAPTER 1 Classical stuff - title to be changed later 1. Positive Definite Kernels To start with something simple and elegant, we choose positive definite kernels which appear at every corner in functional
More informationINTRODUCTION TO REAL ANALYTIC GEOMETRY
INTRODUCTION TO REAL ANALYTIC GEOMETRY KRZYSZTOF KURDYKA 1. Analytic functions in several variables 1.1. Summable families. Let (E, ) be a normed space over the field R or C, dim E
More informationElements of linear algebra
Elements of linear algebra Elements of linear algebra A vector space S is a set (numbers, vectors, functions) which has addition and scalar multiplication defined, so that the linear combination c 1 v
More informationAlgebra Exam Topics. Updated August 2017
Algebra Exam Topics Updated August 2017 Starting Fall 2017, the Masters Algebra Exam will have 14 questions. Of these students will answer the first 8 questions from Topics 1, 2, and 3. They then have
More informationMATH COURSE NOTES - CLASS MEETING # Introduction to PDEs, Fall 2011 Professor: Jared Speck
MATH 8.52 COURSE NOTES - CLASS MEETING # 6 8.52 Introduction to PDEs, Fall 20 Professor: Jared Speck Class Meeting # 6: Laplace s and Poisson s Equations We will now study the Laplace and Poisson equations
More informationRecall that any inner product space V has an associated norm defined by
Hilbert Spaces Recall that any inner product space V has an associated norm defined by v = v v. Thus an inner product space can be viewed as a special kind of normed vector space. In particular every inner
More informationUNBOUNDED OPERATORS ON HILBERT SPACES. Let X and Y be normed linear spaces, and suppose A : X Y is a linear map.
UNBOUNDED OPERATORS ON HILBERT SPACES EFTON PARK Let X and Y be normed linear spaces, and suppose A : X Y is a linear map. Define { } Ax A op = sup x : x 0 = { Ax : x 1} = { Ax : x = 1} If A
More informationALGEBRAIC PROPERTIES OF OPERATOR ROOTS OF POLYNOMIALS
ALGEBRAIC PROPERTIES OF OPERATOR ROOTS OF POLYNOMIALS TRIEU LE Abstract. Properties of m-selfadjoint and m-isometric operators have been investigated by several researchers. Particularly interesting to
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More informationProblem 1A. Calculus. Problem 3A. Real analysis. f(x) = 0 x = 0.
Problem A. Calculus Find the length of the spiral given in polar coordinates by r = e θ, < θ 0. Solution: The length is 0 θ= ds where by Pythagoras ds = dr 2 + (rdθ) 2 = dθ 2e θ, so the length is 0 2e
More informationFunctional Analysis I
Functional Analysis I Course Notes by Stefan Richter Transcribed and Annotated by Gregory Zitelli Polar Decomposition Definition. An operator W B(H) is called a partial isometry if W x = X for all x (ker
More informationSpectral theory for linear operators on L 1 or C(K) spaces
Spectral theory for linear operators on L 1 or C(K) spaces Ian Doust, Florence Lancien, and Gilles Lancien Abstract It is known that on a Hilbert space, the sum of a real scalar-type operator and a commuting
More informationarxiv: v4 [math.fa] 19 Apr 2010
SEMIPROJECTIVITY OF UNIVERSAL C*-ALGEBRAS GENERATED BY ALGEBRAIC ELEMENTS TATIANA SHULMAN arxiv:0810.2497v4 [math.fa] 19 Apr 2010 Abstract. Let p be a polynomial in one variable whose roots either all
More information1. Partial Fraction Expansion All the polynomials in this note are assumed to be complex polynomials.
Partial Fraction Expansion All the polynomials in this note are assumed to be complex polynomials A rational function / is a quotient of two polynomials P, Q with 0 By Fundamental Theorem of algebra, =
More informationCompact symetric bilinear forms
Compact symetric bilinear forms Mihai Mathematics Department UC Santa Barbara IWOTA 2006 IWOTA 2006 Compact forms [1] joint work with: J. Danciger (Stanford) S. Garcia (Pomona
More informationContractivity of the H -calculus and Blaschke products
Contractivity of the H -calculus and Blaschke products Christoph Kriegler and Lutz Weis Abstract. It is well known that a densely defined operator A on a Hilbert space is accretive if and only if A has
More informationNovember 18, 2013 ANALYTIC FUNCTIONAL CALCULUS
November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral
More informationSolutions to Complex Analysis Prelims Ben Strasser
Solutions to Complex Analysis Prelims Ben Strasser In preparation for the complex analysis prelim, I typed up solutions to some old exams. This document includes complete solutions to both exams in 23,
More informationMULTICENTRIC CALCULUS AND THE RIESZ PROJECTION
JOURNAL OF NUMERICAL ANALYSIS AND APPROXIMATION THEORY J. Numer. Anal. Approx. Theory, vol. 44 (2015) no. 2, pp. 127 145 ictp.acad.ro/jnaat MULTICENTRIC CALCULUS AND THE RIESZ PROJECTION DIANA APETREI
More informationChapter 2: Linear Independence and Bases
MATH20300: Linear Algebra 2 (2016 Chapter 2: Linear Independence and Bases 1 Linear Combinations and Spans Example 11 Consider the vector v (1, 1 R 2 What is the smallest subspace of (the real vector space
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationComplex Analysis Important Concepts
Complex Analysis Important Concepts Travis Askham April 1, 2012 Contents 1 Complex Differentiation 2 1.1 Definition and Characterization.............................. 2 1.2 Examples..........................................
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationKey words. GMRES method, convergence bounds, worst-case GMRES, ideal GMRES, field of values
THE FIELD OF VALUES BOUNDS ON IDEAL GMRES JÖRG LIESEN AND PETR TICHÝ 27.03.2018) Abstract. A widely known result of Elman, and its improvements due to Starke, Eiermann and Ernst, gives a bound on the worst-case
More informationCONSERVATIVE DILATIONS OF DISSIPATIVE N-D SYSTEMS: THE COMMUTATIVE AND NON-COMMUTATIVE SETTINGS. Dmitry S. Kalyuzhnyĭ-Verbovetzkiĭ
CONSERVATIVE DILATIONS OF DISSIPATIVE N-D SYSTEMS: THE COMMUTATIVE AND NON-COMMUTATIVE SETTINGS Dmitry S. Kalyuzhnyĭ-Verbovetzkiĭ Department of Mathematics Ben-Gurion University of the Negev P.O. Box 653
More informationMidterm solutions. (50 points) 2 (10 points) 3 (10 points) 4 (10 points) 5 (10 points)
Midterm solutions Advanced Linear Algebra (Math 340) Instructor: Jarod Alper April 26, 2017 Name: } {{ } Read all of the following information before starting the exam: You may not consult any outside
More informationMath 220A - Fall Final Exam Solutions
Math 22A - Fall 216 - Final Exam Solutions Problem 1. Let f be an entire function and let n 2. Show that there exists an entire function g with g n = f if and only if the orders of all zeroes of f are
More informationSOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION. Dagmar Medková
29 Kragujevac J. Math. 31 (2008) 29 42. SOLUTION OF THE DIRICHLET PROBLEM WITH L p BOUNDARY CONDITION Dagmar Medková Czech Technical University, Faculty of Mechanical Engineering, Department of Technical
More informationTrotter s product formula for projections
Trotter s product formula for projections Máté Matolcsi, Roman Shvydkoy February, 2002 Abstract The aim of this paper is to examine the convergence of Trotter s product formula when one of the C 0-semigroups
More informationNumerical Range in C*-Algebras
Journal of Mathematical Extension Vol. 6, No. 2, (2012), 91-98 Numerical Range in C*-Algebras M. T. Heydari Yasouj University Abstract. Let A be a C*-algebra with unit 1 and let S be the state space of
More informationA Characterization of Group Generators on Hilbert Spaces and the H -calculus
Semigroup Forum Vol. 66 (2003) 288 304 c 2002 Springer-Verlag New York Inc. DOI: 10.1007/s002330010165 ESEACH ATICLE A Characterization of Group Generators on Hilbert Spaces and the H -calculus Markus
More informationINTEGRATION WORKSHOP 2003 COMPLEX ANALYSIS EXERCISES
INTEGRATION WORKSHOP 23 COMPLEX ANALYSIS EXERCISES DOUGLAS ULMER 1. Meromorphic functions on the Riemann sphere It s often useful to allow functions to take the value. This exercise outlines one way to
More informationThe second dual of a C*-algebra
The second dual of a C*-algebra The main goal is to expose a short proof of the result of Z Takeda, [Proc Japan Acad 30, (1954), 90 95], that the second dual of a C*-algebra is, in effect, the von Neumann
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
More informationA von Neumann Wold decomposition of two-isometries
Acta Sci. Math. (Szeged) 70 (2004), 715 726 A von Neumann Wold decomposition of two-isometries Anders Olofsson Communicated by L. Kérchy Abstract. We present a von Neumann Wold decomposition of a twoisometric
More informationUNIFORM APPROXIMATION BY b HARMONIC FUNCTIONS
UNIFORM APPROXIMATION BY b HARMONIC FUNCTIONS JOHN T. ANDERSON Abstract. The Mergelyan and Ahlfors-Beurling estimates for the Cauchy transform give quantitative information on uniform approximation by
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More information. Then g is holomorphic and bounded in U. So z 0 is a removable singularity of g. Since f(z) = w 0 + 1
Now we describe the behavior of f near an isolated singularity of each kind. We will always assume that z 0 is a singularity of f, and f is holomorphic on D(z 0, r) \ {z 0 }. Theorem 4.2.. z 0 is a removable
More informationCOMMUTATIVITY OF OPERATORS
COMMUTATIVITY OF OPERATORS MASARU NAGISA, MAKOTO UEDA, AND SHUHEI WADA Abstract. For two bounded positive linear operators a, b on a Hilbert space, we give conditions which imply the commutativity of a,
More informationComplex Analysis Qual Sheet
Complex Analysis Qual Sheet Robert Won Tricks and traps. traps. Basically all complex analysis qualifying exams are collections of tricks and - Jim Agler Useful facts. e z = 2. sin z = n=0 3. cos z = z
More informationOPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS
PORTUGALIAE MATHEMATICA Vol. 59 Fasc. 2 2002 Nova Série OPTIMAL CONTROL AND STRANGE TERM FOR A STOKES PROBLEM IN PERFORATED DOMAINS J. Saint Jean Paulin and H. Zoubairi Abstract: We study a problem of
More informationGENERATORS WITH INTERIOR DEGENERACY ON SPACES OF L 2 TYPE
Electronic Journal of Differential Equations, Vol. 22 (22), No. 89, pp. 3. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GENERATORS WITH INTERIOR
More informationReal Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi
Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationSPECTRAL THEORY EVAN JENKINS
SPECTRAL THEORY EVAN JENKINS Abstract. These are notes from two lectures given in MATH 27200, Basic Functional Analysis, at the University of Chicago in March 2010. The proof of the spectral theorem for
More information1. The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions.
Complex Analysis Qualifying Examination 1 The COMPLEX PLANE AND ELEMENTARY FUNCTIONS: Complex numbers; stereographic projection; simple and multiple connectivity, elementary functions 2 ANALYTIC FUNCTIONS:
More informationMATH 5640: Fourier Series
MATH 564: Fourier Series Hung Phan, UMass Lowell September, 8 Power Series A power series in the variable x is a series of the form a + a x + a x + = where the coefficients a, a,... are real or complex
More informationUNCERTAINTY PRINCIPLES FOR THE FOCK SPACE
UNCERTAINTY PRINCIPLES FOR THE FOCK SPACE KEHE ZHU ABSTRACT. We prove several versions of the uncertainty principle for the Fock space F 2 in the complex plane. In particular, for any unit vector f in
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More informationMH 7500 THEOREMS. (iii) A = A; (iv) A B = A B. Theorem 5. If {A α : α Λ} is any collection of subsets of a space X, then
MH 7500 THEOREMS Definition. A topological space is an ordered pair (X, T ), where X is a set and T is a collection of subsets of X such that (i) T and X T ; (ii) U V T whenever U, V T ; (iii) U T whenever
More informationSalah Mecheri. Communicated by S. L. Troyanski
Serdica Math J 26 (2000), 119-126 ON MINIMIZING S (AX XB) p p Salah Mecheri Communicated by S L Troyanski Abstract In this paper, we minimize the map F p (X)= S (AX XB) p p, where the pair (A, B) has the
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
More informationSHARPENED FORMS OF A VON NEUMANN
MATH. SCAND. 102 (2008), 265 282 SHARPENED FORMS OF A VON NEUMANN INEQUALITY FOR ρ-contractions GILLES CASSIER NICOLAE SUCIU Abstract The purpose of this paper is to give some sharpened forms of the von
More informationGeometry and the Kato square root problem
Geometry and the Kato square root problem Lashi Bandara Centre for Mathematics and its Applications Australian National University 7 June 2013 Geometric Analysis Seminar University of Wollongong Lashi
More informationTWO REMARKS ABOUT NILPOTENT OPERATORS OF ORDER TWO
TWO REMARKS ABOUT NILPOTENT OPERATORS OF ORDER TWO STEPHAN RAMON GARCIA, BOB LUTZ, AND DAN TIMOTIN Abstract. We present two novel results about Hilbert space operators which are nilpotent of order two.
More informationHille-Yosida Theorem and some Applications
Hille-Yosida Theorem and some Applications Apratim De Supervisor: Professor Gheorghe Moroșanu Submitted to: Department of Mathematics and its Applications Central European University Budapest, Hungary
More informationSEMICROSSED PRODUCTS OF THE DISK ALGEBRA
SEMICROSSED PRODUCTS OF THE DISK ALGEBRA KENNETH R. DAVIDSON AND ELIAS G. KATSOULIS Abstract. If α is the endomorphism of the disk algebra, A(D), induced by composition with a finite Blaschke product b,
More informationarxiv: v1 [math.na] 26 Jul 2013
Decay properties for functions of matrices over C -algebras Michele Benzi a,1, Paola Boito b arxiv:137.7119v1 [math.na] 26 Jul 213 a Department of Mathematics and Computer Science, Emory University, Atlanta,
More information5.3 The Upper Half Plane
Remark. Combining Schwarz Lemma with the map g α, we can obtain some inequalities of analytic maps f : D D. For example, if z D and w = f(z) D, then the composition h := g w f g z satisfies the condition
More informationMeans of unitaries, conjugations, and the Friedrichs operator
J. Math. Anal. Appl. 335 (2007) 941 947 www.elsevier.com/locate/jmaa Means of unitaries, conjugations, and the Friedrichs operator Stephan Ramon Garcia Department of Mathematics, Pomona College, Claremont,
More informationA functional model for commuting pairs of contractions and the symmetrized bidisc
A functional model for commuting pairs of contractions and the symmetrized bidisc Nicholas Young Leeds and Newcastle Universities Lecture 2 The symmetrized bidisc Γ and Γ-contractions St Petersburg, June
More informationObstacle problems and isotonicity
Obstacle problems and isotonicity Thomas I. Seidman Revised version for NA-TMA: NA-D-06-00007R1+ [June 6, 2006] Abstract For variational inequalities of an abstract obstacle type, a comparison principle
More informationAn introduction to some aspects of functional analysis
An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms
More informationTheorem 5.3. Let E/F, E = F (u), be a simple field extension. Then u is algebraic if and only if E/F is finite. In this case, [E : F ] = deg f u.
5. Fields 5.1. Field extensions. Let F E be a subfield of the field E. We also describe this situation by saying that E is an extension field of F, and we write E/F to express this fact. If E/F is a field
More informationMATH 104 : Final Exam
MATH 104 : Final Exam 10 May, 2017 Name: You have 3 hours to answer the questions. You are allowed one page (front and back) worth of notes. The page should not be larger than a standard US letter size.
More informationu xx + u yy = 0. (5.1)
Chapter 5 Laplace Equation The following equation is called Laplace equation in two independent variables x, y: The non-homogeneous problem u xx + u yy =. (5.1) u xx + u yy = F, (5.) where F is a function
More information1 Holomorphic functions
Robert Oeckl CA NOTES 1 15/09/2009 1 1 Holomorphic functions 11 The complex derivative The basic objects of complex analysis are the holomorphic functions These are functions that posses a complex derivative
More information2.2. OPERATOR ALGEBRA 19. If S is a subset of E, then the set
2.2. OPERATOR ALGEBRA 19 2.2 Operator Algebra 2.2.1 Algebra of Operators on a Vector Space A linear operator on a vector space E is a mapping L : E E satisfying the condition u, v E, a R, L(u + v) = L(u)
More informationLecture 1 The complex plane. z ± w z + w.
Lecture 1 The complex plane Exercise 1.1. Show that the modulus obeys the triangle inequality z ± w z + w. This allows us to make the complex plane into a metric space, and thus to introduce topological
More informationSpectral theorems for bounded self-adjoint operators on a Hilbert space
Chapter 10 Spectral theorems for bounded self-adjoint operators on a Hilbert space Let H be a Hilbert space. For a bounded operator A : H H its Hilbert space adjoint is an operator A : H H such that Ax,
More informationPart II Probability and Measure
Part II Probability and Measure Theorems Based on lectures by J. Miller Notes taken by Dexter Chua Michaelmas 2016 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More information1. Introduction Since the pioneering work by Leray [3] in 1934, there have been several studies on solutions of Navier-Stokes equations
Math. Res. Lett. 13 (6, no. 3, 455 461 c International Press 6 NAVIER-STOKES EQUATIONS IN ARBITRARY DOMAINS : THE FUJITA-KATO SCHEME Sylvie Monniaux Abstract. Navier-Stokes equations are investigated in
More informationComplex Analysis Qualifying Exam Solutions
Complex Analysis Qualifying Exam Solutions May, 04 Part.. Let log z be the principal branch of the logarithm defined on G = {z C z (, 0]}. Show that if t > 0, then the equation log z = t has exactly one
More informationReal Analysis Notes. Thomas Goller
Real Analysis Notes Thomas Goller September 4, 2011 Contents 1 Abstract Measure Spaces 2 1.1 Basic Definitions........................... 2 1.2 Measurable Functions........................ 2 1.3 Integration..............................
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationLINEAR FRACTIONAL COMPOSITION OPERATORS ON H 2
J Integral Equations and Operator Theory (988, 5 60 LINEAR FRACTIONAL COMPOSITION OPERATORS ON H 2 CARL C COWEN Abstract If ϕ is an analytic function mapping the unit disk D into itself, the composition
More informationA Brief Introduction to Functional Analysis
A Brief Introduction to Functional Analysis Sungwook Lee Department of Mathematics University of Southern Mississippi sunglee@usm.edu July 5, 2007 Definition 1. An algebra A is a vector space over C with
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationarxiv: v1 [math.cv] 21 Jun 2016
B. N. Khabibullin, N. R. Tamindarova SUBHARMONIC TEST FUNCTIONS AND THE DISTRIBUTION OF ZERO SETS OF HOLOMORPHIC FUNCTIONS arxiv:1606.06714v1 [math.cv] 21 Jun 2016 Abstract. Let m, n 1 are integers and
More informationLocally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem
56 Chapter 7 Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem Recall that C(X) is not a normed linear space when X is not compact. On the other hand we could use semi
More informationFunctional Analysis. Franck Sueur Metric spaces Definitions Completeness Compactness Separability...
Functional Analysis Franck Sueur 2018-2019 Contents 1 Metric spaces 1 1.1 Definitions........................................ 1 1.2 Completeness...................................... 3 1.3 Compactness......................................
More informationSpectral theory for compact operators on Banach spaces
68 Chapter 9 Spectral theory for compact operators on Banach spaces Recall that a subset S of a metric space X is precompact if its closure is compact, or equivalently every sequence contains a Cauchy
More informationMATH 225 Summer 2005 Linear Algebra II Solutions to Assignment 1 Due: Wednesday July 13, 2005
MATH 225 Summer 25 Linear Algebra II Solutions to Assignment 1 Due: Wednesday July 13, 25 Department of Mathematical and Statistical Sciences University of Alberta Question 1. [p 224. #2] The set of all
More information