How Rutishauser may have found the qd and LR algorithms, t
|
|
- Warren Carroll
- 6 years ago
- Views:
Transcription
1 How Rutishauser may have found the qd and LR algorithms, the fore-runners of QR Martin H Gutnecht and Beresford N Parlett Seminar for Applied Mathematics ETH Zurich Department of Mathematics University of California Bereley 3rd Biennial Conference on Numerical Analysis, U Strathclyde, Glasgow, Scotland 3 6 June 9 Martin H Gutnecht and Beresford N Parlett
2 Rutishauser s qd algorithm: early papers Rutishauser (954a, ZAMP): Der Quotienten Differenzen Algorithmus Rutishauser (954b, ZAMP): Anwendungen des Quotienten Differenzen Algorithmus Rutishauser (954c, ArchMath): Ein infinitesimales Analogon zum Quotienten Differenzen Algorithmus Rutishauser (955a, ZAMP): Bestimmung der Eigenwerte und Eigenvetoren einer Matrix mit Hilfe des Quotienten Differenzen Algorithmus Rutihauser (957a, Mitt IAM, ETH): Der Quotienten Differenzen Algorithmus Henrici (958, NBS boo): The Quotient-Difference Algorithm Martin H Gutnecht and Beresford N Parlett
3 Eduard Stiefel s suggestion ( 953) Stiefel s suggestion to Rutishauser: Given A, x, y, use the Schwarz constants (= moments = Marov parameters) to find the eigenvalues of A s : y T A x ( =,,,) () Daniel Bernoulli (73), J König (884): s ν+ s ν λ as ν if λ > λ λ Note: It turned out that for the other eigenvalues, Stiefel s proposal was a bad idea, since the dependence of the EVals from the moments is highly ill-conditioned (Gautschi (968)) Martin H Gutnecht and Beresford N Parlett
4 Moments and their generating function Given: N N matrix A and x, y R N, let f(z) : y,(zi A) x = y, z (I z A) x () f is a rational function of type (N, N), so f( ) = The poles of f are eigenvalues of A f can be expanded into a power series in z : where f(z) = = s z + = s z + s z + s + (3) z3 s = y T A x So, f is the generating function of the moments Martin H Gutnecht and Beresford N Parlett
5 Alternative formulations of the problem Clearly, there are several equivalent problems: Find eigenvalues of A Find poles of generating (rational) function f Find zeros of the denominator polynomial of f (Bernoulli) In theory, the problem had been solved before by Hadamard (89) (his PhD thesis!), de Montessus de Ballore (9/95), Aiten (96/93) But none of them had an efficient algorithm Rutishauser cites Hadamard and Aiten, but never de Montessus de Ballore, who proved the convergence of Padé approximants with fixed denominator degree Martin H Gutnecht and Beresford N Parlett
6 Hadamard s theorem (89) Given the power series of f in z of (3), let H (ν) :, and define the Hanel determinants s ν s ν+ s ν+ H (ν) = s ν+ s ν+ s ν+ s ν+ s ν+ s ν+ ( =,,; ν =,, ) THEOREM [Hadamard (89)] If λ + < Λ < λ, then, as ν, ( ) H (ν) Λ ν ] = const (λ λ ) [ ν + O λ For a simpler proof see Henrici (958) or Henrici (974) Martin H Gutnecht and Beresford N Parlett
7 Hadamard s theorem (89) (cont d) COROLLARY If f has N simple poles, then H (ν) ( =,,N) for large enough ν, and H (ν) N+ = ( ν) If λ > λ + then H (ν+) H (ν) λ λ λ as ν (4) 3 If λ > λ > λ + then q (ν) : H (ν+) H (ν) H (ν) H (ν+) λ as ν (5) Martin H Gutnecht and Beresford N Parlett
8 Aiten s scheme (93) (cont d) Computing, for fixed ν, the Hanel determinants H (ν),, H (ν) N (if nonzero) requires the LU decomposition of the matrix H (ν) N Aiten (96, 93) used what is now called Jacobi identity ( theorem of compound determinants ) ( H (ν) ) (ν ) = H H (ν+) + H (ν ) + H (ν+) (6) It had also been nown to Hadamard, but Aiten used it to build up from the left or from the top the table H () H () H () H () H () H () 3 H (3) H () H () 3 H () 4 Martin H Gutnecht and Beresford N Parlett
9 Rutishauser s qd algorithm (QD-Algorithmus) Rutishauser (954a) new Aiten s wor and refers to (4), H (ν+) H (ν) λ λ λ as ν as the ey to computing non-dominant poles But instead of computing the H (ν) table, he headed directly for recurrences for q (ν) : H (ν+) H (ν) H (ν) H (ν+) and : H (ν) + H (ν) H (ν+) H (ν+) (7) In Rutishauser (954a) he derives the formulas needed for q (ν), and then states recursions for general Martin H Gutnecht and Beresford N Parlett
10 Rutishauser s qd algorithm (cont d) qd table (QD Schema): q () q () q () q (3) q (4) e () e () e () e (3) e (4) + q () q () q () q (3) + e () e () e () e (3) e () N e () N e () N q () N q () N e () q () = q () e () q () + e () = e () + q () Martin H Gutnecht and Beresford N Parlett
11 Rutishauser s qd algorithm (cont d) Rhombus rules (called so by Stiefel, 955) of qd algorithm: For building up the table columnwise from left to right: := e (ν+) + q (ν+) q (ν) q (ν) + := q (ν+) e (ν+) ( =,,) (8) For building up the table row-wise, from top to bottom: q (ν+) := q (ν) + e (ν+) e (ν+) := q (ν+) q (ν) + ( =,,) (9) Recursions (9) are the basis of the progressive qd algorithm (the relevant version) Martin H Gutnecht and Beresford N Parlett
12 Rutishauser s qd algorithm (cont d) In Rutishauser (954a) the correctness of the rhombus rules follows later from the connections to continued fractions (probably Stiefel s argument) Originally, Rutishauser derived them probably from Hadamard s Jacobi identity ( H (ν) ) (ν ) = H H (ν+) + H (ν ) + H (ν+) Henrici (958), who was in contact with Rutishauser, pointed out that one rhombus rules (+) can be derived by combining two applications of this formula, the other ( ) just by using the definitions (9) of q (ν) and The details have been wored out in Parlett (996), a TR entitled What Hadamard missed Martin H Gutnecht and Beresford N Parlett
13 From power series to a continued fractions By a standard operation the given power series (3) in z of f, f(z) = = s z + = s z + s z + s z 3 + can be turned into a continued fraction (which typically converges in a much larger region) We may also write f(z) = s z + s z + + s ν z ν + f ν(z) z ν () and expand the remainder f ν (z) of the power series into a continued fractions In each case, two different types of continued fractions can be used So we get two whole series of continued fractions It turns out that their coefficients are related by the rhombus rules Martin H Gutnecht and Beresford N Parlett
14 Continued fractions: J fractions and S fractions f ν (z) : = s ν+ = zν z+ ( ν f(z) = s z + can be expanded both into a Jacobi fraction or J fraction f ν (z) = s ν z q (ν) q(ν) z q (ν) and into a formal Stieltjes fraction or S fraction ) () q(ν) z q (ν) 3 () f ν (z) = s ν z q(ν) e(ν) z q(ν) e(ν) z (3) The J fraction is the so-called even part of the S fraction obtained by merging two successive terms into one Martin H Gutnecht and Beresford N Parlett
15 Continued fractions: J fractions and S fractions The odd part of the S fraction is another formal J fraction, obtained by merging the two differently chosen successive terms into one, f ν(z) = sν z + q (ν) z q (ν) q(ν) z q (ν) By comparing this J fraction with the one for q(ν) 3 z q (ν) 3 3 (4) f ν+ (z) = zf ν (z) s ν, (5) one recovers Rutishauser s rhombus rules of the qd algorithm This is the nicest derivation of the qd algorithm, but not the original one Rutishauser (954a) indicates that it may have been suggested to him by Stiefel Martin H Gutnecht and Beresford N Parlett
16 Continued fractions, Padé approximations, FOPs The partial sums = convergents = approximants of the continued fractions are confluent rational interpolants of f They are Padé approximants (at ) associated with the moments s +ν ( =,,; ν fixed) of the function f ν (z) defined by f(z) = s z + s z + + s ν z ν + f ν(z) z ν For fixed ν, the denominators of the convergents (= Padé approximants) are (formal) orthogonal polynomials p (ν) (z) They can be arranged in a table that he called p table (Analogous to the Padé table) Martin H Gutnecht and Beresford N Parlett
17 Associated polynomials and their p table p table (P Schema): p () p () p () p (3) p (4) p () p () p () p (3) p () In the last column, p () N q () z p() p () p () p () 3 p () 3 p () 3 p () N p () N (z) := z p() (z) q() p() (z) = p() N p () N p () N = is the minimal polynomial Martin H Gutnecht and Beresford N Parlett
18 Continued fractions, Padé approximations, FOPs (cont d) Rutishauser realized that they are also the Lanczos polynomials of the (nonsymmetric) Lanczos algorithm (Lanczos, 95) for A started with the pair (y, A ν x ) Rutishauser never mentions Padé approximants, but he had no need, since they are just the convergents of the J fractions and S fractions For him, actually only the FOPs in the denominator matter Later, , Householder, Gragg, and Stewart stress the connection to Padé approximants in several papers NB: Hadamard s theorem (89) de Montessus de Ballore s theorem (9/95) Martin H Gutnecht and Beresford N Parlett
19 Associated polynomials: recurrences The p table can be built up from the initial column p (ν) by the left-to-right recurrence p (ν) (z) := zp (ν+) (z) q(ν) p (ν) (z) (6) Rutishauser (954a) derived also a top-to-bottom recurrence p (ν+) (z) := p (ν) (z) p (ν+) (z) (7) and the diagonal 3-term recurrence (with :, p (ν) : ) [ ] p (ν) + (z) := z q (ν) + e(ν) p (ν) (z) + q(ν) p (ν) (z) (8) Martin H Gutnecht and Beresford N Parlett
20 Further relations and applications So, in addition to introducing and investigating the qd algorithm Rutishauser (954a) [rec 5 Aug 953] (954b) [rec 8 Sep 953], (955a) [rec 9 July 954)] explained many connections to other topics and gave many applications; eg, in (954a): the connection to continued fractions, the connection to the Lanczos BIO algorithm, the connection to the CG algorithm, computing partial fraction decompositions of rational fcts In (954b): summation of badly converging series, solving algebraic equations = computing zeros of polynomials, quadratic convergence by using shifts / double shifts Martin H Gutnecht and Beresford N Parlett
21 Further relations and applications (cont d) In (955a): computing EVals by combining Lanczos BIO alg and the progressive qd algorithm, computing EVecs (several new algorithms are suggested), EVals and EVecs of infinite matrices Still missing: tridiagonal matrices (except for computing shifts), LU decomposition of these tridiagonal matrices, LR algorithm Martin H Gutnecht and Beresford N Parlett
22 FOPs and tridiagonal matrices Rutishauser new well (see Rutishauser (953) on the Lanczos BIO algorithm) that associated to the 3-term recurrence (8), T (ν) n = [ ] p (ν) + (z) := z q (ν) + e(ν) p (ν) (z) + q(ν) p (ν) (z) (with fixed ν) there is a nested set of tridiagonal matrices q (ν) q(ν) + q (ν) q(ν) + q (ν) 3 n q(ν) n such that p (ν) n (z) is the characteristic polynomial of T (ν) n Since he was interested in the limit of the zeros of p n (ν) ν it was natural to loo at T (ν) n Martin H Gutnecht and Beresford N Parlett n + q(ν) n as
23 LU (LR) decomposition of a tridiagonal matrix Clearly, T (ν) n has the LU decomposition (LR-Zerlegung) T (ν) n = L (ν) n R (ν) n (9) with L (ν) n = n, R (ν) n = q (ν) q (ν) q (ν) 3 q (ν) n Martin H Gutnecht and Beresford N Parlett
24 The LU (LR) transformation At some historic moment in 954, Rutishauser must have realized that his progressive qd algorithm (9) can be interpreted as computing this LU factorization T (ν) n = L (ν) n R (ν) n and then forming + q (ν) R (ν) n L (ν) n = = T (ν+) n q(ν) + q (ν) q(ν) q (ν) 3 n q(ν) n q (ν) n Martin H Gutnecht and Beresford N Parlett
25 The LR transformation (cont d) So, the qd algorithm consists of performing the step T (ν) n = L (ν) n R (ν) n R (ν) n L (ν) n = T (ν+) n called LR transformation, which is a similarity transformation: ( ) T (ν+) n = R (ν) n T (ν) n R (ν) n The liely motivation: Diagonals ( rows ) of qd table correspond to J fractions, FOPS (Lanczos polynomials), and tridiagonal matrices Rhombus rules lead us from one diagonal to the next They are matched by construction with J and S fractions There are corresponding rules for the polynomials Hence, there must be a rule for transforming one tridiagonal matrix into the next Martin H Gutnecht and Beresford N Parlett
26 The LR algorithm LR algorithm: succession of LR transformations (LR steps) Convergence of ( =,,n) as ν means: Convergence of L (ν) n to diagonal matrix as ν, Convergence of T (ν) n to upper bidiagonal matrix as ν, The diagonals of T (ν) n A, and R (ν) n ultimately contain eigenvalues of Generalization to full matrices is immediate, but unimportant Martin H Gutnecht and Beresford N Parlett
27 The LR algorithm: publications The first two publication on the LR algorithm were in French, two two-page notes in the Comptes Rendues: Rutishauser (955e) [séance du 3 janvier 955], Rutishauser/Bauer (955) [séance du 5 avril 955] 956 Rutishauser produced a mimeographed 5-page ETH research report in English, entitled Report on the Solution of Eigenvalue Problems with the LR transformation Two years later it got properly published in a National Bureau of Standards (NBS) boo series (Rutishauser, 958a) In the same issue: Henrici s review paper on the qd algorithm, and Stiefel s paper on ernel polynomials in NLA In 957, Rutishauser included a 5-page appendix on the LR transformation in the (German) boolet that compiled and updated most of his previous wor on qd (Rutishauser, 957a) Martin H Gutnecht and Beresford N Parlett
28 Conclusions The discovery of the qd and the LR algorithms probably evolved in the following steps: Generalizing Aiten s wor qd table / algorithm Considering the corresponding p table (gen Lanczos wor) and finding the diagonal 3-term recurrence for this table Maing the connection to continued fractions and Lanczos polynomials (and as well to many other topics) Maing the connection to tridiagonal matrices Noticing their extremely simple LU decomposition Noticing that qd algorithm = LR algorithm for tridiagonal matrices Generalizing the LR algorithm to full matrices Martin H Gutnecht and Beresford N Parlett
The kernel structure of rectangular Hankel and Toeplitz matrices
The ernel structure of rectangular Hanel and Toeplitz matrices Martin H. Gutnecht ETH Zurich, Seminar for Applied Mathematics ØØÔ»»ÛÛÛºÑ Ø º Ø Þº» Ñ 3rd International Conference on Structured Matrices
More informationIntroduction. Chapter One
Chapter One Introduction The aim of this book is to describe and explain the beautiful mathematical relationships between matrices, moments, orthogonal polynomials, quadrature rules and the Lanczos and
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Second Edition Translated by R. Bartels, W. Gautschi, and C. Witzgall With 35 Illustrations Springer Contents Preface to the Second Edition Preface
More informationIassume you share the view that the rapid
the Top T HEME ARTICLE THE QR ALGORITHM After a brief sketch of the early days of eigenvalue hunting, the author describes the QR Algorithm and its major virtues. The symmetric case brings with it guaranteed
More informationOn the Vorobyev method of moments
On the Vorobyev method of moments Zdeněk Strakoš Charles University in Prague and Czech Academy of Sciences http://www.karlin.mff.cuni.cz/ strakos Conference in honor of Volker Mehrmann Berlin, May 2015
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 informationMatrix Algorithms. Volume II: Eigensystems. G. W. Stewart H1HJ1L. University of Maryland College Park, Maryland
Matrix Algorithms Volume II: Eigensystems G. W. Stewart University of Maryland College Park, Maryland H1HJ1L Society for Industrial and Applied Mathematics Philadelphia CONTENTS Algorithms Preface xv xvii
More informationInternational Journal of Scientific Research and Reviews
Research article Available online www.ijsrr.org ISSN: 2279 0543 International Journal of Scientific Research and Reviews Eigen Values of non Symmetric Matrices Suman *1 and Kumar Naveen 2 1 Research Scholar,
More informationChasing the Bulge. Sebastian Gant 5/19/ The Reduction to Hessenberg Form 3
Chasing the Bulge Sebastian Gant 5/9/207 Contents Precursers and Motivation 2 The Reduction to Hessenberg Form 3 3 The Algorithm 5 4 Concluding Remarks 8 5 References 0 ntroduction n the early days of
More information6.4 Krylov Subspaces and Conjugate Gradients
6.4 Krylov Subspaces and Conjugate Gradients Our original equation is Ax = b. The preconditioned equation is P Ax = P b. When we write P, we never intend that an inverse will be explicitly computed. P
More informationMatching moments and matrix computations
Matching moments and matrix computations Jörg Liesen Technical University of Berlin and Petr Tichý Czech Academy of Sciences and Zdeněk Strakoš Charles University in Prague and Czech Academy of Sciences
More informationSensitivity of Gauss-Christoffel quadrature and sensitivity of Jacobi matrices to small changes of spectral data
Sensitivity of Gauss-Christoffel quadrature and sensitivity of Jacobi matrices to small changes of spectral data Zdeněk Strakoš Academy of Sciences and Charles University, Prague http://www.cs.cas.cz/
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 informationLinear Least-Squares Data Fitting
CHAPTER 6 Linear Least-Squares Data Fitting 61 Introduction Recall that in chapter 3 we were discussing linear systems of equations, written in shorthand in the form Ax = b In chapter 3, we just considered
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationMiller Objectives Alignment Math
Miller Objectives Alignment Math 1050 1 College Algebra Course Objectives Spring Semester 2016 1. Use algebraic methods to solve a variety of problems involving exponential, logarithmic, polynomial, and
More informationChapter 6. Algebraic eigenvalue problems Introduction Introduction 113. Das also war des Pudels Kern!
6.0. Introduction 113 Chapter 6 Algebraic eigenvalue problems Das also war des Pudels Kern! GOETHE. 6.0. Introduction Determination of eigenvalues and eigenvectors of matrices is one of the most important
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationThe Lanczos and conjugate gradient algorithms
The Lanczos and conjugate gradient algorithms Gérard MEURANT October, 2008 1 The Lanczos algorithm 2 The Lanczos algorithm in finite precision 3 The nonsymmetric Lanczos algorithm 4 The Golub Kahan bidiagonalization
More informationBounds for Eigenvalues of Tridiagonal Symmetric Matrices Computed. by the LR Method. By Gene H. Golub
Bounds for Eigenvalues of Tridiagonal Symmetric Matrices Computed by the LR Method By Gene H. Golub 1. Introduction. In recent years, a number of methods have been proposed for finding the eigenvalues
More information1 Number Systems and Errors 1
Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........
More informationMoments, Model Reduction and Nonlinearity in Solving Linear Algebraic Problems
Moments, Model Reduction and Nonlinearity in Solving Linear Algebraic Problems Zdeněk Strakoš Charles University, Prague http://www.karlin.mff.cuni.cz/ strakos 16th ILAS Meeting, Pisa, June 2010. Thanks
More informationNUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING
NUMERICAL COMPUTATION IN SCIENCE AND ENGINEERING C. Pozrikidis University of California, San Diego New York Oxford OXFORD UNIVERSITY PRESS 1998 CONTENTS Preface ix Pseudocode Language Commands xi 1 Numerical
More informationPADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS
ELLIPTIC INTEGRABLE SYSTEMS PADÉ INTERPOLATION TABLE AND BIORTHOGONAL RATIONAL FUNCTIONS A.S. ZHEDANOV Abstract. We study recurrence relations and biorthogonality properties for polynomials and rational
More informationBlock Krylov Space Solvers: a Survey
Seminar for Applied Mathematics ETH Zurich Nagoya University 8 Dec. 2005 Partly joint work with Thomas Schmelzer, Oxford University Systems with multiple RHSs Given is a nonsingular linear system with
More informationDirect methods for symmetric eigenvalue problems
Direct methods for symmetric eigenvalue problems, PhD McMaster University School of Computational Engineering and Science February 4, 2008 1 Theoretical background Posing the question Perturbation theory
More informationBeresford N. Parlett Λ. April 22, Early History of Eigenvalue Computations 1. 2 The LU and QR Algorithms 3. 3 The Symmetric Case 8
The QR Algorithm Beresford N. Parlett Λ April 22, 2002 Contents 1 Early History of Eigenvalue Computations 1 2 The LU and QR Algorithms 3 3 The Symmetric Case 8 4 The Discovery of the Algorithms 9 Λ Mathematics
More informationBy J H McCabe. to provide the convergents of the simple continued fraction expansions of 5.
An Observation on the Padé Table for + x and the simple continued fraction expansions for, and 5. By J H McCabe. Abstract. The main staircase sequence Padé approximants for + x are shown to yield the convergents
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis Volume 38, pp 184-201, 2011 Copyright 2011, ISSN 1068-9613 ETNA http://etnamathentedu POSITIVITY OF DLV AND MDLVS ALGORITHMS FOR COMPUTING SINGULAR VALUES
More informationZero localization: from 17th-century algebra to challenges of today
Zero localization: from 7th-century algebra to challenges of today UC Berkeley & TU Berlin Bremen summer school July 203 Descartes rule of signs Theorem [Descartes]. The number of positive zeros of a real
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 informationTHE RELATION BETWEEN THE QR AND LR ALGORITHMS
SIAM J. MATRIX ANAL. APPL. c 1998 Society for Industrial and Applied Mathematics Vol. 19, No. 2, pp. 551 555, April 1998 017 THE RELATION BETWEEN THE QR AND LR ALGORITHMS HONGGUO XU Abstract. For an Hermitian
More informationKrylov Space Solvers
Seminar for Applied Mathematics ETH Zurich International Symposium on Frontiers of Computational Science Nagoya, 12/13 Dec. 2005 Sparse Matrices Large sparse linear systems of equations or large sparse
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Algorithms Notes for 2016-10-31 There are several flavors of symmetric eigenvalue solvers for which there is no equivalent (stable) nonsymmetric solver. We discuss four algorithmic ideas: the workhorse
More informationBlock Bidiagonal Decomposition and Least Squares Problems
Block Bidiagonal Decomposition and Least Squares Problems Åke Björck Department of Mathematics Linköping University Perspectives in Numerical Analysis, Helsinki, May 27 29, 2008 Outline Bidiagonal Decomposition
More informationCALCULATION OF GAUSS-KRONROD QUADRATURE RULES. 1. Introduction A(2n+ 1)-point Gauss-Kronrod integration rule for the integral.
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 1133 1145 S 0025-5718(97)00861-2 CALCULATION OF GAUSS-KRONROD QUADRATURE RULES DIRK P. LAURIE Abstract. The Jacobi matrix of the (2n+1)-point
More informationFrancis s Algorithm. David S. Watkins
Francis s Algorithm David S Watkins Abstract John Francis s implicitly shifted QR algorithm turned the problem of matrix eigenvalue computation from difficult to routine almost overnight about fifty years
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 informationIntroduction to Applied Linear Algebra with MATLAB
Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation
More informationINTRODUCTION TO PADÉ APPROXIMANTS. E. B. Saff Center for Constructive Approximation
INTRODUCTION TO PADÉ APPROXIMANTS E. B. Saff Center for Constructive Approximation H. Padé (1863-1953) Student of Hermite Histhesiswon French Academy of Sciences Prize C. Hermite (1822-1901) Used Padé
More informationA new algorithm for computing the Geronimus transformation with large shifts
Noname manuscript No (will be inserted by the editor) A new algorithm for computing the Geronimus transformation with large shifts M I Bueno A Deaño E Tavernetti Received: date / Accepted: date Abstract
More informationREORTHOGONALIZATION FOR THE GOLUB KAHAN LANCZOS BIDIAGONAL REDUCTION: PART I SINGULAR VALUES
REORTHOGONALIZATION FOR THE GOLUB KAHAN LANCZOS BIDIAGONAL REDUCTION: PART I SINGULAR VALUES JESSE L. BARLOW Department of Computer Science and Engineering, The Pennsylvania State University, University
More information9.5. Polynomial and Rational Inequalities. Objectives. Solve quadratic inequalities. Solve polynomial inequalities of degree 3 or greater.
Chapter 9 Section 5 9.5 Polynomial and Rational Inequalities Objectives 1 3 Solve quadratic inequalities. Solve polynomial inequalities of degree 3 or greater. Solve rational inequalities. Objective 1
More informationIntroduction to Numerical Analysis
J. Stoer R. Bulirsch Introduction to Numerical Analysis Translated by R. Bartels, W. Gautschi, and C. Witzgall Springer Science+Business Media, LLC J. Stoer R. Bulirsch Institut fiir Angewandte Mathematik
More informationThe spectral decomposition of near-toeplitz tridiagonal matrices
Issue 4, Volume 7, 2013 115 The spectral decomposition of near-toeplitz tridiagonal matrices Nuo Shen, Zhaolin Jiang and Juan Li Abstract Some properties of near-toeplitz tridiagonal matrices with specific
More informationFinding Zeros of a Polynomial By the Q-D Algorithm
A. C. DOWNING, JR., Editor Finding Zeros of a Polynomial By the Q-D Algorithm P. ~-IENRICI University of California,* Los Angeles, and Eidgenhssisehe Technisehe Hochschule, Zurich, Switzerland BRUCE O.
More informationThe QR Algorithm. Marco Latini. February 26, 2004
The QR Algorithm Marco Latini February 26, 2004 Contents 1 Introduction 2 2 The Power and Inverse Power method 2 2.1 The Power Method............................... 2 2.2 The Inverse power method...........................
More information1 Extrapolation: A Hint of Things to Come
Notes for 2017-03-24 1 Extrapolation: A Hint of Things to Come Stationary iterations are simple. Methods like Jacobi or Gauss-Seidel are easy to program, and it s (relatively) easy to analyze their convergence.
More informationLinear algebra & Numerical Analysis
Linear algebra & Numerical Analysis Eigenvalues and Eigenvectors Marta Jarošová http://homel.vsb.cz/~dom033/ Outline Methods computing all eigenvalues Characteristic polynomial Jacobi method for symmetric
More information6: Polynomials and Polynomial Functions
6: Polynomials and Polynomial Functions 6-1: Polynomial Functions Okay you know what a variable is A term is a product of constants and powers of variables (for example: x ; 5xy ) For now, let's restrict
More informationA Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations
A Method for Constructing Diagonally Dominant Preconditioners based on Jacobi Rotations Jin Yun Yuan Plamen Y. Yalamov Abstract A method is presented to make a given matrix strictly diagonally dominant
More informationIndex. for generalized eigenvalue problem, butterfly form, 211
Index ad hoc shifts, 165 aggressive early deflation, 205 207 algebraic multiplicity, 35 algebraic Riccati equation, 100 Arnoldi process, 372 block, 418 Hamiltonian skew symmetric, 420 implicitly restarted,
More informationLanczos tridiagonalization, Krylov subspaces and the problem of moments
Lanczos tridiagonalization, Krylov subspaces and the problem of moments Zdeněk Strakoš Institute of Computer Science AS CR, Prague http://www.cs.cas.cz/ strakos Numerical Linear Algebra in Signals and
More informationDecomposition of a recursive family of polynomials
Decomposition of a recursive family of polynomials Andrej Dujella and Ivica Gusić Abstract We describe decomposition of polynomials f n := f n,b,a defined by f 0 := B, f 1 (x := x, f n+1 (x = xf n (x af
More information3.3 Dividing Polynomials. Copyright Cengage Learning. All rights reserved.
3.3 Dividing Polynomials Copyright Cengage Learning. All rights reserved. Objectives Long Division of Polynomials Synthetic Division The Remainder and Factor Theorems 2 Dividing Polynomials In this section
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 informationNEW APPROACH TO THE RATIONAL INTERPOLATION PROBLEM
Luminy paper - 1990 1 NEW APPROACH TO THE RATIONAL INTERPOLATION PROBLEM Marc Van Barel, A. Bultheel Abstract We shall reformulate the classical Newton-Padé rational interpolation problem (NPRIP) to take
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 informationSylvester Matrix and GCD for Several Univariate Polynomials
Sylvester Matrix and GCD for Several Univariate Polynomials Manuela Wiesinger-Widi Doctoral Program Computational Mathematics Johannes Kepler University Linz 4040 Linz, Austria manuela.wiesinger@dk-compmath.jku.at
More informationFEM and sparse linear system solving
FEM & sparse linear system solving, Lecture 9, Nov 19, 2017 1/36 Lecture 9, Nov 17, 2017: Krylov space methods http://people.inf.ethz.ch/arbenz/fem17 Peter Arbenz Computer Science Department, ETH Zürich
More informationM. VAN BAREL Department of Computing Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium
MATRIX RATIONAL INTERPOLATION WITH POLES AS INTERPOLATION POINTS M. VAN BAREL Department of Computing Science, K.U.Leuven, Celestijnenlaan 200A, B-3001 Heverlee, Belgium B. BECKERMANN Institut für Angewandte
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 16: Reduction to Hessenberg and Tridiagonal Forms; Rayleigh Quotient Iteration Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationAlg II Syllabus (First Semester)
Alg II Syllabus (First Semester) Unit 1: Solving linear equations and inequalities Lesson 01: Solving linear equations Lesson 02: Solving linear inequalities (See Calculator Appendix A and associated video.)
More informationIn this section again we shall assume that the matrix A is m m, real and symmetric.
84 3. The QR algorithm without shifts See Chapter 28 of the textbook In this section again we shall assume that the matrix A is m m, real and symmetric. 3.1. Simultaneous Iterations algorithm Suppose we
More informationWAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME)
WAITING FOR A BAT TO FLY BY (IN POLYNOMIAL TIME ITAI BENJAMINI, GADY KOZMA, LÁSZLÓ LOVÁSZ, DAN ROMIK, AND GÁBOR TARDOS Abstract. We observe returns of a simple random wal on a finite graph to a fixed node,
More informationPartial Fractions. Combining fractions over a common denominator is a familiar operation from algebra: 2 x 3 + 3
Partial Fractions Combining fractions over a common denominator is a familiar operation from algebra: x 3 + 3 x + x + 3x 7 () x 3 3x + x 3 From the standpoint of integration, the left side of Equation
More informationScientific Computing with Case Studies SIAM Press, Lecture Notes for Unit VII Sparse Matrix
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit VII Sparse Matrix Computations Part 1: Direct Methods Dianne P. O Leary c 2008
More informationIterative Methods for Linear Systems of Equations
Iterative Methods for Linear Systems of Equations Projection methods (3) ITMAN PhD-course DTU 20-10-08 till 24-10-08 Martin van Gijzen 1 Delft University of Technology Overview day 4 Bi-Lanczos method
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.4: Krylov Subspace Methods 2 / 51 Krylov Subspace Methods In this chapter we give
More informationAnalysis of Carleman Representation of Analytical Recursions
JOURNAL OF MATHEMATICAL ANALYSIS AND APPLICATIONS 224, 8190 1998 ARTICLE NO. AY985986 Analysis of Carleman Representation of Analytical Recursions G. Berolaio Department of Physics, Bar-Ilan Uniersity,
More informationTHE QR METHOD A = Q 1 R 1
THE QR METHOD Given a square matrix A, form its QR factorization, as Then define A = Q 1 R 1 A 2 = R 1 Q 1 Continue this process: for k 1(withA 1 = A), A k = Q k R k A k+1 = R k Q k Then the sequence {A
More informationOn rational approximation of algebraic functions. Julius Borcea. Rikard Bøgvad & Boris Shapiro
On rational approximation of algebraic functions http://arxiv.org/abs/math.ca/0409353 Julius Borcea joint work with Rikard Bøgvad & Boris Shapiro 1. Padé approximation: short overview 2. A scheme of rational
More informationA QR-decomposition of block tridiagonal matrices generated by the block Lanczos process
1 A QR-decomposition of block tridiagonal matrices generated by the block Lanczos process Thomas Schmelzer Martin H Gutknecht Oxford University Computing Laboratory Seminar for Applied Mathematics, ETH
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationStructured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries
Structured Krylov Subspace Methods for Eigenproblems with Spectral Symmetries Fakultät für Mathematik TU Chemnitz, Germany Peter Benner benner@mathematik.tu-chemnitz.de joint work with Heike Faßbender
More informationMS&E 318 (CME 338) Large-Scale Numerical Optimization
Stanford University, Management Science & Engineering (and ICME MS&E 38 (CME 338 Large-Scale Numerical Optimization Course description Instructor: Michael Saunders Spring 28 Notes : Review The course teaches
More informationApplied Numerical Analysis
Applied Numerical Analysis Using MATLAB Second Edition Laurene V. Fausett Texas A&M University-Commerce PEARSON Prentice Hall Upper Saddle River, NJ 07458 Contents Preface xi 1 Foundations 1 1.1 Introductory
More informationComputational Methods. Eigenvalues and Singular Values
Computational Methods Eigenvalues and Singular Values Manfred Huber 2010 1 Eigenvalues and Singular Values Eigenvalues and singular values describe important aspects of transformations and of data relations
More informationThe Lanczos Algorithm and Hankel Matrix Factorization
The Lanczos Algorithm and Hankel Matrix Factorization Daniel L Boley Department of Computer Science University of Minnesota Minneapolis, Minnesota Tong J Lee and Franklin T Luk School of Electrical Engineering
More informationPakistan Urdu School- Kingdom of Bahrain. Curriculum Implementation Plan for MATHEMATICS Grade XI Sc.
Pakistan Urdu School- Kingdom of Bahrain No. Curriculum Implementation Plan for MATHEMATICS Grade XI Sc. Unit 9 Fundamentals of Trigonometry 1 st Term 1 April/4 th Fundamental Identities 2 23-04-2018 Values
More informationLecture # 11 The Power Method for Eigenvalues Part II. The power method find the largest (in magnitude) eigenvalue of. A R n n.
Lecture # 11 The Power Method for Eigenvalues Part II The power method find the largest (in magnitude) eigenvalue of It makes two assumptions. 1. A is diagonalizable. That is, A R n n. A = XΛX 1 for some
More informationESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM
ESTIMATES OF THE TRACE OF THE INVERSE OF A SYMMETRIC MATRIX USING THE MODIFIED CHEBYSHEV ALGORITHM GÉRARD MEURANT In memory of Gene H. Golub Abstract. In this paper we study how to compute an estimate
More informationAlexander Ostrowski
Ostrowski p. 1/3 Alexander Ostrowski 1893 1986 Walter Gautschi wxg@cs.purdue.edu Purdue University Ostrowski p. 2/3 Collected Mathematical Papers Volume 1 Determinants Linear Algebra Algebraic Equations
More informationThe Singular Value Decomposition
The Singular Value Decomposition An Important topic in NLA Radu Tiberiu Trîmbiţaş Babeş-Bolyai University February 23, 2009 Radu Tiberiu Trîmbiţaş ( Babeş-Bolyai University)The Singular Value Decomposition
More informationResearch Article On Maslanka s Representation for the Riemann Zeta Function
International Mathematics and Mathematical Sciences Volume 200, Article ID 7447, 9 pages doi:0.55/200/7447 Research Article On Maslana s Representation for the Riemann Zeta Function Luis Báez-Duarte Departamento
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 informationThe purpose of the second Lanczos paper on the subject was \to adopt the general principles of the previous investigation to the specic demands that a
THE UNSYMMETRIC LANCZOS ALGORITHMS AND THEIR RELATIONS TO PADE APPROXIMATION, CONTINUED FRACTIONS, AND THE QD ALGORITHM MARTIN H. GUTKNECHT Abstract. First, several algorithms based on the unsymmetric
More informationRepresentation of doubly infinite matrices as non-commutative Laurent series
Spec. Matrices 217; 5:25 257 Research Article Open Access María Ivonne Arenas-Herrera and Luis Verde-Star* Representation of doubly infinite matrices as non-commutative Laurent series https://doi.org/1.1515/spma-217-18
More informationTikhonov Regularization of Large Symmetric Problems
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 2000; 00:1 11 [Version: 2000/03/22 v1.0] Tihonov Regularization of Large Symmetric Problems D. Calvetti 1, L. Reichel 2 and A. Shuibi
More informationnonlinear simultaneous equations of type (1)
Module 5 : Solving Nonlinear Algebraic Equations Section 1 : Introduction 1 Introduction Consider set of nonlinear simultaneous equations of type -------(1) -------(2) where and represents a function vector.
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 95-222, 2002. Copyright 2002,. ISSN 068-963. ETNA ASYMPTOTICS FOR QUADRATIC HERMITE-PADÉ POLYNOMIALS ASSOCIATED WITH THE EXPONENTIAL FUNCTION
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 informationNumerical Analysis Fall. Gauss Elimination
Numerical Analysis 2015 Fall Gauss Elimination Solving systems m g g m m g x x x k k k k k k k k k 3 2 1 3 2 1 3 3 3 2 3 2 2 2 1 0 0 Graphical Method For small sets of simultaneous equations, graphing
More informationTEACHING NUMERICAL LINEAR ALGEBRA AT THE UNDERGRADUATE LEVEL by Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University
TEACHING NUMERICAL LINEAR ALGEBRA AT THE UNDERGRADUATE LEVEL by Biswa Nath Datta Department of Mathematical Sciences Northern Illinois University DeKalb, IL 60115 E-mail: dattab@math.niu.edu What is Numerical
More informationA Modified Method for Reconstructing Periodic Jacobi Matrices
mathematics of computation volume 42. number 165 january 1984, pages 143 15 A Modified Method for Reconstructing Periodic Jacobi Matrices By Daniel Boley and Gene H. Golub* Abstract. In this note, we discuss
More informationNon-stationary extremal eigenvalue approximations in iterative solutions of linear systems and estimators for relative error
on-stationary extremal eigenvalue approximations in iterative solutions of linear systems and estimators for relative error Divya Anand Subba and Murugesan Venatapathi* Supercomputer Education and Research
More informationIntegration of Rational Functions by Partial Fractions
Title Integration of Rational Functions by MATH 1700 MATH 1700 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 2 / 11 Rational functions A rational function is one of the form where P and Q are
More informationMatrices, Moments and Quadrature, cont d
Jim Lambers CME 335 Spring Quarter 2010-11 Lecture 4 Notes Matrices, Moments and Quadrature, cont d Estimation of the Regularization Parameter Consider the least squares problem of finding x such that
More informationFoundations of Discrete Mathematics
Foundations of Discrete Mathematics Chapter 0 By Dr. Dalia M. Gil, Ph.D. Statement Statement is an ordinary English statement of fact. It has a subject, a verb, and a predicate. It can be assigned a true
More information