Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: Eigenvalues works very slow
Replies: 3   Last Post: Oct 25, 2012 11:35 PM

 Messages: [ Previous | Next ]
 Jure Posts: 10 Registered: 9/25/12
Eigenvalues works very slow
Posted: Oct 24, 2012 3:30 AM

Hello,
I'm trying to calculate eigenvalues of different sized matrices (from 10x10 to 1000x1000) using mathematica's built-in function - Eigenvalues. For smaller matrices it works ok, but for larger matrices it's just too slow. I tried writing the code in another system and it finds eigenvalues very quickly (in a second or two at most) even for big matrices. I'm not sure whether am I doing something wrong or is the speed difference really this big.
My code (you only need to change d to change the size of matrix):
d = 25;
q = Table[Table[1/2*Sqrt[i + j + 1]*KroneckerDelta[Abs[i - j], 1], {j, 0, d-1}],{i, 0, d - 1}];
h0 = Table[Table[(i + 1/2)*KroneckerDelta[i, j], {j, 0, d - 1}], {i, 0,d-1}];
lambda = 1/2;
q4 = lambda*q.q.q.q;
N[Eigenvalues[q4 + h0]]

Thanks for help.

Date Subject Author
10/24/12 Jure
10/25/12 David Bailey
10/25/12 Bill Rowe
10/25/12 Michael Weyrauch