Functions of a Matrix: Theory and Computation
|
|
- Lily Nichols
- 5 years ago
- Views:
Transcription
1 Functions of a Matrix: Theory and Computation Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk Landscapes in Mathematical Science, University of Bath, November
2 Outline 1 Definition of f(a) 2 Motivation and MATLAB 3 e A and its Frechét derivative 4 A 1/2 : Modified Newton Methods MIMS Nick Higham Functions of a Matrix 2 / 45
3 Function of a Matrix f : C n n C n n for an underlying scalar function f. These are not matrix functions: trace(a), det(a). The adjugate (or adjoint) matrix. Transfer function f(t) = B(tI A) 1 C. sin A = (sin a ij ). These are matrix functions: e A = I + A + A2 2! +. log(i + A) = A A2 2 + A3 3 A 1, A 1/2. +, ρ(a) < 1. MIMS Nick Higham Functions of a Matrix 3 / 45
4 Multiplicity of Definitions There have been proposed in the literature since 1880 eight distinct definitions of a matric function, by Weyr, Sylvester and Buchheim, Giorgi, Cartan, Fantappiè, Cipolla, Schwerdtfeger and Richter. R. F. Rinehart, The Equivalence of Definitions of a Matric Function, Amer. Math. Monthly (1955) MIMS Nick Higham Functions of a Matrix 4 / 45
5 Jordan Canonical Form Z 1 AZ = J = diag(j 1,...,J p ), J }{{} k = m k m k λ k 1 λ k λ k Definition f(a) = Zf(J)Z 1 = Zdiag(f(J k ))Z 1, f(λ k ) f f (m k 1) )(λ k ) (λ k )... (m k 1)! f(j k ) = f(λ k ) f (λ k ) f(λ k ) MIMS Nick Higham Functions of a Matrix 5 / 45
6 The Formula for f(j k ) Write J k = λ k I + E k C m k m k. For mk = 3 we have E k = , Ek 2 = , Ek 3 = Assume f has Taylor expansion f(t) = f(λ k ) + f (λ k )(t λ k ) + + f (j) (λ k )(t λ k ) j j! +. Then f(j k ) = f(λ k )I + f (λ k )E k + + f (m k 1) (λ k )E m k 1 k. (m k 1)! MIMS Nick Higham Functions of a Matrix 6 / 45
7 Interpolation Definition (Sylvester, 1883; Buchheim, 1886) Distinct e vals λ 1,...,λ s, n i = max size of Jordan blocks for λ i. Then f(a) = r(a), where r is unique Hermite interpolating poly of degree < s i=1 n i satisfying r (j) (λ i ) = f (j) (λ i ), j = 0: n i 1, i = 1: s. MIMS Nick Higham Functions of a Matrix 7 / 45
8 Interpolation Definition (Sylvester, 1883; Buchheim, 1886) Distinct e vals λ 1,...,λ s, n i = max size of Jordan blocks for λ i. Then f(a) = r(a), where r is unique Hermite interpolating poly of degree < s i=1 n i satisfying r (j) (λ i ) = f (j) (λ i ), j = 0: n i 1, i = 1: s. [ ] 2 2 Example. Let f(t) = t 1/2, A =, λ(a) = {1, 4}. 1 3 Taking +ve square roots, r(t) = f(1) t f(4) t = 1 (t + 2). 3 [ A 1/2 = r(a) = 1 3 (A + 2I) = 1 3 ]. MIMS Nick Higham Functions of a Matrix 7 / 45
9 Cauchy Integral Theorem Definition f(a) = 1 f(z)(zi A) 1 dz, 2πi Γ where f is analytic on and inside a closed contour Γ that encloses λ(a). MIMS Nick Higham Functions of a Matrix 8 / 45
10 Equivalence of Definitions Theorem The three definitions are equivalent, modulo analyticity assumption for Cauchy. Interpolation: for basic properties. JCF: for solving matrix equations (e.g., X 2 = A, e X = A). For evaluation (normal A). Cauchy: various uses. MIMS Nick Higham Functions of a Matrix 9 / 45
11 Equivalence of Definitions Theorem The three definitions are equivalent, modulo analyticity assumption for Cauchy. Interpolation: for basic properties. JCF: for solving matrix equations (e.g., X 2 = A, e X = A). For evaluation (normal A). Cauchy: various uses. For computation: Use the definitions (with care). Schur decomposition for general f. Methods specific to particular f. MIMS Nick Higham Functions of a Matrix 9 / 45
12 Outline 1 Definition of f(a) 2 Motivation and MATLAB 3 e A and its Frechét derivative 4 A 1/2 : Modified Newton Methods MIMS Nick Higham Functions of a Matrix 10 / 45
13 Toolbox of Matrix Functions Want to have techniques for evaluating interesting f at matrix args as well as scalar args. Example: has solution d 2 y dt 2 + Ay = 0, y(0) = y 0, y (0) = y 0 y(t) = cos( At)y 0 + ( A ) 1 sin( At)y 0, where A is any square root of A. MATLAB hasexpm,logm,sqrtm,funm. MIMS Nick Higham Functions of a Matrix 11 / 45
14 Classic MATLAB < M A T L A B > Version of 01/10/84 HELP is available <> help Type HELP followed by INTRO (To get started) NEWS (recent revisions) ABS ANS ATAN BASE CHAR CHOL CHOP CLEA COND CONJ COS DET DIAG DIAR DISP EDIT EIG ELSE END EPS EXEC EXIT EXP EYE FILE FLOP FLPS FOR FUN HESS HILB IF IMAG INV KRON LINE LOAD LOG LONG LU MACR MAGI NORM ONES ORTH PINV PLOT POLY PRIN PROD QR RAND RANK RCON RAT REAL RETU RREF ROOT ROUN SAVE SCHU SHOR SEMI SIN SIZE SQRT STOP SUM SVD TRIL TRIU USER WHAT WHIL WHO WHY < > ( ) =., ; \ / + - * : MIMS Nick Higham Functions of a Matrix 12 / 45
15 Classic MATLAB <> help fun FUN For matrix arguments X, the functions SIN, COS, ATAN, SQRT, LOG, EXP and X**p are computed using eigenvalues D and eigenvectors V. If <V,D> = EIG(X) then f(x) = V*f(D)/V. This method may give inaccurate results if V is badly conditioned. Some idea of the accuracy can be obtained by comparing X**1 with X. For vector arguments, the function is applied to each component. The availability of [FUN] in early versions of MATLAB quite possibly contributed to the system s technical and commercial success. Cleve Moler (2003) MIMS Nick Higham Functions of a Matrix 13 / 45
16 Outline 1 Definition of f(a) 2 Motivation and MATLAB 3 e A and its Frechét derivative 4 A 1/2 : Modified Newton Methods MIMS Nick Higham Functions of a Matrix 14 / 45
17 Matrix Exponential Large literature. Early exposition in Frazer, Duncan & Collar. Elementary Matrices and Some Applications to Dynamics and Differential Equations. CUP, Moler & Van Loan. Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later, SIAM Rev., 45 (2003). Over 500 citations on ISI Citation Index. MIMS Nick Higham Functions of a Matrix 15 / 45
18 Scaling and Squaring Method B A/2 s so B 1 r m (B) = [m/m] Padé approximant to e B X = r m (B) 2s e A Used byexpm in MATLAB. Originates with Lawson (1967). Ward (1977): algorithm, with rounding error analysis and a posteriori error bound. Moler & Van Loan (1978): give backward error analysis covering truncation error in Padé approximants, allowing choice of s and m. MIMS Nick Higham Functions of a Matrix 16 / 45
19 Padé Approximations r m to e x r m (x) = p m (x)/q m (x) known explicitly: p m (x) = m j=0 (2m j)!m! x j (2m)!(m j)! j! and q m (x) = p m ( x). Error satisfies e x r m (x) = ( 1) m (m!) 2 (2m)!(2m + 1)! x2m+1 +O(x 2m+2 ). MIMS Nick Higham Functions of a Matrix 17 / 45
20 Analysis Let e A r m (A) = I + G = e H and assume G < 1. Then H = log(i + G) G j /j = log(1 G ). j=1 Hence r m (A) = e A e H = e A+H. Rewrite as r m (A/2 s ) 2s = e A+E, where E = 2 s H satisfies E 2 s log(1 G ). MIMS Nick Higham Functions of a Matrix 18 / 45
21 Result Theorem Let e 2 sa r m (2 s A) = I + G, where G < 1. Then the Padé approximant r m satisfies r m (2 s A) 2s = e A+E, where E A log(1 G ). 2 s A Need to bound G given m and 2 s A. Then E / A bounded in terms of m, s and A. Now select best (s, m) pair for given A. MIMS Nick Higham Functions of a Matrix 19 / 45
22 Bounding G ρ(x) := e x r m (x) 1 = i=2m+1 c i x i converges absolutely for x < min{ t : q m (t) = 0 } =: ν m. Hence, with θ := 2 s A < ν m, G = ρ(2 s A) i=2m+1 Thus E / A log(1 f(θ))/θ). c i θ i =: f(θ). ( ) If only A known, ( ) is optimal bound on G. Moler & Van Loan (1978) bound less sharp; Dieci & Papini (2000) bound a different error. MIMS Nick Higham Functions of a Matrix 20 / 45
23 Summary of Computations Solve bound = unit roundoff by summing 150 terms of series in 250 digit arithmetic. Work out cost of evaluating r m (B) and of the squaring phase. Minimize cost subject to retaining numerical stability in evaluation of r m (B). Precision m θ m IEEE single IEEE double IEEE quad MIMS Nick Higham Functions of a Matrix 21 / 45
24 Scaling and Squaring Algorithm for e A Algorithm 1 (H, 2005; MATLAB 7.2, Mathematica 5.1) 1 for m = [ ] 2 if A 1 θ m 3 X = r m (A). 4 quit 5 end 6 end 7 A A/2 s with s 0 minimal s.t. A/2 s 1 θ 13 = A 2 = A 2, A 4 = A 2 2, A 6 = A 2 A 4 U = A [ A 6 (b 13 A 6 + b 11 A 4 + b 9 A 2 ) + b 7 A 6 + b 5 A 4 + b 3 A 2 + b 1 I ] 10 V = A 6 (b 12 A 6 + b 10 A 4 + b 8 A 2 ) + b 6 A 6 + b 4 A 4 + b 2 A 2 + b 0 I 11 Solve ( U + V)r 13 = U + V for r X = r s 13 by repeated squaring. MIMS Nick Higham Functions of a Matrix 22 / 45
25 Comparison with Existing Algorithms Method m max 2 s A Alg Ward (1977) [θ 8 = 1.5] Old MATLABexpm [θ 6 = 0.54] Sidje (1998) A 1 > 1: Alg 1 requires 1 2 fewer mat mults than Ward, 2 3 fewer thanexpm. A 1 (2, 2.1): Alg 1 Ward expm Sidje mults MIMS Nick Higham Functions of a Matrix 23 / 45
26 Squaring Phase The bound A 2 fl(a 2 ) γ n A 2, γ n = nu 1 nu. shows the dangers in matrix squaring. Open question: are errors in squaring phase consistent with conditioning of the problem? Our choice of parameters uses 1 5 fewer matrix squarings than previous algorithms. Hence has potential accuracy advantages. MIMS Nick Higham Functions of a Matrix 24 / 45
27 Numerical Experiment 70 test matrices, dimension Evaluated 1-norm relative error X X 1 / X 1. Notation: expm: Alg 1 (MATLAB 7.2). old_expm: MATLAB 7.1. funm: MATLAB 7. padm: Sidje. e A+E e A 2 cond(a) = lim max. ǫ 0 E 2 ǫ A 2 ǫ e A 2 MIMS Nick Higham Functions of a Matrix 25 / 45
28 Different S&S Codes andfunm old_expm 10 6 padm funm 10 8 expm (Alg 1) cond*u MIMS Nick Higham Functions of a Matrix 26 / 45
29 Performance Profiles Dolan & Moré (2002) propose a new type of performance profile. For the given set of solvers and test problems, plot x-axis: α y-axis: probability that solver has error within factor α of smallest error over all solvers on the test set. MIMS Nick Higham Functions of a Matrix 27 / 45
30 Performance Profile expm (Alg 1) padm 0.5 p 0.4 funm old_expm α MIMS Nick Higham Functions of a Matrix 28 / 45
31 Frechét Derivative Fréchet derivative of f : C n n C n n at X C n n A linear mapping L : C n n C n n s.t. for all E C n n f(x + E) f(x) L(X, E) = o( E ). Example For f(x) = X 2 we have so L(X, E) = XE + EX. f(x + E) f(x) = XE + EX + E 2, MIMS Nick Higham Functions of a Matrix 29 / 45
32 Frechét Derivative of e A L(A, E) = 1 L(A) := max L(A, E) / E. 0 e A(1 s) Ee As ds. Condition no. of the exponential: κ exp (A) = L(A) A. e A L(A) L(A, I) = e A κ exp (A) A. Theorem If A C n n is normal then in the 2-norm, κ exp (A) = A 2. If A R n n is a nonnegative scalar multiple of a stochastic matrix then in the -norm, κ exp (A) = A. MIMS Nick Higham Functions of a Matrix 30 / 45
33 Quadrature Repeated trap rule: 1 f(t) dt 1 ( 1f 0 m ( f 1 + f f m 1 + 1f 2 m), f i := f(i/m) gives L(A, E) 1 1 m 2 ea E + m 1 i=1 ea(1 i/m) Ee Ai/m + ). 1 2 EeA Requires e A/m, e 2A/m,..., e (m 1)A/m, e A. Lemma Consider R 1 (A, E) = p i=1 w i e A(1 t i) Ee At i and denote its m-times repeated form by R m (A, E). If Q s = R 1 (2 s A, 2 s E), Q i 1 = e 2 ia Q i + Q i e 2 ia, i = s : 1 : 1, then Q 0 = R 2 s(a, E). MIMS Nick Higham Functions of a Matrix 31 / 45
34 Scaling and Squaring Recall L x 2(A, E) = AE + EA. Applying chain rule to e A = (e A/2 ) 2 gives ( L exp (A, E) = L x 2 e A/2, L exp (A/2, E/2) ) = e A/2 L exp (A/2, E/2) + L exp (A/2, E/2)e A/2. Recurrence for L 0 = L exp (A, E): L s = L exp (2 s A, 2 s E), L i 1 = e 2 ia L i + L i e 2 ia, i = s : 1 : 1. MIMS Nick Higham Functions of a Matrix 32 / 45
35 Quadrature Algorithm Algorithm (Kenney & Laub, 1998; Mathias, 1993) Approximates L = L exp (A, E) via repeated Simpson rule; order of magnitude estimate. 1 B = A/2 s with s 0 minimal s.t. A/2 s 1 1/2 2 X = e B 3 X = e B/2 4 Q s = 2 s (XE + 4 XE X + EX)/6 5 for i = s: 1: 1 6 if i < s, X = e 2 ia, end 7 Q i 1 = XQ i + Q i X 8 end 9 L = Q 0 MIMS Nick Higham Functions of a Matrix 33 / 45
36 Kronecker Formula L(A, E) = 1 0 e A(1 s) Ee As ds. MIMS Nick Higham Functions of a Matrix 34 / 45
37 Kronecker Formula L(A, E) = 1 A B = (a ij B) C n2 n 2. A B = A I n + I m B. e A B = e A e B. vec(axb) = (B T A)vec(X). 0 e A(1 s) Ee As ds. MIMS Nick Higham Functions of a Matrix 34 / 45
38 Kronecker Formula L(A, E) = 1 A B = (a ij B) C n2 n 2. A B = A I n + I m B. e A B = e A e B. vec(axb) = (B T A)vec(X). Theorem vec(l(a, E)) = K(A)vec(E), where 0 e A(1 s) Ee As ds. K(A) = 1 2 (eat e A )τ ( 1 2 [AT ( A)] ) C n2 n 2, τ(x) = tanh(x)/x and 1 2 AT ( A) < π/2 assumed. MIMS Nick Higham Functions of a Matrix 34 / 45
39 Approximating the Frechét Derivative [ A T ( A) θi ] vec(e) = [ A T I I A θi ] vec(e) = vec(ea AE θi) = vec(e(a θ/2 I) (A + θ/2 I)E). Obtain Padé approximant τ(x) = tanh(x)/x r m (x) = m (x/β i 1) 1 (x/α i 1) by truncating τ(x) = i=1 1 x 2 /(1 3) x 2 /(3 5) 1 + x 2 /((2k 1) (2k + 1)) 1 +. MIMS Nick Higham Functions of a Matrix 35 / 45
40 Algorithm (Fréchet derivative; Kenney & Laub, 1998) Evaluates L = L exp (A, E) using scaling and squaring and [8/8] Padé approximant to τ. 1 B = A/2 s with s 0 minimal s.t. A/2 s G 0 = 2 s E 3 for i = 1: 8 4 Solve (I + B/β i ) G i + G i (I B/β i ) = (I + B/α i )G i 1 + G i 1 (I B/α i ). 5 end 6 X = e B 7 L s = (G m X + XG m )/2 8 for i = s: 1: 1 9 if i < s, X = e 2 ia, end 10 L i 1 = XL i + L i X 11 end 12 L = L 0
41 Outline 1 Definition of f(a) 2 Motivation and MATLAB 3 e A and its Frechét derivative 4 A 1/2 : Modified Newton Methods MIMS Nick Higham Functions of a Matrix 37 / 45
42 Matrix Square Root X is a square root of A C n n X 2 = A. Number of square roots may be zero, finite or infinite. Definition For A with no eigenvalues on R = {x R : x 0} the principal square root A 1/2 is unique square root X with spectrum in open right half-plane. MIMS Nick Higham Functions of a Matrix 38 / 45
43 Newton s Method for Square Root Apply Newton to F(X) = X 2 A = 0: X 0 given, } Solve X k E k + E k X k = A Xk 2 k = 0, 1, 2,... X k+1 = X k + E k Modified Newton iteration: freeze Fréchet derivative at X 0 : } Solve X 0 E k + E k X 0 = A Xk 2 k = 0, 1, 2,..., X k+1 = X k + E k X 0 diagonal cheap to solve for E k. MIMS Nick Higham Functions of a Matrix 39 / 45
44 Visser Iteration Set X 0 = (2α) 1 I in modified Newton: Visser iteration (1937) X k+1 = X k + α(a X 2 k ), X 0 = (2α) 1 I. Stationary iteration. Richardson iteration. Linear convergence. Choice of α? MIMS Nick Higham Functions of a Matrix 40 / 45
45 Visser History X k+1 = X k + α(a X 2 k ), X 0 = (2α) 1 I. Used with α = 1/2 by Visser (1937) to show positive operator on Hilbert space has a positive square root. Likewise in functional analysis texts, e.g. Riesz & Sz.-Nagy (1956). Enables proof of existence of A 1/2 without using spectral theorem. Used computationally by Liebl (1965), Babuška, Práger & Vitásek (1966), Späth (1966), Duke (1969), Elsner (1970). Elsner proves cgce for A C n n with real, positive ei vals if 0 < α ρ(a) 1/2. MIMS Nick Higham Functions of a Matrix 41 / 45
46 Visser Transformations Let X k = θy k, β = θα and à = θ 2 A. Set β = 1/2: Y k+1 = Y k + β(ã Y 2 k ), Y 0 = 1 2β I. Y k+1 = Y k (à Y 2 k ), Y 0 = I. With à I C and Y k = I P k : Q k = P k /2: P k+1 = 1 2 (C + P2 k), P 0 = 0. Q k+1 = Q 2 k + C 4, Q 0 = 0. MIMS Nick Higham Functions of a Matrix 42 / 45
47 Visser Convergence X k+1 = X k + α(a Xk 2 ), X 0 = (2α) 1 I. Theorem (H, 2006) Let A C n n and α > 0. If Λ(I 4α 2 A) lies in the cardioid D = { 2z z 2 : z C, z < 1 } then A 1/2 exists and X k A 1/2 linearly MIMS Nick Higham Functions of a Matrix 43 / 45
48 Example A R spd with a ii = i 2, a ij = 0.1, i j. Aim for rel residual < nu in IEEE DP arithmetic. Pulay iteration D = diag(a): θ = 0.191, 9 iters. Visser iteration α = (hand optimized), 245 iters MIMS Nick Higham Functions of a Matrix 44 / 45
49 In Conclusion Many applications of f(a), e.g. control theory, computer graphics, theoretical physics. Need better understanding of conditioning of f(a). Can we exploit structure, e.g. A matrix automorphism group or Jordan or Lie algebra? Krylov methods needed for large, sparse A. How to use Cauchy integral computationally (H & Trefethen)? MIMS Nick Higham Functions of a Matrix 45 / 45
50 References I D. A. Bini, N. J. Higham, and B. Meini. Algorithms for the matrix pth root. Numer. Algorithms, 39(4): , C.-H. Guo and N. J. Higham. A Schur Newton method for the matrix pth root and its inverse. SIAM J. Matrix Anal. Appl., 28(3): , N. J. Higham. Functions of a Matrix: Theory and Computation. Book in preparation. MIMS Nick Higham Functions of a Matrix 44 / 45
51 References II N. J. Higham. The scaling and squaring method for the matrix exponential revisited. SIAM J. Matrix Anal. Appl., 26(4): , C. S. Kenney and A. J. Laub. A Schur Fréchet algorithm for computing the logarithm and exponential of a matrix. SIAM J. Matrix Anal. Appl., 19(3): , R. Mathias. Evaluating the Frechet derivative of the matrix exponential. Numer. Math., 63(2): , MIMS Nick Higham Functions of a Matrix 45 / 45
Lecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
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 informationResearch CharlieMatters
Research CharlieMatters Van Loan and the Matrix Exponential February 25, 2009 Nick Nick Higham Director Schoolof ofresearch Mathematics The School University of Mathematics of Manchester http://www.maths.manchester.ac.uk/~higham
More informationThe Scaling and Squaring Method for the Matrix Exponential Revisited. Nicholas J. Higham
The Scaling and Squaring Method for the Matrix Exponential Revisited Nicholas J. Higham 2005 MIMS EPrint: 2006.394 Manchester Institute for Mathematical Sciences School of Mathematics The University of
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 informationImproved Inverse Scaling and Squaring Algorithms for the Matrix Logarithm. Al-Mohy, Awad H. and Higham, Nicholas J. MIMS EPrint: 2011.
Improved Inverse Scaling and Squaring Algorithms for the Matrix Logarithm Al-Mohy, Awad H. and Higham, Nicholas J. 2011 MIMS EPrint: 2011.83 Manchester Institute for Mathematical Sciences School of Mathematics
More informationFunctions of Matrices. Nicholas J. Higham. November MIMS EPrint: Manchester Institute for Mathematical Sciences School of Mathematics
Functions of Matrices Nicholas J. Higham November 2005 MIMS EPrint: 2005.21 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available from: And
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationComputing the pth Roots of a Matrix. with Repeated Eigenvalues
Applied Mathematical Sciences, Vol. 5, 2011, no. 53, 2645-2661 Computing the pth Roots of a Matrix with Repeated Eigenvalues Amir Sadeghi 1, Ahmad Izani Md. Ismail and Azhana Ahmad School of Mathematical
More informationFunctions of Matrices and
Functions of Matrices and Nearest Research Correlation MattersMatrices February 25, 2009 Nick Higham School Nick Higham of Mathematics The Director University of Research of Manchester Nicholas.J.Higham@manchester.ac.uk
More informationTesting matrix function algorithms using identities. Deadman, Edvin and Higham, Nicholas J. MIMS EPrint:
Testing matrix function algorithms using identities Deadman, Edvin and Higham, Nicholas J. 2014 MIMS EPrint: 2014.13 Manchester Institute for Mathematical Sciences School of Mathematics The University
More informationKrylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms
Krylov Subspace Methods for the Evaluation of Matrix Functions. Applications and Algorithms 2. First Results and Algorithms Michael Eiermann Institut für Numerische Mathematik und Optimierung Technische
More informationMatrix Functions: A Short Course. Higham, Nicholas J. and Lijing, Lin. MIMS EPrint:
Matrix Functions: A Short Course Higham, Nicholas J. and Lijing, Lin 2013 MIMS EPrint: 2013.73 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports
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 informationChapter 16 Numerical Linear Algebra
16.1 Sets of Linear Equations Chapter 16 Numerical Linear Algebra MATLAB was developed to handle problems involving matrices and vectors in an efficient way. One of the most basic problems of this type
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 informationThe Complex Step Approximation to the Fréchet Derivative of a Matrix Function. Al-Mohy, Awad H. and Higham, Nicholas J. MIMS EPrint: 2009.
The Complex Step Approximation to the Fréchet Derivative of a Matrix Function Al-Mohy, Awad H. and Higham, Nicholas J. 2009 MIMS EPrint: 2009.31 Manchester Institute for Mathematical Sciences School of
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
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 informationThe Leja method revisited: backward error analysis for the matrix exponential
The Leja method revisited: backward error analysis for the matrix exponential M. Caliari 1, P. Kandolf 2, A. Ostermann 2, and S. Rainer 2 1 Dipartimento di Informatica, Università di Verona, Italy 2 Institut
More informationTaylor s Theorem for Matrix Functions with Applications to Condition Number Estimation. Deadman, Edvin and Relton, Samuel. MIMS EPrint: 2015.
Taylor s Theorem for Matrix Functions with Applications to Condition Number Estimation Deadman, Edvin and Relton, Samuel 215 MIMS EPrint: 215.27 Manchester Institute for Mathematical Sciences School of
More informationComputing Real Logarithm of a Real Matrix
International Journal of Algebra, Vol 2, 2008, no 3, 131-142 Computing Real Logarithm of a Real Matrix Nagwa Sherif and Ehab Morsy 1 Department of Mathematics, Faculty of Science Suez Canal University,
More informationA more accurate Briggs method for the logarithm
Numer Algor (2012) 59:393 402 DOI 10.1007/s11075-011-9496-z ORIGINAL PAPER A more accurate Briggs method for the logarithm Awad H. Al-Mohy Received: 25 May 2011 / Accepted: 15 August 2011 / Published online:
More informationMath 307 Learning Goals. March 23, 2010
Math 307 Learning Goals March 23, 2010 Course Description The course presents core concepts of linear algebra by focusing on applications in Science and Engineering. Examples of applications from recent
More informationComputing the matrix cosine
Numerical Algorithms 34: 13 26, 2003. 2003 Kluwer Academic Publishers. Printed in the Netherlands. Computing the matrix cosine Nicholas J. Higham and Matthew I. Smith Department of Mathematics, University
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 informationSolving Burnup Equations in Serpent:
Solving Burnup Equations in Serpent: Matrix Exponential Method CRAM Maria Pusa September 15th, 2011 Outline Burnup equations and matrix exponential solution Characteristics of burnup matrices Established
More informationDense LU factorization and its error analysis
Dense LU factorization and its error analysis Laura Grigori INRIA and LJLL, UPMC February 2016 Plan Basis of floating point arithmetic and stability analysis Notation, results, proofs taken from [N.J.Higham,
More informationA Cholesky LR algorithm for the positive definite symmetric diagonal-plus-semiseparable eigenproblem
A Cholesky LR algorithm for the positive definite symmetric diagonal-plus-semiseparable eigenproblem Bor Plestenjak Department of Mathematics University of Ljubljana Slovenia Ellen Van Camp and Marc Van
More information3. Array and Matrix Operations
3. Array and Matrix Operations Almost anything you learned about in your linear algebra classmatlab has a command to do. Here is a brief summary of the most useful ones for physics. In MATLAB matrices
More informationComputing Matrix Functions. Higham, Nicholas J. and Al-Mohy, Awad H. MIMS EPrint:
Computing Matrix Functions Higham, Nicholas J. and Al-Mohy, Awad H. 2010 MIMS EPrint: 2010.18 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationMATH 5524 MATRIX THEORY Problem Set 4
MATH 5524 MATRIX THEORY Problem Set 4 Posted Tuesday 28 March 217. Due Tuesday 4 April 217. [Corrected 3 April 217.] [Late work is due on Wednesday 5 April.] Complete any four problems, 25 points each.
More informationNumerical methods for matrix functions
Numerical methods for matrix functions SF2524 - Matrix Computations for Large-scale Systems Lecture 13 Numerical methods for matrix functions 1 / 26 Reading material Lecture notes online Numerical methods
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2018 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationPositive Denite Matrix. Ya Yan Lu 1. Department of Mathematics. City University of Hong Kong. Kowloon, Hong Kong. Abstract
Computing the Logarithm of a Symmetric Positive Denite Matrix Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A numerical method for computing the logarithm
More informationN. J. Higham. Numerical Analysis Report No January Manchester Centre for Computational Mathematics Numerical Analysis Reports
UMIST ISSN 1360-1725 A New sqrtm for Matlab N. J. Higham Numerical Analysis Report No. 336 January 1999 Manchester Centre for Computational Mathematics Numerical Analysis Reports DEPARTMENTS OF MATHEMATICS
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
More informationJim Lambers MAT 610 Summer Session Lecture 2 Notes
Jim Lambers MAT 610 Summer Session 2009-10 Lecture 2 Notes These notes correspond to Sections 2.2-2.4 in the text. Vector Norms Given vectors x and y of length one, which are simply scalars x and y, the
More informationMath 307 Learning Goals
Math 307 Learning Goals May 14, 2018 Chapter 1 Linear Equations 1.1 Solving Linear Equations Write a system of linear equations using matrix notation. Use Gaussian elimination to bring a system of linear
More informationThe Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix
The Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix Chun-Yueh Chiang Center for General Education, National Formosa University, Huwei 632, Taiwan. Matthew M. Lin 2, Department of
More informationOn-the-fly backward error estimate for matrix exponential approximation by Taylor algorithm
On-the-fly backward error estimate for matrix exponential approximation by Taylor algorithm M. Caliari a,, F. Zivcovich b a Department of Computer Science, University of Verona, Italy b Department of Mathematics,
More informationCayley-Hamilton Theorem
Cayley-Hamilton Theorem Massoud Malek In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n Let A be an n n matrix Although det (λ I n A
More informationSome Formulas for the Principal Matrix pth Root
Int. J. Contemp. Math. Sciences Vol. 9 014 no. 3 141-15 HIKARI Ltd www.m-hiari.com http://dx.doi.org/10.1988/ijcms.014.4110 Some Formulas for the Principal Matrix pth Root R. Ben Taher Y. El Khatabi and
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationRehabilitating Research Correlations, Matters Avoiding Inversion, and Extracting Roots
Rehabilitating Research Correlations, Matters Avoiding Inversion, and Extracting Roots February 25, 2009 Nick Nick Higham Higham Director School of of Research Mathematics The University of Manchester
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationEigenvalues, Eigenvectors and the Jordan Form
EE/ME 701: Advanced Linear Systems Eigenvalues, Eigenvectors and the Jordan Form Contents 1 Introduction 3 1.1 Review of basic facts about eigenvectors and eigenvalues..... 3 1.1.1 Looking at eigenvalues
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 informationSome Useful Results in the Theory of Matrix Functions
Some Useful Results in the Theory of Matrix Functions Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Matrix Analysis and Applications,
More informationSolving Systems of Fuzzy Differential Equation
International Mathematical Forum, Vol. 6, 2011, no. 42, 2087-2100 Solving Systems of Fuzzy Differential Equation Amir Sadeghi 1, Ahmad Izani Md. Ismail and Ali F. Jameel School of Mathematical Sciences,
More informationAutomatica, 33(9): , September 1997.
A Parallel Algorithm for Principal nth Roots of Matrices C. K. Koc and M. _ Inceoglu Abstract An iterative algorithm for computing the principal nth root of a positive denite matrix is presented. The algorithm
More informationThe Leja method revisited: backward error analysis for the matrix exponential
The Leja method revisited: backward error analysis for the matrix exponential M. Caliari 1, P. Kandolf 2, A. Ostermann 2, and S. Rainer 2 1 Dipartimento di Informatica, Università di Verona, Italy 2 Institut
More informationMatrix functions that preserve the strong Perron- Frobenius property
Electronic Journal of Linear Algebra Volume 30 Volume 30 (2015) Article 18 2015 Matrix functions that preserve the strong Perron- Frobenius property Pietro Paparella University of Washington, pietrop@uw.edu
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 informationANSWERS (5 points) Let A be a 2 2 matrix such that A =. Compute A. 2
MATH 7- Final Exam Sample Problems Spring 7 ANSWERS ) ) ). 5 points) Let A be a matrix such that A =. Compute A. ) A = A ) = ) = ). 5 points) State ) the definition of norm, ) the Cauchy-Schwartz inequality
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 informationMATRIX ARITHMETIC-GEOMETRIC MEAN AND THE COMPUTATION OF THE LOGARITHM
SIAM J. MATRIX ANAL. APPL. Vol. 37, No., pp. 719 743 c 016 Society for Industrial and Applied Mathematics MATRIX ARITHMETIC-GEOMETRIC MEAN AND THE COMPUTATION OF THE LOGARITHM JOÃO R. CARDOSO AND RUI RALHA
More informationMATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N.
MATH 3795 Lecture 13. Numerical Solution of Nonlinear Equations in R N. Dmitriy Leykekhman Fall 2008 Goals Learn about different methods for the solution of F (x) = 0, their advantages and disadvantages.
More informationNovember 18, 2013 ANALYTIC FUNCTIONAL CALCULUS
November 8, 203 ANALYTIC FUNCTIONAL CALCULUS RODICA D. COSTIN Contents. The spectral projection theorem. Functional calculus 2.. The spectral projection theorem for self-adjoint matrices 2.2. The spectral
More 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 informationA Schur Parlett Algorithm for Computing Matrix Functions. Davies, Philip I. and Higham, Nicholas J. MIMS EPrint:
A Schur Parlett Algorithm for Computing Matrix Functions Davies, Philip I. and Higham, Nicholas J. 2003 MIMS EPrint: 2006.13 Manchester Institute for Mathematical Sciences School of Mathematics The University
More informationUsing MATLAB. Linear Algebra
Using MATLAB in Linear Algebra Edward Neuman Department of Mathematics Southern Illinois University at Carbondale One of the nice features of MATLAB is its ease of computations with vectors and matrices.
More informationLinear Algebraic Equations
Linear Algebraic Equations 1 Fundamentals Consider the set of linear algebraic equations n a ij x i b i represented by Ax b j with [A b ] [A b] and (1a) r(a) rank of A (1b) Then Axb has a solution iff
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 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 informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 12: Gauss for Linear Systems Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationEigenvalue Problems and Singular Value Decomposition
Eigenvalue Problems and Singular Value Decomposition Sanzheng Qiao Department of Computing and Software McMaster University August, 2012 Outline 1 Eigenvalue Problems 2 Singular Value Decomposition 3 Software
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 informationStable iterations for the matrix square root
Numerical Algorithms 5 (997) 7 4 7 Stable iterations for the matrix square root Nicholas J. Higham Department of Mathematics, University of Manchester, Manchester M3 9PL, England E-mail: higham@ma.man.ac.u
More informationLAPACK-Style Codes for Pivoted Cholesky and QR Updating. Hammarling, Sven and Higham, Nicholas J. and Lucas, Craig. MIMS EPrint: 2006.
LAPACK-Style Codes for Pivoted Cholesky and QR Updating Hammarling, Sven and Higham, Nicholas J. and Lucas, Craig 2007 MIMS EPrint: 2006.385 Manchester Institute for Mathematical Sciences School of Mathematics
More informationPerturbations of functions of diagonalizable matrices
Electronic Journal of Linear Algebra Volume 20 Volume 20 (2010) Article 22 2010 Perturbations of functions of diagonalizable matrices Michael I. Gil gilmi@bezeqint.net Follow this and additional works
More informationScaled and squared subdiagonal Padé approximation for the matrix exponential. Güttel, Stefan and Yuji, Nakatsukasa. MIMS EPrint: 2015.
Scaled and squared subdiagonal Padé approximation for the matrix exponential Güttel, Stefan and Yuji, Nakatsukasa 215 MIMS EPrint: 215.46 Manchester Institute for Mathematical Sciences School of Mathematics
More informationUniversity of Maryland Department of Computer Science TR-5009 University of Maryland Institute for Advanced Computer Studies TR April 2012
University of Maryland Department of Computer Science TR-5009 University of Maryland Institute for Advanced Computer Studies TR-202-07 April 202 LYAPUNOV INVERSE ITERATION FOR COMPUTING A FEW RIGHTMOST
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationChapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors
Chapter 7 Canonical Forms 7.1 Eigenvalues and Eigenvectors Definition 7.1.1. Let V be a vector space over the field F and let T be a linear operator on V. An eigenvalue of T is a scalar λ F such that there
More informationLAPACK-Style Codes for Pivoted Cholesky and QR Updating
LAPACK-Style Codes for Pivoted Cholesky and QR Updating Sven Hammarling 1, Nicholas J. Higham 2, and Craig Lucas 3 1 NAG Ltd.,Wilkinson House, Jordan Hill Road, Oxford, OX2 8DR, England, sven@nag.co.uk,
More informationBasic Elements of Linear Algebra
A Basic Review of Linear Algebra Nick West nickwest@stanfordedu September 16, 2010 Part I Basic Elements of Linear Algebra Although the subject of linear algebra is much broader than just vectors and matrices,
More informationOn the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination
On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination J.M. Peña 1 Introduction Gaussian elimination (GE) with a given pivoting strategy, for nonsingular matrices
More informationExponential integrators and functions of the matrix exponential
Exponential integrators and functions of the matrix exponential Paul Matthews, Stephen Cox, Hala Ashi and Linda Cummings School of Mathematical Sciences, University of Nottingham, UK Introduction to exponential
More informationComputational Economics and Finance
Computational Economics and Finance Part II: Linear Equations Spring 2016 Outline Back Substitution, LU and other decomposi- Direct methods: tions Error analysis and condition numbers Iterative methods:
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 informationComputing the Action of the Matrix Exponential, with an Application to Exponential Integrators. Al-Mohy, Awad H. and Higham, Nicholas J.
Computing the Action of the Matrix Exponential, with an Application to Exponential Integrators Al-Mohy, Awad H. and Higham, Nicholas J. 2010 MIMS EPrint: 2010.30 Manchester Institute for Mathematical Sciences
More informationSolving linear equations with Gaussian Elimination (I)
Term Projects Solving linear equations with Gaussian Elimination The QR Algorithm for Symmetric Eigenvalue Problem The QR Algorithm for The SVD Quasi-Newton Methods Solving linear equations with Gaussian
More informationHigher-order Chebyshev Rational Approximation Method (CRAM) and application to Burnup Calculations
Higher-order Chebyshev Rational Approximation Method (CRAM) and application to Burnup Calculations Maria Pusa September 18th 2014 1 Outline Burnup equations and matrix exponential solution Characteristics
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More informationSquare root of a symmetric matrix
Square root of a symmetric matrix Mario Berljafa Stefan Güttel June 2015 Contents 1 Introduction 1 2 Running rkfit 1 3 Evaluating the rational approximant 2 4 Some different choices for the initial poles
More information2 DIECI AND PAPINI how to choose k; on one hand, we need k large enough because Pade approximations are accurate only if A=2 k is suciently small (see
PADE APPROXIMATION FOR THE EXPONENTIAL OF A BLOCK TRIANGULAR MATRIX. LUCA DIECI AND ALESSANDRA PAPINI Abstract. In this work we obtain improved error bounds for Pade approximations to e A when A is block
More informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More information13-2 Text: 28-30; AB: 1.3.3, 3.2.3, 3.4.2, 3.5, 3.6.2; GvL Eigen2
The QR algorithm The most common method for solving small (dense) eigenvalue problems. The basic algorithm: QR without shifts 1. Until Convergence Do: 2. Compute the QR factorization A = QR 3. Set A :=
More informationG1110 & 852G1 Numerical Linear Algebra
The University of Sussex Department of Mathematics G & 85G Numerical Linear Algebra Lecture Notes Autumn Term Kerstin Hesse (w aw S w a w w (w aw H(wa = (w aw + w Figure : Geometric explanation of the
More informationOn a quadratic matrix equation associated with an M-matrix
Article Submitted to IMA Journal of Numerical Analysis On a quadratic matrix equation associated with an M-matrix Chun-Hua Guo Department of Mathematics and Statistics, University of Regina, Regina, SK
More informationLinear Algebra. Matrices Operations. Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0.
Matrices Operations Linear Algebra Consider, for example, a system of equations such as x + 2y z + 4w = 0, 3x 4y + 2z 6w = 0, x 3y 2z + w = 0 The rectangular array 1 2 1 4 3 4 2 6 1 3 2 1 in which the
More informationFinal Examination. CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University
Final Examination CS 205A: Mathematical Methods for Robotics, Vision, and Graphics (Fall 2013), Stanford University The exam runs for 3 hours. The exam contains eight problems. You must complete the first
More informationCompound matrices and some classical inequalities
Compound matrices and some classical inequalities Tin-Yau Tam Mathematics & Statistics Auburn University Dec. 3, 04 We discuss some elegant proofs of several classical inequalities of matrices by using
More information