Abstract:This thesis is concerned with two classical topics in matrix computations : The QR algorithm for solving nonsymmetric eigenvalue problems and the computation of matrix exponentials for two types of structured matrices. We focus on the performance in the former topic and on accuracy in the latter one. For computing all eigenvalues of a non-Hermitian matrix, the QR algorithm which iteratively computes a Schur decomposition of the matrix is the method of choice. We present a new parallel implementation of the multishift QR algorithm targeting distributed memory architectures. Starting from recent developments of the parallel multishift QR algorithm, we propose a number of algorithmic and implementation improvements. Guidelines concerning several important tunable algorithmic parameters are also provided. Numerous computational experiments confirm that our new implementation significantly outperforms previous parallel implementations of the QR algorithm. The computation of the exponential of a square matrix is also an important task in matrix computations. For a general dense matrix, the scaling and squaring method coupled with Pade approximation is the most popular approach. However, for an essentially nonnegative matrix (a real square matrix with nonnegative off-diagonal entries), truncated Taylor series rather than Pade approximation is preferred to achieve componentwise accuracy in the matrix exponential. We propose a method which efficiently computes all entries of the exponential of an essentially nonnegative matrix to high relative accuracy. Truncation and rounding error bounds, as well as numerical experiments demonstrate the efficiency and accuracy of our method. When the matrix is banded, the entries of its matrix exponential decay exponentially away from the main diagonal. We analyze the decay property for the exponentials of several classes of doubly-infinite skew-Hermitian matrices. Then finite section methods based on the decay property are established. We also propose a repeated doubling strategy which works well even when a priori error estimates are pessimistic or not easy to compute. Finally, numerical experiments are presented to illustrate the effectiveness of the finite section method.

Convergence Characteristic of Single Direction Shrink QR Algorithm in the Singular Value Decomposition

Global Convergence of Hessenberg Shifted QR I: Exact Arithmetic

A Quadratically Convergent QR-like Method Without Shifts for the Hermitian Eigenvalue Problem

A Convergence Analysis of the Inexact Rayleigh Quotient Iteration and Simplified Jacobi-Davidson Method for the Large Hermitian Matrix Eigenproblem

A Parallelizable Method for Computing the Singular Subspaces of a General Matrix

A randomized algorithm for the QR decomposition-based approximate SVD

Linear Convergence of an Algorithm for Largest Singular Value of a Nonnegative Rectangular Tensor

The Singular Value Decomposition, Applications and Beyond.

Dense and Structured Matrix Computations —the Parallel QR Algorithm and Matrix Exponentials

Convergence of a Ulm-like Method for Square Inverse Singular Value Problems with Multiple and Zero Singular Values.

The Convergence of Krylov Subspace Methods for Large Unsymmetric Linear Systems

The Generalized Triangular Decomposition.

The Joint Bidiagonalization Method for Large GSVD Computations in Finite Precision.

Perturbation Analysis of Randomized SVD and its Applications to High-dimensional Statistics

Singular Value Decomposition For Unitary Symmetric Matrix

An Implicitly Restarted Refined Bidiagonalization Lanczos Method for Computing a Partial Singular Value Decomposition.

A Power Method for Computing Singular Value Decomposition

Effective Methods of QR-Decompositions of Square Complex Matrices by Fast Discrete Signal-Induced Heap Transforms

Implementing Communication-Optimal Parallel and Sequential QR Factorizations

High Order Numerical Methods To Approximate The Singular Value Decomposition

Towards a more robust algorithm for computing the restricted singular value decomposition