>>> numpy.linalg.eig(numpy.random.rand(3,3))
You'll (almost always) get a wrong answer back from numpy. Something like:
(array([ 1.72167898, -0.07251007, -0.07251007]),
array([[ 0.47908847, 0.72095163, 0.72095163],
[ 0.56659142, -0.46403504, -0.46403504],
[ 0.67040914, 0.01361572, 0.01361572]]))
The return value should be complex (unless rand() just happens to return something symmetric).
It really needs to either throw an exception, or preferably for this function, just go ahead and return something complex, like the numpy.dft functions do.
On the other hand it, would be nice to stick with plain doubles if the output isn't complex, but I'd rather get the right output all the time than get the minimal type that will handle the output.
This is with beta 1.
Incidentally, I tried logging into the Trac here: http://projects.scipy.org/scipy/scipy
to file a bug, but it wouldn't let me in under the account I've been using for a while now. Is the login system broken? Were passwords reset or something?
--bb
------------------------------------------------------------------------- Using Tomcat but need to do more? Need to support web services, security? Get stuff done quickly with pre-integrated technology to make your job easier Download IBM WebSphere Application Server v.1.0.1 based on Apache Geronimo http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&dat=121642
_______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion
