Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Skip to main content

A note on the recursive calculation of dominant singular subspaces

  • Published:
Numerical Algorithms Aims and scope Submit manuscript

Abstract

A recursive procedure for computing an approximation of the left and right dominant singular subspaces of a given matrix is proposed in [1]. The method is particularly suited for matrices with many more rows than columns. The procedure consists of a few steps. In one of these steps a Householder transformation is multiplied to an upper triangular matrix. The following step consists in recomputing an upper triangular matrix from the latter product. In [1] it is said that the latter step is accomplished in O(k3) operations, where k is the order of the triangular matrix. In this short note we show that this step can be accomplished in O(k2) operations.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save

Springer+ Basic
$34.99 /Month
  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime
Subscribe now

Buy Now

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  1. Y. Chahlaoui, K. Gallivan and P. Van Dooren, Recursive calculation of dominant singular subspaces, SIAM J. Matrix Anal. Appl. 25(2) (2004) 445–463.

    Google Scholar 

  2. T.F. Chan, An improved algorithm for computing the singular value decomposition, ACM Trans. Math. Software 8(1) (1982) 72–83.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Nicola Mastronardi.

Additional information

Communicated by P. Van Dooren

This research was partially supported by MIUR, grant number 2002014121 (first author) and by the Research Council K.U.Leuven, project OT/00/16 (SLAP: Structured Linear Algebra Package), by the Fund for Scientific Research–Flanders (Belgium), projects G.0078.01 (SMA: Structured Matrices and their Applications), G.0176.02 (ANCILA: Asymptotic aNalysis of the Convergence behavior of Iterative methods in numerical Linear Algebra), G.0184.02 (CORFU: Constructive study of Orthogonal Functions) and G.0455.0 (RHPH: Riemann–Hilbert problems, random matrices and Padé–Hermite approximation), and by the Belgian Programme on Interuniversity Poles of Attraction, initiated by the Belgian State, Prime Minister’s Office for Science, Technology and Culture, project IUAP V-22 (Dynamical Systems and Control: Computation, Identification & Modelling) (second and third author). The scientific responsibility rests with the authors.

AMS subject classification

15A15, 15A09, 15A23

Rights and permissions

Reprints and permissions

About this article

Cite this article

Mastronardi, N., Van Barel, M. & Vandebril, R. A note on the recursive calculation of dominant singular subspaces. Numer Algor 38, 237–242 (2005). https://doi.org/10.1007/s11075-004-4338-x

Download citation

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11075-004-4338-x

Keywords