Localizando os menores autovetores de matriz esparsa grande, 100x mais lenta no SciPy do que no Octave
Estou tentando calcular alguns (5-500) autovetores correspondentes aos menores autovalores de matrizes esparsas quadradas simétricas grandes (até 30000x30000) com menos de 0,1% dos valores diferentes de zero. Atualmente, estou usando scipy.sparse.linalg.eigsh no modo shift-invert (sigma = 0,0),...