[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

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

reply via email to

[Prev in Thread] Current Thread [Next in Thread]