Я имею большую NxN плотную симметричную матрицу и хочу, чтобы собственные векторы соответствовали k наибольшим собственным значениям. Каков наилучший способ их найти (предпочтительно используя numpy, но, возможно, в общем случае, используя blas / atlas / lapack, если это единственный путь)? Обычно N намного больше, чем k (скажем, N & gt; 5000, k & lt; 10).

Кажется, что Numpy имеет функции для нахождения k наибольших собственных значений, если моя стартовая матрица разрежена.