A Kogbetliantz-type algorithm for the hyperbolic SVD

03/14/2020
by   Vedran Novaković, et al.
0

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.

READ FULL TEXT
POST COMMENT

Comments

There are no comments yet.

Authors

page 1

page 2

page 3

page 4

05/15/2020

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 ...
07/19/2019

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 ...
02/16/2022

Vectorization of the Jacobi-type singular value decomposition method

The eigenvalue decomposition (EVD) of (a batch of) Hermitian matrices of...
10/13/2020

Projection techniques to update the truncated SVD of evolving matrices

This paper considers the problem of updating the rank-k truncated Singul...
09/05/2021

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

Recent studies have demonstrated improved skill in numerical weather pre...
07/04/2017

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...
08/31/2019

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.