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.
