MathGroup Archive 2009

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

Search the Archive

Re: Eigenvalues broken?

  • To: mathgroup at
  • Subject: [mg102837] Re: [mg102827] Eigenvalues broken?
  • From: Vadim Oganesyan <vadim.oganesyan at>
  • Date: Sat, 29 Aug 2009 06:30:35 -0400 (EDT)
  • References: <> <h50sjj$6pl$>

Windows XP Pro x64

Interestingly, 32 bit laptop executes this code flawlessly.

On Fri, Aug 28, 2009 at 9:47 AM, Sseziwa Mukasa <mukasa at> wrote:

> 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.
> Regards,
> Sseziwa

  • Prev by Date: Re: Eigenvalues broken?
  • Next by Date: Re: simplex representation
  • Previous by thread: Re: Eigenvalues broken?
  • Next by thread: Re: Re: Eigenvalues broken?