Fast numerical computation of eigenvalues

asked 2018-12-01 12:10:10 -0600

I am trying to calculate the eigenvalues of a large (say, n=1000 to 10000) sparse Hermitian square matrix using SAGE

numpy.linalg.eigvalsh(MyMatrix)

takes very long. I noticed it utilizes only a single core of my CPU.

How would one go about speeding the calculation? Specifically, I'm looking for a solution using parallel computation, or maybe something which is "more compiled".

Thank you.

edit retag flag offensive close merge delete