Fast numerical computation of eigenvalues

asked 2018-12-01 12:10:10 -0500

SageMathematician gravatar image

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

Comments

I'm having a similar issue with Sage notebooks on my AMD machine! Maybe OpenBLAS is being forced to run in single-threaded mode? I wish I knew how to change this...

droberts01 gravatar imagedroberts01 ( 2020-09-02 00:17:54 -0500 )edit