find eigenvalues of a very large sparse matrix

57 visualizaciones (últimos 30 días)
Zeyu Zhang
Zeyu Zhang el 11 de Sept. de 2018
Respondida: Andrew Knyazev el 21 de Sept. de 2018
I need to find the eigenvalues of a very large sparse matrix (5digit x 5digit dimension). I tried to use 'eig' command but apparently, I don't have enough memory in my RAM ( the return message is always run of out memory). But it is worth to mention that I did this with all the zero entries filled with zeros. This is because when the zero entries are squeezed out (in sparse matrix form), it took forever to run the 'eig' command due to unknown reason so I have to give up. I also tried to do it on GPU, but it was the same situation (run out of memory) when zeros are squeezed out, and GPU can't handle matrix in parse form so the other option is not available. I wonder if it is possible to partition the large matrix into blocks of small segments, do some computation separately first then combine the results from individual parts and somehow obtain the final result from here. Is that a mathematical viable approach? If yes how to do it then? Thanks.

Respuestas (2)

M
M el 11 de Sept. de 2018
  2 comentarios
Zeyu Zhang
Zeyu Zhang el 19 de Sept. de 2018
Editada: Zeyu Zhang el 19 de Sept. de 2018
Hi, I checked eigs and it seems to work. However for example d = eigs(Afun,1500,6,'smallestabs') returns the eigenvalues with smallest magnitude. What if I want it to return the most negative eigenvalues?
M
M el 19 de Sept. de 2018
Are all eigenvalues real? If yes you can use 'smallestreal' option.

Iniciar sesión para comentar.


Andrew Knyazev
Andrew Knyazev el 21 de Sept. de 2018
If the matrix is real symmetric or Hermitian, you may also want to try https://www.mathworks.com/matlabcentral/fileexchange/48-lobpcg-m

Categorías

Más información sobre Linear Algebra en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by