How does one diagonalize a large sparse symmetric matrix to get the eigenvalues and the eigenvectors?
The problem is the matrix could be very large (though it is sparse), at most $2500\times 2500$. Is there a good algorithm to do that and most importantly, one that I can implement it into my own code? Thanks a lot!