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

Question

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

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

Answer 1

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