Hi, 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. Regards Bruce ---- Original message ----

Date: Thu, 3 Jun 2004 14:18:56 +0300 From: Janne Sinkkonen Subject: [Numpy-discussion] RandomArray.random() To: numpy-discussion@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 variates.

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.

-- Janne

------------------------------------------------------- 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@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion