Robust level-3 BLAS Inverse Iteration from the Hessenberg Matrix

01/13/2021
by   Angelika Schwarz, et al.
0

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.

READ FULL TEXT

Authors

page 1

page 2

page 3

page 4

06/24/2021

A sparse approximate inverse for triangular matrices based on Jacobi iteration

In this paper, we propose a sparse approximate inverse for triangular ma...
05/13/2022

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...
05/20/2020

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...
10/19/2020

Probabilistic Linear Solvers for Machine Learning

Linear systems are the bedrock of virtually all numerical computation. M...
01/06/2022

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

Computing effective eigenvalues for neutron transport often requires a f...
05/31/2022

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.