[Numpy-discussion] Numeric eigenvectors problem on Linux

[Nadav Horesh]

On a linux system (RH9) with python 2.3.3 and Numeric 23.1, I get the right results. Maybe you should just to upgrade to Numeric23.1.

  Nadav.

-----Original Message-----
From:	Alexandre Gillet [mailto:gillet at scripps.edu]
Sent:	Wed 04-Feb-04 00:43
To:	numpy-discussion at lists.sourceforge.net
Cc:	Michel Sanner
Subject:	[Numpy-discussion] Numeric eigenvectors problem on Linux
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






-------------------------------------------------------
The SF.Net email is sponsored by EclipseCon 2004
Premiere Conference on Open Tools Development and Integration
See the breadth of Eclipse activity. February 3-5 in Anaheim, CA.
http://www.eclipsecon.org/osdn
_______________________________________________
Numpy-discussion mailing list
Numpy-discussion at lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/numpy-discussion