Nonlinear Eigenvalue Problems: An Introduction
|
|
- Mervyn Glenn
- 5 years ago
- Views:
Transcription
1 Nonlinear Eigenvalue Problems: An Introduction Cedric Effenberger Seminar for Applied Mathematics ETH Zurich Pro*Doc Workshop Disentis, August 18 21, 2010 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 1 / 37
2 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 2 / 37 Overview 1 Definition 2 Applications 3 Solution 4 Solution of polynomial / rational EVPs
3 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 3 / 37 Linear Eigenvalue Problems Standard eigenvalue problem Given A C n n, we seek eigenvalues λ C, and eigenvectors x C n \ {0} such that Ax = λx. Generalized eigenvalue problem Given A, B C n n, we seek eigenvalues λ C, and eigenvectors x C n \ {0} such that Ax = λbx.
4 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 4 / 37 Generalization to nonlinear EVPs Both problems can be reformulated as T (λ)x = 0, where T (λ) = { A λi, A λb, for the standard EVP, for the generalized EVP. T (λ) = λ d A d + + λa 1 + A 0 polynomial EVP T (λ) = f 1 (λ)a f r (λ)a r nonlinear EVP A j C n n constant matrices, f j : D C analytic, D C open, connected
5 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 5 / 37 Free vibrations of mechanical systems with damping Equations of motion Mü + C u + Ku = 0 M: mass matrix C: damping matrix K: stiffness matrix System of second-order, linear, homogeneous ODEs Eigenfrequencies ω of the mechanical system are eigenvalues of Quadratic eigenvalue problem (ω 2 M + ωc + K)u = 0.
6 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 6 / 37 Free vibrations of string with elastically attached mass vibrating string u = λu in Ω = (0, 1) left end clamped u(0) = 0 mass elastically attached to right end u (1) + αλ λ α u(1) = 0 λ-dependent boundary condition! Discretization with n piecewise linear FE leads to [ K + αλ λ α e nen T ] λm x = 0.
7 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 7 / 37 Delay eigenvalue problems Linear delay differential equation with r discrete delays u(t) = A 0 u(t) + A 1 u(t τ 1 ) + + A r u(t τ r ) models influences which take effect only after some time famous example: Hot shower problem stability analysis involves nonlinear EVPs Delay eigenvalue problem ( λi + A 0 + e τ 1λ A e τr λ A r )x = 0
8 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 8 / 37 Electronic bandstructure computations photonic crystal: lattice of mixed dielectric media control electromagnetic waves by designing the crystal such that it inhibits their propagation complete photonic band gap: frequency range with no propagation of electromagnetic waves of any polarization travelling in any direction
9 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 9 / 37 Electronic bandstructure computations II 2D crystal: periodic in x- and y-direction; homogeneous in z-direction consider only electromagnetic waves propagating in the xy-plane
10 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 10 / 37 Electronic bandstructure computations III Time-harmonic modes of electromagnetic wave (E, H) can be decomposed into transverse electric (TE) polarized modes (E x, E y, 0, 0, 0, H z ), transverse magnetic (TM) polarized modes (0, 0, E z, H x, H y, 0). For TM polarized modes, the macroscopic Maxwell Equations reduce to a scalar equation for E z, r: spatial variable ω: frequency ε: relative permittivity E z = ω 2 ε(r, ω)e z. ω-dependent material parameter!
11 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 11 / 37 Electronic bandstructure computations IV Bloch ansatz: E z = e ik r u(r) ( + ik) ( + ik)u(r) = ω 2 ε(r, ω)u(r) k: wave vector u: periodic function on lattice assumption: lattice consists of 2 materials, one of which is air { ε 1 = 1, r Ω 1 Ω = Ω 1 Ω 2, ε(r, ω) = ε 2 (ω), r Ω 2 discretization using discontinuous Galerkin with p-enhancement (Engström & Wang, 2010) [ G ω 2 M 1 ω 2 ε 2 (ω)m 2 ] u = 0
12 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 12 / 37 Electronic bandstructure computations V nonlinearity caused by ω-dependency of ε 2 popular model for permittivity: Lorentz model K ξ k ε 2 (ω) = α + η k ω 2 iγ k ω k=1 rational EVP ω 2π
13 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 13 / 37 Spectrum of nonlinear EVPs Definition Let T : D C n n be holomorphic. We call ρ(t ) := {λ D : T (λ) is invertible} the resolvent set of T, σ(t ) := D \ ρ(t ) the spectrum of T. Theorem Either the resolvent set of T is empty or the spectrum of T consists of isolated eigenvalues. The number of eigenvalues in σ(t ) may be infinite!
14 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 14 / 37 Newton-based methods Two possibilities to apply Newton s method: 1 det T (λ) = 0 (scalar equation in C) 2 T (λ)x = 0 (vector equation in C n ) We focus on the second case. x, λ n + 1 unknowns T (λ)x = 0 n constraints We have to add one additional constraint. v H x = 1
15 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 15 / 37 Nonlinear inverse iteration F (x, λ) := [ ] T (λ)x v H x 1 Application of Newton s method ] [ ] [ˆx 0 =! x T (λ)ˆx + (ˆλ λ)t F (x, λ) + DF (x, λ) = (λ)x ˆλ λ v H ˆx 1 yields ˆx = (ˆλ λ)t (λ) 1 T (λ)x 1 ˆλ = λ v H T (λ) 1 T (λ)x.
16 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 16 / 37 Nonlinear inverse iteration II Algorithm Input: approximate eigenpair (x 0, λ 0 ) with v H x 0 = 1 for j = 0, 1,... until convergence do solve T (λ j ) x j+1 = T (λ j )x j for x j+1 update λ j+1 := λ j v H x j v H x j+1 normalize x j+1 := 1 v H u j+1 x j+1 end for local quadratic convergence to simple eigenvalues main computational work lies in the solution of the linear system
17 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 17 / 37 Reducing the computational work Replacing T (λ j ) 1 by T (σ) 1 for a fixed σ leads to misconvergence! But... ˆx = (ˆλ λ)t (λ) 1 T (λ)x = x T (λ) 1 T (λ)x (ˆλ λ)t (λ) 1 T (λ)x = x T (λ) 1[ T (λ) + (ˆλ λ)t (λ) ] x = x T (λ) 1 T (ˆλ)x + O ( ˆλ λ 2) Now T (σ) 1 can be used in place of T (λ) 1 safely. Even inexact solution of the linear system is possible.
18 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 18 / 37 Residual inverse iteration Algorithm Input: approximate eigenpair (x 0, σ) for j = 0, 1,... until convergence do solve v H T (σ) 1 T (λ j+1 )x j = 0 for λ j+1 solve T (σ) x = T (λ j+1 )x j for x update x j+1 := x j + x normalize x j+1 := 1 v H x j+1 x j+1 end for local linear convergence to simple eigenvalues convergence speed dependent on distance from σ to eigenvalue
19 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 19 / 37 Convergence string with elastically attached mass smallest magnitude eigenvalue using residual inverse iteration initial guess: σ = 0, x discrete version of function u(z) = z rel. error iterations it. λ
20 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 20 / 37 Difficulties no global convergence can only handle one eigenvalue at a time cannot handle multiple eigenvalues consecutive runs may converge to the same eigenvalue Linear eigensolvers (ARPACK, Jacobi-Davidson) exploit linear independence of eigenvectors to prevent reconvergence. Example (loss of linear independence for NLEVPs) The eigenvalues 1 and 2 of ([ ] [ ] [ ]) λ + λ x = share the same eigenvector [ 1 1 ].
21 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 21 / 37 Several eigenvalues in the linear case Ax = λx Let (x 1, λ 1 ), (x 2, λ 2 ) be two eigenpairs: Ax 1 = λ 1 x 1, Ax 2 = λ 2 x 2. The above equations can be merged: A [ [ ] ] [ ] λ1 x 1 x 2 = x1 x 2 }{{} λ 2 =:X }{{} =:Λ Hence, we have AX = X Λ. X spans an invariant subspace of A.
22 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 22 / 37 Generalization to the nonlinear case f 1 (λ)a 1 x + + f r (λ)a r x = 0 Assume r = 2 and let (x 1, λ 1 ), (x 2, λ 2 ) be two eigenpairs. A 1 x 1 f 1 (λ 1 ) + A 2 x 1 f 2 (λ 1 ) = 0, A 1 x 2 f 1 (λ 2 ) + A 2 x 2 f 2 (λ 2 ) = 0. These equations can again be merged [ ] [ ] f A 1 x1 x 1 (λ 1 ) [ ] [ f 2 + A f 1 (λ 2 ) 2 x1 x 2 (λ 1 ) 2 ] = 0 f 2 (λ 2 ) With X = [ x 1 x 2 ] and Λ = diag (λ1, λ 2 ) as before A 1 Xf 1 (Λ) + A 2 Xf 2 (Λ) = 0.
23 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 23 / 37 Invariant pairs Definition (X, Λ) C n m C m m is called an invariant pair if A 1 Xf 1 (Λ) + + A r Xf r (Λ) = 0. need to exclude degenerate situations, such as X = 0 linear case: require X to have full column rank not suitable for the nonlinear case: [ ] [ ] X =, Λ = is a perfectly reasonable invariant pair of ([ ] [ ] [ ]) λ + λ x =
24 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 24 / 37 Minimality Definition X = An invariant pair (X, Λ) is called minimal if X X Λ. X Λ l 1 has full column rank for some integer l. The smallest such l is called the minimality index of (X, Λ). [ ] [ 1 1 1, Λ = ] [ ] 1 1 X = 1 1 X Λ
25 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 25 / 37 Minimal invariant pairs (X, Λ) is a minimal invariant pair with minimality index 2. Theorem Let (X, Λ) be a minimal invariant pair of a nonlinear EVP T (λ)x = 0. Then every eigenvalue of Λ is an eigenvalue of T.
26 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 26 / 37 Computation of invariant pairs Apply Newton s method to T(X, Λ) := A 1 Xf 1 (Λ) + + A r Xf r (Λ) = 0. X, Λ n m + m 2 unknowns T(X, Λ) = 0 n m constraints We have to impose m 2 additional constraints. N(X, Λ) := V H Theorem X X Λ. X Λ l 1 I = 0, V suitably chosen Let (X, Λ) be a minimal invariant pair of T. The Fréchet derivative DF of F := [ T N ] at (X, Λ) is invertible iff the multiplicities of Λ s eigenvalues match those of T.
27 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 27 / 37 Ongoing work: Continuation of invariant pairs Joint project with Daniel Kressner (SAM, ETHZ) and Wolf-Jürgen Beyn (University of Bielefeld) Parameter-dependent nonlinear EVP T (λ, s)x = 0 Goal: track several eigenvalues as the parameter s varies Idea: use invariant pair ( X (s), Λ(s) ) at parameter value s as an initial guess for invariant pair at s + s apply Newton to obtain invariant pair ( X (s + s), Λ(s + s) ) at parameter value s + s repeat
28 Part II: Solution of polynomial / rational EVPs Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 28 / 37
29 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 29 / 37 Linearization Polynomial EVPs can be solved by transforming them into linear ones. Example: (λ 2 M + λc + K)x = 0 introduce auxiliary variable y = λx rewrite as λmy + Cy + Kx = 0 This can be written as a generalized linear EVP with the same eigenvalues: [ ] [ ] [ ] [ ] C K y M 0 y = λ. I 0 x 0 I x linearizations not unique could also rewrite as λmy + λcx + Kx = 0 [ ] [ ] [ ] [ ] K 0 x C M x = λ 0 I y I 0 y
30 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 30 / 37 Linearization II There are entire vector spaces of linearizations (Mackey et al., 2006)! choice of linearization depends on underlying application structure preservation of interest Example The eigenvalues of an even polynomial eigenvalue problem (λ 2 M + λc + K)x = 0, M = M T, C = C T, K = K T occur in pairs (λ, λ). If M is invertible, the symmetric / skew-symmetric linearization [ ] [ ] [ ] [ ] K 0 x C M x = λ 0 M y M 0 y preserves this property.
31 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 31 / 37 Linearization III linearization of polynomial EVPs of degree d > 2 analogously by introducing the auxiliary vector x λx. λ d 1 x Conclusion 1: polynomial EVP of size n linearized EVP of size dn polynomial EVP has exactly dn eigenvalues (counting multiplicities) Conclusion 2: The extended vectors above will always be linearly independent, even though the eigenvectors x themselves may not.
32 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 32 / 37 Solution of rational EVPs most straightforward way: multiplication by the common denominator of all rational terms polynomial EVP degree of resulting polynomial EVP may be very large sometimes smarter ways to solve rational EVPs: Ongoing joint project with Daniel Kressner and Christian Engström (both SAM, ETHZ) [ G ω 2 M 1 ω 2( ξ ] α + )M η ω 2 2 u = 0 iγω M 1 + M 2 is a splitting of the total mass matrix of the problem.
33 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 33 / 37 Rational EVPs in electronic bandstructure computations by suitable ordering of nodes: [ 0 M 2 = R ˇM2] n n with a positive definite ˇM 2 R m m, m n by Cholesky decomposition: M 2 = F T F, F R m n Idea (Bai & Su, 2008): write rational term(s) in transfer function form ω 2 ξ η ω 2 iγω = ξ + bt (A ωe) 1 b
34 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 34 / 37 Conversion to transfer function form by partial fraction expansion p(ω) q(ω) = ζ + σ σ k ω ρ 1 ω ρ k ( ω = ζ + ρ ) 1 1 ( ω + + ρ ) k 1 σ 1 σ 1 σ k σ k = ζ + b T (A ωe) 1 b, where 1 b =., A = 1 ρ 1 σ 1... ρ k σ k, E = 1 σ σ k
35 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 35 / 37 Linearization of rational EVPs EVP with rational term in transfer function form: Ĝu b T (A ωe) 1 bf T F u = ω 2 ˆMu, where Ĝ = G + ξm 2, ˆM = M 1 + αm 2 rearrange: b T (A ωe)bf T F = b T (A ωe) 1 b F T F = B T (A ωe) 1 B where B = b F, A = A I, E = E I
36 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 36 / 37 Linearization of rational EVPs II Ĝu B T (A ωe) 1 Bu = ω 2 ˆMu Introducing the auxiliary variables v := (A ωe) 1 Bu, we obtain the linearized EVP Ĝ B T u ˆM u B A v = ω E v I w I w w := ωu, problem size standard lin. rational lin. p #DoF size comp. time size comp. time s s s s s s s
37 Cedric Effenberger (SAM, ETHZ) NLEVPs: An Introduction Pro*Doc Disentis 37 / 37 Summary introduced nonlinear EVPs showed several sample applications discussed Newton-based technique for determining one or several eigenpairs demonstrated linearization techniques for polynomial and rational EVPs outlined two ongoing projects
Definite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationPolynomial Jacobi Davidson Method for Large/Sparse Eigenvalue Problems
Polynomial Jacobi Davidson Method for Large/Sparse Eigenvalue Problems Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan April 28, 2011 T.M. Huang (Taiwan Normal Univ.)
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 informationResearch Matters. February 25, The Nonlinear Eigenvalue Problem. Nick Higham. Part III. Director of Research School of Mathematics
Research Matters February 25, 2009 The Nonlinear Eigenvalue Problem Nick Higham Part III Director of Research School of Mathematics Françoise Tisseur School of Mathematics The University of Manchester
More informationSOLVING RATIONAL EIGENVALUE PROBLEMS VIA LINEARIZATION
SOLVNG RATONAL EGENVALUE PROBLEMS VA LNEARZATON YANGFENG SU AND ZHAOJUN BA Abstract Rational eigenvalue problem is an emerging class of nonlinear eigenvalue problems arising from a variety of physical
More informationEigenvalues and Eigenvectors
Contents Eigenvalues and Eigenvectors. Basic Concepts. Applications of Eigenvalues and Eigenvectors 8.3 Repeated Eigenvalues and Symmetric Matrices 3.4 Numerical Determination of Eigenvalues and Eigenvectors
More informationRational Krylov methods for linear and nonlinear eigenvalue problems
Rational Krylov methods for linear and nonlinear eigenvalue problems Mele Giampaolo mele@mail.dm.unipi.it University of Pisa 7 March 2014 Outline Arnoldi (and its variants) for linear eigenproblems Rational
More informationRobust successive computation of eigenpairs for nonlinear eigenvalue problems
Robust successive computation of eigenpairs for nonlinear eigenvalue problems C. Effenberger July 21, 2012 Newton-based methods are well-established techniques for solving nonlinear eigenvalue problems.
More informationYORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #1. July 11, 2013 Solutions
YORK UNIVERSITY Faculty of Science Department of Mathematics and Statistics MATH 222 3. M Test # July, 23 Solutions. For each statement indicate whether it is always TRUE or sometimes FALSE. Note: For
More informationRecent Advances in the Numerical Solution of Quadratic Eigenvalue Problems
Recent Advances in the Numerical Solution of Quadratic Eigenvalue Problems Françoise Tisseur School of Mathematics The University of Manchester ftisseur@ma.man.ac.uk http://www.ma.man.ac.uk/~ftisseur/
More informationScaling, Sensitivity and Stability in Numerical Solution of the Quadratic Eigenvalue Problem
Scaling, Sensitivity and Stability in Numerical Solution of the Quadratic Eigenvalue Problem Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/
More informationNumerical Methods I Eigenvalue Problems
Numerical Methods I Eigenvalue Problems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 2nd, 2014 A. Donev (Courant Institute) Lecture
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More information1. Introduction. We consider nonlinear eigenvalue problems of the form
ROBUST SUCCESSIVE COMPUTATION OF EIGENPAIRS FOR NONLINEAR EIGENVALUE PROBLEMS C. EFFENBERGER Abstract. Newton-based methods are well-established techniques for solving nonlinear eigenvalue problems. If
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.4.1
More informationAlgorithms for Solving the Polynomial Eigenvalue Problem
Algorithms for Solving the Polynomial Eigenvalue Problem Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with D. Steven Mackey
More informationSolution of eigenvalue problems. Subspace iteration, The symmetric Lanczos algorithm. Harmonic Ritz values, Jacobi-Davidson s method
Solution of eigenvalue problems Introduction motivation Projection methods for eigenvalue problems Subspace iteration, The symmetric Lanczos algorithm Nonsymmetric Lanczos procedure; Implicit restarts
More informationEigenvalues and Eigenvectors
Eigenvalues and Eigenvectors week -2 Fall 26 Eigenvalues and eigenvectors The most simple linear transformation from R n to R n may be the transformation of the form: T (x,,, x n ) (λ x, λ 2,, λ n x n
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 informationLecture 3: Inexact inverse iteration with preconditioning
Lecture 3: Department of Mathematical Sciences CLAPDE, Durham, July 2008 Joint work with M. Freitag (Bath), and M. Robbé & M. Sadkane (Brest) 1 Introduction 2 Preconditioned GMRES for Inverse Power Method
More informationOn Dominant Poles and Model Reduction of Second Order Time-Delay Systems
On Dominant Poles and Model Reduction of Second Order Time-Delay Systems Maryam Saadvandi Joint work with: Prof. Karl Meerbergen and Dr. Elias Jarlebring Department of Computer Science, KULeuven ModRed
More informationSeminar 6: COUPLED HARMONIC OSCILLATORS
Seminar 6: COUPLED HARMONIC OSCILLATORS 1. Lagrangian Equations of Motion Let consider a system consisting of two harmonic oscillators that are coupled together. As a model, we will use two particles attached
More informationA computer-assisted Band-Gap Proof for 3D Photonic Crystals
A computer-assisted Band-Gap Proof for 3D Photonic Crystals Henning Behnke Vu Hoang 2 Michael Plum 2 ) TU Clausthal, Institute for Mathematics, 38678 Clausthal, Germany 2) Univ. of Karlsruhe (TH), Faculty
More informationMultiparameter eigenvalue problem as a structured eigenproblem
Multiparameter eigenvalue problem as a structured eigenproblem Bor Plestenjak Department of Mathematics University of Ljubljana This is joint work with M Hochstenbach Będlewo, 2932007 1/28 Overview Introduction
More informationACCURATE SOLUTIONS OF POLYNOMIAL EIGENVALUE PROBLEMS
ACCURATE SOLUTIONS OF POLYNOMIAL EIGENVALUE PROBLEMS YILING YOU, JOSE ISRAEL RODRIGUEZ, AND LEK-HENG LIM Abstract. Quadratic eigenvalue problems (QEP) and more generally polynomial eigenvalue problems
More informationChapter 5 Eigenvalues and Eigenvectors
Chapter 5 Eigenvalues and Eigenvectors Outline 5.1 Eigenvalues and Eigenvectors 5.2 Diagonalization 5.3 Complex Vector Spaces 2 5.1 Eigenvalues and Eigenvectors Eigenvalue and Eigenvector If A is a n n
More informationA Jacobi Davidson Method for Nonlinear Eigenproblems
A Jacobi Davidson Method for Nonlinear Eigenproblems Heinrich Voss Section of Mathematics, Hamburg University of Technology, D 21071 Hamburg voss @ tu-harburg.de http://www.tu-harburg.de/mat/hp/voss Abstract.
More informationAdvanced Vibrations. Elements of Analytical Dynamics. By: H. Ahmadian Lecture One
Advanced Vibrations Lecture One Elements of Analytical Dynamics By: H. Ahmadian ahmadian@iust.ac.ir Elements of Analytical Dynamics Newton's laws were formulated for a single particle Can be extended to
More informationFAME-matlab Package: Fast Algorithm for Maxwell Equations Tsung-Ming Huang
FAME-matlab Package: Fast Algorithm for Maxwell Equations Tsung-Ming Huang Modelling, Simulation and Analysis of Nonlinear Optics, NUK, September, 4-8, 2015 1 2 FAME group Wen-Wei Lin Department of Applied
More informationIterative projection methods for sparse nonlinear eigenvalue problems
Iterative projection methods for sparse nonlinear eigenvalue problems Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Iterative projection
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 information3 Constitutive Relations: Macroscopic Properties of Matter
EECS 53 Lecture 3 c Kamal Sarabandi Fall 21 All rights reserved 3 Constitutive Relations: Macroscopic Properties of Matter As shown previously, out of the four Maxwell s equations only the Faraday s and
More informationEIGENVALUE PROBLEMS. EIGENVALUE PROBLEMS p. 1/4
EIGENVALUE PROBLEMS EIGENVALUE PROBLEMS p. 1/4 EIGENVALUE PROBLEMS p. 2/4 Eigenvalues and eigenvectors Let A C n n. Suppose Ax = λx, x 0, then x is a (right) eigenvector of A, corresponding to the eigenvalue
More informationNonlinear Eigenvalue Problems
115 Nonlinear Eigenvalue Problems Heinrich Voss Hamburg University of Technology 115.1 Basic Properties........................................ 115-2 115.2 Analytic matrix functions.............................
More informationModeling of Thermoelastic Damping in MEMS Resonators
Modeling of Thermoelastic Damping in MEMS Resonators T. Koyama a, D. Bindel b, S. Govindjee a a Dept. of Civil Engineering b Computer Science Division 1 University of California, Berkeley MEMS Resonators
More informationNumerical 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 information1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det
What is the determinant of the following matrix? 3 4 3 4 3 4 4 3 A 0 B 8 C 55 D 0 E 60 If det a a a 3 b b b 3 c c c 3 = 4, then det a a 4a 3 a b b 4b 3 b c c c 3 c = A 8 B 6 C 4 D E 3 Let A be an n n matrix
More informationAn Algorithm for. Nick Higham. Françoise Tisseur. Director of Research School of Mathematics.
An Algorithm for the Research Complete Matters Solution of Quadratic February Eigenvalue 25, 2009 Problems Nick Higham Françoise Tisseur Director of Research School of Mathematics The School University
More informationMATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian.
MATH 304 Linear Algebra Lecture 10: Linear independence. Wronskian. Spanning set Let S be a subset of a vector space V. Definition. The span of the set S is the smallest subspace W V that contains S. If
More informationA numerical method for polynomial eigenvalue problems using contour integral
A numerical method for polynomial eigenvalue problems using contour integral Junko Asakura a Tetsuya Sakurai b Hiroto Tadano b Tsutomu Ikegami c Kinji Kimura d a Graduate School of Systems and Information
More informationarxiv: v1 [math-ph] 4 Mar 2008
arxiv:0803.0456v1 [math-ph] 4 Mar 2008 On Spectral Analysis of an Operator Pencil with Applications to Wave Propagation in Periodic and Frequency Dependent Materials Christian Engström and Markus Richter
More informationEcon Slides from Lecture 7
Econ 205 Sobel Econ 205 - Slides from Lecture 7 Joel Sobel August 31, 2010 Linear Algebra: Main Theory A linear combination of a collection of vectors {x 1,..., x k } is a vector of the form k λ ix i for
More information20D - Homework Assignment 5
Brian Bowers TA for Hui Sun MATH D Homework Assignment 5 November 8, 3 D - Homework Assignment 5 First, I present the list of all matrix row operations. We use combinations of these steps to row reduce
More informationMatrix stabilization using differential equations.
Matrix stabilization using differential equations. Nicola Guglielmi Universitá dell Aquila and Gran Sasso Science Institute, Italia NUMOC-2017 Roma, 19 23 June, 2017 Inspired by a joint work with Christian
More informationLocal convergence analysis of several inexact Newton-type algorithms for general nonlinear eigenvalue problems. Fei Xue and Daniel B.
Local convergence analysis of several inexact Newton-type algorithms for general nonlinear eigenvalue problems Fei Xue and Daniel B. Szyld Report 11-08-09 August 2011 This report is available in the World
More informationWave Propagation in Uniaxial Media. Reflection and Transmission at Interfaces
Lecture 5: Crystal Optics Outline 1 Homogeneous, Anisotropic Media 2 Crystals 3 Plane Waves in Anisotropic Media 4 Wave Propagation in Uniaxial Media 5 Reflection and Transmission at Interfaces Christoph
More informationEigenvalue Problems CHAPTER 1 : PRELIMINARIES
Eigenvalue Problems CHAPTER 1 : PRELIMINARIES Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Preliminaries Eigenvalue problems 2012 1 / 14
More informationComputing Transfer Function Dominant Poles of Large Second-Order Systems
Computing Transfer Function Dominant Poles of Large Second-Order Systems Joost Rommes Mathematical Institute Utrecht University rommes@math.uu.nl http://www.math.uu.nl/people/rommes joint work with Nelson
More informationReduction of nonlinear eigenproblems with JD
Reduction of nonlinear eigenproblems with JD Henk van der Vorst H.A.vanderVorst@math.uu.nl Mathematical Institute Utrecht University July 13, 2005, SIAM Annual New Orleans p.1/16 Outline Polynomial Eigenproblems
More informationVariational Principles for Nonlinear Eigenvalue Problems
Variational Principles for Nonlinear Eigenvalue Problems Heinrich Voss voss@tuhh.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Variational Principles for Nonlinear EVPs
More informationPreconditioned inverse iteration and shift-invert Arnoldi method
Preconditioned inverse iteration and shift-invert Arnoldi method Melina Freitag Department of Mathematical Sciences University of Bath CSC Seminar Max-Planck-Institute for Dynamics of Complex Technical
More informationMath 215 HW #9 Solutions
Math 5 HW #9 Solutions. Problem 4.4.. If A is a 5 by 5 matrix with all a ij, then det A. Volumes or the big formula or pivots should give some upper bound on the determinant. Answer: Let v i be the ith
More informationthat determines x up to a complex scalar of modulus 1, in the real case ±1. Another condition to normalize x is by requesting that
Chapter 3 Newton methods 3. Linear and nonlinear eigenvalue problems When solving linear eigenvalue problems we want to find values λ C such that λi A is singular. Here A F n n is a given real or complex
More informationCHAPTER 5. Basic Iterative Methods
Basic Iterative Methods CHAPTER 5 Solve Ax = f where A is large and sparse (and nonsingular. Let A be split as A = M N in which M is nonsingular, and solving systems of the form Mz = r is much easier than
More informationKnowledge Discovery and Data Mining 1 (VO) ( )
Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory
More informationThe Jordan Normal Form and its Applications
The and its Applications Jeremy IMPACT Brigham Young University A square matrix A is a linear operator on {R, C} n. A is diagonalizable if and only if it has n linearly independent eigenvectors. What happens
More informationSolving Polynomial Eigenproblems by Linearization
Solving Polynomial Eigenproblems by Linearization Nick Higham School of Mathematics University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with D. Steven Mackey and Françoise
More informationKU Leuven Department of Computer Science
Backward error of polynomial eigenvalue problems solved by linearization of Lagrange interpolants Piers W. Lawrence Robert M. Corless Report TW 655, September 214 KU Leuven Department of Computer Science
More informationA Jacobi-Davidson method for two real parameter nonlinear eigenvalue problems arising from delay differential equations
A Jacobi-Davidson method for two real parameter nonlinear eigenvalue problems arising from delay differential equations Heinrich Voss voss@tuhh.de Joint work with Karl Meerbergen (KU Leuven) and Christian
More informationThe quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying
I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem QEP is to find scalars λ and nonzero vectors u satisfying where Qλx = 0, 1.1 Qλ = λ 2 M + λd + K, M, D and K are given
More informationAlgebra Workshops 10 and 11
Algebra Workshops 1 and 11 Suggestion: For Workshop 1 please do questions 2,3 and 14. For the other questions, it s best to wait till the material is covered in lectures. Bilinear and Quadratic Forms on
More informationJUST THE MATHS UNIT NUMBER 9.8. MATRICES 8 (Characteristic properties) & (Similarity transformations) A.J.Hobson
JUST THE MATHS UNIT NUMBER 9.8 MATRICES 8 (Characteristic properties) & (Similarity transformations) by A.J.Hobson 9.8. Properties of eigenvalues and eigenvectors 9.8. Similar matrices 9.8.3 Exercises
More informationDimension. Eigenvalue and eigenvector
Dimension. Eigenvalue and eigenvector Math 112, week 9 Goals: Bases, dimension, rank-nullity theorem. Eigenvalue and eigenvector. Suggested Textbook Readings: Sections 4.5, 4.6, 5.1, 5.2 Week 9: Dimension,
More informationSolving the Polynomial Eigenvalue Problem by Linearization
Solving the Polynomial Eigenvalue Problem by Linearization Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with Ren-Cang Li,
More informationChapters 5 & 6: Theory Review: Solutions Math 308 F Spring 2015
Chapters 5 & 6: Theory Review: Solutions Math 308 F Spring 205. If A is a 3 3 triangular matrix, explain why det(a) is equal to the product of entries on the diagonal. If A is a lower triangular or diagonal
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationMATH 583A REVIEW SESSION #1
MATH 583A REVIEW SESSION #1 BOJAN DURICKOVIC 1. Vector Spaces Very quick review of the basic linear algebra concepts (see any linear algebra textbook): (finite dimensional) vector space (or linear space),
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 informationMathematical Methods wk 2: Linear Operators
John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
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 informationAA242B: MECHANICAL VIBRATIONS
AA242B: MECHANICAL VIBRATIONS 1 / 50 AA242B: MECHANICAL VIBRATIONS Undamped Vibrations of n-dof Systems These slides are based on the recommended textbook: M. Géradin and D. Rixen, Mechanical Vibrations:
More informationMATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,
MATH 205 HOMEWORK #3 OFFICIAL SOLUTION Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. a F = R, V = R 3, b F = R or C, V = F 2, T = T = 9 4 4 8 3 4 16 8 7 0 1
More informationCHAPTER 5 : NONLINEAR EIGENVALUE PROBLEMS
CHAPTER 5 : NONLINEAR EIGENVALUE PROBLEMS Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Mathematics TUHH Heinrich Voss Nonlinear eigenvalue problems Eigenvalue problems
More informationParallelizing large scale time domain electromagnetic inverse problem
Parallelizing large scale time domain electromagnetic inverse problems Eldad Haber with: D. Oldenburg & R. Shekhtman + Emory University, Atlanta, GA + The University of British Columbia, Vancouver, BC,
More informationDefinition (T -invariant subspace) Example. Example
Eigenvalues, Eigenvectors, Similarity, and Diagonalization We now turn our attention to linear transformations of the form T : V V. To better understand the effect of T on the vector space V, we begin
More informationEigenvalues and Eigenvectors. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB
Eigenvalues and Eigenvectors Consider the equation A x = λ x, where A is an nxn matrix. We call x (must be non-zero) an eigenvector of A if this equation can be solved for some value of λ. We call λ an
More informationA Jacobi Davidson-type projection method for nonlinear eigenvalue problems
A Jacobi Davidson-type projection method for nonlinear eigenvalue problems Timo Betce and Heinrich Voss Technical University of Hamburg-Harburg, Department of Mathematics, Schwarzenbergstrasse 95, D-21073
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
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 informationStudy guide - Math 220
Study guide - Math 220 November 28, 2012 1 Exam I 1.1 Linear Equations An equation is linear, if in the form y + p(t)y = q(t). Introducing the integrating factor µ(t) = e p(t)dt the solutions is then in
More informationEigenvector error bound and perturbation for nonlinear eigenvalue problems
Eigenvector error bound and perturbation for nonlinear eigenvalue problems Yuji Nakatsukasa School of Mathematics University of Tokyo Joint work with Françoise Tisseur Workshop on Nonlinear Eigenvalue
More informationMATH 56A: STOCHASTIC PROCESSES CHAPTER 0
MATH 56A: STOCHASTIC PROCESSES CHAPTER 0 0. Chapter 0 I reviewed basic properties of linear differential equations in one variable. I still need to do the theory for several variables. 0.1. linear differential
More informationJUST THE MATHS SLIDES NUMBER 9.6. MATRICES 6 (Eigenvalues and eigenvectors) A.J.Hobson
JUST THE MATHS SLIDES NUMBER 96 MATRICES 6 (Eigenvalues and eigenvectors) by AJHobson 96 The statement of the problem 962 The solution of the problem UNIT 96 - MATRICES 6 EIGENVALUES AND EIGENVECTORS 96
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationLinear Algebra M1 - FIB. Contents: 5. Matrices, systems of linear equations and determinants 6. Vector space 7. Linear maps 8.
Linear Algebra M1 - FIB Contents: 5 Matrices, systems of linear equations and determinants 6 Vector space 7 Linear maps 8 Diagonalization Anna de Mier Montserrat Maureso Dept Matemàtica Aplicada II Translation:
More informationABSTRACT. Name: Daniil R. Sarkissian Department: Mathematical Sciences
ABSTRACT Name: Daniil R. Sarkissian Department: Mathematical Sciences Title: Theory and Computations of Partial Eigenvalue and Eigenstructure Assignment Problems in Matrix Second-Order and Distributed-Parameter
More informationJust solve this: macroscopic Maxwell s equations
ust solve this: macroscopic Maxwell s equations Faraday: E = B t constitutive equations (here, linear media): Ampere: H = D t + (nonzero frequency) Gauss: D = ρ B = 0 D = ε E B = µ H electric permittivity
More informationA comparison of solvers for quadratic eigenvalue problems from combustion
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS [Version: 2002/09/18 v1.01] A comparison of solvers for quadratic eigenvalue problems from combustion C. Sensiau 1, F. Nicoud 2,, M. van Gijzen 3,
More informationThe matrix sign function
The matrix sign function 1 Re x > 0, sign(x) = 1 Re x < 0, undefined Re x = 0. Suppose the Jordan form of A is reblocked as A = [ V 1 V 2 ] [ J 1 J 2 ] [ V 1 V 2 ] 1, where J 1 contains all eigenvalues
More informationCS 143 Linear Algebra Review
CS 143 Linear Algebra Review Stefan Roth September 29, 2003 Introductory Remarks This review does not aim at mathematical rigor very much, but instead at ease of understanding and conciseness. Please see
More informationThe Behavioral Approach to Systems Theory
The Behavioral Approach to Systems Theory... and DAEs Patrick Kürschner July 21, 2010 1/31 Patrick Kürschner The Behavioral Approach to Systems Theory Outline 1 Introduction 2 Mathematical Models 3 Linear
More information1. Reminder: E-Dynamics in homogenous media and at interfaces
0. Introduction 1. Reminder: E-Dynamics in homogenous media and at interfaces 2. Photonic Crystals 2.1 Introduction 2.2 1D Photonic Crystals 2.3 2D and 3D Photonic Crystals 2.4 Numerical Methods 2.4.1
More informationPutzer s Algorithm. Norman Lebovitz. September 8, 2016
Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),
More informationSolution of Matrix Eigenvalue Problem
Outlines October 12, 2004 Outlines Part I: Review of Previous Lecture Part II: Review of Previous Lecture Outlines Part I: Review of Previous Lecture Part II: Standard Matrix Eigenvalue Problem Other Forms
More informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
More informationContour Integral Method for the Simulation of Accelerator Cavities
Contour Integral Method for the Simulation of Accelerator Cavities V. Pham-Xuan, W. Ackermann and H. De Gersem Institut für Theorie Elektromagnetischer Felder DESY meeting (14.11.2017) November 14, 2017
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationEigenvalues, Eigenvectors, and an Intro to PCA
Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Changing Basis We ve talked so far about re-writing our data using a new set of variables, or a new basis.
More information