Abstract
In this paper, we consider the product eigenvalue problem for a wide class of structured matrices containing the well-known Vandermonde and Cauchy matrices. A periodic qd-type reduction method is developed for computing eigenvalues of products of these rectangular matrices so that no subtraction of like-signed numbers occurs. Consequently, all the eigenvalues of such a product are computed to high relative accuracy in a preferable complexity. Error analysis and numerical experiments are provided to confirm the high relative accuracy.
Similar content being viewed by others
References
Alfa, A.S., Xue, J., Ye, Q.: Accurate computation of the smallest eigenvalue of a diagonally dominant M-matrix. Math. Comput. 71, 217–236 (2002)
Barlow, J., Demmel, J.: Computing accurate eigensystems of scaled diagonally dominant matrices. SIAM J. Numer. Anal. 27, 762–791 (1990)
Bevilacqua, R., Bozzo, E., Corso, G.M.D.: qd-Type methods for quasiseparable matrices. SIAM J. Matrix Anal. Appl. 32, 722–747 (2011)
Chu, D., Lathauwer, L.D., Moor, B.D.: A QR-type reduction for computing the SVD of a general matrix product/quotient. Numer. Math. 95, 101–121 (2003)
Demmel, J.: Accurate singular value decompositions of structured matrices. SIAM J. Matrix Anal. Appl. 21, 562–580 (1999)
Demmel, J., Gragg, W.: On computing accurate singular values and eigenvalues of acyclic matrices. Linear Algebra Appl. 185, 203–218 (1993)
Demmel, J., Ioana, D., Holtz, O., Koev, P.: Accurate and efficient expression evaluation and linear algebra. Acta Numer. 17, 87–145 (2008)
Demmel, J., Kahan, W.: Accurate singular values of bidiagonal matrices. SIAM J. Sci. Stat. Comput. 11, 873–912 (1990)
Dhilion, I.S., Parlett, B.N.: Orthogonal eigenvectors and relative gaps. SIAM J. Matrix Anal. Appl. 25, 858–899 (2004)
Dopico, F.M., Koev, P.: Accurate symmetric rank revealing and eigendecompositions of symmetric structured matrices. SIAM J. Matrix Anal. Appl. 28, 1126–1156 (2006)
Dopico, F.M., Koev, P., Molera, J.M.: Implicit standard Jacobi Givens high relative accuracy. Numer. Math. 113, 519–553 (2009)
Dopico, F.M., Molera, J.M., Moro, J.: An orthogonal high relative accuracy algorithm for the symmetric eigenproblem. SIAM J. Matrix Anal. Appl. 25, 301–351 (2003)
Fernando, K., Parlett, B.: Accurate signular values and differential qd algorithms. Numer. Math. 67, 191–229 (1994)
Golub, G.H., Sølna, K., Dooren, P.V.: Computing the SVD of a general matrix product/quotient. SIAM J. Matrix Anal. Appl. 22, 1–19 (2000)
Heath, M.T., Laub, A.J., Paige, C.C., Ward, R.C.: Computing the singular value decomposition of a product of two matrices. SIAM J. Sci. Stat. Comput. 7, 1147–1159 (1986)
Huang, R.: Rank structure properties of rectangular matrices admitting bidiagonal-type factorizations. Linear Algebra Appl. 465, 1–14 (2015)
Huang, R., Chu, D.L.: Relative perturbation analysis for eigenvalues and singular values of totally nonpositive matrices. SIAM J. Matrix Anal. Appl. 36, 476–495 (2015)
Kaufman, L.: The LZ-algorithm to solve the generalized eigenvalue problem. SIAM J. Numer. Anal. 11, 997–1024 (1974)
Koev, P.: Accurate eigenvalues and SVDs of totally nonnegative matrices. SIAM J. Matrix Anal. Appl. 27, 1–23 (2005)
Koev, P., Dopico, F.M.: Accurate eigenvalues of certain sign regular matrices. Linear Algebra Appl. 424, 435–447 (2007)
Moler, C.B., Stewart, G.W.: An algorithm for generalized matrix eigenvalue problems. SIAM J. Numer. Anal. 10, 241–256 (1973)
Peláez, M.J., Moro, J.: Accurate factorization and eigenvalues algorithms for symmetric DSTU and TSC matrices. SIAM J. Matrix Anal. Appl. 28, 1173–1198 (2006)
Pan, V.Y.: Structured Matrices and Polynomials: Unified Superfast Algorithms. Birkhäuser Boston Inc., Boston, MA (2001)
Pan, V.Y.: How bad are Vandermonde matrices. SIAM J. Matrix Anal. Appl. 27, 676–694 (2016)
Parlett, B.: The new qd algorithms. Acta Numer. 17, 459–491 (1995)
Phillips, G.M.: Interpolation and Approximation by Polynomials. Springer, New York (2003)
Press, W.H., Teukolsky, S.A., Vetterling, W.T., Flannery, B.P.: Numerical Recipes: The Art of Scientific Computing, 3rd edn. Cambridge University Press, Cambridge (2007)
Stewart, D.: A new algorithm for the SVD of a long product of matrices and the stability of products. Electron. Trans. Numer. Anal. 5, 29–47 (1997)
Watkins, D.S.: Product eigenvalue problems. SIAM Rev. 47, 3–40 (2005)
Watkins, D.S.: The Matrix Eigenvalue Problem: GR and KrylovSubsp ace Methods. SIAM, Philadelphia (2007)
Acknowledgements
The author would like to thank the Associate Editor and the anonymous referees for their valuable comments and suggestions which have helped to improve the overall presentation of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
Research supported by the National Natural Science Foundation of China (Grant No. 11471279), the Natural Science Foundation for Distinguished Young Scholars of Hunan Province (Grant No. 2017JJ1025) and the Research Foundation of Education Bureau of Hunan Province (Grant No. 14B178).
Rights and permissions
About this article
Cite this article
Huang, R. A Periodic qd-Type Reduction for Computing Eigenvalues of Structured Matrix Products to High Relative Accuracy. J Sci Comput 75, 1229–1261 (2018). https://doi.org/10.1007/s10915-017-0584-7
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10915-017-0584-7