We discuss some potential advantages of the orthogonal symmetric-diagonal reduction in two main versions of the Schur-QR method for symmetric positive definite generalized eigenvalue problems. We also advise and use the appropriate reductions as preprocessing on the solvers, mainly the Cholesky-QR method, of the considered problems. We discuss numerical stability of the methods via providing upper bound for backward error of the computed eigenpairs and via investigating two kinds of scaled residual errors. We also propose and apply two kinds of symmetrizing which improve the stability and the performance of the methods. Numerical experiments show that the implemented versions of the Schur-QR method and the preprocessed versions of the Cholesky-QR method are usually more stable than the Cholesky-QR method.
Ahmadnasab, M. (2023). Symmetric-diagonal reductions as preprocessing for symmetric positive definite generalized eigenvalue solvers. Journal of Mathematical Modeling, 11(2), 301-322. doi: 10.22124/jmm.2023.23734.2120
MLA
Morad Ahmadnasab. "Symmetric-diagonal reductions as preprocessing for symmetric positive definite generalized eigenvalue solvers". Journal of Mathematical Modeling, 11, 2, 2023, 301-322. doi: 10.22124/jmm.2023.23734.2120
HARVARD
Ahmadnasab, M. (2023). 'Symmetric-diagonal reductions as preprocessing for symmetric positive definite generalized eigenvalue solvers', Journal of Mathematical Modeling, 11(2), pp. 301-322. doi: 10.22124/jmm.2023.23734.2120
VANCOUVER
Ahmadnasab, M. Symmetric-diagonal reductions as preprocessing for symmetric positive definite generalized eigenvalue solvers. Journal of Mathematical Modeling, 2023; 11(2): 301-322. doi: 10.22124/jmm.2023.23734.2120