MATLAB: Optimizing calculation of eigenvectors and eigenvalues

eigeigenvalueeigenvectoroptimizationoptimize

I have a quadratic matrix A with a size of about 2000×2000. I want to calculate its eigenvectors and eigenvalues. The eigenvalues must be in a vector (not a diagonal matrix as usual). I use the following code:
eigenvals = abs(real(eig(A)));
[eigenvecs, ~] = eig(A);
eigenvecs = real(eigenvecs);
The problem is that it takes a lot of time (about 15 seconds) and I need to repeat this process a lot of times in a loop. So I tried to optimize it and changed it to this:
[eigenvecs, eigenvals]=eig(A);
eigenvecs=real(eigenvecs);
eigenvals=abs(real(diag(eigenvals)));
This runs about 5 seconds faster. The calculated eigenvectors are the same for both code snippets, however the eigenvalues are not and differ a bit. How could I further optimize the code in a way that it delivers equivalent results to the first code snippet?
Thank you very much in advance!

Best Answer

Hi Guiseppi,
there is nothing guaranteed about the order of the eigenvalues that are produced by eig. And if you do the first method, which is basically
lambda = eig(A)
% and
[v, ~] = eig(A)
the order of eigenvalues may not be the same in each case, meaning that the order of eigenvalues and the order of rows of the eigenvector matrix may not match up. In fact for some examples I tried, they do not match up, which is similar to your experience. This means the first method does not work.
If you would like a vector rather than a diagonal matrix for lambda, there is the option
[v lambda] = eig(A,'vector')
which only calls eig once and so is faster than the first method, in addition to being correct.
Related Question