Собственное значение Probem

Question

Я пытаюсь диагоналировать матрицу в джулии, но сама матрица слишком велика для расчета. Но я могу эффективно рассчитать матричное векторное произведение. Я ищу решение, чтобы я мог использовать матричное векторное произведение в eigs и решать с помощью метода lanzcos. Есть ли пакет, который дает мне вывод в виде матричного векторного произведения?

Может кто-нибудь помочь здесь?

Answer 1

Вы можете использовать ARPACK на основе eigs функция, которая включена в Julia 0.6.2 или вы можете использовать чистую реализацию Julia связанного метода в https://github.com/JuliaMath/IterativeSolvers.jl. Последний не основан на Lanczos, но я думаю, что это не имеет значения, пока он дает вам собственные значения.