from scipy.sparse import identity
from scipy.sparse.linalg import eigs
A = identity(10, format='csc')
A.setdiag(list(range(1, 11)))
eigs(A, 3, sigma=0)
Вычисление N наименьших собственных значений разреженной матрицы в Python