Работает ли Eigen :: EigenSolver над объектом класса SparseMatrix?

Question

Работает ли Eigen::EigenSolver на объекте класса SparseMatrix? В частности, я использую RcppEigen.

Во-вторых, какой алгоритм используется? Это тот же QR-алгоритм, что и в EISPACK, LAPACK, Wilkinson (1965) ... и т. Д.

Answer 1

1. EigenSolver работает только с плотными матрицами, а не с разреженными матрицами.
2. Да, EigenSolver использует алгоритм QR. Это очень похоже на то, что в EISPACK; документация гласит: Реализация адаптирована из JAMA (общественное достояние). Их код основан на EISPACK. Я считаю, что LAPACK использует несколько иной вариант алгоритма QR с несколькими сменами.