ISSN:
1070-5325

Keywords:
complex symmetric matrices
;
singular value decomposition
;
iterative methods
;
Engineering
;
Numerical Methods and Modeling

Source:
Wiley InterScience Backfile Collection 1832-2000

Topics:
Mathematics

Notes:
We present an algorithm for the approximation of the dominant singular values and corresponding right and left singular vectors of a complex symmetric matrix. The method is based on two short-term recurrences first proposed by Saunders, Simon and Yip [24] for a non-Hermitian linear system solver. With symmetric matrices, the recurrence can be modified so as to generate a tridiagonal symmetric matrix from which the original triplets can be approximated. The recurrence formally resembles the Lanczos method, in spite of substantial differences which make usual convergence results inapplicable. Implementation aspects are discussed, such as re-orthogonalization and the use of alternative representation matrices. The method is very efficient over existing approaches which do not exploit the symmetry of the problem. Numerical experiments on application problems validate the analysis, while showing satisfactory results, especially on dense matrices. © 1997 by John Wiley & Sons, Ltd.

Additional Material:
4 Ill.

Type of Medium:
Electronic Resource

Permalink