Nonlinear Eigenvalue Problems: Theory and Applications
|
|
- Leo West
- 6 years ago
- Views:
Transcription
1 Nonlinear Eigenvalue Problems: Theory and Applications David Bindel 1 Department of Computer Science Cornell University 12 January Joint work with Amanda Hood Berkeley 1 / 37
2 The Computational Science & Engineering Picture Application Analysis Computation MEMS Smart grids Networks Systems Linear algebra Approximation theory Symmetry + structure Optimization HPC / cloud Simulators Solvers Frameworks Berkeley 2 / 37
3 Why eigenvalues? The polynomial connection The optimization connection The approximation connection The Fourier/quadrature/special function connection The dynamics connection Berkeley 3 / 37
4 Why nonlinear eigenvalues? y Ay = y + By + Ky = y Ay By(t 1) = T (d/dt)y = y=e λt v (λi A)v = y=e λt v (λ 2 I + λb + K)v = y=e λt v (λi A Be λ )v = y=e λt v T (λ)v = Higher-order ODEs Dynamic element formulations Delay differential equations Boundary integral equation eigenproblems Radiation boundary conditions Berkeley 4 / 37
5 Big and little y + By + Ky = (λ 2 I + λb + K)u = v=y [ ] v + y ( v=λu λi + [ B K I [ B K I ] [ ] v = y ]) [ ] v = u Tradeoff: more variables = more linear. Berkeley 5 / 37
6 The big picture T (λ)v =, v. where T : Ω C n n analytic, Ω C simply connected Regularity: det(t ) Nonlinear spectrum: Λ(T ) = {z Ω : T (z) singular}. What do we want? Qualitative information (e.g. no eigenvalues in right half plane) Error bounds on computed/estimated eigenvalues Assurances that we know all the eigenvalues in some region Berkeley 6 / 37
7 Standard solver strategies A common approach: T (λ)v =, v. 1 Approximate T locally (linear, rational, etc) 2 Solve approximate problem. 3 Repeat as needed. How should we choose solver parameters? What about global behavior? What can we trust? What might we miss? Berkeley 7 / 37
8 Analyticity to the rescue T (λ)v =, v. where T : Ω C n n analytic, Ω C simply connected Regularity: det(t ) Nonlinear spectrum: Λ(T ) = {z Ω : T (z) singular}. Goal: Use analyticity to compare and to count Berkeley 8 / 37
9 Winding and Cauchy s argument principle Γ Ω W Γ (f ) = 1 f (z) 2πi Γ f (z) dz = # zeros # poles Berkeley 9 / 37
10 Winding, Rouché, and Gohberg-Sigal f f + g Ω Γ s Γ Analytic f, g : Ω C T, E : Ω C n n Winding # 1 f (z) 2πi Γ f (z) dz tr ( 1 2πi Γ T (z) 1 T (z) dz ) Theorem Rouché (1862): Gohberg-Sigal (1971): g < f on Γ = T 1 E < 1 on Γ = same # zeros of f, f + g same # eigs of T, T + E Berkeley 1 / 37
11 Comparing NEPs Suppose T, E : Ω C n n analytic Γ Ω a simple closed contour T (z) + se(z) nonsingular s [, 1], z Γ Then T and T + E have the same number of eigenvalues inside Γ. Pf: Constant winding number around Γ. Berkeley 11 / 37
12 Nonlinear pseudospectra Λ ɛ (T ) {z Ω : T (z) 1 > ɛ 1 } Berkeley 12 / 37
13 Pseudospectral comparison U ɛ Ω Ω ɛ E analytic, E(z) < ɛ on Ω ɛ. Then Λ(T + E) Ω ɛ Λ ɛ (T ) Ω ɛ Also, if U ɛ a component of Λ ɛ and Ū ɛ Ω ɛ, then Λ(T + E) U ɛ = Λ(T ) U ɛ Berkeley 13 / 37
14 Pseudospectral comparison U ɛ Ω Ω ɛ Most useful when T is linear Even then, can be expensive to compute! What about related tools? Berkeley 14 / 37
15 The Gershgorin picture (linear case) A = D + F, D = diag(d i ), ρ i = j f ij ρ 1 ρ 2 ρ 3 d 1 d 2 d 3 Berkeley 15 / 37
16 Gershgorin (+ɛ) Write A = D + F, D = diag(d 1,..., d n ). Gershgorin disks are: G i = z C : z d i f ij. j Useful facts: Pf: Spectrum of A lies in m i=1 G i i I G i disjoint from other disks = contains I eigenvalues. A zi strictly diagonally dominant outside m i=1 G i. Eigenvalues of D sf, s 1, are continuous. Berkeley 16 / 37
17 Nonlinear Gershgorin Write T (z) = D(z) + F (z). Gershgorin regions are G i = z C : d i(z) f ij (z). j Useful facts: Spectrum of T lies in m i=1 G i Bdd connected component of m i=1 G i strictly in Ω = same number of eigs of D and T in component = at least one eig per component of G i involved Pf: Strict diag dominance test + continuity of eigs Berkeley 17 / 37
18 Mini-example: Counting contributions G 1 2 G 1 1 G 1 2 z 1 T (z) = z Berkeley 18 / 37
19 Mini-example: Domain boundaries T (z) = [ ] z.2 z z 1 Ω = C (, ] G 1 1 G 1 2 det(d(z)) =( z.1 i.99) ( z.1 + i.99) det(t (z)) =( z +.1 i.99) D has two eigenvalues in Ω; T hides both eigenvalues behind a branch cut. ( z i.99) Berkeley 19 / 37
20 Example I: Hadeler T (z) = (e z 1)B + z 2 A αi, A, B R Berkeley 2 / 37
21 Comparison to simplified problem Bauer-Fike idea: apply a similarity! T (z) = (e z 1)B + z 2 A αi T (z) = U T T (z)u = (e z 1)D B + z 2 I αe = D(z) αe G i = {z : β i (e z 1) + z 2 < ρ i }. Berkeley 21 / 37
22 Gershgorin regions Berkeley 22 / 37
23 A different comparison Approximate e z 1 by a Chebyshev interpolant: T (z) = (e z 1)B + z 2 A αi T (z) = q(z)b + z 2 A α T (z) = T (z) + r(z)b Linearize T and transform both: T (z) D C zi T (z) D C zi + r(z)e Restrict to Ω ɛ = {z : r(z) < ɛ}: G i Ĝ i = {z : z µ i < ρ i ɛ}, ρ i = j e ij Berkeley 23 / 37
24 Spectrum of T Berkeley 24 / 37
25 Ĝ i for ɛ < Berkeley 25 / 37
26 Ĝ i for ɛ = Berkeley 25 / 37
27 Ĝ i for ɛ = Berkeley 25 / 37
28 Example II: Resonance problem V a b ψ() = ( d 2 ) dx 2 + V λ ψ = on (, b), ψ (b) = i λψ(b), Berkeley 26 / 37
29 Reduction via shooting V R a (λ) a b R ab (λ) R a (λ) R ab (λ) ψ() =, [ ] ψ() ψ = () [ ] ψ(b) ψ = (b) [ ] ψ(a) ψ, (a) [ ] ψ(b) ψ, (b) ψ (b) = i λψ(b) Berkeley 27 / 37
30 Reduction via shooting First-order form: du dx = [ ] 1 u, where u(x) V λ [ ] ψ(x) ψ. (x) On region (c, d) where V is constant: ( [ ]) 1 u(d) = R cd (λ)u(c), R cd (λ) = exp (d c) V λ Reduce resonance problem to 6D NEP: R a (λ) I T (λ)u all [ R ab (λ) I u() ] [ ] 1 i u(a) =. u(b) λ 1 Berkeley 28 / 37
31 Expansion via rational approximation Consider the equation ([ A B C D ] λi ) [ ] u = v Partial Gaussian elimination gives the spectral Schur complement ( A λi B(D λi ) 1 C ) u = Idea: Given T (λ) = A λi F (λ), find a rational approximation F (λ) B(D λi ) 1 C. Berkeley 29 / 37
32 Expansion via rational approximation u() u(a) u(b) =. Berkeley 3 / 37
33 Analyzing the expanded system ˆT (z) is a Schur complement in K zm So Λ( ˆT ) is easy to compute. Or: think T (z) is a Schur complement in K zm + E(z) Compare ˆT (z) to T (z) or compare K zm + E(z) to K zm Berkeley 31 / 37
34 Analyzing the expanded system Q: Can we find all eigs in a region not missing anything? Concrete plan (ɛ = 1 8 ) T = shooting system ˆT = rational approximation Find region D with boundary Γ s.t. D Ω ɛ (i.e. T ˆT < ɛ on D) Γ does not intersect Λ ɛ (T ) = Same eigenvalue counts for T, ˆT = Eigs of ˆT in components of Λ ɛ (T ) Converse holds if D Ω ɛ/2 Can refine eigs of ˆT in D via Newton. Berkeley 32 / 37
35 6 8 Resonance approximation Berkeley 33 / 37
36 8 1 Resonance approximation Berkeley 33 / 37
37 Example III: Bounding dynamics DDE model (of a laser with feedback) u = Au(t) + Bu(t 1). Solution with u(+) = u and u(t) = φ(t) on [ 1, ): u(t) = Ψ(t)u + 1 Ψ(t s)bφ(s t) ds where Ψ(t) is a fundamental matrix for shock solutions. Berkeley 34 / 37
38 Representing Ψ(t) Associated nonlinear eigenvalue problem involves T (z) = zi A Be z Assume all eigs in left half plane; via inverse Laplace transform, Ψ(t) = 1 T (z) 1 e zt dz. 2πi for contour Γ right of the spectrum. Γ Berkeley 35 / 37
39 Bounding pseudospectra = bounding dynamics Bound fundamental solution Ψ via Ψ(t) = 1 T (z) 1 e zt dz 2πi 1 2π Γ Γ T (z) 1 e zt dz. Berkeley 36 / 37
40 For more Localization theorems for nonlinear eigenvalues. David Bindel and Amanda Hood, SIAM Review 57(4), Dec 215 Pseudospectral bounds on transient growth for higher order and constant delay differential equations. Amanda Hood and David Bindel, submitted to SIMAX arxiv: , Nov 216 Berkeley 37 / 37
41 Trailers! Application Analysis Computation Berkeley 37 / 37
42 Spectral Topic Modeling A B T C B Berkeley 37 / 37
43 Music of the Microspheres Berkeley 37 / 37
44 Fast Fingerprints for Power Systems Berkeley 37 / 37
45 Response Surfaces for Global Optimization Berkeley 37 / 37
46 x x Graph Densities of States Berkeley 37 / 37
47 Fin Berkeley 37 / 37
Numerical Analysis for Nonlinear Eigenvalue Problems
Numerical Analysis for Nonlinear Eigenvalue Problems D. Bindel Department of Computer Science Cornell University 14 Sep 29 Outline Nonlinear eigenvalue problems Resonances via nonlinear eigenproblems Sensitivity
More informationBounds and Error Estimates for Nonlinear Eigenvalue Problems
Bounds and Error Estimates for Nonlinear Eigenvalue Problems D. Bindel Courant Institute for Mathematical Sciences New York Univerity 8 Oct 28 Outline Resonances via nonlinear eigenproblems Sensitivity
More informationLOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS
LOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS DAVID BINDEL AND AMANDA HOOD Abstract. Let T : Ω C n n be a matrix-valued function that is analytic on some simplyconnected domain Ω C. A point λ
More informationLOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS
LOCALIZATION THEOREMS FOR NONLINEAR EIGENVALUE PROBLEMS DAVID BINDEL AND AMANDA HOOD Abstract. Let T : Ω C n n be a matrix-valued function that is analytic on some simplyconnected domain Ω C. A point λ
More information1 Sylvester equations
1 Sylvester equations Notes for 2016-11-02 The Sylvester equation (or the special case of the Lyapunov equation) is a matrix equation of the form AX + XB = C where A R m m, B R n n, B R m n, are known,
More informationEIGENVALUE PROBLEMS. Background on eigenvalues/ eigenvectors / decompositions. Perturbation analysis, condition numbers..
EIGENVALUE PROBLEMS Background on eigenvalues/ eigenvectors / decompositions Perturbation analysis, condition numbers.. Power method The QR algorithm Practical QR algorithms: use of Hessenberg form and
More informationResearch Matters. February 25, The Nonlinear Eigenvalue. Director of Research School of Mathematics
Research Matters February 25, 2009 The Nonlinear Eigenvalue Nick Problem: HighamPart I Director of Research School of Mathematics Françoise Tisseur School of Mathematics The University of Manchester Woudschoten
More informationEigenvalues and Eigenvectors
/88 Chia-Ping Chen Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Eigenvalue Problem /88 Eigenvalue Equation By definition, the eigenvalue equation for matrix
More informationA Tale of Two Eigenvalue Problems
A Tale of Two Eigenvalue Problems Music of the Microspheres + Hearing the Shape of a Graph David Bindel Department of Computer Science Cornell University CAM Visit Talk, 28 Feb 214 Brown bag 1 / 22 The
More informationStructure-preserving model reduction for MEMS
Structure-preserving model reduction for MEMS David Bindel Department of Computer Science Cornell University SIAM CSE Meeting, 1 Mar 2011 Collaborators Sunil Bhave Emmanuel Quévy Zhaojun Bai Tsuyoshi Koyama
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More informationEigenvalue Problems in Surface Acoustic Wave Filter Simulations
Eigenvalue Problems in Surface Acoustic Wave Filter Simulations S. Zaglmayr, J. Schöberl, U. Langer FWF Start Project Y-192 3D hp-finite Elements: Fast Solvers and Adaptivity JKU + RICAM, Linz in collaboration
More informationLecture 4 Eigenvalue problems
Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.2: Fundamentals 2 / 31 Eigenvalues and Eigenvectors Eigenvalues and eigenvectors of
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 informationDiagonalization of the Coupled-Mode System.
Diagonalization of the Coupled-Mode System. Marina Chugunova joint work with Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaborators: Mason A. Porter, California Institute
More informationCollocation approximation of the monodromy operator of periodic, linear DDEs
p. Collocation approximation of the monodromy operator of periodic, linear DDEs Ed Bueler 1, Victoria Averina 2, and Eric Butcher 3 13 July, 2004 SIAM Annual Meeting 2004, Portland 1=Dept. of Math. Sci.,
More informationOn rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro
On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational
More informationMATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2.
MATH 304 Linear Algebra Lecture 23: Diagonalization. Review for Test 2. Diagonalization Let L be a linear operator on a finite-dimensional vector space V. Then the following conditions are equivalent:
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 informationQualifying Exam Complex Analysis (Math 530) January 2019
Qualifying Exam Complex Analysis (Math 53) January 219 1. Let D be a domain. A function f : D C is antiholomorphic if for every z D the limit f(z + h) f(z) lim h h exists. Write f(z) = f(x + iy) = u(x,
More informationCOMPLEX ANALYSIS Spring 2014
COMPLEX ANALYSIS Spring 24 Homework 4 Solutions Exercise Do and hand in exercise, Chapter 3, p. 4. Solution. The exercise states: Show that if a
More informationNotes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.
Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.. Tuesday, June 2, 202 Combinatorics is the study of finite structures that
More informationEvans function review
Evans function review Part I: History, construction, properties and applications Strathclyde, April 18th 2005 Simon J.A. Malham http://www.ma.hw.ac.uk/ simonm/talks/ Acknowledgments T.J. Bridges, C.K.R.T.
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationEigenvalue and Eigenvector Problems
Eigenvalue and Eigenvector Problems An attempt to introduce eigenproblems Radu Trîmbiţaş Babeş-Bolyai University April 8, 2009 Radu Trîmbiţaş ( Babeş-Bolyai University) Eigenvalue and Eigenvector Problems
More informationPart IB. Further Analysis. Year
Year 2004 2003 2002 2001 10 2004 2/I/4E Let τ be the topology on N consisting of the empty set and all sets X N such that N \ X is finite. Let σ be the usual topology on R, and let ρ be the topology on
More informationConstrained Nonlinear Optimization Algorithms
Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu Institute for Mathematics and its Applications University of Minnesota August 4, 2016
More informationLecture notes: Applied linear algebra Part 2. Version 1
Lecture notes: Applied linear algebra Part 2. Version 1 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 First, some exercises: xercise 0.1 (2 Points) Another least
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationACM 104. Homework Set 5 Solutions. February 21, 2001
ACM 04 Homework Set 5 Solutions February, 00 Franklin Chapter 4, Problem 4, page 0 Let A be an n n non-hermitian matrix Suppose that A has distinct eigenvalues λ,, λ n Show that A has the eigenvalues λ,,
More informationNonlinear Eigenvalue Problems: Interpolatory Algorithms and Transient Dynamics
Nonlinear Eigenvalue Problems: Interpolatory Algorithms and Transient Dynamics Mark Embree Virginia Tech with Serkan Gugercin Jonathan Baker, Michael Brennan, Alexander Grimm Virginia Tech SIAM Conference
More informationScalable kernel methods and their use in black-box optimization
with derivatives Scalable kernel methods and their use in black-box optimization David Eriksson Center for Applied Mathematics Cornell University dme65@cornell.edu November 9, 2018 1 2 3 4 1/37 with derivatives
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationMath 185 Fall 2015, Sample Final Exam Solutions
Math 185 Fall 2015, Sample Final Exam Solutions Nikhil Srivastava December 12, 2015 1. True or false: (a) If f is analytic in the annulus A = {z : 1 < z < 2} then there exist functions g and h such that
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 informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationMTH 3102 Complex Variables Final Exam May 1, :30pm-5:30pm, Skurla Hall, Room 106
Name (Last name, First name): MTH 02 omplex Variables Final Exam May, 207 :0pm-5:0pm, Skurla Hall, Room 06 Exam Instructions: You have hour & 50 minutes to complete the exam. There are a total of problems.
More informationBarycentric interpolation via the AAA algorithm
Barycentric interpolation via the AAA algorithm Steven Elsworth Stefan Güttel November 217 Contents 1 Introduction 1 2 A simple scalar example 1 3 Solving a nonlinear eigenproblem 2 4 References 6 1 Introduction
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationSOLUTIONS TO SELECTED PROBLEMS FROM ASSIGNMENTS 3, 4
SOLUTIONS TO SELECTED POBLEMS FOM ASSIGNMENTS 3, 4 Problem 5 from Assignment 3 Statement. Let be an n-dimensional bounded domain with smooth boundary. Show that the eigenvalues of the Laplacian on with
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More informationChapter 12 Solving secular equations
Chapter 12 Solving secular equations Gérard MEURANT January-February, 2012 1 Examples of Secular Equations 2 Secular equation solvers 3 Numerical experiments Examples of secular equations What are secular
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationSample Solutions of Assignment 9 for MAT3270B
Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically
More informationEigenvalue problems and optimization
Notes for 2016-04-27 Seeking structure For the past three weeks, we have discussed rather general-purpose optimization methods for nonlinear equation solving and optimization. In practice, of course, we
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH College of Informatics and Electronics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MS425 SEMESTER: Autumn 25/6 MODULE TITLE: Applied Analysis DURATION OF EXAMINATION:
More informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationAnnouncements Wednesday, November 7
Announcements Wednesday, November 7 The third midterm is on Friday, November 6 That is one week from this Friday The exam covers 45, 5, 52 53, 6, 62, 64, 65 (through today s material) WeBWorK 6, 62 are
More informationx x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)
Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)
More informationEvaluation of Chebyshev pseudospectral methods for third order differential equations
Numerical Algorithms 16 (1997) 255 281 255 Evaluation of Chebyshev pseudospectral methods for third order differential equations Rosemary Renaut a and Yi Su b a Department of Mathematics, Arizona State
More informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationConsidering our result for the sum and product of analytic functions, this means that for (a 0, a 1,..., a N ) C N+1, the polynomial.
Lecture 3 Usual complex functions MATH-GA 245.00 Complex Variables Polynomials. Construction f : z z is analytic on all of C since its real and imaginary parts satisfy the Cauchy-Riemann relations and
More informationEigenvalues, Resonance Poles, and Damping in MEMS
Eigenvalues,, and Damping in D. Bindel Computer Science Division Department of EECS University of California, Berkeley Matrix computations seminar, 19 Apr 26 Outline 1 2 3 4 Outline 1 2 3 4 Resonating
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More informationLecture 5 Rational functions and partial fraction expansion
EE 102 spring 2001-2002 Handout #10 Lecture 5 Rational functions and partial fraction expansion (review of) polynomials rational functions pole-zero plots partial fraction expansion repeated poles nonproper
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 informationNonlinear Programming Algorithms Handout
Nonlinear Programming Algorithms Handout Michael C. Ferris Computer Sciences Department University of Wisconsin Madison, Wisconsin 5376 September 9 1 Eigenvalues The eigenvalues of a matrix A C n n are
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 informationStability of flow past a confined cylinder
Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder
More informationSolutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011
Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 8, Solutions to Problems 5 are omitted since they involve topics no longer covered on the Comprehensive
More informationMATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5
MATH 722, COMPLEX ANALYSIS, SPRING 2009 PART 5.. The Arzela-Ascoli Theorem.. The Riemann mapping theorem Let X be a metric space, and let F be a family of continuous complex-valued functions on X. We have
More informationParameter-Dependent Eigencomputations and MEMS Applications
Parameter-Dependent Eigencomputations and MEMS Applications David Bindel UC Berkeley, CS Division Parameter-Dependent Eigencomputations and MEMS Applications p.1/36 Outline Some applications Mathematical
More informationUnderstanding graphs through spectral densities
Understanding graphs through spectral densities David Bindel Department of Computer Science Cornell University SCAN seminar, 29 Feb 216 SCAN 1 / 1 Can One Hear the Shape of a Drum? 2 u = λu on Ω, u = on
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
More informationMath Homework 2
Math 73 Homework Due: September 8, 6 Suppose that f is holomorphic in a region Ω, ie an open connected set Prove that in any of the following cases (a) R(f) is constant; (b) I(f) is constant; (c) f is
More informationCyclic Reduction History and Applications
ETH Zurich and Stanford University Abstract. We discuss the method of Cyclic Reduction for solving special systems of linear equations that arise when discretizing partial differential equations. In connection
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationNumerical and semi-analytical structure-preserving model reduction for MEMS
Numerical and semi-analytical structure-preserving model reduction for MEMS David Bindel ENUMATH 07, 10 Sep 2007 Collaborators Tsuyoshi Koyama Sanjay Govindjee Sunil Bhave Emmanuel Quévy Zhaojun Bai Resonant
More informationMATH 115A: SAMPLE FINAL SOLUTIONS
MATH A: SAMPLE FINAL SOLUTIONS JOE HUGHES. Let V be the set of all functions f : R R such that f( x) = f(x) for all x R. Show that V is a vector space over R under the usual addition and scalar multiplication
More informationKrylov methods for the computation of matrix functions
Krylov methods for the computation of matrix functions Jitse Niesen (University of Leeds) in collaboration with Will Wright (Melbourne University) Heriot-Watt University, March 2010 Outline Definition
More informationSecond Midterm Exam Name: Practice Problems March 10, 2015
Math 160 1. Treibergs Second Midterm Exam Name: Practice Problems March 10, 015 1. Determine the singular points of the function and state why the function is analytic everywhere else: z 1 fz) = z + 1)z
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationFinal Exam December 20, 2011
Final Exam December 20, 2011 Math 420 - Ordinary Differential Equations No credit will be given for answers without mathematical or logical justification. Simplify answers as much as possible. Leave solutions
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationQuadratic Matrix Polynomials
Research Triangularization Matters of Quadratic Matrix Polynomials February 25, 2009 Nick Françoise Higham Tisseur Director School of of Research Mathematics The University of Manchester School of Mathematics
More informationEvaluation of integrals
Evaluation of certain contour integrals: Type I Type I: Integrals of the form 2π F (cos θ, sin θ) dθ If we take z = e iθ, then cos θ = 1 (z + 1 ), sin θ = 1 (z 1 dz ) and dθ = 2 z 2i z iz. Substituting
More information6 Linear Equation. 6.1 Equation with constant coefficients
6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ
More informationPart IB. Complex Analysis. Year
Part IB Complex Analysis Year 2018 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2018 Paper 1, Section I 2A Complex Analysis or Complex Methods 7 (a) Show that w = log(z) is a conformal
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 informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationClassical iterative methods for linear systems
Classical iterative methods for linear systems Ed Bueler MATH 615 Numerical Analysis of Differential Equations 27 February 1 March, 2017 Ed Bueler (MATH 615 NADEs) Classical iterative methods for linear
More informationThen x 1,..., x n is a basis as desired. Indeed, it suffices to verify that it spans V, since n = dim(v ). We may write any v V as r
Practice final solutions. I did not include definitions which you can find in Axler or in the course notes. These solutions are on the terse side, but would be acceptable in the final. However, if you
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationMath 126 Winter 2012: Rough lecture notes
Math 126 Winter 2012: Rough lecture notes Alex Barnett, Dartmouth College March 1, 2012 1 Introduction Numerical mathematics is at the intersection of analysis (devising and proving theorems), computation
More informationPoisson Image Editing
Poisson Image Editing Sébastien Boisgérault, Mines ParisTech, under CC BY-NC-SA 4.0 January 21, 2016 Contents Introduction 1 Modelling the Problem 2 Harmonic Functions 5 The irichlet Problem 8 Harmonic/Fourier
More informationEigenvalues and Eigenvectors
Sec. 6.1 Eigenvalues and Eigenvectors Linear transformations L : V V that go from a vector space to itself are often called linear operators. Many linear operators can be understood geometrically by identifying
More informationDiagonalization of Matrices
LECTURE 4 Diagonalization of Matrices Recall that a diagonal matrix is a square n n matrix with non-zero entries only along the diagonal from the upper left to the lower right (the main diagonal) Diagonal
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
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 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 information