Menemukan vektor eigen terkecil dari matriks jarang besar, lebih dari 100x lebih lambat di SciPy daripada di Octave
Saya mencoba untuk menghitung beberapa (5-500) vektor eigen yang sesuai dengan nilai eigen terkecil dari matriks persegi jarang simetris (hingga 30000x30000) dengan kurang dari 0,1% dari nilai yang bukan nol. Saat ini saya menggunakan scipy.sparse.linalg.eigsh dalam mode shift-invert (sigma =...