[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Bug-gsl] bug in 'gsl_eigen_symmv'
From: |
Brian Gough |
Subject: |
Re: [Bug-gsl] bug in 'gsl_eigen_symmv' |
Date: |
Mon, 15 Aug 2011 16:37:34 +0100 |
User-agent: |
Wanderlust/2.15.6 (Almost Unreal) Emacs/23.3 Mule/6.0 (HANACHIRUSATO) |
At Sat, 6 Aug 2011 01:08:06 +0200,
Sazonov, Vasily (address@hidden) wrote:
> I suppose, I found a bug in 'gsl_eigen_symmv' function. I have examined it on
> matrix:
>
> 7.255197 0.000000
> -2.666667 17.177062
>
I noticed that your matrix is not symmetrical, so gsl_eigen_symmv is
not going to give the result you expect - it assumes a symmetrical
matrix.
The lower triangular part is used, so you are actually computing the
eigenvalues of [7.255197, -2.666667; -2.666667, 17.177062], and the
correct results for that are 6.5839, 17.8484 as given.
--
Brian Gough