eigenvectors from numpy.eig not orthogonal

Question

my problem is the following: using scipy.linalg.eig to get eigenvectors and eigenvalues i see that all my eigenvalues have multiplicity 1 yet when i run the code below it doesn\

Answer 1

Why should they be orthogonal? Your matrix

a=matM(60000)

is far from being symmetric,

abs(a-a.T).max() -> 2.16

with

abs(a).max() -> 1.08

so I wouldn't necessarily expect orthogonal eigenvectors. Is it possibile that the function matM or the data matA or matB is wrong?