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.

On the finite section method for computing exponentials of doubly-infinite skew-Hermitian matrices

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

High-Performance Computation of the Exponential of a Large Sparse Matrix

Efficient scaling and squaring method for the matrix exponential

Computation of the exponential function of matrices by a formula without oscillatory integrals on infinite intervals

Preconditioned Hermitian and Skew-Hermitian Splitting Methods for Non-Hermitian Positive Semidefinite Linear Systems

A Skew-Symmetric Lanczos Bidiagonalization Method for Computing Several Extremal Eigenpairs of a Large Skew-Symmetric Matrix

Concerning the stability of exponential systems and Fourier matrices

On Analysis Of Projection Methods For Rational Function Approximation To The Matrix Exponential

Computing exponentials of essentially non-negative matrices entrywise to high relative accuracy.

On the Exponentiation of Interval Matrices

Shift-splitting Iteration Methods for a Class of Large Sparse Linear Matrix Equations

An Interesting Method For The Exponentials For Some Special Matrices

The eigenvalue decomposition of normal matrices by the decomposition of the skew-symmetric part with applications to orthogonal matrices

Matrix‐less methods for the spectral approximation of large non‐Hermitian Toeplitz matrices: A concise theoretical analysis and a numerical study

Fast Non-Hermitian Toeplitz Eigenvalue Computations, Joining Matrixless Algorithms and FDE Approximation Matrices

Low-synchronization Arnoldi Methods for the Matrix Exponential with Application to Exponential Integrators

Asymmetric Hermitian and skew-Hermitian splitting methods for positive definite linear systems

Skew-Symmetric Matrix Decompositions on Shared-Memory Architectures

Explicit formulas for the inverses of Toeplitz matrices, with applications

Local approximations of inverse block Toeplitz matrices and Baxter-type theorems for long-memory processes