A Task-Based Algorithm for Reordering the Eigenvalues of a Matrix in Real Schur Form

Representations for the image-kernel (p,q)-inverses of block matrices in rings

Georgian Mathematical Journal ◽

10.1515/gmj-2018-0042 ◽

2020 ◽

Vol 27 (2) ◽

pp. 297-305

Author(s):

Dijana Mosić

Keyword(s):

Block Matrix ◽

Block Matrices ◽

Schur Form

AbstractWe present the conditions for a block matrix of a ring to have the image-kernel{(p,q)}-inverse in the generalized Banachiewicz–Schur form. We give representations for the image-kernel inverses of the sum and the product of two block matrices. Some characterizations of the image-kernel{(p,q)}-inverse in a ring with involution are investigated too.

Download Full-text

Algorithm 1019: A Task-based Multi-shift QR/QZ Algorithm with Aggressive Early Deflation

ACM Transactions on Mathematical Software ◽

10.1145/3495005 ◽

2022 ◽

Vol 48 (1) ◽

pp. 1-36

Author(s):

Mirko Myllykoski

Keyword(s):

Eigenvalue Problems ◽

Critical Path ◽

Hessenberg Matrix ◽

Single Node ◽

Standard Case ◽

Qr Algorithm ◽

Schur Form ◽

Tightly Coupled ◽

Qz Algorithm ◽

New Ideas

The QR algorithm is one of the three phases in the process of computing the eigenvalues and the eigenvectors of a dense nonsymmetric matrix. This paper describes a task-based QR algorithm for reducing an upper Hessenberg matrix to real Schur form. The task-based algorithm also supports generalized eigenvalue problems (QZ algorithm) but this paper concentrates on the standard case. The task-based algorithm adopts previous algorithmic improvements, such as tightly-coupled multi-shifts and Aggressive Early Deflation (AED) , and also incorporates several new ideas that significantly improve the performance. This includes, but is not limited to, the elimination of several synchronization points, the dynamic merging of previously separate computational steps, the shortening and the prioritization of the critical path, and experimental GPU support. The task-based implementation is demonstrated to be multiple times faster than multi-threaded LAPACK and ScaLAPACK in both single-node and multi-node configurations on two different machines based on Intel and AMD CPUs. The implementation is built on top of the StarPU runtime system and is part of the open-source StarNEig library.

Download Full-text