Mailman 3 December 2008 - NumPy-Discussion

Unexpected MaskedArray behavior
by Ryan May 17 Dec '08

17 Dec '08

Hi, I just noticed the following and I was kind of surprised: >>>a = ma.MaskedArray([1,2,3,4,5], mask=[False,True,True,False,False]) >>>b = a*5 >>>b masked_array(data = [5 -- -- 20 25], mask = [False True True False False], fill_value=999999) >>>b.data array([ 5, 10, 15, 20, 25]) I was expecting that the underlying data wouldn't get modified while masked. Is this actual behavior expected? Ryan -- Ryan May Graduate Research Assistant School of Meteorology University of Oklahoma

3 4

genloadtxt : last call
by Pierre GM 17 Dec '08

17 Dec '08

All, Here's the latest version of genloadtxt, with some recent corrections. With just a couple of tweaking, we end up with some decent speed: it's still slower than np.loadtxt, but only 15% so according to the test at the end of the package. And so, now what ? Should I put the module in numpy.lib.io ? Elsewhere ? Thx for any comment and suggestions.

5 10

Win64 build?
by Klaus Noekel 17 Dec '08

17 Dec '08

Dear all, I would like to use numpy under Windows Vista 64-bit, but I am scared a bit by compiling a 64bit build myself. Is there an installer for Win64 somewhere, or are there any plans for one? Thanks for any advice, including directions for building myself (sigh), if that should be the only way. Klaus Noekel

4 4

Re: [Numpy-discussion] Efficient removal of duplicates
by Daran Rife 16 Dec '08

16 Dec '08

How about a solution inspired by recipe 18.1 in the Python Cookbook, 2nd Ed: import numpy as np a = [(x0,y0), (x1,y1), ...] l = a.tolist() l.sort() unique = [x for i, x in enumerate(l) if not i or x != b[l-1]] a_unique = np.asarray(unique) Performance of this approach should be highly scalable. Daran -- Hi, I the following problem: I have a relatively long array of points [(x0,y0), (x1,y1), ...]. Apparently, I have some duplicate entries, which prevents the Delaunay triangulation algorithm from completing its task. Question, is there an efficent way, of getting rid of the duplicate entries? All I can think of involves loops. Thanks and regards, Hanno

7 10

Mersenne twister seeds
by Michael Gilbert 15 Dec '08

15 Dec '08

According to wikipedia [1], some common Mersenne twister algorithms use a linear congruential gradient (LCG) to generate seeds. LCGs have been known to produce poor random numbers. Does numpy's Mersenne twister do this? And if so, is this potentially a problem? http://en.wikipedia.org/wiki/Linear_congruential_generator

4 3

Concatenating Arrays to make Views
by Benjamin Haynor 15 Dec '08

15 Dec '08

Hi, I was wondering if I can concatenate 3 arrays, where the result will be a view of the original three arrays, instead of a copy of the data. For example, suppose I write the following import numpy as n a = n.array([[1,2],[3,4]]) b = n.array([[5,6],[7,8]]) c = n.array([[9,10],[11,12]]) c = n.r_[a,b] Now c = : [[1,2], [3,4], [5,6], [7,8], [9,10], [11,12]] I was hoping to get an array, such that, when I change d, a, b, and c will also change appropriately. Any ideas? - Ben bhaynor(a)hotmail.com _________________________________________________________________ Send e-mail anywhere. No map, no compass. http://windowslive.com/Explore/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_anywhe…

3 2

Efficient removal of duplicates
by Hanno Klemm 15 Dec '08

15 Dec '08

Hi, I the following problem: I have a relatively long array of points [(x0,y0), (x1,y1), ...]. Apparently, I have some duplicate entries, which prevents the Delaunay triangulation algorithm from completing its task. Question, is there an efficent way, of getting rid of the duplicate entries? All I can think of involves loops. Thanks and regards, Hanno -- Hanno Klemm klemm(a)phys.ethz.ch

4 3

ANNOUNCE: ETS 3.1.0 released!
by Dave Peterson 14 Dec '08

14 Dec '08

I'm pleased to announce that the Enthought Tool Suite (ETS) 3.1.0 has been tagged, released, and uploaded to PyPi[1]! Both source distributions (.tar.gz) and binary (.egg) for Windows have been built and uploaded to PyPi. You can update an existing ETS install to v3.1.0 like so: easy_install -U ETS==3.1.0 What is ETS? ------------------ The Enthought Tool Suite (ETS) is a collection of projects developed by members of the OSS community, including Enthought employees, which we use every day to construct custom scientific applications. It includes a wide variety of components, including: * an extensible application framework * application building blocks * 2-D and 3-D graphics libraries * scientific and math libraries * developer tools The cornerstone on which these tools rest is the Traits project, which provides explicit type declarations in Python; its features include initialization, validation, delegation, notification, and visualization of typed attributes. More information is available for all these packages from the Enthought Tool Suite development home page: http://code.enthought.com/projects/index.php -- Dave

4 5

Plot directive in numpy docs
by Gael Varoquaux 13 Dec '08

13 Dec '08

Hi, What is the guideline on using the plot directive in the numpy docs? It can make some examples much easier to understand, but on the other hand, it can clutter the docstrings. In addition, I am not sure how our documentation pipeline deals with it. Cheers, Gaël

3 2

bugfix for np.random logseries and hypergeometric are verified
by josef.pktd＠gmail.com 12 Dec '08

12 Dec '08

Hi, Now that I managed to compile numpy, I tried out the bugfixes in ticket:921 and ticket:923. For both, checking the results gives now correct results. I attached the scripts and results to the tickets. Both distributions will be tested again in scipy.stats, once I remove the skip and add the arguments that previously resulted in wrong numbers. Can someone committ these, or can I commit them myself? (This will be the maximum that I touch C, changing one word and one inequality) Josef

2 1