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

05/15/2020
by   Vedran Novaković, et al.
0

In this paper a vectorized algorithm for simultaneously computing up to eight singular value decompositions (SVDs, each of the form A=UΣ V^∗) of real or complex matrices of order two is proposed. The algorithm extends to a batch of matrices of an arbitrary length n, that arises, for example, in the annihilation part of the parallel Kogbetliantz algorithm for the SVD of a square matrix of order 2n. The SVD algorithm for a single matrix of order two is derived first. It scales, in most instances error-free, the input matrix A such that its singular values Σ_ii cannot overflow whenever its elements are finite, and then computes the URV factorization of the scaled matrix, followed by the SVD of a non-negative upper-triangular middle factor. A vector-friendly data layout for the batch is then introduced, where the same-indexed elements of each of the input and the output matrices form vectors, and the algorithm's steps over such vectors are described. The vectorized approach is then shown to be about three times faster than processing each matrix in isolation, while slightly improving accuracy over the straightforward method for the 2× 2 SVD.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
03/23/2022

Phase Factors in Singular Value Decomposition and Schmidt Decomposition

In singular value decomposition (SVD) of a complex matrix A, the singula...
research
02/16/2022

Vectorization of a thread-parallel Jacobi singular value decomposition method

The eigenvalue decomposition (EVD) of (a batch of) Hermitian matrices of...
research
03/14/2020

A Kogbetliantz-type algorithm for the hyperbolic SVD

In this paper a two-sided, parallel Kogbetliantz-type algorithm for the ...
research
02/22/2023

Singular value decomposition based matrix surgery

This paper aims to develop a simple procedure to reduce and control the ...
research
09/28/2010

Use of multiple singular value decompositions to analyze complex intracellular calcium ion signals

We compare calcium ion signaling (Ca^2+) between two exposures; the data...
research
03/10/2020

Error Estimation for Sketched SVD via the Bootstrap

In order to compute fast approximations to the singular value decomposit...
research
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...

Please sign up or login with your details

Forgot password? Click here to reset