rcpp - Eigen::EigenSolver 是否适用于 SparseMatrix 类的对象？

Question

Eigen::EigenSolver对类对象起作用吗SparseMatrix？特别是，我正在使用 RcppEigen。

其次，使用哪种算法？它是否与 EISPACK、LAPACK、Wilkinson (1965) 等中的 QR 算法相同？

Answer 1

1. EigenSolver 仅适用于密集矩阵，不适用于稀疏矩阵。
2. 是的，EigenSolver 使用 QR 算法。它与 EISPACK 中的非常相似；文档说：该实现改编自 JAMA（公共领域）。他们的代码基于 EISPACK。我相信 LAPACK 使用了一个略有不同的 QR 算法变体，具有多个班次。