[Numpy-discussion] Numeric eigenvectors problem on Linux

Alexandre Gillet gillet at scripps.edu
Tue Feb 3 14:44:00 CST 2004


Hi,
I would like to know if anybody as experience any problem in their 
result using eigenvector function on Linux

The same calculation return the correct value on Windows,SUN,SGI.
All the Numeric were built without out modifying the setup.py file. I 
just run python setup.py install  (python2.3, Numeric 23.0)

The test script we run is the following:
import Numeric
import LinearAlgebra
e = Numeric.array([[ 0.6875,-0.3125,0.3125],
                    [-0.3125 , 0.6875 , 0.3125],
                    [ 0.3125 , 0.3125 , 0.6875]])
evals, evecs = LinearAlgebra.eigenvectors(e)
print "----------------------------------"
print "evecs",evecs


on Linux we  get:
evals [ 1.      0.0625  1.    ]
evecs [[ 0.81649658 -0.40824829  0.40824829]
	 [-0.57735027 -0.57735027  0.57735027]
  	[-0.26840956  0.80201256  0.533603  ]]

On all others platforms we get:
evals [ 1.      0.0625  1.    ]
evecs [[ 0.81649658 -0.40824829  0.40824829]
  	[-0.57735027 -0.57735027  0.57735027]
	 [-0.22645541  0.79259392  0.56613852]]


Any idea? Is it a bug in Numeric?

Thanks
Alex


-- 
   o Alexandre Gillet               email: gillet at scripps.edu
  /  The Scripps Research Institute,
o   Dept. Molecular Biology,  MB-5,
  \  10550  North Torrey Pines Road,
   o La Jolla,  CA 92037-1000,  USA.
  /  tel: (858) 784-2053
o   fax: (858) 784-2860








More information about the Numpy-discussion mailing list