Re: [Numpy-discussion] Producing a Histogram When Bins Are Known

2009-11-26 Thread josef . pktd
On Fri, Nov 27, 2009 at 12:16 AM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: I have a list that already has the frequencies from 0 to 255. However, I'd like to make a histogram  that has say 32 bins whose ranges are 0-7, 8-15, ... 248-255. Is it possible? If they have equal sized (width)

Re: [Numpy-discussion] Help making better use of numpy array functions

2009-11-25 Thread josef . pktd
On Wed, Nov 25, 2009 at 4:13 PM, mdekauwe mdeka...@gmail.com wrote: I tried redoing the internal logic for example by using the where function but I can't seem to work out how to match up the logic. For example (note slightly different from above). If I change the main loop to lst =

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-24 Thread josef . pktd
On Sun, Nov 22, 2009 at 1:40 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sun, Nov 22, 2009 at 6:28 AM, josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 2:35 AM, David Cournapeau courn...@gmail.com wrote: On Wed, Nov 18, 2009 at 6:10 PM, David Cournapeau courn...@gmail.com

[Numpy-discussion] recompiling everything ?

2009-11-23 Thread josef . pktd
I have several scikits with c extension compiled against an older trunk version of 1.4 but didn't recompile them with the numpy trunk version that I currently use. scikits.timeseries just crashed on me taking several hours worth of examples with it. (scipy crashed during testing, so I rebuilt it

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-22 Thread josef . pktd
On Sun, Nov 22, 2009 at 2:35 AM, David Cournapeau courn...@gmail.com wrote: On Wed, Nov 18, 2009 at 6:10 PM, David Cournapeau courn...@gmail.com wrote: On Wed, Nov 18, 2009 at 12:38 AM,  josef.p...@gmail.com wrote: Now numpy builds without problems. When I run the tests I get 16 failures

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-22 Thread josef . pktd
On Sun, Nov 22, 2009 at 10:01 AM, David Cournapeau courn...@gmail.com wrote: On Sun, Nov 22, 2009 at 11:45 PM, Charles R Harris charlesr.har...@gmail.com wrote: Might be nice to print out the actual values of np.spacing and np.nextafter here. Yes, I should add some utilities to print those

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-22 Thread josef . pktd
On Sun, Nov 22, 2009 at 10:25 AM, David Cournapeau courn...@gmail.com wrote: On Mon, Nov 23, 2009 at 12:14 AM,  josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 10:01 AM, David Cournapeau courn...@gmail.com wrote: On Sun, Nov 22, 2009 at 11:45 PM, Charles R Harris

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-22 Thread josef . pktd
On Sun, Nov 22, 2009 at 10:37 AM, josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 10:25 AM, David Cournapeau courn...@gmail.com wrote: On Mon, Nov 23, 2009 at 12:14 AM,  josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 10:01 AM, David Cournapeau courn...@gmail.com wrote: On Sun, Nov

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-22 Thread josef . pktd
On Sun, Nov 22, 2009 at 10:47 AM, josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 10:37 AM,  josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at 10:25 AM, David Cournapeau courn...@gmail.com wrote: On Mon, Nov 23, 2009 at 12:14 AM,  josef.p...@gmail.com wrote: On Sun, Nov 22, 2009 at

Re: [Numpy-discussion] Correlation filter

2009-11-22 Thread josef . pktd
On Fri, Nov 20, 2009 at 4:27 PM, josef.p...@gmail.com wrote: On Fri, Nov 20, 2009 at 2:28 PM, Keith Goodman kwgood...@gmail.com wrote: On Fri, Nov 20, 2009 at 11:17 AM,  josef.p...@gmail.com wrote: On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman kwgood...@gmail.com wrote: def corr3(x, y):    

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 10:51 AM, Keith Goodman kwgood...@gmail.com wrote: I have a short 1d array x and a large 2d array y. I'd like to locate the places in the y array that are most like (correlated to) the x array. My first attempt, corr1, is too slow. My second attempt, corr2, is faster

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman kwgood...@gmail.com wrote: On Fri, Nov 20, 2009 at 8:53 AM,  josef.p...@gmail.com wrote: scipy.signal.correlate  would be fast, but it will not be easy to subtract the correct moving mean. Subtracting a standard moving mean would subtract

Re: [Numpy-discussion] Correlation filter

2009-11-20 Thread josef . pktd
On Fri, Nov 20, 2009 at 2:28 PM, Keith Goodman kwgood...@gmail.com wrote: On Fri, Nov 20, 2009 at 11:17 AM,  josef.p...@gmail.com wrote: On Fri, Nov 20, 2009 at 1:51 PM, Keith Goodman kwgood...@gmail.com wrote: def corr3(x, y):    x = x - x.mean()    x /= x.std()    nx = x.size    one =

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-17 Thread josef . pktd
On Mon, Nov 16, 2009 at 10:52 PM, David Cournapeau courn...@gmail.com wrote: On Tue, Nov 17, 2009 at 3:33 AM,  josef.p...@gmail.com wrote: Now, the numpy build runs for a while then breaks while building umath. Any ideas? The behavior of distutils with config files is mysterious, I gave up

Re: [Numpy-discussion] Strange inversion in shape with some slices

2009-11-17 Thread josef . pktd
On Tue, Nov 17, 2009 at 11:24 AM, Jean-Baptiste Rudant boogalo...@yahoo.fr wrote: Hello, I think there's something strange with shape when a slice is given by an array. import numpy as N my_array = N.ones((2, 3, 6)) ind = N.arange(4) #you hope to find (3, 4) print my_array[0, :, ind].shape

Re: [Numpy-discussion] failure building trunk with mingw

2009-11-16 Thread josef . pktd
On Mon, Nov 16, 2009 at 12:48 AM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: josef.p...@gmail.com wrote: Are there new changes to the configuration needed?   mingw 3.4.5, WindowsXP, Python 2.5.2 Python.h is not picked up anymore: _configtest.c:1:20: Python.h: No such file or

[Numpy-discussion] failure building trunk with mingw

2009-11-15 Thread josef . pktd
Are there new changes to the configuration needed? mingw 3.4.5, WindowsXP, Python 2.5.2 Python.h is not picked up anymore: _configtest.c:1:20: Python.h: No such file or directory Josef C:\Josef\_progs\Subversion\numpy-trunksetup.py bdist Running from numpy source directory. F2PY Version

Re: [Numpy-discussion] Initial implementation of histogram_discrete()

2009-11-14 Thread josef . pktd
On Sat, Nov 14, 2009 at 6:53 AM, Priit Laes pl...@plaes.org wrote: Ühel kenal päeval, R, 2009-11-13 kell 13:36, kirjutas Ernest Adrogué: 13/11/09 @ 09:41 (+0200), thus spake Priit Laes: Does anyone have a scenario where one would actually have both negative and positive numbers (integers) in

Re: [Numpy-discussion] Initial implementation of histogram_discrete()

2009-11-14 Thread josef . pktd
On Sat, Nov 14, 2009 at 7:10 AM, josef.p...@gmail.com wrote: On Sat, Nov 14, 2009 at 6:53 AM, Priit Laes pl...@plaes.org wrote: Ühel kenal päeval, R, 2009-11-13 kell 13:36, kirjutas Ernest Adrogué: 13/11/09 @ 09:41 (+0200), thus spake Priit Laes: Does anyone have a scenario where one would

Re: [Numpy-discussion] Initial implementation of histogram_discrete()

2009-11-12 Thread josef . pktd
On Thu, Nov 12, 2009 at 6:00 PM, Priit Laes pl...@plaes.org wrote: Hey, I cooked up an initial implementation for one-dimensional histogram_discrete(). Example: import numpy numpy.histogram_discrete([-1, 9, 9, 0, 3, 5, 3]) array([1, 1, 0, 0, 2, 0, 1, 0, 0, 0, 2])

Re: [Numpy-discussion] Initial implementation of histogram_discrete()

2009-11-12 Thread josef . pktd
On Thu, Nov 12, 2009 at 6:16 PM, David Warde-Farley d...@cs.toronto.edu wrote: On 12-Nov-09, at 6:09 PM, josef.p...@gmail.com wrote: On Thu, Nov 12, 2009 at 6:00 PM, Priit Laes pl...@plaes.org wrote: Hey, I cooked up an initial implementation for one-dimensional histogram_discrete().

Re: [Numpy-discussion] Initial implementation of histogram_discrete()

2009-11-12 Thread josef . pktd
On Thu, Nov 12, 2009 at 6:31 PM, Robert Kern robert.k...@gmail.com wrote: On Thu, Nov 12, 2009 at 17:25,  josef.p...@gmail.com wrote: If I could make a related wish, I wish np.bincount to take also a 2d array as weights. It does. If it does, I don't manage to make it work, and the new docs

Re: [Numpy-discussion] Multiple Regression

2009-11-12 Thread josef . pktd
On Thu, Nov 12, 2009 at 6:44 PM, Robert Kern robert.k...@gmail.com wrote: On Thu, Nov 12, 2009 at 17:38, Alexey Tigarev alexey.tiga...@gmail.com wrote: Hi All! I have implemented multiple regression in a following way: def multipleRegression(x, y):     Perform linear regression using

Re: [Numpy-discussion] finding close together points.

2009-11-10 Thread josef . pktd
On Tue, Nov 10, 2009 at 7:48 PM, Anne Archibald peridot.face...@gmail.com wrote: 2009/11/10 Christopher Barker chris.bar...@noaa.gov: Hi all, I have a bunch of points in 2-d space, and I need to find out which pairs of points are within a certain distance of one-another (regular old

Re: [Numpy-discussion] Use-case for np.choose

2009-11-09 Thread josef . pktd
On Mon, Nov 9, 2009 at 7:54 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: May I infer from the sudden silence that I finally have it? I think so, I assume that the result of broadcasting is unique, I haven't seen an example yet where broadcasting would depend on the sequence in which it is

Re: [Numpy-discussion] Use-case for np.choose

2009-11-09 Thread josef . pktd
On Mon, Nov 9, 2009 at 7:59 PM, josef.p...@gmail.com wrote: On Mon, Nov 9, 2009 at 7:54 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: May I infer from the sudden silence that I finally have it? I think so, I assume that the result of broadcasting is unique, I haven't seen an example

Re: [Numpy-discussion] Use-case for np.choose

2009-11-08 Thread josef . pktd
On Sun, Nov 8, 2009 at 5:00 AM, David Goldsmith d.l.goldsm...@gmail.com wrote: On Sun, Nov 8, 2009 at 12:57 AM, Anne Archibald peridot.face...@gmail.com wrote: 2009/11/8 David Goldsmith d.l.goldsm...@gmail.com: On Sat, Nov 7, 2009 at 11:59 PM, Anne Archibald peridot.face...@gmail.com

Re: [Numpy-discussion] Use-case for np.choose

2009-11-08 Thread josef . pktd
On Sun, Nov 8, 2009 at 10:03 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: OK, now I'm trying to wrap my brain around broadcasting in choose when both `a` *and* `choices` need to be (non-trivially) broadcast in order to arrive at a common shape, e.g.: c=np.arange(4).reshape((2,1,2)) #

Re: [Numpy-discussion] Use-case for np.choose

2009-11-08 Thread josef . pktd
On Sun, Nov 8, 2009 at 10:40 PM, Anne Archibald peridot.face...@gmail.com wrote: As Josef said, this is not correct. I think the key point of confusion is this: Do not pass choose two arrays. Pass it one array and a *list* of arrays. The fact that choices can be an array is a quirk we

Re: [Numpy-discussion] Use-case for np.choose

2009-11-08 Thread josef . pktd
On Sun, Nov 8, 2009 at 10:57 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: On Sun, Nov 8, 2009 at 7:23 PM, josef.p...@gmail.com wrote: On Sun, Nov 8, 2009 at 10:03 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: OK, now I'm trying to wrap my brain around broadcasting in choose when

Re: [Numpy-discussion] How to get rid of the loop?

2009-11-07 Thread josef . pktd
On Sat, Nov 7, 2009 at 1:51 PM, Stas K stanc...@gmail.com wrote: Can I get rid of the loop in this example? And what is the fastest way to get v in the example? ar = array([1,2,3]) for a in ar:    for b in ar:        v = a**2+b**2 ar[:,None]**2 + ar**2 array([[ 2, 5, 10], [ 5, 8,

Re: [Numpy-discussion] Use-case for np.choose

2009-11-07 Thread josef . pktd
On Sat, Nov 7, 2009 at 7:53 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: Thanks, Anne. On Sat, Nov 7, 2009 at 1:32 PM, Anne Archibald peridot.face...@gmail.com wrote: 2009/11/7 David Goldsmith d.l.goldsm...@gmail.com: snip Also, my experimenting suggests that the index array

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-05 Thread josef . pktd
On Thu, Nov 5, 2009 at 6:36 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Nov 5, 2009 at 4:26 PM, David Warde-Farley d...@cs.toronto.edu wrote: On 5-Nov-09, at 4:54 PM, David Goldsmith wrote: Interesting thread, which leaves me wondering two things: is it documented

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-05 Thread josef . pktd
On Thu, Nov 5, 2009 at 9:23 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Nov 5, 2009 at 7:04 PM, josef.p...@gmail.com wrote: On Thu, Nov 5, 2009 at 6:36 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Nov 5, 2009 at 4:26 PM, David Warde-Farley

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-05 Thread josef . pktd
On Thu, Nov 5, 2009 at 10:42 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: On Thu, Nov 5, 2009 at 3:26 PM, David Warde-Farley d...@cs.toronto.edu wrote: On 5-Nov-09, at 4:54 PM, David Goldsmith wrote: Interesting thread, which leaves me wondering two things: is it documented

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-05 Thread josef . pktd
On Thu, Nov 5, 2009 at 10:28 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Nov 5, 2009 at 7:53 PM, josef.p...@gmail.com wrote: On Thu, Nov 5, 2009 at 9:23 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Thu, Nov 5, 2009 at 7:04 PM, josef.p...@gmail.com wrote:

Re: [Numpy-discussion] converting discrete data to unique integers

2009-11-04 Thread josef . pktd
On Wed, Nov 4, 2009 at 3:57 PM, David Warde-Farley d...@cs.toronto.edu wrote: Thanks Alan and Robert, I probably should have mentioned that I was interested in obtaining the corresponding integer for each value in the array d, in which case the dictionary bit works but would require a further

Re: [Numpy-discussion] Ignorance question

2009-11-01 Thread josef . pktd
On Sun, Nov 1, 2009 at 9:58 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: I Googled scipy brownian and the top hit was the doc for numpy.random.wald, but said doc has a tone that suggests there are more sophisticated ways to generate a random Brownian signal? Or is wald indeed SotA? 

Re: [Numpy-discussion] Ignorance question

2009-11-01 Thread josef . pktd
On Sun, Nov 1, 2009 at 10:26 PM, josef.p...@gmail.com wrote: On Sun, Nov 1, 2009 at 9:58 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: I Googled scipy brownian and the top hit was the doc for numpy.random.wald, but said doc has a tone that suggests there are more sophisticated ways to

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-01 Thread josef . pktd
On Sun, Nov 1, 2009 at 10:20 PM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: Thomas Robitaille wrote: Hi, I'm trying to generate random 64-bit integer values for integers and floats using Numpy, within the entire range of valid values for that type. To generate random 32-bit floats,

Re: [Numpy-discussion] Random int64 and float64 numbers

2009-11-01 Thread josef . pktd
On Sun, Nov 1, 2009 at 10:55 PM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: josef.p...@gmail.com wrote: array([ Inf,  Inf,  Inf,  Inf,  Inf,  Inf,  Inf,  Inf,  Inf,  Inf]) might actually be the right answer if you want a uniform distribution on the real line. Does it make sense

Re: [Numpy-discussion] Using matplotlib's prctile on masked arrays

2009-10-28 Thread josef . pktd
On Wed, Oct 28, 2009 at 9:52 AM, Gökhan Sever gokhanse...@gmail.com wrote: On Tue, Oct 27, 2009 at 12:23 PM, Pierre GM pgmdevl...@gmail.com wrote: On Oct 27, 2009, at 7:56 AM, Gökhan Sever wrote: Unfortunately, matplotlib.mlab's prctile cannot handle this division: Actually, the

Re: [Numpy-discussion] TypeError when calling numpy.kaiser()

2009-10-21 Thread josef . pktd
On Sun, Oct 18, 2009 at 3:11 PM, Jeffrey McGee jeffamc...@gmail.com wrote: Howdy, I'm having trouble getting the kaiser window to work. Anytime I try to call numpy.kaiser(), it throws an exception. Here's the output when I run the example code from

Re: [Numpy-discussion] GSOC 2010

2009-10-21 Thread josef . pktd
On Wed, Oct 21, 2009 at 3:02 PM, Charles R Harris charlesr.har...@gmail.com wrote: Hi All, I don't feel that numpy/scipy did as well in GSOC 2009 as it could have.  I think this was mostly due to lack of preparation on our part, we weren't ready when the students started showing up on the

Re: [Numpy-discussion] GSOC 2010

2009-10-21 Thread josef . pktd
On Wed, Oct 21, 2009 at 3:23 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Wed, Oct 21, 2009 at 1:11 PM, josef.p...@gmail.com wrote: On Wed, Oct 21, 2009 at 3:02 PM, Charles R Harris charlesr.har...@gmail.com wrote: Hi All, I don't feel that numpy/scipy did as well in GSOC

Re: [Numpy-discussion] Optimized sum of squares

2009-10-20 Thread josef . pktd
On Sun, Oct 18, 2009 at 6:06 AM, Gary Ruben gru...@bigpond.net.au wrote: Hi Gaël, If you've got a 1D array/vector called a, I think the normal idiom is np.dot(a,a) For the more general case, I think np.tensordot(a, a, axes=something_else) should do it, where you should be able to figure

Re: [Numpy-discussion] Optimized sum of squares

2009-10-20 Thread josef . pktd
On Tue, Oct 20, 2009 at 3:09 PM, Anne Archibald peridot.face...@gmail.com wrote: 2009/10/20  josef.p...@gmail.com: On Sun, Oct 18, 2009 at 6:06 AM, Gary Ruben gru...@bigpond.net.au wrote: Hi Gaël, If you've got a 1D array/vector called a, I think the normal idiom is np.dot(a,a) For the

Re: [Numpy-discussion] Optimized sum of squares

2009-10-18 Thread josef . pktd
On Sun, Oct 18, 2009 at 12:06 PM, Skipper Seabold jsseab...@gmail.com wrote: On Sun, Oct 18, 2009 at 8:09 AM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Sun, Oct 18, 2009 at 09:06:15PM +1100, Gary Ruben wrote: Hi Gaël, If you've got a 1D array/vector called a, I think the normal

Re: [Numpy-discussion] vectorized version of logsumexp? (from scipy.maxentropy)

2009-10-17 Thread josef . pktd
On Sat, Oct 17, 2009 at 2:02 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Oct 17, 2009 at 11:54 AM, josef.p...@gmail.com wrote: On Sat, Oct 17, 2009 at 1:20 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Oct 17, 2009 at 9:36 AM, per freem

Re: [Numpy-discussion] vectorized version of logsumexp? (from scipy.maxentropy)

2009-10-17 Thread josef . pktd
On Sat, Oct 17, 2009 at 7:46 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Oct 17, 2009 at 5:27 PM, josef.p...@gmail.com wrote: On Sat, Oct 17, 2009 at 2:02 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Oct 17, 2009 at 11:54 AM, josef.p...@gmail.com

Re: [Numpy-discussion] vectorized version of logsumexp? (from scipy.maxentropy)

2009-10-17 Thread josef . pktd
On Sat, Oct 17, 2009 at 8:00 PM, josef.p...@gmail.com wrote: On Sat, Oct 17, 2009 at 7:46 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Oct 17, 2009 at 5:27 PM, josef.p...@gmail.com wrote: On Sat, Oct 17, 2009 at 2:02 PM, Charles R Harris charlesr.har...@gmail.com wrote:

Re: [Numpy-discussion] extension questions: f2py and cython

2009-10-15 Thread josef . pktd
On Thu, Oct 15, 2009 at 8:53 AM, Robin robi...@gmail.com wrote: Hi, I have another question about distributing a Python extension which uses f2py wrapped code. Ideally I'd like to keep pure Python/Numpy alternatives and just use fortran version if available - but I think that should be OK.

Re: [Numpy-discussion] Google Groups archive?

2009-10-15 Thread josef . pktd
On Fri, Oct 16, 2009 at 12:14 AM, David Warde-Farley d...@cs.toronto.edu wrote: Does anyone know what happened to the Google Groups archive of this list? when I try to access it, I see: Cannot find numpy-discussion The group named numpy-discussion has been removed because it violated

Re: [Numpy-discussion] finding nonzero elements in list

2009-10-12 Thread josef . pktd
On Mon, Oct 12, 2009 at 10:18 AM, per freem perfr...@gmail.com wrote: hi all, i'm trying to find nonzero elements in an array, as follows: a = array([[1, 0],       [1, 1],       [1, 1],       [0, 1]]) i want to find all elements that are [1,1]. i tried: nonzero(a == [1,0]) but i cannot

Re: [Numpy-discussion] tostring() for array rows

2009-10-09 Thread josef . pktd
On Wed, Oct 7, 2009 at 1:20 PM, Christopher Barker chris.bar...@noaa.gov wrote: josef.p...@gmail.com wrote: I wanted to avoid the python loop and thought creating the view will be faster with large arrays. But for this I need to know the memory length of a row of arbitrary types for the

Re: [Numpy-discussion] Tabular data package

2009-10-06 Thread josef . pktd
On Mon, Oct 5, 2009 at 5:22 PM, Elaine Angelino elaine.angel...@gmail.com wrote: Hi there, We are writing to announce the release of Tabular, a package of Python modules for working with tabular data. Tabular is a package of Python modules for working with tabular data. Its main object is

Re: [Numpy-discussion] numpy-discussion in google groups ?

2009-10-05 Thread josef . pktd
On Mon, Oct 5, 2009 at 9:00 AM, Paul Rudin p...@rudin.co.uk wrote: denis bzowy denis-bz...@t-online.de writes: Folks,   http://groups.google.com/group/numpy-discussion - The group named numpy-discussion has been removed because it violated Google's Terms Of Service however scipy-user

[Numpy-discussion] poly class question

2009-10-02 Thread josef . pktd
Is there a way in numpy (or scipy) to get an infinite expansion for the inverse of a polynomial (for a finite number of terms) np.poly1d([ -0.8, 1])**(-1) application for example the MA representation of an AR(1) and fractional powers np.poly1d([ -1, 1])**0.5 this is useful for fractionally

Re: [Numpy-discussion] fixed_pt some progress and a question

2009-09-29 Thread josef . pktd
On Tue, Sep 29, 2009 at 10:22 AM, Neal Becker ndbeck...@gmail.com wrote: This doesn't work either:    def as_double (self):        import math        def _as_double_1 (x):            return math.ldexp (x, -self.frac_bits)        vecfunc = np.vectorize (_as_double_1, otypes=[np.float])      

Re: [Numpy-discussion] making the distinction between -0.0 and 0.0..

2009-09-29 Thread josef . pktd
On Tue, Sep 29, 2009 at 6:48 PM, Christopher Barker chris.bar...@noaa.gov wrote: I'm assuming it's a bug that was fixed somewhere in between? It works on my 2.5, on a PPC: In [10]: struct.pack('d', -0.0) Out[10]: '\x80\x00\x00\x00\x00\x00\x00\x00' In [11]: struct.pack('d', -0.0) Out[11]:

Re: [Numpy-discussion] Tuple outer product?

2009-09-25 Thread josef . pktd
On Fri, Sep 25, 2009 at 2:53 PM, Keith Goodman kwgood...@gmail.com wrote: On Fri, Sep 25, 2009 at 11:05 AM, Mads Ipsen m...@comxnet.dk wrote: Sure, and you could do the same thing using two nested for loops. But its bound to be slow when the arrays become large (and they will). The question is

Re: [Numpy-discussion] simple indexing question

2009-09-23 Thread josef . pktd
On Wed, Sep 23, 2009 at 11:12 AM, Neal Becker ndbeck...@gmail.com wrote: Robert Cimrman wrote: Neal Becker wrote: I have an array: In [12]: a Out[12]: array([[0, 1, 2, 3, 4],        [5, 6, 7, 8, 9]]) And a selection array: In [13]: b Out[13]: array([1, 1, 1, 1, 1]) I want a

Re: [Numpy-discussion] numpy docstring sphinx pre-processors

2009-09-21 Thread josef . pktd
On Mon, Sep 21, 2009 at 1:08 PM, Robert Kern robert.k...@gmail.com wrote: On Mon, Sep 21, 2009 at 12:03, Elaine Angelino elaine.angel...@gmail.com wrote: Hi there -- I have been working on a small Python package whose central data object comes from Numpy (the record array object). I would

[Numpy-discussion] how can we concatenate structured arrays ?

2009-09-16 Thread josef . pktd
I have two structured arrays of different types. How can I horizontally concatenate the two arrays? Is there a direct way, or do I need to start from scratch? nobs = 10 testdata = np.random.randint(3, size=(nobs,4)).view([('a',int),('b',int),('c',int),('d',int)]) testdatacont = np.random.normal(

Re: [Numpy-discussion] how can we concatenate structured arrays ?

2009-09-16 Thread josef . pktd
On Wed, Sep 16, 2009 at 4:25 AM, josef.p...@gmail.com wrote: I have two structured arrays of different types. How can I horizontally concatenate the two arrays? Is there a direct way, or do I need to start from scratch? nobs = 10 testdata = np.random.randint(3,

Re: [Numpy-discussion] how can we concatenate structured arrays ?

2009-09-16 Thread josef . pktd
On Wed, Sep 16, 2009 at 4:42 AM, Pierre GM pgmdevl...@gmail.com wrote: On Sep 16, 2009, at 4:25 AM, josef.p...@gmail.com wrote: I have two structured arrays of different types. How can I horizontally concatenate the two arrays? Is there a direct way, or do I need to start from scratch?

Re: [Numpy-discussion] `missing` argument in genfromtxt only a string?

2009-09-15 Thread josef . pktd
On Tue, Sep 15, 2009 at 10:44 AM, Skipper Seabold jsseab...@gmail.com wrote: On Tue, Sep 15, 2009 at 9:43 AM, Bruce Southey bsout...@gmail.com wrote: On 09/14/2009 09:31 PM, Skipper Seabold wrote: On Mon, Sep 14, 2009 at 9:59 PM, Pierre GMpgmdevl...@gmail.com  wrote: [snip] OK, I see the

Re: [Numpy-discussion] Why is the truth value of ndarray not simply size0 ?

2009-09-13 Thread josef . pktd
On Sun, Sep 13, 2009 at 9:05 AM, Alan G Isaac ais...@american.edu wrote: On 9/13/2009 7:46 AM, Robert wrote: 2 ways seem to be consistently Pythonic and logical: size 0; or any(a) (*); and the later option may be more 'numerical'. Well, *there's* the problem. As a user I have felt more than

Re: [Numpy-discussion] Dot product performance on python 2.6 (windows)

2009-09-11 Thread josef . pktd
On Fri, Sep 11, 2009 at 5:25 AM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: V. Armando Solé wrote: David Cournapeau wrote: V. Armando Solé wrote: Hello, It seems to point towards a packaging problem. In python 2.5, I can do: import numpy.core._dotblas as dotblas

Re: [Numpy-discussion] Row-wise dot product?

2009-09-07 Thread josef . pktd
On Mon, Sep 7, 2009 at 10:09 AM, Hans-Andreas Engeleng...@deshaw.com wrote: From: T J tjhnson at gmail.com Is there a better way to achieve the following, perhaps without the python for loop? x.shape (1,3) y.shape (1,3) z = empty(len(x)) for i in range(1): ...    z[i] =

Re: [Numpy-discussion] Behavior from a change in dtype?

2009-09-07 Thread josef . pktd
On Mon, Sep 7, 2009 at 6:36 PM, Skipper Seaboldjsseab...@gmail.com wrote: Hello all, I ran into a problem with some of my older code (since figured out the user error).  However, in trying to give a simple example that replicates the problem I was having, I ran into this. In [19]: a =

Re: [Numpy-discussion] Behavior from a change in dtype?

2009-09-07 Thread josef . pktd
On Mon, Sep 7, 2009 at 8:01 PM, Skipper Seaboldjsseab...@gmail.com wrote: On Mon, Sep 7, 2009 at 7:35 PM, josef.p...@gmail.com wrote: On Mon, Sep 7, 2009 at 6:36 PM, Skipper Seaboldjsseab...@gmail.com wrote: Hello all, I ran into a problem with some of my older code (since figured out the

Re: [Numpy-discussion] help creating a reversed cumulative histogram

2009-09-03 Thread josef . pktd
On Thu, Sep 3, 2009 at 9:23 AM, Tim Michelsentimmichel...@gmx-topmail.de wrote: Hello, I have checked the snippets you proposed. It does what I wanted to achieve. Obviously, I had to substract the values as Robert demonstrated. This could also be perceived from the figure I posted. I

Re: [Numpy-discussion] help creating a reversed cumulative histogram

2009-09-03 Thread josef . pktd
On Thu, Sep 3, 2009 at 12:58 PM, Tim Michelsentimmichel...@gmx-topmail.de wrote: My first stop is usually wikipedia: [...] Thanks. So I I'known that I have to call the beast a empirical inverse survival function, Robert would also have foundit easier to help. Anyway, step by step... In

Re: [Numpy-discussion] help creating a reversed cumulative histogram

2009-09-02 Thread josef . pktd
On Wed, Sep 2, 2009 at 7:26 PM, Robert Kernrobert.k...@gmail.com wrote: On Wed, Sep 2, 2009 at 18:15, Tim Michelsentimmichel...@gmx-topmail.de wrote: Hello fellow numy users, I posted some questions on histograms recently [1, 2] but still couldn't find  a solution. I am trying to create a

Re: [Numpy-discussion] What type should / return in python 3k when applied to two integer types?

2009-08-28 Thread josef . pktd
On Fri, Aug 28, 2009 at 9:55 AM, Pauli Virtanenp...@iki.fi wrote: Fri, 28 Aug 2009 09:46:39 -0400, Neal Becker kirjoitti: Robert Kern wrote: On Thu, Aug 27, 2009 at 14:22, Christopher Barkerchris.bar...@noaa.gov wrote: By the way -- is there something about py3k that changes all this? Or

Re: [Numpy-discussion] What type should / return in python 3k when applied to two integer types?

2009-08-28 Thread josef . pktd
On Fri, Aug 28, 2009 at 10:46 AM, Neal Beckerndbeck...@gmail.com wrote: Charles R Harris wrote: On Fri, Aug 28, 2009 at 8:08 AM, josef.p...@gmail.com wrote: On Fri, Aug 28, 2009 at 9:55 AM, Pauli Virtanenp...@iki.fi wrote: Fri, 28 Aug 2009 09:46:39 -0400, Neal Becker kirjoitti: Robert

Re: [Numpy-discussion] What type should / return in python 3k when applied to two integer types?

2009-08-28 Thread josef . pktd
On Fri, Aug 28, 2009 at 12:46 PM, Charles R Harrischarlesr.har...@gmail.com wrote: On Fri, Aug 28, 2009 at 8:08 AM, josef.p...@gmail.com wrote: On Fri, Aug 28, 2009 at 9:55 AM, Pauli Virtanenp...@iki.fi wrote: Fri, 28 Aug 2009 09:46:39 -0400, Neal Becker kirjoitti: Robert Kern wrote:

Re: [Numpy-discussion] histogram: sum up values in each bin

2009-08-27 Thread josef . pktd
On Thu, Aug 27, 2009 at 8:23 AM, alexander bakerbaker.alexan...@gmail.com wrote: Here is an example, this does something a extra at the end but shows how the bins can be used. Regards Alex Baker. from scipy.stats import norm r = norm.rvs(size=1) import numpy as np p, bins =

Re: [Numpy-discussion] histogram: sum up values in each bin

2009-08-27 Thread josef . pktd
On Thu, Aug 27, 2009 at 9:23 AM, Vincent Schutsc...@sarvision.nl wrote: Tim Michelsen wrote: Hello, I need some advice on histograms. If I interpret the documentation [1, 2] for numpy.histogram correctly, the result of the function is a count of the occurences sorted into each bin. (n,

Re: [Numpy-discussion] histogram: sum up values in each bin

2009-08-27 Thread josef . pktd
On Thu, Aug 27, 2009 at 12:49 PM, Tim Michelsentimmichel...@gmx-topmail.de wrote: Tim, do you mean, that you want to apply other functions, e.g. mean or variance, to the original values but calculated per bin? Sorry that I forgot to add this. Shame. I would like to apply these mathematical

Re: [Numpy-discussion] What type should / return in python 3k when applied to two integer types?

2009-08-27 Thread josef . pktd
On Thu, Aug 27, 2009 at 3:57 PM, Charles R Harrischarlesr.har...@gmail.com wrote: On Thu, Aug 27, 2009 at 1:46 PM, Robert Kern robert.k...@gmail.com wrote: On Thu, Aug 27, 2009 at 12:43, Charles R Harrischarlesr.har...@gmail.com wrote: On Thu, Aug 27, 2009 at 1:27 PM, Robert Kern

Re: [Numpy-discussion] histogram: sum up values in each bin

2009-08-27 Thread josef . pktd
On Thu, Aug 27, 2009 at 1:27 PM, josef.p...@gmail.com wrote: On Thu, Aug 27, 2009 at 12:49 PM, Tim Michelsentimmichel...@gmx-topmail.de wrote: Tim, do you mean, that you want to apply other functions, e.g. mean or variance, to the original values but calculated per bin? Sorry that I forgot to

Re: [Numpy-discussion] mirr test correctly fails for given input.

2009-08-26 Thread josef . pktd
On Wed, Aug 26, 2009 at 10:08 AM, Skipper Seaboldjsseab...@gmail.com wrote: On Wed, Aug 26, 2009 at 1:45 AM, josef.p...@gmail.com wrote: On Tue, Aug 25, 2009 at 11:38 PM, Charles R Harrischarlesr.har...@gmail.com wrote: So is it a bug in the test or a bug in the implementation? The problem is

Re: [Numpy-discussion] mirr test correctly fails for given input.

2009-08-26 Thread josef . pktd
On Wed, Aug 26, 2009 at 11:25 AM, josef.p...@gmail.com wrote: On Wed, Aug 26, 2009 at 10:08 AM, Skipper Seaboldjsseab...@gmail.com wrote: On Wed, Aug 26, 2009 at 1:45 AM, josef.p...@gmail.com wrote: On Tue, Aug 25, 2009 at 11:38 PM, Charles R Harrischarlesr.har...@gmail.com wrote: So is it a

Re: [Numpy-discussion] mirr test correctly fails for given input.

2009-08-25 Thread josef . pktd
On Tue, Aug 25, 2009 at 11:38 PM, Charles R Harrischarlesr.har...@gmail.com wrote: So is it a bug in the test or a bug in the implementation? The problem is that the slice values[1:] when values =  [-12,39000,3,21000,37000,46000] contains no negative number and a nan is returned. This

Re: [Numpy-discussion] itemsize() doesn't work

2009-08-20 Thread josef . pktd
2009/8/20 Stéfan van der Walt ste...@sun.ac.za: 2009/8/20 Dr. Phillip M. Feldman pfeld...@verizon.net: I've been reading the online NumPy tutorial at the following URL: http://numpy.scipy.org/numpydoc/numpy-10.html When I try the following example, I get an error message: In [1]:

Re: [Numpy-discussion] unique rows of array

2009-08-18 Thread josef . pktd
On Tue, Aug 18, 2009 at 2:01 AM, Maria Liukisliu...@usc.edu wrote: Josef, Many thanks for the example! It should become an official NumPy recipe :) Thanks again, Masha liu...@usc.edu Actually, there is also an implementation of unique rows in scipy.stats._support. It

Re: [Numpy-discussion] Strange crash in numpy.linalg.lstsq.

2009-08-17 Thread josef . pktd
On Mon, Aug 17, 2009 at 1:42 PM, Jonathan Taylorjonathan.tay...@utoronto.ca wrote: Hi, I am getting a strange crash in numpy.linalg.lstsq.  I have put the code that causes the crash along with two data files on my website at: http://www.cs.toronto.edu/~jtaylor/crash/ I would be interested

Re: [Numpy-discussion] unique rows of array

2009-08-17 Thread josef . pktd
On Tue, Aug 18, 2009 at 12:30 AM, Maria Liukisliu...@usc.edu wrote: Hello everybody, While re-implementing some Matlab code in Python, I've run into a problem of finding a NumPy function analogous to the Matlab's unique(array, 'rows') to get unique rows of an array. Searching the web, I've

Re: [Numpy-discussion] unique rows of array

2009-08-17 Thread josef . pktd
On Tue, Aug 18, 2009 at 12:59 AM, Maria Liukisliu...@usc.edu wrote: On Aug 17, 2009, at 9:51 PM, Charles R Harris wrote: On Mon, Aug 17, 2009 at 10:30 PM, Maria Liukis liu...@usc.edu wrote: Hello everybody, While re-implementing some Matlab code in Python, I've run into a problem of

Re: [Numpy-discussion] unique rows of array

2009-08-17 Thread josef . pktd
On Tue, Aug 18, 2009 at 1:03 AM, josef.p...@gmail.com wrote: On Tue, Aug 18, 2009 at 12:59 AM, Maria Liukisliu...@usc.edu wrote: On Aug 17, 2009, at 9:51 PM, Charles R Harris wrote: On Mon, Aug 17, 2009 at 10:30 PM, Maria Liukis liu...@usc.edu wrote: Hello everybody, While re-implementing

Re: [Numpy-discussion] Howto create a record array from arrays without copying their data

2009-08-12 Thread josef . pktd
On Wed, Aug 12, 2009 at 11:28 AM, Ryan Mayrma...@gmail.com wrote: On Wed, Aug 12, 2009 at 10:22 AM, Ralph Heinkel ra...@dont-mind.de wrote: Hi, I'm creating (actually calculating) a set of very large 1-d arrays (vectors), which I would like to assemble into a record array so I can access

[Numpy-discussion] is this a chararray decode bug

2009-08-12 Thread josef . pktd
(copied from the lengthy unicode thread in scipy-dev, so it doesn't get lost) this looks like a bug ? or is it a known limitation that chararrays cannot be 0-d b0= np.array(u'\xe9','U1').view(np.chararray) print b0.encode('cp1252') Traceback (most recent call last): File pyshell#47, line 1,

Re: [Numpy-discussion] add axis to results of reduction (mean, min, ...)

2009-08-10 Thread josef . pktd
On Mon, Aug 10, 2009 at 11:55 AM, Keith Goodmankwgood...@gmail.com wrote: On Thu, Aug 6, 2009 at 9:07 AM, Robert Kernrobert.k...@gmail.com wrote: On Thu, Aug 6, 2009 at 11:03, Keith Goodmankwgood...@gmail.com wrote: On Thu, Aug 6, 2009 at 8:55 AM, josef.p...@gmail.com wrote: What's the best

Re: [Numpy-discussion] speeding up getting a subset of a data array

2009-08-10 Thread josef . pktd
On Mon, Aug 10, 2009 at 8:52 PM, Brennan Williamsbrennan.willi...@visualreservoir.com wrote: Hi No doubt asked many times before so apologies I'm pulling a subset array out of a data array where I have a list of the indices I want (could be an array rather than a list actually - I have

Re: [Numpy-discussion] [SciPy-User] Sanity checklist for those attending the SciPy'09 tutorials

2009-08-09 Thread josef . pktd
On Sun, Aug 9, 2009 at 5:12 AM, Fernando Perezfperez@gmail.com wrote: Hi all, [ sorry for spamming the list, but even though I sent this to all the email addresses I have on file for tutorial attendees, I know I am missing a few, so I hope they see this message. ] In order to make your

Re: [Numpy-discussion] Power distribution

2009-08-07 Thread josef . pktd
On Fri, Aug 7, 2009 at 5:25 PM, Andrew Hawrylukhawr...@novachem.com wrote: Hmm ... good point. It appears to give a probability distribution proportional to x**(a-1), but I see no good reason why the domain should be limited to [0,1]. def test(a):    nums =

Re: [Numpy-discussion] Power distribution

2009-08-07 Thread josef . pktd
On Fri, Aug 7, 2009 at 5:42 PM, josef.p...@gmail.com wrote: On Fri, Aug 7, 2009 at 5:25 PM, Andrew Hawrylukhawr...@novachem.com wrote: Hmm ... good point. It appears to give a probability distribution proportional to x**(a-1), but I see no good reason why the domain should be limited to [0,1].

Re: [Numpy-discussion] Power distribution

2009-08-07 Thread josef . pktd
On Fri, Aug 7, 2009 at 6:13 PM, josef.p...@gmail.com wrote: On Fri, Aug 7, 2009 at 5:42 PM, josef.p...@gmail.com wrote: On Fri, Aug 7, 2009 at 5:25 PM, Andrew Hawrylukhawr...@novachem.com wrote: Hmm ... good point. It appears to give a probability distribution proportional to x**(a-1), but I

<    5   6   7   8   9   10   11   12   >