A Kogbetliantz-type algorithm for the hyperbolic SVD

by   Vedran Novaković, et al.

In this paper a two-sided, parallel Kogbetliantz-type algorithm for the hyperbolic singular value decomposition (HSVD) of real and complex square matrices is developed, with a single assumption that the input matrix, of order n, admits such a decomposition into the product of a unitary, a non-negative diagonal, and a J-unitary matrix, where J is a given diagonal matrix of positive and negative signs. When J=± I, the proposed algorithm computes the ordinary SVD. The paper's most important contribution – a derivation of formulas for the HSVD of 2× 2 matrices – is presented first, followed by the details of their implementation in floating-point arithmetic. Next, the effects of the hyperbolic transformations on the columns of the iteration matrix are discussed. These effects then guide a redesign of the dynamic pivot ordering, being already a well-established pivot strategy for the ordinary Kogbetliantz algorithm, for the general, n× n HSVD. A heuristic but sound convergence criterion is then proposed, which contributes to high accuracy demonstrated in the numerical testing results. Such a J-Kogbetliantz algorithm as presented here is intrinsically slow, but is nevertheless usable for matrices of small orders.



There are no comments yet.


page 1

page 2

page 3

page 4


Batched computation of the singular value decompositions of order two by the AVX-512 vectorization

In this paper a vectorized algorithm for simultaneously computing up to ...

The LAPW method with eigendecomposition based on the Hari–Zimmermann generalized hyperbolic SVD

In this paper we propose an accurate, highly parallel algorithm for the ...

Vectorization of the Jacobi-type singular value decomposition method

The eigenvalue decomposition (EVD) of (a batch of) Hermitian matrices of...

Projection techniques to update the truncated SVD of evolving matrices

This paper considers the problem of updating the rank-k truncated Singul...

Efficient computation of matrix-vector products with full observation weighting matrices in data assimilation

Recent studies have demonstrated improved skill in numerical weather pre...

The Normalized Singular Value Decomposition of Non-Symmetric Matrices Using Givens fast Rotations

In this paper we introduce the algorithm and the fixed point hardware to...

Implicit Hari–Zimmermann algorithm for the generalized SVD on the GPUs

A parallel, blocked, one-sided Hari–Zimmermann algorithm for the general...
This week in AI

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