RE: [Numpy-discussion] Numeric eigenvectors problem on Linux

[Nadav H. <na...@vi...>]

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:gi...@sc...]
Sent:	Wed 04-Feb-04 00:43
To:	num...@li...
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=20
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=20
just run python setup.py install  (python2.3, Numeric 23.0)

The test script we run is the following:
import Numeric
import LinearAlgebra
e =3D Numeric.array([[ 0.6875,-0.3125,0.3125],
                    [-0.3125 , 0.6875 , 0.3125],
                    [ 0.3125 , 0.3125 , 0.6875]])
evals, evecs =3D 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


--=20
   o Alexandre Gillet               email: gi...@sc...
  /  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
Num...@li...
https://lists.sourceforge.net/lists/listinfo/numpy-discussion