I know it is late; but I came across a similar issue. You could use the symmetry of the matrix, something like: from scipy.sparse.linalg import eigsh and eigsh(A, k=1, which='LM', return_eigenvectors=False). np.linalg.norm(A, ord=2) is painfully slow.