Premium
GR decompositions and their relations to Cholesky‐like factorizations
Author(s) -
Benner Peter,
Penke Carolin
Publication year - 2021
Publication title -
pamm
Language(s) - English
Resource type - Journals
ISSN - 1617-7061
DOI - 10.1002/pamm.202000065
Subject(s) - cholesky decomposition , qr decomposition , mathematics , matrix decomposition , factorization , skew symmetric matrix , matrix (chemical analysis) , scalar (mathematics) , symplectic geometry , pure mathematics , algebra over a field , symmetric matrix , eigenvalues and eigenvectors , algorithm , square matrix , geometry , physics , materials science , quantum mechanics , composite material
Abstract For a given matrix, we are interested in computing GR decompositions A = GR , where G is an isometry with respect to given scalar products. The orthogonal QR decomposition is the representative for the Euclidian scalar product. For a signature matrix, a respective factorization is given as the hyperbolic QR decomposition. Considering a skew‐symmetric matrix leads to the symplectic QR decomposition. The standard approach for computing GR decompositions is based on the successive elimination of subdiagonal matrix entries. For the hyperbolic and symplectic case, this approach does in general not lead to a satisfying numerical accuracy. An alternative approach computes the QR decomposition via a Cholesky factorization, but also has bad stability. It is improved by repeating the procedure a second time. In the same way, the hyperbolic and the symplectic QR decomposition are related to the LDL T and a skew‐symmetric Cholesky‐like factorization. We show that methods exploiting this connection can provide better numerical stability than elimination‐based approaches.