MathGroup Archive 2009

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: Eigenvalues broken?

  • To: mathgroup at
  • Subject: [mg102836] Re: [mg102827] Eigenvalues broken?
  • From: Sseziwa Mukasa <mukasa at>
  • Date: Sat, 29 Aug 2009 06:30:24 -0400 (EDT)
  • References: <> <h50sjj$6pl$> <>

On Aug 28, 2009, at 12:47 AM, vadim.oganesyan at 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.



  • Prev by Date: Re: random variable
  • Next by Date: Re: Eigenvalues broken?
  • Previous by thread: Eigenvalues broken?
  • Next by thread: Re: Eigenvalues broken?