[Numpy-discussion] Problem with numpy.linalg.eig?
David Cournapeau
david@ar.media.kyoto-u.ac...
Mon Nov 12 22:26:56 CST 2007
Geoffrey Zhu wrote:
> On Nov 12, 2007 12:37 PM, Keith Goodman <kwgoodman@gmail.com> wrote:
>> On Nov 12, 2007 10:10 AM, Peter Creasey <p.e.creasey.00@googlemail.com> wrote:
>>> The following code calling numpy v1.0.4 fails to terminate on my machine,
>>> which was not the case with v1.0.3.1
>>>
>>> from numpy import arange, float64
>>> from numpy.linalg import eig
>>> a = arange(13*13, dtype = float64)
>>> a.shape = (13,13)
>>> a = a%17
>>> eig(a)
>> It sounds like the same problem that was reported in this thread:
>>
>> http://thread.gmane.org/gmane.comp.python.numeric.general/17456/focus=17465
>
> The code hangs on my machine too. In the thread you mentioned above, I
> wrote that using the EGG instead of MSI appears to fix the
> numpy.test() problem, but maybe it just somehow hides it.
When you use the MSI, can you always reproduce the problem ? As I said
previously, it is hard to know for sure without being able to reproduce
the bug on our own workstation, but if this is a problem between fortran
and C argument passing, then the result can be pretty random since the
problem reduced to a pointer pointing at a wrong address (crash, wrong
value, etc...).
cheers,
David
More information about the Numpy-discussion
mailing list