southey at uiuc.edu
Wed Jun 9 08:06:11 CDT 2004
I would be prepared to at least look at it with numarray if you are prepared to
share the scripts or rather the relevant part of them. On the surface, it would
appear that this is a 32 vs 64 bit problem in defining the type precision.
---- Original message ----
>Date: Thu, 3 Jun 2004 14:18:56 +0300
>From: Janne Sinkkonen <janne.sinkkonen at hut.fi>
>Subject: [Numpy-discussion] RandomArray.random()
>To: numpy-discussion at lists.sourceforge.net
>Related to a recent discussion on random numbers:
>People using RandomArray.random() on 64-bit architectures should also be aware
>that sometimes the function returns exactly zero or one, which in turn causes
>problems in many other routines which generate e.g. multinomial or dirichlet
>I'm not absolutely sure that the problem persists in the newest version of
>Numeric, and I have not tested it on numarray. Anyway, I have seen it both in
>Alphas and now later in an AMD Opteron machine - over the years.
>Unfortunately, as we do not have had time to dig any deeper around here, we
>have just used wrapper with a test and a loop.
>This SF.Net email is sponsored by the new InstallShield X.
>From Windows to Linux, servers to mobile, InstallShield X is the one
>installation-authoring solution that does it all. Learn more and
>evaluate today! http://www.installshield.com/Dev2Dev/0504
>Numpy-discussion mailing list
>Numpy-discussion at lists.sourceforge.net
More information about the Numpy-discussion