Re: [Numpy-discussion] One question about the numpy.linalg.eig() routine

2012-04-03 Thread Hongbin Zhang
Dears, Though it might sounds strange, but the eigenvectors of my 2X2 matrix is rather different if I get it calculated in a loop over many other similar matrices: for instance: matrix: [[ 0.6000+0.j -1.97537668-0.09386068j] [-1.97537668+0.09386068j -0.6000+0.j]] eigen

Re: [Numpy-discussion] One question about the numpy.linalg.eig() routine

2012-04-03 Thread Hongbin Zhang
Hej Val, Thank you very much for your replies. Yes, I know that both eigenvectors are correct while they are indeed related to each other by unitary transformations (unitary matrices). Actually, what I am trying to do is to evaluate the Berry phase which is closely related to the gauge chosen

[Numpy-discussion] One question about the numpy.linalg.eig() routine

2012-04-02 Thread Hongbin Zhang
Dear Python-users, I am currently very confused about the Scipy routine to obtain the eigenvectors of a complex matrix.In attached you find two files to diagonalize a 2X2 complex Hermitian matrix, however, on my computer, If I run python, I got: [[ 0.80322132+0.j 0.59500941+0.02827207j]