Adjacency uint8 редкие собственные значения

Question

У меня есть исходная и целевая матрица, которую я создаю матрицей смежности, используя разреженную функцию.

W = sparse(w1(:,1),w1(:,2),1,n,n); 

my n приблизительно равен 540 000 (это подмножество данных, которые у меня есть). Создание самой матрицы довольно легко вычислительно. Однако мне нужно получить собственные значения указанной матрицы. Точнее

rho = max(abs(eigs(W))); 

Эта одна строка кода берет навсегда. Немного поучав в этом, мне показалось, что получение собственных значений разреженной матрицы является дорогостоящим вычислительным. Я не могу использовать полную команду, например, uint8 (full (W)), так как в итоге я превысил максимальный размер матрицы/емкость (130 ГБ памяти). Однако, похоже, что если формат был uint8, я мог бы создать указанную полную матрицу смежности.

Любая помощь приветствуется. Мне кажется, что мне нужно одно из следующего:

• Альтернативный способ построения матрицы смежности непосредственно в виде полного uint8;
• Способ преобразования матрицы смежности, которую я создал для полного uint8; или
• Альтернативный способ расчета этих собственных

Answer 1

С W = sparse(w1(:,1),w1(:,2),1,n,n); вы создаете матрицу, которая держит только 0 или 1. В то время как тип данных uint8 не представляется возможным, вы можете создать sparselogical матрицу вместо этого. Таким образом, можно избежать использования типов данных 64-битных

W = sparse(w1(:,1),w1(:,2),true,n,n); 
W = uint8(full(W)); 

Answer 2

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

rho = abs(eigs(W,1));