How to Detect Definite Hermitian Pairs
|
|
- Melvin Anderson
- 5 years ago
- Views:
Transcription
1 How to Detect Definite Hermitian Pairs Françoise Tisseur School of Mathematics The University of Manchester Joint work with Chun-Hua Guo and Nick Higham 17th Scottish Computational Mathematics Symposium, Edinburgh, September 10, 2008
2 Hermitian Positive Definite Matrix A C n n is Hermitian positive definite if x Ax > 0 for all x 0. In the eigenproblem Ax = λx, real, positive eigenvalues, unitary matrix of eigenvectors. Testing for definiteness: attempt a Cholesky decomposition, A = R R. Françoise Tisseur Definite Pairs 2 / 28
3 Generalized Eigenproblem Ax = λbx with A and B Hermitian and B positive definite. Equivalent to the standard eigenproblem Hy B 1/2 AB 1/2 y = λy. All eigenvalues real. A and B are simultaneously diagonalizable. These properties extend to definite pairs (A, B) defined by γ(a, B) := min z C n (z Az) 2 + (z Bz) 2 > 0. z 2 =1 γ is the Crawford number. Françoise Tisseur Definite Pairs 3 / 28
4 Definite Pairs Theorem Let (A, B) be a definite Hermitian pair, and for t R let A(t) = A cos t B sin t, B(t) = A sin t + B cos t. Then there is a t [0, 2π) s.t. B(t) > 0 and γ(a, B) = λ min (B(t)). Note that the eigenproblem A(t)x = µb(t)x share the same e vecs as Ax = λbx and their e vals are related by λ = µ cos t + sint cos t µ sin t. Françoise Tisseur Definite Pairs 4 / 28
5 Why Test for Definiteness? Both theoretical and computational advantages accrue in associated eigenproblem Ax = λbx: eigenvalues are real, A and B simultaneously diagonalizable, if t known s.t. B(t) = A sin t + B cos t > 0 then definiteness of B(t) is exploited to compute e vals. Detect hyperbolicity of quadratic matrix polynomials. Allow CG iterations for nonsymmetric saddle point linear systems. Françoise Tisseur Definite Pairs 5 / 28
6 Hyperbolic Quadratics Q(λ) = λ 2 A + λb + C, A, B, C C n n is hyperbolic if A, B, C Hermitian, A > 0, and (x Bx) 2 > 4(x Ax)(x Cx) for all nonzero x C n. Eigenvalues are all real and semisimple. Q > 0 Q indef. Q < 0 Q indef. Q > 0 + λ 2n λ n+1 λ n λ 1 Q is hyperbolic iff Q(µ) < 0 for some µ R. Françoise Tisseur Definite Pairs 6 / 28
7 Testing for Hyperbolicity Theorem Q(λ) = λ 2 A + λb + C with A, B and C Hermitian and A > 0 is hyperbolic iff the pair (A, B) is definite, where [ ] [ ] C 0 B A A =, B =. 0 A A 0 Proof. (A, B) definite iff αa+βb > 0, α 2 +β 2 = 1. For α 0 we have [ I β αa + βb = I ] [ α αc βb β 2 A 0 ] [ ] α I 0 0 I 0 αa β I I. α So αa + βb is congruent to α diag( Q(µ), A), µ = β/α. Françoise Tisseur Definite Pairs 7 / 28
8 Solving Hyperbolic Q(λ) = λ 2 A + λb + C Let µ s.t. Q(µ) < 0. Hence Q µ (λ) Q(λ + µ) = λ 2 A + λ(b + 2µA) + C + µb + µ 2 A = λ 2 A µ + λb µ + C µ, with C µ = Q(µ) < 0 and A µ = A > 0. The pencils [ ] Aµ 0 λ + 0 C µ [ ] Bµ C µ, λ C µ 0 [ ] [ ] 0 Aµ Aµ 0 + A µ B µ 0 C µ are both Hermitian definite linearizations of Q µ (they share the same e vals as Q). Françoise Tisseur Definite Pairs 8 / 28
9 Linear Systems in Saddle Point Problems Involve matrices of the form [ A B A = B C where A = A T R n n, A > 0 and C = C T R m m, C 0. A is indefinite: it has n positive e vals and rank(c + BA 1 B T ) negative e vals. Theorem (Liesen & Parlett 08) [ ] In 0 If the symmetric pair (A, J ) with J = is definite 0 I m and µ R s.t. A µj > 0 is known then there exists a well-define CG method for solving linear systems with J A. ], Françoise Tisseur Definite Pairs 9 / 28
10 How can we efficiently test whether (A, B) is definite, i.e., whether γ(a, B)>0? γ(a, B) := min (z Az) 2 + (z Bz) 2 z C n z 2 =1 = min z (A + ib)z z C n z 2 =1 ( ) = max max 0 t 2π λ min (A cos t + B sin t), 0. One-dimensional global optimization problem. Unfortunately, not a convex (or concave) maximization problem. Françoise Tisseur Definite Pairs 10 / 28
11 Numerical Methods for Testing Definiteness J-orthogonal Jacobi algorithm [Veselic, 1993]. Bisection algorithm [Higham, T & Van Dooren, 2002]. Level set algorithm [Higham, T & Van Dooren, 2002]. Chebfun package. Crawford & Moon algorithm [1983], improved [Guo, Higham & T, 2008]. Françoise Tisseur Definite Pairs 11 / 28
12 Finding Global Maximum Our aim is to compute max λ ( ) min(a cos t + B sin t) = maxλ min M(z), 0 t 2π z =1 where M(z) = (z 1 C + zc )/2 and C = A + ib Given ξ, is ξ = λ min (M(z)) for some z on unit circle? Françoise Tisseur Definite Pairs 12 / 28
13 Finding Global Maximum Our aim is to compute max λ ( ) min(a cos t + B sin t) = maxλ min M(z), 0 t 2π z =1 where M(z) = (z 1 C + zc )/2 and C = A + ib Given ξ, is ξ = λ min (M(z)) for some z on unit circle? For nonzero z, det(m(z) ξi) = 0 det(c 2ξzI + z 2 C ) = 0. For each e vals z j of C 2ξzI + z 2 C on the unit circle test if λ min (M(z j )) = ξ. Françoise Tisseur Definite Pairs 12 / 28
14 Bisection Algorithm Given Hermitian pair (A, B), determine whether (A, B) is definite. a = max{λ min (A), λ max (A),λ min (B), λ max (B)}. b = σ n ([ A B ]) if a > 0, quit (pair is definite), end if b 0, quit with a = b = 0 (pair is indefinite), end while b a > tol ξ = (a + b)/2 Compute eigenvalues z j of Q(z). If λ min (M(z j )) = ξ for some z j on unit circle a = ξ if a > 0,quit (pair is definite), end else b = ξ if b 0, quit with a = b = 0 (pair is indefinite), end end end Françoise Tisseur Definite Pairs 13 / 28
15 Bisection Algorithm, cont. Produces a bracket for γ(a, B) shrinking to zero. Does not compute t [0, 2π) such that B(t) > 0. Each iteration requires solution of a QEP and up to 2n Hermitian eigenproblems so is expensive. Françoise Tisseur Definite Pairs 14 / 28
16 Level Set Approach E vals of B(t) = A cos t + B sin t are piecewise analytic functions of t. Theorem (A, B) is definite if B(t) has exactly n consecutive increasing zero crossings followed by n consecutive decreasing zero crossings. The zero crossings are determined by the unit modulus e vals of Q(z) = z 2 C + C, where C = A + ib. The sign of t λ j( B(t) ) tk determines if λ j ( B(t) ) increases or decreases at each zero crossing. Françoise Tisseur Definite Pairs 15 / 28
17 Example A and B are random, 5 5, real and symmetric. Plot eigenvalues of A cosθ + B sin θ for θ [ π,π]: Eigenvalues of Acosθ+Bsinθ θ Françoise Tisseur Definite Pairs 16 / 28
18 Level Set Algorithm Given Hermitian pair (A, B), determine whether (A, B) is definite. C = A + ib Compute the eigenvalues of z 2 C + C. If there are 2n eigenvalues of unit modulus Compute the derivatives (using associated e vecs). If there are n consecutive increasing and n consecutive decreasing zero crossings The pair is definite; return end end The pair is not definite. Return an interval of t such that B(t) > 0. Françoise Tisseur Definite Pairs 17 / 28
19 Chebfun ( ) γ(a, B) = max max λ min(a cos t + B sin t), 0. 0 t 2π Use new Chebfun toolbox to solve the 1-variable minimization problem. f min(eig(cos(t)*a+sin(t)*b)); fc = chebfun(f,[0 2*pi]); [maxf,t opt] = max(fc); test = (maxf > 0) gamma = max(max(maxf,0)) Françoise Tisseur Definite Pairs 18 / 28
20 Algorithm of Crawford and Moon Bisection-like algorithm. Received little attention in the literature. PDFIND: Algorithm 646, ACM TOMS, Lack of clarity in the derivation and statement of alg and in the explanation of its properties. Our contributions: Give clearer derivation with no assumptions on definiteness of (A, B). Prove convergence for all (A, B) definite or not. Show that PDFIND is numerically unreliable. Formulate an improved algorithm. Françoise Tisseur Definite Pairs 19 / 28
21 Notations For a Hermitian pencil (A, B) define f(x) = x (A + ib)x x (A + ib)x, x Cn, x (A + ib)x 0. f(x) lies on the unit circle. For a, b of unit modulus, arc[a, b] is the shorter arc on the unit circle connecting a and b. θ[a, b] is the length of the arc [a, b]. When a = b, define arc[a, b] = π. Françoise Tisseur Definite Pairs 20 / 28
22 Range of f(x) = x (A + ib)x/ x (A + ib)x Property 1: When (A, B) is definite the range of f denoted by arc[ α, β] is a closed arc on the unit circle of length θ[ α, β] < π. Property 2 If c = e it and B(t) = A sin t + B cos t 0 then for any x 0 with x B(t)x 0 and x (A + ib)x 0, we have π θ[f(x), c] π/2. Property 3 If (A, B) is definite and c = cos t + i sin t is the midpoint of arc[ α, β] then B( t) > 0. Françoise Tisseur Definite Pairs 21 / 28
23 Pictorial Explanation b a c a b Françoise Tisseur Definite Pairs 22 / 28
24 Crawford and Moon algorithm Input: A and B; starting arc [a, b]. Output: t [0, 2π) when (A, B) definite. (Main loop) 1 c = (a + b)/ a + b = sin t + i cos t 2 if B(t) > 0, quit (pair definite), end 3 Find unit norm vector x s.t. x B(t)x 0. 4 d = f(x) 5 if d is in arc[ a, b] and d a 6 b = d, goto line 1 7 elseif d is in arc[a, b] and d b 8 a = d, goto line 1 9 else 10 quit (pair is indefinite) 11 end Françoise Tisseur Definite Pairs 23 / 28
25 Convergence Lemma On the kth step, θ[a k, b k ] π(1 2 k ), regardless of whether the pair (A, B) is definite or not. Unless alg. converges in finite number of steps, arc[a k, b k ] [π(1 2 k ),π] k. Hence sequence converges to π. (A, B) declared numerically indefinite if no determination is made in k steps and 2 k u (unit roundoff). Françoise Tisseur Definite Pairs 24 / 28
26 An Improved Crawford and Moon algorithm New starting phase; no assumption on definiteness of (A, B). (Main loop) 1 c = (a + b)/ a + b = sin t + i cos t 2 if B(t) > 0, quit (pair definite), end 3 Find unit norm vector x s.t. x B(t)x 0. 4 if x (A + ib)x = 0, quit (pair is indefinite), end 5 d = f(x) 6 if d is in arc[ a, b] and d a 7 b = d, goto line 1 8 elseif d is in arc[a, b] and d b 9 a = d, goto line 1 10 else 11 quit (pair is indefinite) 12 end Françoise Tisseur Definite Pairs 25 / 28
27 Testing Definiteness, Negative Curvature Is B(t) = A sin t + B cos t > 0? If not how to compute x 0 s.t. x B(t)x 0? Definiteness test: run Cholesky alg. until it completes or encounter a negative or zero pivot. Direction of negative curvature: at step k of Cholesky, [ ] k n k ( k n k R C = 11 ( ) k 0 0 R 11 R 12 + n k 0 S k R 12 If s (k) 11[ 0 factorization breaks down and R 1 x = 11 R ] 12 e I 1 satisfies x B(t)x 0. can suffers from numerical instability. Françoise Tisseur Definite Pairs 26 / 28 ),
28 Unreliability of Cholesky Approach: Example R = [ ] , A = R R + e 3 e T 4 + e 4 e T 3, B = diag(0, 1, 1, 1). γ(a, B) = 0.75 so the pair is safely definite. Crawford and Moon alg. attempt to Cholesky factorize B(π/2) = A which is indefinite. In exact arithmetic, x B(pi/2)x = 0. In floating point arithmetic, the computed x is s.t. fl( x B(π/2)x) = > 0, and the alg. breaks down. Françoise Tisseur Definite Pairs 27 / 28
29 Concluding Remarks Determining whether a given pair (A, B) is definite or not is not straightforward. Both theoretical and computational advantages accrue t such that A sin t + B cos t > 0 is known. The (modified) Crawford and Moon algorithm is more efficient than other numerical methods for testing definiteness. Work in progress: numerically stable way of computing the direction of negative curvature. For papers and Eprints, Françoise Tisseur Definite Pairs 28 / 28
Detecting a definite Hermitian pair and a hyperbolic or elliptic quadratic eigenvalue problem, and associated nearness problems
Linear Algebra and its Applications 351 352 (2002) 455 474 www.elsevier.com/locate/laa Detecting a definite Hermitian pair and a hyperbolic or elliptic quadratic eigenvalue problem, and associated nearness
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 informationHermitian Matrix Polynomials with Real Eigenvalues of Definite Type. Part I: Classification. Al-Ammari, Maha and Tisseur, Francoise
Hermitian Matrix Polynomials with Real Eigenvalues of Definite Type. Part I: Classification Al-Ammari, Maha and Tisseur, Francoise 2010 MIMS EPrint: 2010.9 Manchester Institute for Mathematical Sciences
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 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 informationLinear Algebra and its Applications
Linear Algebra and its Applications 436 (2012) 3954 3973 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Hermitian matrix polynomials
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 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 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 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 informationAvailable online at ScienceDirect. Procedia Engineering 100 (2015 ) 56 63
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 100 (2015 ) 56 63 25th DAAAM International Symposium on Intelligent Manufacturing and Automation, DAAAM 2014 Definite Quadratic
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 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 informationEigenvalue Problems. Eigenvalue problems occur in many areas of science and engineering, such as structural analysis
Eigenvalue Problems Eigenvalue problems occur in many areas of science and engineering, such as structural analysis Eigenvalues also important in analyzing numerical methods Theory and algorithms apply
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 4 Eigenvalue Problems Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationDefinite Matrix Polynomials and their Linearization by Definite Pencils. Higham, Nicholas J. and Mackey, D. Steven and Tisseur, Françoise
Definite Matrix Polynomials and their Linearization by Definite Pencils Higham Nicholas J and Mackey D Steven and Tisseur Françoise 2007 MIMS EPrint: 200797 Manchester Institute for Mathematical Sciences
More informationNick Higham. Director of Research School of Mathematics
Exploiting Research Tropical Matters Algebra in Numerical February 25, Linear 2009 Algebra Nick Higham Françoise Tisseur Director of Research School of Mathematics The School University of Mathematics
More informationAPPLIED NUMERICAL LINEAR ALGEBRA
APPLIED NUMERICAL LINEAR ALGEBRA James W. Demmel University of California Berkeley, California Society for Industrial and Applied Mathematics Philadelphia Contents Preface 1 Introduction 1 1.1 Basic Notation
More informationin Numerical Linear Algebra
Exploiting ResearchTropical MattersAlgebra in Numerical Linear Algebra February 25, 2009 Nick Françoise Higham Tisseur Director Schoolof ofresearch Mathematics The School University of Mathematics of Manchester
More informationA conjugate gradient-based algorithm for large-scale quadratic programming problem with one quadratic constraint
Noname manuscript No. (will be inserted by the editor) A conjugate gradient-based algorithm for large-scale quadratic programming problem with one quadratic constraint A. Taati, M. Salahi the date of receipt
More informationAnalysis of the Cholesky Method with Iterative Refinement for Solving the Symmetric Definite Generalized Eigenproblem
Analysis of the Cholesky Method with Iterative Refinement for Solving the Symmetric Definite Generalized Eigenproblem Davies, Philip I. and Higham, Nicholas J. and Tisseur, Françoise 2001 MIMS EPrint:
More informationSolving large scale eigenvalue problems
arge scale eigenvalue problems, Lecture 5, March 23, 2016 1/30 Lecture 5, March 23, 2016: The QR algorithm II http://people.inf.ethz.ch/arbenz/ewp/ Peter Arbenz Computer Science Department, ETH Zürich
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 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 informationLecture 10: Eigenvectors and eigenvalues (Numerical Recipes, Chapter 11)
Lecture 1: Eigenvectors and eigenvalues (Numerical Recipes, Chapter 11) The eigenvalue problem, Ax= λ x, occurs in many, many contexts: classical mechanics, quantum mechanics, optics 22 Eigenvectors and
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 informationSparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices
Sparsity-Preserving Difference of Positive Semidefinite Matrix Representation of Indefinite Matrices Jaehyun Park June 1 2016 Abstract We consider the problem of writing an arbitrary symmetric matrix as
More informationBackward Error of Polynomial Eigenproblems Solved by Linearization. Higham, Nicholas J. and Li, Ren-Cang and Tisseur, Françoise. MIMS EPrint: 2006.
Backward Error of Polynomial Eigenproblems Solved by Linearization Higham, Nicholas J and Li, Ren-Cang and Tisseur, Françoise 2007 MIMS EPrint: 2006137 Manchester Institute for Mathematical Sciences School
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 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 informationON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS
ON SYLVESTER S LAW OF INERTIA FOR NONLINEAR EIGENVALUE PROBLEMS ALEKSANDRA KOSTIĆ AND HEINRICH VOSS Key words. eigenvalue, variational characterization, principle, Sylvester s law of inertia AMS subject
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 informationAn Algorithm for the Complete Solution of Quadratic Eigenvalue Problems. Hammarling, Sven and Munro, Christopher J. and Tisseur, Francoise
An Algorithm for the Complete Solution of Quadratic Eigenvalue Problems Hammarling, Sven and Munro, Christopher J. and Tisseur, Francoise 2011 MIMS EPrint: 2011.86 Manchester Institute for Mathematical
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 informationEIGENVALUE PROBLEMS (EVP)
EIGENVALUE PROBLEMS (EVP) (Golub & Van Loan: Chaps 7-8; Watkins: Chaps 5-7) X.-W Chang and C. C. Paige PART I. EVP THEORY EIGENVALUES AND EIGENVECTORS Let A C n n. Suppose Ax = λx with x 0, then x is a
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lectures 13 15, Tuesday 8 th and Friday 11 th November 016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE
More informationNumerical Algorithms. IE 496 Lecture 20
Numerical Algorithms IE 496 Lecture 20 Reading for This Lecture Primary Miller and Boxer, Pages 124-128 Forsythe and Mohler, Sections 1 and 2 Numerical Algorithms Numerical Analysis So far, we have looked
More informationMath 489AB Exercises for Chapter 2 Fall Section 2.3
Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial
More informationSymmetric Matrices and Eigendecomposition
Symmetric Matrices and Eigendecomposition Robert M. Freund January, 2014 c 2014 Massachusetts Institute of Technology. All rights reserved. 1 2 1 Symmetric Matrices and Convexity of Quadratic Functions
More informationConjugate Gradient Method
Conjugate Gradient Method direct and indirect methods positive definite linear systems Krylov sequence spectral analysis of Krylov sequence preconditioning Prof. S. Boyd, EE364b, Stanford University Three
More information1 Positive and completely positive linear maps
Part 5 Functions Matrices We study functions on matrices related to the Löwner (positive semidefinite) ordering on positive semidefinite matrices that A B means A B is positive semi-definite. Note that
More informationMICHIEL E. HOCHSTENBACH
VARIATIONS ON HARMONIC RAYLEIGH RITZ FOR STANDARD AND GENERALIZED EIGENPROBLEMS MICHIEL E. HOCHSTENBACH Abstract. We present several variations on the harmonic Rayleigh Ritz method. First, we introduce
More informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
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 informationNonlinear palindromic eigenvalue problems and their numerical solution
Nonlinear palindromic eigenvalue problems and their numerical solution TU Berlin DFG Research Center Institut für Mathematik MATHEON IN MEMORIAM RALPH BYERS Polynomial eigenvalue problems k P(λ) x = (
More informationSolving large scale eigenvalue problems
arge scale eigenvalue problems, Lecture 4, March 14, 2018 1/41 Lecture 4, March 14, 2018: The QR algorithm http://people.inf.ethz.ch/arbenz/ewp/ Peter Arbenz Computer Science Department, ETH Zürich E-mail:
More information5.3 The Power Method Approximation of the Eigenvalue of Largest Module
192 5 Approximation of Eigenvalues and Eigenvectors 5.3 The Power Method The power method is very good at approximating the extremal eigenvalues of the matrix, that is, the eigenvalues having largest and
More informationPrincipal pivot transforms of quasidefinite matrices and semidefinite Lagrangian subspaces. Federico Poloni and Nataša Strabi.
Principal pivot transforms of quasidefinite matrices and semidefinite Lagrangian subspaces Federico Poloni and Nataša Strabi October 2015 MIMS EPrint: 2015.92 Manchester Institute for Mathematical Sciences
More informationHOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)
HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe
More informationTwo hours. To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER. 29 May :45 11:45
Two hours MATH20602 To be provided by Examinations Office: Mathematical Formula Tables. THE UNIVERSITY OF MANCHESTER NUMERICAL ANALYSIS 1 29 May 2015 9:45 11:45 Answer THREE of the FOUR questions. If more
More information5.3.5 The eigenvalues are 3, 2, 3 (i.e., the diagonal entries of D) with corresponding eigenvalues. Null(A 3I) = Null( ), 0 0
535 The eigenvalues are 3,, 3 (ie, the diagonal entries of D) with corresponding eigenvalues,, 538 The matrix is upper triangular so the eigenvalues are simply the diagonal entries, namely 3, 3 The corresponding
More informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
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 informationECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues. Prof. Peter Bermel January 23, 2013
ECE 595, Section 10 Numerical Simulations Lecture 7: Optimization and Eigenvalues Prof. Peter Bermel January 23, 2013 Outline Recap from Friday Optimization Methods Brent s Method Golden Section Search
More informationMath Homework 8 (selected problems)
Math 102 - Homework 8 (selected problems) David Lipshutz Problem 1. (Strang, 5.5: #14) In the list below, which classes of matrices contain A and which contain B? 1 1 1 1 A 0 0 1 0 0 0 0 1 and B 1 1 1
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
More informationMath 489AB Exercises for Chapter 1 Fall Section 1.0
Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T
More informationHOMOGENEOUS JACOBI DAVIDSON. 1. Introduction. We study a homogeneous Jacobi Davidson variant for the polynomial eigenproblem
HOMOGENEOUS JACOBI DAVIDSON MICHIEL E. HOCHSTENBACH AND YVAN NOTAY Abstract. We study a homogeneous variant of the Jacobi Davidson method for the generalized and polynomial eigenvalue problem. Special
More informationNumerical Linear Algebra And Its Applications
Numerical Linear Algebra And Its Applications Xiao-Qing JIN 1 Yi-Min WEI 2 August 29, 2008 1 Department of Mathematics, University of Macau, Macau, P. R. China. 2 Department of Mathematics, Fudan University,
More informationUniversity of Colorado Denver Department of Mathematical and Statistical Sciences Applied Linear Algebra Ph.D. Preliminary Exam January 23, 2015
University of Colorado Denver Department of Mathematical and Statistical Sciences Applied Linear Algebra PhD Preliminary Exam January 23, 2015 Name: Exam Rules: This exam lasts 4 hours and consists of
More information7. Symmetric Matrices and Quadratic Forms
Linear Algebra 7. Symmetric Matrices and Quadratic Forms CSIE NCU 1 7. Symmetric Matrices and Quadratic Forms 7.1 Diagonalization of symmetric matrices 2 7.2 Quadratic forms.. 9 7.4 The singular value
More informationStructured Mapping Problems for Matrices Associated with Scalar Products Part I: Lie and Jordan Algebras
Structured Mapping Problems for Matrices Associated with Scalar Products Part I: Lie and Jordan Algebras Mackey, D. Steven and Mackey, Niloufer and Tisseur, Françoise 27 MIMS EPrint: 26.44 Manchester Institute
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationThe Singular Value Decomposition
The Singular Value Decomposition Philippe B. Laval KSU Fall 2015 Philippe B. Laval (KSU) SVD Fall 2015 1 / 13 Review of Key Concepts We review some key definitions and results about matrices that will
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 informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationBasic Concepts in Matrix Algebra
Basic Concepts in Matrix Algebra An column array of p elements is called a vector of dimension p and is written as x p 1 = x 1 x 2. x p. The transpose of the column vector x p 1 is row vector x = [x 1
More information, b = 0. (2) 1 2 The eigenvectors of A corresponding to the eigenvalues λ 1 = 1, λ 2 = 3 are
Quadratic forms We consider the quadratic function f : R 2 R defined by f(x) = 2 xt Ax b T x with x = (x, x 2 ) T, () where A R 2 2 is symmetric and b R 2. We will see that, depending on the eigenvalues
More informationComputing the Action of the Matrix Exponential
Computing the Action of the Matrix Exponential Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with Awad H. Al-Mohy 16th ILAS
More informationComputing Matrix Functions by Iteration: Convergence, Stability and the Role of Padé Approximants
Computing Matrix Functions by Iteration: Convergence, Stability and the Role of Padé Approximants Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/
More informationThere are six more problems on the next two pages
Math 435 bg & bu: Topics in linear algebra Summer 25 Final exam Wed., 8/3/5. Justify all your work to receive full credit. Name:. Let A 3 2 5 Find a permutation matrix P, a lower triangular matrix L with
More informationComputable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization
Computable error bounds for nonlinear eigenvalue problems allowing for a minmax characterization Heinrich Voss voss@tuhh.de Joint work with Kemal Yildiztekin Hamburg University of Technology Institute
More informationANALYSIS OF STRUCTURED POLYNOMIAL EIGENVALUE PROBLEMS
ANALYSIS OF STRUCTURED POLYNOMIAL EIGENVALUE PROBLEMS A thesis submitted to the University of Manchester for the degree of Doctor of Philosophy in the Faculty of Engineering and Physical Sciences 2011
More informationLINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12,
LINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12, 2000 74 6 Summary Here we summarize the most important information about theoretical and numerical linear algebra. MORALS OF THE STORY: I. Theoretically
More informationJ-Orthogonal Matrices: Properties and Generation
SIAM REVIEW Vol. 45,No. 3,pp. 504 519 c 2003 Society for Industrial and Applied Mathematics J-Orthogonal Matrices: Properties and Generation Nicholas J. Higham Abstract. A real, square matrix Q is J-orthogonal
More informationCHARACTERIZATIONS. is pd/psd. Possible for all pd/psd matrices! Generating a pd/psd matrix: Choose any B Mn, then
LECTURE 6: POSITIVE DEFINITE MATRICES Definition: A Hermitian matrix A Mn is positive definite (pd) if x Ax > 0 x C n,x 0 A is positive semidefinite (psd) if x Ax 0. Definition: A Mn is negative (semi)definite
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationMath (P)refresher Lecture 8: Unconstrained Optimization
Math (P)refresher Lecture 8: Unconstrained Optimization September 2006 Today s Topics : Quadratic Forms Definiteness of Quadratic Forms Maxima and Minima in R n First Order Conditions Second Order Conditions
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 informationPositive Definite Matrix
1/29 Chia-Ping Chen Professor Department of Computer Science and Engineering National Sun Yat-sen University Linear Algebra Positive Definite, Negative Definite, Indefinite 2/29 Pure Quadratic Function
More informationSolving Symmetric Semi-definite (ill-conditioned) Generalized Eigenvalue Problems
Solving Symmetric Semi-definite (ill-conditioned) Generalized Eigenvalue Problems Zhaojun Bai University of California, Davis Berkeley, August 19, 2016 Symmetric definite generalized eigenvalue problem
More information1 Overview. 2 A Characterization of Convex Functions. 2.1 First-order Taylor approximation. AM 221: Advanced Optimization Spring 2016
AM 221: Advanced Optimization Spring 2016 Prof. Yaron Singer Lecture 8 February 22nd 1 Overview In the previous lecture we saw characterizations of optimality in linear optimization, and we reviewed the
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Eigenvalue Problems; Similarity Transformations Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical Analysis I 1 / 18 Eigenvalue
More informationLecture Notes 6: Dynamic Equations Part B: Second and Higher-Order Linear Difference Equations in One Variable
University of Warwick, EC9A0 Maths for Economists Peter J. Hammond 1 of 56 Lecture Notes 6: Dynamic Equations Part B: Second and Higher-Order Linear Difference Equations in One Variable Peter J. Hammond
More informationHomework 1 Elena Davidson (B) (C) (D) (E) (F) (G) (H) (I)
CS 106 Spring 2004 Homework 1 Elena Davidson 8 April 2004 Problem 1.1 Let B be a 4 4 matrix to which we apply the following operations: 1. double column 1, 2. halve row 3, 3. add row 3 to row 1, 4. interchange
More informationTropical roots as approximations to eigenvalues of matrix polynomials. Noferini, Vanni and Sharify, Meisam and Tisseur, Francoise
Tropical roots as approximations to eigenvalues of matrix polynomials Noferini, Vanni and Sharify, Meisam and Tisseur, Francoise 2014 MIMS EPrint: 2014.16 Manchester Institute for Mathematical Sciences
More informationIs there a Small Skew Cayley Transform with Zero Diagonal?
Is there a Small Skew Cayley Transform with Zero Diagonal? Abstract The eigenvectors of an Hermitian matrix H are the columns of some complex unitary matrix Q. For any diagonal unitary matrix Ω the columns
More informationLinear Algebra And Its Applications Chapter 6. Positive Definite Matrix
Linear Algebra And Its Applications Chapter 6. Positive Definite Matrix KAIST wit Lab 2012. 07. 10 남성호 Introduction The signs of the eigenvalues can be important. The signs can also be related to the minima,
More informationOrthogonalization and least squares methods
Chapter 3 Orthogonalization and least squares methods 31 QR-factorization (QR-decomposition) 311 Householder transformation Definition 311 A complex m n-matrix R = [r ij is called an upper (lower) triangular
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 informationAlgorithms and Perturbation Theory for Matrix Eigenvalue Problems and the SVD
Algorithms and Perturbation Theory for Matrix Eigenvalue Problems and the SVD Yuji Nakatsukasa PhD dissertation University of California, Davis Supervisor: Roland Freund Householder 2014 2/28 Acknowledgment
More informationMATH 5720: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 2018
MATH 57: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 18 1 Global and Local Optima Let a function f : S R be defined on a set S R n Definition 1 (minimizers and maximizers) (i) x S
More information1. Select the unique answer (choice) for each problem. Write only the answer.
MATH 5 Practice Problem Set Spring 7. Select the unique answer (choice) for each problem. Write only the answer. () Determine all the values of a for which the system has infinitely many solutions: x +
More informationPage 712. Lecture 42: Rotations and Orbital Angular Momentum in Two Dimensions Date Revised: 2009/02/04 Date Given: 2009/02/04
Page 71 Lecture 4: Rotations and Orbital Angular Momentum in Two Dimensions Date Revised: 009/0/04 Date Given: 009/0/04 Plan of Attack Section 14.1 Rotations and Orbital Angular Momentum: Plan of Attack
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 informationElements of Floating-point Arithmetic
Elements of Floating-point Arithmetic Sanzheng Qiao Department of Computing and Software McMaster University July, 2012 Outline 1 Floating-point Numbers Representations IEEE Floating-point Standards Underflow
More informationHere each term has degree 2 (the sum of exponents is 2 for all summands). A quadratic form of three variables looks as
Reading [SB], Ch. 16.1-16.3, p. 375-393 1 Quadratic Forms A quadratic function f : R R has the form f(x) = a x. Generalization of this notion to two variables is the quadratic form Q(x 1, x ) = a 11 x
More informationDiagonalization of Matrix
of Matrix King Saud University August 29, 2018 of Matrix Table of contents 1 2 of Matrix Definition If A M n (R) and λ R. We say that λ is an eigenvalue of the matrix A if there is X R n \ {0} such that
More informationChapter 6. Eigenvalues. Josef Leydold Mathematical Methods WS 2018/19 6 Eigenvalues 1 / 45
Chapter 6 Eigenvalues Josef Leydold Mathematical Methods WS 2018/19 6 Eigenvalues 1 / 45 Closed Leontief Model In a closed Leontief input-output-model consumption and production coincide, i.e. V x = x
More information