From: Nils W. <nw...@ia...> - 2006-08-18 07:23:33
|
Bill Baxter wrote: > If you do this: > >>> 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 > Num...@li... > https://lists.sourceforge.net/lists/listinfo/numpy-discussion > AFAIK this problem is fixed. http://projects.scipy.org/scipy/numpy/ticket/215 I have no problem wrt the Trac system. Nils |