You are here: MIMS > EPrints
MIMS EPrints

2006.69: $J$-Orthogonal Matrices: Properties and Generation

2006.69: Nicholas J. Higham (2003) $J$-Orthogonal Matrices: Properties and Generation. SIAM Review, 45 (3). pp. 504-519. ISSN 0036-1445

Full text available as:

PDF - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
191 Kb

DOI: 10.1137/S0036144502414930


A real, square matrix $Q$ is $J$-orthogonal if $Q^TJQ = J$, where the signature matrix $J = \diag(\pm 1)$. $J$-orthogonal matrices arise in the analysis and numerical solution of various matrix problems involving indefinite inner products, including, in particular, the downdating of Cholesky factorizations. We present techniques and tools useful in the analysis, application and construction of these matrices, giving a self-contained treatment that provides new insights. First, we define and explore the properties of the exchange operator, which maps $J$-orthogonal matrices to orthogonal matrices and vice versa. Then we show how the exchange operator can be used to obtain a hyperbolic CS decomposition of a $J$-orthogonal matrix directly from the usual CS decomposition of an orthogonal matrix. We employ the decomposition to derive an algorithm for constructing random $J$-orthogonal matrices with specified norm and condition number. We also give a short proof of the fact that $J$-orthogonal matrices are optimally scaled under two-sided diagonal scalings. We introduce the indefinite polar decomposition and investigate two iterations for computing the $J$-orthogonal polar factor: a Newton iteration involving only matrix inversion and a Schulz iteration involving only matrix \mult. We show that these iterations can be used to $J$-orthogonalize a matrix that is not too far from being $J$-orthogonal.

Item Type:Article
Uncontrolled Keywords:$J$-orthogonal matrix, exchange operator, gyration operator, sweep operator, principal pivot transform, hyperbolic CS decomposition, two-sided scaling, indefinite least squares problem, hyperbolic QR factorization, indefinite polar decomposition, Newton's method, Schulz iteration
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2006.69
Deposited By:Nick Higham
Deposited On:02 May 2006

Download Statistics: last 4 weeks
Repository Staff Only: edit this item