Robust level-3 BLAS Inverse Iteration from the Hessenberg Matrix

by   Angelika Schwarz, et al.

Inverse iteration is known to be an effective method for computing eigenvectors corresponding to simple and well-separated eigenvalues. In the non-symmetric case, the solution of shifted Hessenberg systems is a central step. Existing inverse iteration solvers approach the solution of the shifted Hessenberg systems with either RQ or LU factorizations and, once factored, solve the corresponding systems. This approach has limited level-3 BLAS potential since distinct shifts have distinct factorizations. This paper rearranges the RQ approach such that data shared between distinct shifts is exposed. Thereby the backward substitution with the triangular R factor can be expressed mostly with matrix-matrix multiplications (level-3 BLAS). The resulting algorithm computes eigenvectors in a tiled, overflow-free, and task-parallel fashion. The numerical experiments show that the new algorithm outperforms existing inverse iteration solvers for the computation of both real and complex eigenvectors.



page 1

page 2

page 3

page 4


A sparse approximate inverse for triangular matrices based on Jacobi iteration

In this paper, we propose a sparse approximate inverse for triangular ma...

Global Convergence of Hessenberg Shifted QR III: Approximate Ritz Values via Shifted Inverse Iteration

We give a self-contained randomized algorithm based on shifted inverse i...

Single-step triangular splitting iteration method for a class of complex symmetric linear system

In this work, a single-step triangular splitting (SSTS) iteration method...

Probabilistic Linear Solvers for Machine Learning

Linear systems are the bedrock of virtually all numerical computation. M...

A low-rank power iteration scheme for neutron transport criticality problems

Computing effective eigenvalues for neutron transport often requires a f...

On Forward Kinematics of a 3SPR Parallel Manipulator

In this paper, a new numerical method to solve the forward kinematics (F...
This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.