[Date Index] [Thread Index] [Author Index]
Re: Eigenvalues broken?
On Aug 28, 2009, at 12:47 AM, vadim.oganesyan at gmail.com wrote: > I have been doing full diagonalization of sparse symmetric real > matrices (Hamiltonians) of different sizes and just noticed that > Eigenvalues (and Eigensystem) break once the matrices are 256by256 > (maybe smaller, but that's where i first saw it). > > i first noticed the problem when checking the quality of the solutions > by comparing eigenvalues with expectation of the matrix in each of > eigenvectors. > > even more dramatically, simply rerunning Eigenvalues twice on the same > matrix produces somewhat different results!!! > > If I chose to get only few eigenvales (thus forcing it to use Arnoldi > method) the problem goes away. > > I was recently told that the version of Lapack used by Mathemica was > recompiled by people at Wolfram for added functionality -- did they > break it too? > My trusty C++ code has no problems with these matrices.... > > thoughts? suggestions? What platform are you using? I've had problems with the numeric linear algebra routines on OS X, the same code seems to work on Windows. Regards, Sseziwa