3
Виявлення найменших власних векторів великої розрідженої матриці, на 100 разів повільніше в SciPy, ніж в Октаві
Я намагаюся обчислити декілька (5-500) власних векторів, що відповідають найменшим власним значенням великих симетричних квадратних розріджених матриць (до 30000x30000), причому менше 0,1% значень є ненульовими. В даний час я використовую scipy.sparse.linalg.eigsh у режимі зсуву-інвертування (сигма = 0,0), що я з'ясував, використовуючи різні пости в темі, є кращим рішенням. Однак для …