Eigenvalue problem i GSL?

Gerard Jungman jungman@lanl.gov
Wed Mar 29 13:36:00 GMT 2000


Øystein O Johansen wrote:

> I'm trying to calculate the eigenvalues from this matrix:
> 
>    | -3.34   1.596 |
> a =|               |
>    | -1.596  3.512 |
> 
> and if I haven't done anything wrong, gsl_eigen_jacobi returns:
> 
> -3.69351 and 3.86551
> 
> The prblem is that this is the wrong answer, according to Matlab, and
> my hp-48 calculator, which both gets:    -0.2946 and 0.3118

The Jacobi method is appropriate for symmetric matrices, which
yours is not. I have no idea what it tells you if the matrix
is not symmetric.

I see that there is nothing yet in the documentation to
indicate this. Perhaps soon there will be...


--
G. Jungman


More information about the Gsl-discuss mailing list