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.

Extension of accurate numerical algorithms for matrix multiplication based on error-free transformation

Improving the numerical stability of fast matrix multiplication

Acceleration of complex matrix multiplication using arbitrary precision floating-point arithmetic

Tuning Technique for Multiple Precision Dense Matrix Multiplication using Prediction of Computational Time

Acceleration of Approximate Matrix Multiplications on GPUs

Performance Enhancement of the Ozaki Scheme on Integer Matrix Multiplication Unit

Highly Accurate Doubling Algorithms for M -Matrix Algebraic Riccati Equations

Acceleration of Multiple Precision Matrix Multiplication using Ozaki scheme

Computing the matrix fractional power with the double exponential formula

Acceleration of multiple precision matrix multiplication based on multi-component floating-point arithmetic using AVX2

Efficient scaling and squaring method for the matrix exponential

Floating-Point Multiply-Add with Approximate Normalization for Low-Cost Matrix Engines

Accelerated Multiple Precision Matrix Multiplication using Strassen's Algorithm and Winograd's Variant

Precision-Aware Iterative Algorithms Based on Group-Shared Exponents of Floating-Point Numbers

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

Permutation-Based Approximate Multiplier with High Accuracy.

Matrix-by-matrix multiplication algorithm with $O(N^2log_2N)$ computational complexity for variable precision arithmetic

A method for accelerating low precision operations by sparse matrix multiplication

Roundoff error analysis of the double exponential formula-based method for the matrix sign function

Efficient numerical algorithms for multi-precision and multi-accuracy calculation of the error functions and Dawson integral with complex arguments

Adaptive multiplication of $\mathcal{H}^2$-matrices with block-relative error control