Re: [Numpy-discussion] in the NA discussion, what can we agree on?

2011-11-04 Thread Benjamin Root
On Fri, Nov 4, 2011 at 1:59 PM, Pauli Virtanen p...@iki.fi wrote: For shorthand, we can refer to the above choices with the nomenclature shorthand ::= propagation destructivity payload_type propagation ::= P | N destructivity ::= d | n | s payload_type ::= S | E | C I really

Re: [Numpy-discussion] in the NA discussion, what can we agree on?

2011-11-04 Thread Benjamin Root
On Fri, Nov 4, 2011 at 10:33 PM, T J tjhn...@gmail.com wrote: On Fri, Nov 4, 2011 at 8:03 PM, Nathaniel Smith n...@pobox.com wrote: On Fri, Nov 4, 2011 at 7:43 PM, T J tjhn...@gmail.com wrote: On Fri, Nov 4, 2011 at 6:31 PM, Pauli Virtanen p...@iki.fi wrote: An acid test for proposed

Re: [Numpy-discussion] what is the point of dx for np.gradient()?

2011-11-05 Thread Benjamin Root
On Fri, Nov 4, 2011 at 1:20 PM, Benjamin Root ben.r...@ou.edu wrote: For np.gradient(), one can specify a sample distance for each axis to apply to the gradient. But, all this does is just divides the gradient by the sample distance. I could easily do that myself with the output from

Re: [Numpy-discussion] Reading automatically all the parameters from a file

2011-11-30 Thread Benjamin Root
On Wednesday, November 30, 2011, Robert Kern robert.k...@gmail.com wrote: On Wed, Nov 30, 2011 at 11:09, Giovanni Plantageneto g.plantagen...@gmail.com wrote: Dear all, I have a simple question. I would like to have all the parameters of a model written in a configuration file (text), and I

Re: [Numpy-discussion] numpy.array() of mixed integers and strings can truncate data

2011-12-01 Thread Benjamin Root
On Thursday, December 1, 2011, Thouis Jones thouis.jo...@curie.fr wrote: On Thu, Dec 1, 2011 at 15:47, Pierre Haessig pierre.haes...@crans.org wrote: Le 01/12/2011 14:52, Thouis (Ray) Jones a écrit : Is this expected behavior? np.array([-345,4,2,'ABC']) array(['-34', '4', '2', 'ABC'],

Re: [Numpy-discussion] build numpy matrix out of smaller matrix

2011-12-01 Thread Benjamin Root
On Thu, Dec 1, 2011 at 10:52 AM, jonasr jonas.rueb...@web.de wrote: Hi, is there any possibility to define a numpy matrix, via a smaller given matrix, i.e. in matlab i can do this like a=[1 2 ; 3 4 ] A=[a a ; a a ] so that i finally get A=[ [1,2,1,2] [3,4,3,4] [1,2,1,2]

Re: [Numpy-discussion] Fast Reading of ASCII files

2011-12-14 Thread Benjamin Root
On Wed, Dec 14, 2011 at 1:22 PM, Ralf Gommers ralf.gomm...@googlemail.comwrote: On Wed, Dec 14, 2011 at 4:11 PM, Bruce Southey bsout...@gmail.com wrote: ** On 12/14/2011 01:03 AM, Chris Barker wrote: On Tue, Dec 13, 2011 at 1:21 PM, Ralf Gommers ralf.gomm...@googlemail.com wrote:

Re: [Numpy-discussion] Problem installing NumPy with Python 3.2.2/MacOS X 10.7.2

2011-12-16 Thread Benjamin Root
On Friday, December 16, 2011, McNicol, Adam amcni...@longroad.ac.uk wrote: Hi There, I am very new to numpy and have really only started investigating it as one of my students needs some functionality from matplotlib. I have managed to install everything under Windows for work in class but I

Re: [Numpy-discussion] find location of maximum values

2011-12-20 Thread Benjamin Root
On Tuesday, December 20, 2011, questions anon questions.a...@gmail.com wrote: ok thanks, a quick try at using it resulted in: IndexError: index out of bounds but I may need to do abit more investigating to understand how it works. thanks The assumption is that these arrays are all the same

Re: [Numpy-discussion] find location of maximum values

2012-01-09 Thread Benjamin Root
On Monday, January 9, 2012, questions anon questions.a...@gmail.com wrote: thanks for the responses. Unfortunately they are not matching shapes print TSFC.shape, TIME.shape, LAT.shape, LON.shape (721, 106, 193) (721,) (106,) (193,) So I still receive index out of bounds error:

Re: [Numpy-discussion] Calculating density based on distance

2012-01-14 Thread Benjamin Root
On Saturday, January 14, 2012, Thiago Franco de Moraes totonixs...@gmail.com wrote: Hi all, I have the following problem: Given a array with dimension Nx3, where N is generally greater than 1.000.000, for each item in this array I have to calculate its density, Where its density is the

Re: [Numpy-discussion] np.zeros(2, 'S') returns empty strings.

2012-01-14 Thread Benjamin Root
On Sat, Jan 14, 2012 at 4:12 PM, Charles R Harris charlesr.har...@gmail.com wrote: This sort of makes sense, but is it the 'correct' behavior? In [20]: zeros(2, 'S') Out[20]: array(['', ''], dtype='|S1') It might be more consistent to return '0' instead, as in In [3]: zeros(2,

Re: [Numpy-discussion] np.zeros(2, 'S') returns empty strings.

2012-01-14 Thread Benjamin Root
On Sat, Jan 14, 2012 at 4:16 PM, Benjamin Root ben.r...@ou.edu wrote: On Sat, Jan 14, 2012 at 4:12 PM, Charles R Harris charlesr.har...@gmail.com wrote: This sort of makes sense, but is it the 'correct' behavior? In [20]: zeros(2, 'S') Out[20]: array(['', ''], dtype='|S1

Re: [Numpy-discussion] Negative indexing.

2012-01-16 Thread Benjamin Root
On Sun, Jan 15, 2012 at 6:54 AM, David Cournapeau courn...@gmail.comwrote: On Sat, Jan 14, 2012 at 11:53 PM, Nathan Faggian nathan.fagg...@gmail.com wrote: Hi, I am finding it less than useful to have the negative index wrapping on nd-arrays. Here is a short example: import numpy as

Re: [Numpy-discussion] Negative indexing.

2012-01-16 Thread Benjamin Root
On Mon, Jan 16, 2012 at 3:24 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Jan 14, 2012 at 4:53 PM, Nathan Faggian nathan.fagg...@gmail.comwrote: Hi, I am finding it less than useful to have the negative index wrapping on nd-arrays. Here is a short example: import

Re: [Numpy-discussion] Negative indexing.

2012-01-16 Thread Benjamin Root
On Mon, Jan 16, 2012 at 3:30 PM, Benjamin Root ben.r...@ou.edu wrote: On Mon, Jan 16, 2012 at 3:24 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Jan 14, 2012 at 4:53 PM, Nathan Faggian nathan.fagg...@gmail.com wrote: Hi, I am finding it less than useful to have

Re: [Numpy-discussion] getting position index from array

2012-01-19 Thread Benjamin Root
On Thursday, January 19, 2012, Ruby Stevenson ruby...@gmail.com wrote: hi, all I am a newbie on numpy ... I am trying to figure out, given an array, how to get back position value based on some conditions. Say, array([1, 0, 0, 0 1], and I want to get a list of indices where it is none-zero,

Re: [Numpy-discussion] Cross-covariance function

2012-01-27 Thread Benjamin Root
On Friday, January 27, 2012, Pierre Haessig pierre.haes...@crans.org wrote: Le 26/01/2012 19:19, josef.p...@gmail.com a écrit : The discussion had this reversed, numpy matches the behavior of MATLAB, while R (statistics) only returns the cross covariance part as proposed. I would also say

Re: [Numpy-discussion] preferred way of testing empty arrays

2012-01-27 Thread Benjamin Root
On Fri, Jan 27, 2012 at 2:46 PM, Bartosz Telenczuk b.telenc...@biologie.hu-berlin.de wrote: I have been using numpy for several years and I am very impressed with its flexibility. However, there is one problem that has always bothered me. Quite often I need to test consistently whether a

Re: [Numpy-discussion] NetCDF4/numpy question

2012-01-27 Thread Benjamin Root
On Fri, Jan 27, 2012 at 4:37 PM, Howard how...@renci.org wrote: I have found, in using tricontourf, that in the mapping from data values to color values, the range of the data seems to include even the data from the masked triangles. This causes the data to be either monochromatic or

Re: [Numpy-discussion] numpy all unexpected result (generator)

2012-01-31 Thread Benjamin Root
On Tuesday, January 31, 2012, Alan G Isaac alan.is...@gmail.com wrote: On 1/31/2012 8:26 AM, Neal Becker wrote: I was just bitten by this unexpected behavior: In [24]: all ([i 0 for i in xrange (10)]) Out[24]: False In [25]: all (i 0 for i in xrange (10)) Out[25]: True Turns out: In

Re: [Numpy-discussion] numpy all unexpected result (generator)

2012-01-31 Thread Benjamin Root
On Tue, Jan 31, 2012 at 9:18 AM, Robert Kern robert.k...@gmail.com wrote: On Tue, Jan 31, 2012 at 15:13, Benjamin Root ben.r...@ou.edu wrote: Is np.all() using np.array() or np.asanyarray()? If the latter, I would expect it to return a numpy array from a generator. Why would you expect

Re: [Numpy-discussion] numpy all unexpected result (generator)

2012-01-31 Thread Benjamin Root
On Tue, Jan 31, 2012 at 10:05 AM, Olivier Delalleau sh...@keba.be wrote: Le 31 janvier 2012 10:50, Robert Kern robert.k...@gmail.com a écrit : On Tue, Jan 31, 2012 at 15:35, Benjamin Root ben.r...@ou.edu wrote: On Tue, Jan 31, 2012 at 9:18 AM, Robert Kern robert.k...@gmail.com wrote

Re: [Numpy-discussion] Broadcasting doesn't work with divide after tile

2012-02-01 Thread Benjamin Root
2012/2/1 martin großhauser mgroszhau...@gmail.com 2012/2/1 martin großhauser mgroszhau...@gmail.com: Hello, when I try in my script to divide a masked array by a scalar I get an error. The instruction is: sppa = sp / 100. sp is a masked array with ndim = 3. error is:

[Numpy-discussion] autocorrelation computation performance : use of np.correlate

2012-02-01 Thread Benjamin Root
On Wednesday, February 1, 2012, Pierre Haessig pierre.haes...@crans.org wrote: Hi, [I'm not sure whether this discussion belongs to numpy-discussion or scipy-dev] In day to day time series analysis I regularly need to look at the data autocorrelation (acorr or acf depending on the software

Re: [Numpy-discussion] fast method to to count a particular value in a large matrix

2012-02-04 Thread Benjamin Root
On Saturday, February 4, 2012, Naresh Pai n...@uark.edu wrote: I am somewhat new to Python (been coding with Matlab mostly). I am trying to simplify (and expedite) a piece of code that is currently a bottleneck in a larger code. I have a large array (7000 rows x 4500 columns) titled say, abc,

Re: [Numpy-discussion] Logical indexing and higher-dimensional arrays.

2012-02-07 Thread Benjamin Root
On Tue, Feb 7, 2012 at 11:11 AM, Jordi Gutiérrez Hermoso jord...@octave.org wrote: Consider the following. Is this a bug? Thanks, - Jordi G. H. --- #!/usr/bin/python import numpy as np x = np.reshape(np.random.uniform(size=2*3*4), [2,3,4])

Re: [Numpy-discussion] numpy.arange() error?

2012-02-09 Thread Benjamin Root
On Thursday, February 9, 2012, Sturla Molden stu...@molden.no wrote: Den 9. feb. 2012 kl. 22:44 skrev eat e.antero.ta...@gmail.com: Maybe this issue is raised also earlier, but wouldn't it be more consistent to let arange operate only with integers (like Python's range) and let linspace

Re: [Numpy-discussion] Migrating issues to GitHub

2012-02-11 Thread Benjamin Root
On Saturday, February 11, 2012, Travis Oliphant tra...@continuum.io wrote: How to people feel about moving the issue tracking for NumPy to Github? It looks like they have improved their issue tracking quite a bit and the workflow and integration with commits looks quite good from what I can see.

Re: [Numpy-discussion] [SciPy-Dev] Discussion with Guido van Rossum and (hopefully) core python-dev on scientific Python and Python3

2012-02-13 Thread Benjamin Root
On Monday, February 13, 2012, Aaron Meurer asmeu...@gmail.com wrote: I'd like the ability to make in (i.e., __contains__) return something other than a bool. Also, the ability to make the x y z syntax would be useful. It's been suggested that the ability to override the boolean operators

Re: [Numpy-discussion] Change in scalar upcasting rules for 1.6.x?

2012-02-13 Thread Benjamin Root
On Monday, February 13, 2012, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 13, 2012 at 9:04 PM, Travis Oliphant tra...@continuum.io wrote: I disagree with your assessment of the subscript operator, but I'm sure we will have plenty of time to discuss that. I don't think it's

Re: [Numpy-discussion] Change in scalar upcasting rules for 1.6.x?

2012-02-13 Thread Benjamin Root
On Monday, February 13, 2012, Travis Oliphant tra...@continuum.io wrote: On Feb 13, 2012, at 10:14 PM, Charles R Harris wrote: On Mon, Feb 13, 2012 at 9:04 PM, Travis Oliphant tra...@continuum.io wrote: I disagree with your assessment of the subscript operator, but I'm sure we will have

Re: [Numpy-discussion] Typecasting changes from 1.5.1 to 1.6.1

2012-02-14 Thread Benjamin Root
On Tuesday, February 14, 2012, Travis Oliphant tra...@continuum.io wrote: Here is the code I used to determine the coercion table of types. I first used *all* of the numeric_ops, narrowed it down to those with 2 inputs and 1 output, and then determined the run-time coercion table. Then, I

Re: [Numpy-discussion] Typecasting changes from 1.5.1 to 1.6.1

2012-02-14 Thread Benjamin Root
On Tuesday, February 14, 2012, Travis Oliphant tra...@continuum.io wrote: As you can see there were changes in each release. Most of these were minor prior to the change from 1.5.1 to 1.6.1. I am still reviewing the changes from 1.5.1 to 1.6.1.At first blush, it looks like there are a lot

Re: [Numpy-discussion] Numpy governance update - was: Updated differences between 1.5.1 to 1.6.1

2012-02-14 Thread Benjamin Root
On Tuesday, February 14, 2012, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Tue, Feb 14, 2012 at 1:54 PM, Travis Oliphant tra...@continuum.io wrote: There is a mailing list for numfocus that you can sign up for if you would like to be part of those discussions. Let me know if you

[Numpy-discussion] Numpy governance update - was: Updated differences between 1.5.1 to 1.6.1

2012-02-14 Thread Benjamin Root
On Tuesday, February 14, 2012, Travis Oliphant tra...@continuum.io wrote: I have to agree with Mathew here, to a point. There has been discussions of these groups before, but I don't recall any announcement of this group. Of course, now that it has been announced, maybe a link to it should be

[Numpy-discussion] Implicit conversion of python datetime to numpy datetime64?

2012-02-14 Thread Benjamin Root
Just a thought I had. Right now, I can pass a list of python ints or floats into np.array() and get a numpy array with a sensible dtype. Is there any reason why we can't do the same for python's datetime? Right now, it is very easy for me to make a list comprehension of datetime objects using

Re: [Numpy-discussion] Implicit conversion of python datetime to numpy datetime64?

2012-02-14 Thread Benjamin Root
On Tuesday, February 14, 2012, Mark Wiebe mwwi...@gmail.com wrote: On Tue, Feb 14, 2012 at 8:17 PM, Benjamin Root ben.r...@ou.edu wrote: Just a thought I had. Right now, I can pass a list of python ints or floats into np.array() and get a numpy array with a sensible dtype. Is there any reason

Re: [Numpy-discussion] recarray field access asymmetry

2012-02-15 Thread Benjamin Root
On Wednesday, February 15, 2012, David Gowers (kampu) 00a...@gmail.com wrote: Hi all, This email is about the difference, given a recarray 'arr', between A) arr.foo.x[0] and B) arr.foo[0].x Specifically, form A returns the 0-th x value, whereas form B raises AttributeError:

Re: [Numpy-discussion] recarray field access asymmetry

2012-02-15 Thread Benjamin Root
On Wednesday, February 15, 2012, Benjamin Root ben.r...@ou.edu wrote: On Wednesday, February 15, 2012, David Gowers (kampu) 00a...@gmail.com wrote: Hi all, This email is about the difference, given a recarray 'arr', between A) arr.foo.x[0] and B) arr.foo[0].x Specifically

Re: [Numpy-discussion] Implicit conversion of python datetime to numpy datetime64?

2012-02-15 Thread Benjamin Root
On Tuesday, February 14, 2012, Mark Wiebe mwwi...@gmail.com wrote: On Tue, Feb 14, 2012 at 9:37 PM, Benjamin Root ben.r...@ou.edu wrote: On Tuesday, February 14, 2012, Mark Wiebe mwwi...@gmail.com wrote: On Tue, Feb 14, 2012 at 8:17 PM, Benjamin Root ben.r...@ou.edu wrote: Just a thought I

Re: [Numpy-discussion] Implicit conversion of python datetime to numpy datetime64?

2012-02-15 Thread Benjamin Root
On Wed, Feb 15, 2012 at 8:29 AM, Benjamin Root ben.r...@ou.edu wrote: On Tuesday, February 14, 2012, Mark Wiebe mwwi...@gmail.com wrote: On Tue, Feb 14, 2012 at 9:37 PM, Benjamin Root ben.r...@ou.edu wrote: On Tuesday, February 14, 2012, Mark Wiebe mwwi...@gmail.com wrote: On Tue, Feb

Re: [Numpy-discussion] Numpy governance update

2012-02-15 Thread Benjamin Root
On Wed, Feb 15, 2012 at 2:03 PM, Alan G Isaac alan.is...@gmail.com wrote: On 2/15/2012 2:46 PM, Benjamin Root wrote: The NA discussion is the perfect example where a governance structure would help resolve disputes. How? I'm not seeing it. Who would have behaved differently and why

Re: [Numpy-discussion] Numpy governance update

2012-02-15 Thread Benjamin Root
On Wed, Feb 15, 2012 at 4:08 PM, T J tjhn...@gmail.com wrote: On Wed, Feb 15, 2012 at 12:45 PM, Alan G Isaac alan.is...@gmail.comwrote: for the core developers. The right way to produce a governance structure is to make concrete proposals and show how these proposals are in the interest of

Re: [Numpy-discussion] Numpy governance update

2012-02-15 Thread Benjamin Root
...@gmail.com wrote: Hi, On Wed, Feb 15, 2012 at 11:46 AM, Benjamin Root ben.r...@ou.edu mailto:ben.r...@ou.edu wrote: On Wed, Feb 15, 2012 at 1:32 PM, Alan G Isaac alan.is...@gmail.com mailto:alan.is...@gmail.com wrote

Re: [Numpy-discussion] Numpy governance update

2012-02-15 Thread Benjamin Root
On Wednesday, February 15, 2012, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Wed, Feb 15, 2012 at 4:27 PM, Dag Sverre Seljebotn d.s.seljeb...@astro.uio.no wrote: On 02/15/2012 02:24 PM, Mark Wiebe wrote: There certainly is governance now, it's just informal. It's a combination of

Re: [Numpy-discussion] Numpy governance update

2012-02-16 Thread Benjamin Root
On Thu, Feb 16, 2012 at 2:13 PM, Nathaniel Smith n...@pobox.com wrote: On Thu, Feb 16, 2012 at 5:17 PM, Travis Vaught tra...@vaught.net wrote: On Feb 16, 2012, at 10:56 AM, Nathaniel Smith wrote: Travis's proposal is that we go from a large number of self-selecting people putting in

[Numpy-discussion] Proposed Roadmap Overview

2012-02-16 Thread Benjamin Root
On Thursday, February 16, 2012, Warren Weckesser wrote: On Thu, Feb 16, 2012 at 4:39 PM, Travis Oliphant tra...@continuum.iowrote: Mark Wiebe and I have been discussing off and on (as well as talking with Charles) a good way forward to balance two competing desires: * addition of

Re: [Numpy-discussion] Numpy governance update

2012-02-16 Thread Benjamin Root
On Thursday, February 16, 2012, John Hunter wrote: On Thu, Feb 16, 2012 at 7:26 PM, Alan G Isaac alan.is...@gmail.comjavascript:_e({}, 'cvml', 'alan.is...@gmail.com'); wrote: On 2/16/2012 7:22 PM, Matthew Brett wrote: This has not been an encouraging episode in striving for consensus.

Re: [Numpy-discussion] Proposed Roadmap Overview

2012-02-17 Thread Benjamin Root
On Fri, Feb 17, 2012 at 1:00 PM, Christopher Jordan-Squire cjord...@uw.eduwrote: On Fri, Feb 17, 2012 at 10:21 AM, Mark Wiebe mwwi...@gmail.com wrote: On Fri, Feb 17, 2012 at 11:52 AM, Eric Firing efir...@hawaii.edu wrote: On 02/17/2012 05:39 AM, Charles R Harris wrote: On Fri,

Re: [Numpy-discussion] The end of numpy as we know it ?

2012-02-18 Thread Benjamin Root
On Saturday, February 18, 2012, Sturla Molden wrote: Den 18. feb. 2012 kl. 17:12 skrev Alan G Isaac alan.is...@gmail.comjavascript:; : How does stream-lined code written for maintainability (i.e., with helpful comments and tests) become *less* accessible to amateurs?? I think

Re: [Numpy-discussion] Proposed Roadmap Overview

2012-02-18 Thread Benjamin Root
On Sat, Feb 18, 2012 at 2:45 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Feb 18, 2012 at 1:39 PM, Matthew Brett matthew.br...@gmail.comwrote: Hi, On Sat, Feb 18, 2012 at 12:35 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Feb 18, 2012 at 12:21 PM,

Re: [Numpy-discussion] Possible roadmap addendum: building better text file readers

2012-02-23 Thread Benjamin Root
On Thu, Feb 23, 2012 at 3:14 PM, Robert Kern robert.k...@gmail.com wrote: On Thu, Feb 23, 2012 at 21:09, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Thu, Feb 23, 2012 at 04:07:04PM -0500, Wes McKinney wrote: In this last case for example, around 500 MB of RAM is taken up for an

Re: [Numpy-discussion] bincount([], minlength=2) should work right?

2012-02-25 Thread Benjamin Root
On Saturday, February 25, 2012, Alan G Isaac wrote: On 2/25/2012 4:44 PM, James Bergstra wrote: bincount([]) makes no sense, I disagree: http://permalink.gmane.org/gmane.comp.python.numeric.general/42041 but if a minlength argument is provided, then the routine should succeed.

Re: [Numpy-discussion] Floating point close function?

2012-03-03 Thread Benjamin Root
On Saturday, March 3, 2012, Robert Kern robert.k...@gmail.com wrote: On Sat, Mar 3, 2012 at 14:31, Ralf Gommers ralf.gomm...@googlemail.com wrote: On Sat, Mar 3, 2012 at 3:05 PM, Robert Kern robert.k...@gmail.com wrote: On Sat, Mar 3, 2012 at 13:59, Ralf Gommers ralf.gomm...@googlemail.com

Re: [Numpy-discussion] all elements equal

2012-03-05 Thread Benjamin Root
On Mon, Mar 5, 2012 at 1:44 PM, Keith Goodman kwgood...@gmail.com wrote: On Mon, Mar 5, 2012 at 11:36 AM, josef.p...@gmail.com wrote: How about numpy.ptp, to follow this line? I would expect it's single pass, but wouldn't short circuit compared to cython of Keith I[1] a = np.ones(10)

Re: [Numpy-discussion] Missing data again

2012-03-07 Thread Benjamin Root
On Wed, Mar 7, 2012 at 1:26 PM, Nathaniel Smith n...@pobox.com wrote: On Wed, Mar 7, 2012 at 5:17 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Wed, Mar 7, 2012 at 9:35 AM, Pierre Haessig pierre.haes...@crans.org Coming back to Travis proposition bit-pattern approaches to

Re: [Numpy-discussion] use for missing (ignored) data?

2012-03-07 Thread Benjamin Root
On Wednesday, March 7, 2012, Neal Becker ndbeck...@gmail.com wrote: Charles R Harris wrote: On Wed, Mar 7, 2012 at 1:05 PM, Neal Becker ndbeck...@gmail.com wrote: I'm wondering what is the use for the ignored data feature? I can use: A[valid_A_indexes] = whatever to process only the

Re: [Numpy-discussion] use for missing (ignored) data?

2012-03-07 Thread Benjamin Root
On Wednesday, March 7, 2012, Nathaniel Smith n...@pobox.com wrote: On Wed, Mar 7, 2012 at 8:05 PM, Neal Becker ndbeck...@gmail.com wrote: I'm wondering what is the use for the ignored data feature? I can use: A[valid_A_indexes] = whatever to process only the 'non-ignored' portions of A.

Re: [Numpy-discussion] draft enum NEP

2012-03-15 Thread Benjamin Root
On Thursday, March 15, 2012, Nathaniel Smith n...@pobox.com wrote: On Wed, Mar 14, 2012 at 1:44 AM, Mark Wiebe mwwi...@gmail.com wrote: On Fri, Mar 9, 2012 at 8:55 AM, Bryan Van de Ven bry...@continuum.io wrote: Hi all, I have started working on a NEP for adding an enumerated type to NumPy.

Re: [Numpy-discussion] float96 on windows32 is float64?

2012-03-15 Thread Benjamin Root
On Thursday, March 15, 2012, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Mar 15, 2012 at 10:17 PM, David Cournapeau courn...@gmail.com wrote: On Thu, Mar 15, 2012 at 11:10 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, Am I right in thinking that float96 on windows 32

Re: [Numpy-discussion] numpy installation problem?

2012-03-18 Thread Benjamin Root
On Sunday, March 18, 2012, ava...@famaf.unc.edu.ar wrote: Dear list, I am having problems installing matplotlib (from source) and fipy. I had installed numpy from source and it is running well: :~$ python -c import numpy; print numpy.__version__ 1.6.1 After being trying to solve this

Re: [Numpy-discussion] Slice specified axis

2012-04-06 Thread Benjamin Root
On Friday, April 6, 2012, Val Kalatsky wrote: The only slicing short-cut I can think of is the Ellipsis object, but it's not going to help you much here. The alternatives that come to my mind are (1) manipulation of shape directly and (2) building a string and running eval on it. Your

Re: [Numpy-discussion] Slice specified axis

2012-04-09 Thread Benjamin Root
On Mon, Apr 9, 2012 at 12:14 PM, Jonathan T. Niehof jnie...@lanl.govwrote: On 04/06/2012 06:54 AM, Benjamin Root wrote: Take a peek at how np.gradient() does it. It creates a list of None with a length equal to the number of dimensions, and then inserts a slice object in the appropriate

Re: [Numpy-discussion] Why is numpy.abs so much slower on complex64 than complex128 under windows 32-bit?

2012-04-10 Thread Benjamin Root
On Tue, Apr 10, 2012 at 12:57 PM, Francesc Alted franc...@continuum.iowrote: On 4/10/12 9:55 AM, Henry Gomersall wrote: On 10/04/2012 16:36, Francesc Alted wrote: In [10]: timeit c = numpy.complex64(numpy.abs(numpy.complex128(b))) 100 loops, best of 3: 12.3 ms per loop In [11]: timeit

Re: [Numpy-discussion] Slice specified axis

2012-04-10 Thread Benjamin Root
On Tue, Apr 10, 2012 at 12:52 PM, Jonathan T. Niehof jnie...@lanl.govwrote: On 04/09/2012 09:11 PM, Tony Yu wrote: I guess I wasn't reading very carefully and assumed that you meant a list of `slice(None)` instead of a list of `None`. My apologies to Ben...I wasn't being pedantic to be a

Re: [Numpy-discussion] partial computations

2012-04-12 Thread Benjamin Root
On Wed, Apr 11, 2012 at 11:38 PM, santhu kumar mesan...@gmail.com wrote: Hello all, I am trying to optimise a code and want your suggestions. A : - NX3 matrix (coordinates of N points) After performing pairwise distance computations(called pdist) between these points, depending upon a

Re: [Numpy-discussion] What is consensus anyway

2012-04-24 Thread Benjamin Root
On Tue, Apr 24, 2012 at 2:12 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Apr 24, 2012 at 9:25 AM, josef.p...@gmail.com wrote: On Tue, Apr 24, 2012 at 9:43 AM, Pierre Haessig pierre.haes...@crans.org wrote: Hi, Le 24/04/2012 15:14, Charles R Harris a écrit : a)

Re: [Numpy-discussion] What is consensus anyway

2012-04-24 Thread Benjamin Root
On Tue, Apr 24, 2012 at 3:23 PM, Stéfan van der Walt ste...@sun.ac.zawrote: On Tue, Apr 24, 2012 at 11:12 AM, Charles R Harris charlesr.har...@gmail.com wrote: The advantage of nans, I suppose, is that they are in the hardware and so Why are we having a discussion on NAN's in a thread on

Re: [Numpy-discussion] What is consensus anyway

2012-04-24 Thread Benjamin Root
On Tuesday, April 24, 2012, Matthew Brett wrote: Hi, On Tue, Apr 24, 2012 at 2:25 PM, Charles R Harris charlesr.har...@gmail.com javascript:; wrote: 2012/4/24 Stéfan van der Walt ste...@sun.ac.za javascript:; On Tue, Apr 24, 2012 at 11:12 AM, Charles R Harris

Re: [Numpy-discussion] What is consensus anyway

2012-04-25 Thread Benjamin Root
On Wednesday, April 25, 2012, Matthew Brett wrote: Hi, On Wed, Apr 25, 2012 at 2:35 PM, Travis Oliphant tra...@continuum.iojavascript:; wrote: Do you agree that Numpy has not been very successful in recruiting and maintaining new developers compared to its large user-base?

Re: [Numpy-discussion] What is consensus anyway

2012-04-25 Thread Benjamin Root
On Wednesday, April 25, 2012, Travis Oliphant wrote: On Apr 25, 2012, at 7:18 PM, josef.p...@gmail.com javascript:; wrote: Except for the big changes like NA and datetime, I think the debate is pretty boring. The main problem that I see for discussing technical issues is whether

Re: [Numpy-discussion] A crazy masked-array thought

2012-04-27 Thread Benjamin Root
On Fri, Apr 27, 2012 at 6:32 AM, Richard Hattersley rhatters...@gmail.comwrote: I know used a somewhat jokey tone in my original posting, but fundamentally it was a serious question concerning a live topic. So I'm curious about the lack of response. Has this all been covered before? Sorry if

Re: [Numpy-discussion] Issue Tracking

2012-04-30 Thread Benjamin Root
On Monday, April 30, 2012, Travis Oliphant wrote: Hey all, We have been doing some investigation of various approaches to issue tracking. The last time the conversation left this list was with Ralf's current list of preferences as: 1) Redmine 2) Trac 3) Github Since that time,

Re: [Numpy-discussion] record arrays initialization

2012-05-02 Thread Benjamin Root
On Wednesday, May 2, 2012, Stéfan van der Walt wrote: On Wed, May 2, 2012 at 4:46 PM, Kevin Jacobs jac...@bioinformed.comjavascript:; bioinfor...@gmail.com javascript:; wrote: A FLANN implementation should be even faster--perhaps by as much as another factor of two. I guess it

Re: [Numpy-discussion] Issue Tracking

2012-05-05 Thread Benjamin Root
On Saturday, May 5, 2012, Pauli Virtanen wrote: 05.05.2012 22:53, Ralf Gommers kirjoitti: [clip] would be great to get it done by end of June.To Charles' list and Ralf's suggestions, I would add setting up a server that can relay pull requests to the mailing list.

Re: [Numpy-discussion] Issue Tracking

2012-05-05 Thread Benjamin Root
On Saturday, May 5, 2012, Charles R Harris wrote: On Sat, May 5, 2012 at 8:50 PM, Benjamin Root ben.r...@ou.edujavascript:_e({}, 'cvml', 'ben.r...@ou.edu'); wrote: On Saturday, May 5, 2012, Pauli Virtanen wrote: 05.05.2012 22:53, Ralf Gommers kirjoitti: [clip] would be great

[Numpy-discussion] spurious space in printing record arrays?

2012-05-10 Thread Benjamin Root
Just noticed this in the output from printing some numpy record arrays: [[('2008081712', -24, -78.0, 20.10381469727, 45.0, -999.0, 0.0)] [ ('2008081718', -18, -79.584741211, 20.70762939453, 45.0, -999.0, 0.0)] [ ('2008081800', -12, -80.3305175781, 21.10381469727, 45.0,

Re: [Numpy-discussion] ANN: NumPy 1.6.2 release candidate 1

2012-05-11 Thread Benjamin Root
On Sat, May 5, 2012 at 2:15 PM, Ralf Gommers ralf.gomm...@googlemail.comwrote: Hi, I'm pleased to announce the availability of the first release candidate of NumPy 1.6.2. This is a maintenance release. Due to the delay of the NumPy 1.7.0, this release contains far more fixes than a regular

[Numpy-discussion] stable sort for structured dtypes?

2012-05-11 Thread Benjamin Root
Hello all, I need to sort a structured array in a stable manner. I am also sorting only by one of the keys, so I don't think lexsort() is stable in that respect. np.sort() allows for choosing 'mergesort', but it appears to not be implemented for structured arrays. Am I going to have to create

Re: [Numpy-discussion] stable sort for structured dtypes?

2012-05-11 Thread Benjamin Root
On Fri, May 11, 2012 at 12:00 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Fri, May 11, 2012 at 9:01 AM, Benjamin Root ben.r...@ou.edu wrote: Hello all, I need to sort a structured array in a stable manner. I am also sorting only by one of the keys, so I don't think lexsort

Re: [Numpy-discussion] Should arr.diagonal() return a copy or a view? (1.7 compatibility issue)

2012-05-12 Thread Benjamin Root
On Saturday, May 12, 2012, Travis Oliphant wrote: Another approach would be to introduce a method: a.diag(copy=False) and leave a.diagonal() alone. Then, a.diagonal() could be deprecated over 2-3 releases. -Travis +1 Ben Root ___

Re: [Numpy-discussion] Should arr.diagonal() return a copy or a view? (1.7 compatibility issue)

2012-05-16 Thread Benjamin Root
On Wed, May 16, 2012 at 9:55 AM, Nathaniel Smith n...@pobox.com wrote: On Tue, May 15, 2012 at 2:49 PM, Frédéric Bastien no...@nouiz.org wrote: Hi, In fact, I would arg to never change the current behavior, but add the flag for people that want to use it. Why? 1) There is

Re: [Numpy-discussion] python import question

2012-05-18 Thread Benjamin Root
On Friday, May 18, 2012, Chao YUE wrote: Dear all, This is only a small python import question. I think I'm right but just want some confirmation. Previously I have installed numpy 1.5.1. and then I used pip install --upgrade numpy to install numpy 1.6.1 But when I try to import numpy

Re: [Numpy-discussion] Checking for views (was: Should arr.diagonal() return a copy or aview?)

2012-05-24 Thread Benjamin Root
On Thu, May 24, 2012 at 10:56 AM, Jonathan T. Niehof jnie...@lanl.govwrote: On 05/23/2012 05:31 PM, T J wrote: It seems that there are a number of ways to check if an array is a view. Do we have a preferred way in the API that is guaranteed to stay available? Or are all of the various

Re: [Numpy-discussion] indexes in an array where value is greater than 1?

2012-05-25 Thread Benjamin Root
On Fri, May 25, 2012 at 11:17 AM, Chris Withers ch...@simplistix.co.ukwrote: Hi All, I have an array: arrrgh = numpy.zeros(1) A sparse collection of elements will have values greater than zero: arrrgh[] = 2 arrrgh[3453453] =42 The *wrong* way to do this is: for i in

Re: [Numpy-discussion] indexes in an array where value is greater than 1?

2012-05-27 Thread Benjamin Root
On Sunday, May 27, 2012, Chao YUE wrote: for me, np.nonzero() and np.where() both work. It seems they have same function. chao They are not identical. Nonzeros is for indices. The where function is really meant for a different purpose, but special-cases for this call signature. Ben Root

Re: [Numpy-discussion] better error message possible?

2012-06-01 Thread Benjamin Root
On Fri, Jun 1, 2012 at 9:14 AM, Nathaniel Smith n...@pobox.com wrote: On Fri, Jun 1, 2012 at 10:46 AM, Chris Withers ch...@simplistix.co.uk wrote: Hi All, Any reason why this: import numpy numpy.zeros(10)[-123] Traceback (most recent call last): File stdin, line 1, in

Re: [Numpy-discussion] some typestrings not recognized anymore

2012-06-03 Thread Benjamin Root
On Sunday, June 3, 2012, Ralf Gommers wrote: On Sun, Jun 3, 2012 at 4:49 PM, Nathaniel Smith n...@pobox.comjavascript:_e({}, 'cvml', 'n...@pobox.com'); wrote: On Sun, Jun 3, 2012 at 3:28 PM, Ralf Gommers ralf.gomm...@googlemail.com javascript:_e({}, 'cvml',

Re: [Numpy-discussion] 1D array sorting ascending and descending by fields

2012-06-04 Thread Benjamin Root
On Monday, June 4, 2012, Chris Barker wrote: On Mon, Jun 4, 2012 at 11:10 AM, Patrick Redmond plredm...@gmail.comjavascript:; wrote: Here's how I sorted primarily by field 'a' descending and secondarily by field 'b' ascending: could you multiply the numeric field by -1, sort, then put

Re: [Numpy-discussion] 1D array sorting ascending and descending by fields

2012-06-05 Thread Benjamin Root
On Tue, Jun 5, 2012 at 10:49 AM, Nathaniel Smith n...@pobox.com wrote: On Tue, Jun 5, 2012 at 1:17 AM, Benjamin Root ben.r...@ou.edu wrote: On Monday, June 4, 2012, Chris Barker wrote: On Mon, Jun 4, 2012 at 11:10 AM, Patrick Redmond plredm...@gmail.com wrote: Here's how I sorted

[Numpy-discussion] boolean indexing of structured arrays

2012-06-06 Thread Benjamin Root
Not sure if this is a bug or not. I am using a fairly recent master branch. # Setting up... import numpy as np a = np.zeros((10, 1), dtype=[('foo', 'f4'), ('bar', 'f4'), ('spam', 'f4')]) a['foo'] = np.random.random((10, 1)) a['bar'] = np.random.random((10, 1)) a['spam'] =

Re: [Numpy-discussion] Good way to develop numpy as popular choice!

2012-06-21 Thread Benjamin Root
On Thursday, June 21, 2012, Robert Kern wrote: On Thu, Jun 21, 2012 at 7:33 PM, eat e.antero.ta...@gmail.comjavascript:; wrote: Heh, On Thu, Jun 21, 2012 at 6:03 PM, Robert Kern robert.k...@gmail.comjavascript:; wrote: On Thu, Jun 21, 2012 at 3:59 PM, bob tnur

Re: [Numpy-discussion] Good way to develop numpy as popular choice!

2012-06-22 Thread Benjamin Root
On Fri, Jun 22, 2012 at 9:42 AM, eat e.antero.ta...@gmail.com wrote: Hi, On Fri, Jun 22, 2012 at 7:51 AM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Thu, Jun 21, 2012 at 08:59:09PM -0400, Benjamin Root wrote: munkres seems to be a pure python implementation

Re: [Numpy-discussion] Good way to develop numpy as popular choice!

2012-06-22 Thread Benjamin Root
On Fri, Jun 22, 2012 at 10:25 AM, Travis Oliphant tra...@continuum.iowrote: Accessing individual elements of NumPy arrays is slower than accessing individual elements of lists --- around 2.5x-3x slower.NumPy has to do more work to figure out what kind of indexing you are trying to do

Re: [Numpy-discussion] Created NumPy 1.7.x branch

2012-06-25 Thread Benjamin Root
On Mon, Jun 25, 2012 at 1:41 PM, Travis Oliphant tra...@continuum.iowrote: C was famous for bugs due to the lack of function prototypes. This was fixed with C99 and the stricter typing was a great help. Bugs are not due to lack of function prototypes. Bugs are due to mistakes that

Re: [Numpy-discussion] Created NumPy 1.7.x branch

2012-06-26 Thread Benjamin Root
On Tue, Jun 26, 2012 at 12:48 PM, David Cournapeau courn...@gmail.comwrote: On Tue, Jun 26, 2012 at 5:24 PM, Travis Oliphant tra...@continuum.io wrote: Let us note that that problem was due to Travis convincing David to include the Datetime work in the release against David's own best

Re: [Numpy-discussion] Matrix rank default tolerance - is it too low?

2012-06-26 Thread Benjamin Root
On Tuesday, June 26, 2012, Charles R Harris wrote: On Tue, Jun 26, 2012 at 3:42 PM, Matthew Brett matthew.br...@gmail.comwrote: Hi, On Mon, Jun 18, 2012 at 3:50 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, On Sun, Jun 17, 2012 at 7:22 PM, Charles R Harris

Re: [Numpy-discussion] Created NumPy 1.7.x branch

2012-06-26 Thread Benjamin Root
On Tuesday, June 26, 2012, Thouis (Ray) Jones wrote: On Tue, Jun 26, 2012 at 10:11 PM, Jason Grout jason-s...@creativetrax.com javascript:; wrote: On 6/26/12 3:06 PM, Dag Sverre Seljebotn wrote: Something the Sage project does very well is meeting often in person Another thing we have

Re: [Numpy-discussion] import numpy performance

2012-07-02 Thread Benjamin Root
On Mon, Jul 2, 2012 at 4:34 PM, Nathaniel Smith n...@pobox.com wrote: On Mon, Jul 2, 2012 at 8:17 PM, Andrew Dalke da...@dalkescientific.com wrote: In this email I propose a few changes which I think are minor and which don't really affect the external NumPy API but which I think could

  1   2   3   4   5   6   >