Abstract
The two-sided Jacobi method for singular value decomposition (SVD) has the advantage of obtaining singular vectors quickly and accurately. In previous research, fast and accurate implementations of the two-sided Jacobi method have been achieved for real matrices. In this study, we implemented SVD for complex matrices using the two-sided Jacobi method. In SVD, given rectangular matrices can be converted into upper-triangular matrices by conducting QR decomposition as a preprocessing step. Then, the upper-triangular matrices are decomposed. In the case where the given matrices are complex, the upper-triangular matrices are complex. Two implementations are proposed: the first requires SVD for complex matrices, whereas the second requires SVD for only real matrices.






Similar content being viewed by others
Availability of data and materials
The data used in the paper will be available upon request.
References
Kogbetliantz E (1955) Solutions of linear equations by diagonalization of coefficient matrices. Q Appl Math 13:123–132
Rutishauser H (1996) The Jacobi method for real symmetric matrices. Numer Math 9(1):1–10
Drmac Z (1994) Computing the singular and the generalized singular value. Ph.D. University of Hagen, Germany
Drmac Z (1999) A posteriori computation of singular vectors in a preconditioned Jacobi SVD algorithm. IMA J Numer Anal 19:191–213
Drmac Z, Veselic K (2005) New Fast and Accurate Jacobi SVD Algorithm: I. LAPACK Working Note 169. http://www.netlib.org/lapack/lawns/downloads/
Drmac Z, Veselic K (2005) New fast and accurate Jacobi SVD algorithm: II. LAPACK working note 170. http://www.netlib.org/lapack/lawns/downloads/
Linear Algebra PACKage. http://www.netlib.org/lapack/
Takata M, Araki S, Miyamae T, Kimura K, Nakamura Y (2021) On an implementation of the one-sided Jacobi method with high accuracy. IPSJ Trans Math Model Appl 14(3):68–75
Araki S, Aoki M, Takata M, Kimura K, Nakamura Y (2021) On an implementation of the two-sided Jacobi method. IPSJ Trans Math Model Appl 14(1):12–20
Sakurai T, Tadano H (2007) CIRR: a Rayleigh–Ritz type method with counter integral for generalized eigenvalue problems. Hokkaido Math J 36:745–757
Odedo V, Yavuz ME, Costen F, Himeno R, Yokota H (2017) Time reversal technique based on spatiotemporal windows for through wall imaging. IEEE Trans Antennas Propag 65(6):3065–3072
Imakura A, Sakurai T (2018) Block SS-CAA: a complex moment-based parallel nonlinear eigensolver using the block communication-avoiding Arnoldi procedure. Parallel Comput 74:34–48
Rutishauser H (1966) Jacobi method for real symmetric matrices. Numer Math 9(1):1–10
Araki S, Aoki M, Takata M, Kimura K, Nakamura Y (2021) Implementation of two-sided Jacobi method. IPSJ Trans Model Appl 14(1):12–20
Acknowledgements
This work was supported by JST SPRING, Grant Number JPMJSP2115. We would like to thank Editage (www.editage.jp) for English language editing.
Funding
JST SPRING, Grant Number JPMJSP2115.
Author information
Authors and Affiliations
Contributions
MC, MT and KK wrote the main manuscript text. All authors reviewed the manuscript.
Corresponding authors
Ethics declarations
Conflict of interest
The authors declare no competing interests.
Ethics approval
Not applicable.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Chiyonobu, M., Miyamae, T., Takata, M. et al. Singular value decomposition for complex matrices using two-sided Jacobi method. J Supercomput 80, 11719–11740 (2024). https://doi.org/10.1007/s11227-024-05903-6
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11227-024-05903-6