Re: [Bug-gsl] Eigen vector problem

[Patrick Alken]

Can you provide a simple program which demonstrates the problem? You 
didn't even list your A and B matrices

On 05/30/2014 01:47 AM, belwafi kais wrote:
> Hello,
>
>
>
> I try to implement the eigenvector using GSL library. To do this I use
> gsl_eigen_genv function. And I try to compare the output of this function
> with the output of Matbal function ( [c d]=eig(A,B), A and B are identical
> matrix, all of these element is one; A=ones(4,4,)). I get two different
> results:
>
>
>
>
>
> Eigenvector
>
> Eigen value
>
> Matlab result
>
> c =
>
>      -1     1     0     0
>
>       0    -1     0     0
>
>       0     0     1     0
>
>       0     0     0     1
>
>
>
> d =
>
>       1     0     0     0
>
>       0     0     0     0
>
>       0     0   NaN     0
>
>       0     0     0   NaN
>
>
>
> GSL result
>
> c =
>
>   -1    -0    -0    -0
>
>   0.7  -1.1 -0.7  0
>
>   0      0     1      0
>
>   0      0     0      1
>
>
>
> d =
>
>       1     0     0       0
>
>       0     inf   0       0
>
>       0      0   -NaN    0
>
>       0     0     0  -NaN
>
>
>
>
>
> Please can you give me any tips to correct this problem?
>
>
>
> Thank you.