Re: [Numpy-discussion] Emulate left outer join?

2010-02-09 Thread josef . pktd
On Tue, Feb 9, 2010 at 6:52 PM, Robert Kern robert.k...@gmail.com wrote: On Tue, Feb 9, 2010 at 17:47, Ralf Gommers ralf.gomm...@googlemail.com wrote: On Wed, Feb 10, 2010 at 6:02 AM, Robert Kern robert.k...@gmail.com wrote: For some reason, numpy.lib.recfunctions isn't in the

[Numpy-discussion] numpy.polynomial.chebyshev (not) in the docs

2010-02-09 Thread josef . pktd
Similar to the recfunctions, I also don't find the new chebychev polynomials in the docs. Are they linked from any rst file? A search in the online sphinx html docs comes up empty, and http://docs.scipy.org/numpy/docs/numpy-docs/reference/routines.poly.rst/#routines-poly doesn't link to the new

Re: [Numpy-discussion] Emulate left outer join?

2010-02-09 Thread josef . pktd
On Tue, Feb 9, 2010 at 7:06 PM, Pierre GM pgmdevl...@gmail.com wrote: On Feb 9, 2010, at 6:52 PM, Robert Kern wrote: On Tue, Feb 9, 2010 at 17:47, Ralf Gommers ralf.gomm...@googlemail.com wrote: On Wed, Feb 10, 2010 at 6:02 AM, Robert Kern robert.k...@gmail.com wrote: For some reason,

Re: [Numpy-discussion] numpy.polynomial.chebyshev (not) in the docs

2010-02-09 Thread josef . pktd
On Tue, Feb 9, 2010 at 9:30 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: Are you talking about absence in the Wiki or absence in a NumPy executable. They're in the former (I've been editing them), and they're in 1.4.0 of the latter: I have them in numpy 1.4, I see them in the doceditor,

Re: [Numpy-discussion] Building Windows binaries on OS X

2010-02-08 Thread josef . pktd
On Mon, Feb 8, 2010 at 7:14 AM, Ralf Gommers ralf.gomm...@googlemail.com wrote: Hi David and all, I have a few questions on setting up the build environment on OS X for Windows binaries. I have Wine installed with Python 2.5 and 2.6, MakeNsis and MinGW. The first question is what is meant in

Re: [Numpy-discussion] Building Windows binaries on OS X

2010-02-08 Thread josef . pktd
On Mon, Feb 8, 2010 at 9:54 AM, Ralf Gommers ralf.gomm...@googlemail.com wrote: On Mon, Feb 8, 2010 at 10:25 PM, josef.p...@gmail.com wrote: On Mon, Feb 8, 2010 at 7:14 AM, Ralf Gommers ralf.gomm...@googlemail.com wrote: Hi David and all, I have a few questions on setting up the build

[Numpy-discussion] np.expand_dims is addaxis

2010-02-08 Thread josef . pktd
np.expand_dims has a name that I never remember and it's difficult to search for in the help. usage: it adds an axis e.g. after a reduce operation Please ignore, this is a message for Mr. Google Josef ___ NumPy-Discussion mailing list

Re: [Numpy-discussion] Removing datetime support for 1.4.x series ?

2010-02-06 Thread josef . pktd
On Sat, Feb 6, 2010 at 8:07 AM, Francesc Alted fal...@pytables.org wrote: A Saturday 06 February 2010 13:17:22 David Cournapeau escrigué: On Sat, Feb 6, 2010 at 4:07 PM, Travis Oliphant oliph...@enthought.com wrote: I think this plan is the least disruptive and satisfies the concerns of all

Re: [Numpy-discussion] Installed NumPy and MatPlotLib in the Wrong Order. How uninstall MPL?

2010-02-06 Thread josef . pktd
On Sat, Feb 6, 2010 at 9:02 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Sat, Feb 6, 2010 at 6:15 AM, Alan G Isaac ais...@american.edu wrote: You should be able to have Matplotlib in Python 2.5 and in Python 2.6, no problem. But you need to get the correct installer. There are

Re: [Numpy-discussion] u in [u+1]

2010-02-05 Thread josef . pktd
On Fri, Feb 5, 2010 at 8:48 AM, Zachary Pincus zachary.pin...@yale.edu wrote: I'm having some trouble here.  I have a list of numpy arrays.  I want to know if an array 'u' is in the list. Try: any(numpy.all(u == l) for l in array_list) standard caveats about float comparisons apply;

Re: [Numpy-discussion] Installed NumPy and MatPlotLib in the Wrong Order. How uninstall MPL?

2010-02-05 Thread josef . pktd
On Fri, Feb 5, 2010 at 10:37 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: See Subject. I'm working in IDLE in Win7. It seems to me MPL gets stuck in site-packages under C:\Python25. Maybe this is as simple as deleting the entry? What does it mean that MPL gets stuck? what kind of

Re: [Numpy-discussion] Removing datetime support for 1.4.x series ?

2010-02-03 Thread josef . pktd
On Wed, Feb 3, 2010 at 12:16 PM, Christopher Barker chris.bar...@noaa.gov wrote: David Cournapeau wrote: Just so that there is no confusion: it is only about removing it for 1.4.x, not about removing datetime altogether. It seems that datetime in 1.4.x has few users, Of course it has few

Re: [Numpy-discussion] Normalized histogram for data ranges 0 .. 1 returns PDF 1

2010-02-02 Thread josef . pktd
On Tue, Feb 2, 2010 at 8:05 AM, Manos Tsagias tsag...@gmail.com wrote: Hi all,  I'm using numpy.histogram with normed=True with 1D data ranging 0 .. 1. The results return probabilities greater than 1. The trapezoidal integral returns 1, but I'm afraid this is due to the bin assigned values.

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-02 Thread josef . pktd
2010/2/2 Ernest Adrogué eadro...@gmx.net:  2/02/10 @ 00:01 (-0700), thus spake Charles R Harris: On Mon, Feb 1, 2010 at 10:57 PM, josef.p...@gmail.com wrote: On Tue, Feb 2, 2010 at 12:31 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 1, 2010 at 10:02 PM,

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-02 Thread josef . pktd
On Tue, Feb 2, 2010 at 8:53 AM, josef.p...@gmail.com wrote: 2010/2/2 Ernest Adrogué eadro...@gmx.net:  2/02/10 @ 00:01 (-0700), thus spake Charles R Harris: On Mon, Feb 1, 2010 at 10:57 PM, josef.p...@gmail.com wrote: On Tue, Feb 2, 2010 at 12:31 AM, Charles R Harris

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-01 Thread josef . pktd
On Mon, Feb 1, 2010 at 12:09 PM, Keith Goodman kwgood...@gmail.com wrote: 2010/2/1 Ernest Adrogué eadro...@gmx.net: Hello, Consider the following code: for j in range(5):        f = np.bincount(x[y == j]) It fails with MemoryError whenever y == j is all False element-wise. In [96]:

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-01 Thread josef . pktd
On Mon, Feb 1, 2010 at 11:45 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 1, 2010 at 9:36 PM, David Cournapeau courn...@gmail.com wrote: On Tue, Feb 2, 2010 at 1:05 PM,  josef.p...@gmail.com wrote: I think this could be considered as a correct answer, the count of any

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-01 Thread josef . pktd
On Tue, Feb 2, 2010 at 12:31 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 1, 2010 at 10:02 PM, josef.p...@gmail.com wrote: On Mon, Feb 1, 2010 at 11:45 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 1, 2010 at 9:36 PM, David Cournapeau

Re: [Numpy-discussion] np.bincount raises MemoryError when given an empty array

2010-02-01 Thread josef . pktd
On Tue, Feb 2, 2010 at 12:57 AM, josef.p...@gmail.com wrote: On Tue, Feb 2, 2010 at 12:31 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Feb 1, 2010 at 10:02 PM, josef.p...@gmail.com wrote: On Mon, Feb 1, 2010 at 11:45 PM, Charles R Harris charlesr.har...@gmail.com wrote:

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 11:49 AM, David Huard david.hu...@gmail.com wrote: Hi, I have a 4D array with a given shape, but the array is never actually created since it is large and distributed over multiple binary files. Typical usage would be to take slices across the 4D array. I'd like to

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 12:32 PM, David Huard david.hu...@gmail.com wrote: On Fri, Jan 29, 2010 at 12:10 PM,  josef.p...@gmail.com wrote: On Fri, Jan 29, 2010 at 11:49 AM, David Huard david.hu...@gmail.com wrote: Hi, I have a 4D array with a given shape, but the array is never actually

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 12:48 PM, josef.p...@gmail.com wrote: On Fri, Jan 29, 2010 at 12:32 PM, David Huard david.hu...@gmail.com wrote: On Fri, Jan 29, 2010 at 12:10 PM,  josef.p...@gmail.com wrote: On Fri, Jan 29, 2010 at 11:49 AM, David Huard david.hu...@gmail.com wrote: Hi, I have a 4D

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 12:32 PM, David Huard david.hu...@gmail.com wrote: On Fri, Jan 29, 2010 at 12:10 PM,  josef.p...@gmail.com wrote: On Fri, Jan 29, 2010 at 11:49 AM, David Huard david.hu...@gmail.com wrote: Hi, I have a 4D array with a given shape, but the array is never actually

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 1:03 PM, Keith Goodman kwgood...@gmail.com wrote: On Fri, Jan 29, 2010 at 9:53 AM,  josef.p...@gmail.com wrote: I forgot about ellipsis, since I never use them, replace ellipsis by [slice(None)]*ndim or something like this I don't know how to access an ellipsis

Re: [Numpy-discussion] How to get the shape of an array slice without doing it

2010-01-29 Thread josef . pktd
On Fri, Jan 29, 2010 at 2:58 PM, David Huard david.hu...@gmail.com wrote: For the record, here is what I came up with. import numpy as np def expand_ellipsis(index, ndim):    Replace the ellipsis, real or implied, of an index expression by slices.    Parameters    --    index :

Re: [Numpy-discussion] Warning on http://scipy.org/ about binary incompatibility ?

2010-01-28 Thread josef . pktd
On Thu, Jan 28, 2010 at 1:39 AM, David Cournapeau da...@silveregg.co.jp wrote: Charles R Harris wrote: On Wed, Jan 27, 2010 at 6:20 PM, David Cournapeau da...@silveregg.co.jp mailto:da...@silveregg.co.jp wrote:     josef.p...@gmail.com mailto:josef.p...@gmail.com wrote:       Can

Re: [Numpy-discussion] Warning on http://scipy.org/ about binary incompatibility ?

2010-01-28 Thread josef . pktd
On Thu, Jan 28, 2010 at 9:09 AM, Dag Sverre Seljebotn da...@student.matnat.uio.no wrote: josef.p...@gmail.com wrote: On Thu, Jan 28, 2010 at 1:39 AM, David Cournapeau da...@silveregg.co.jp wrote: Charles R Harris wrote: On Wed, Jan 27, 2010 at 6:20 PM, David Cournapeau

[Numpy-discussion] Warning on http://scipy.org/ about binary incompatibility ?

2010-01-27 Thread josef . pktd
Can we/someone add a warning on the front page http://scipy.org/ (maybe under news for numpy download) about incompatibility of the binaries on sourceforge of scipy =0.7.1 with numpy 1.4.0 ? It would avoid that users have to find it out for themselves and reduce questions on the mailing list. An

Re: [Numpy-discussion] Warning on http://scipy.org/ about binary incompatibility ?

2010-01-27 Thread josef . pktd
On Wed, Jan 27, 2010 at 8:20 PM, David Cournapeau da...@silveregg.co.jp wrote: josef.p...@gmail.com wrote: Can we/someone add a warning on the front page http://scipy.org/ (maybe under news for numpy download) about incompatibility of the binaries on sourceforge of scipy =0.7.1 with numpy

Re: [Numpy-discussion] [ANN] New open source project for labeled arrays

2010-01-27 Thread josef . pktd
On Wed, Jan 27, 2010 at 10:13 PM, Pierre GM pgmdevl...@gmail.com wrote: On Jan 27, 2010, at 9:10 PM, Keith Goodman wrote: I recently opened sourced one of my packages. It is a labeled array that I call larry. A two-dimensional larry, for example, contains a 2d NumPy array with labels on each

Re: [Numpy-discussion] Warning on http://scipy.org/ about binary incompatibility ?

2010-01-27 Thread josef . pktd
On Wed, Jan 27, 2010 at 8:38 PM, David Cournapeau da...@silveregg.co.jp wrote: josef.p...@gmail.com wrote: If you are also able to provide new scipy binaries, then at least the combination would be usable without intermittent import errors and crashes. The problem is that the new scipy

Re: [Numpy-discussion] should numpy.equal raise an exception instead of returning NotImplemented?

2010-01-26 Thread josef . pktd
2010/1/26 Ernest Adrogué eadro...@gmx.net: Hi, Do you think it is sensible for np.equal to return a NotImplemented object when is given an array of variable length dtype? Consider this code: x = np.array(['xyz','zyx']) np.where(np.equal(x, 'zyx'), [0,0], [1,1]) the last line returns

Re: [Numpy-discussion] is shuffle needlessly slow?

2010-01-26 Thread josef . pktd
On Tue, Jan 26, 2010 at 2:00 PM, Alan G Isaac ais...@american.edu wrote: Is this a fair test? I expected shuffle to be much faster (no array creation). Alan Isaac import timeit setup = ... import numpy as np ... prng = np.random.RandomState() ... N = 10**5 ... indexes = np.arange(N)

Re: [Numpy-discussion] Module Index for numpy?

2010-01-25 Thread josef . pktd
On Mon, Jan 25, 2010 at 11:28 AM, denis denis-bz...@t-online.de wrote: On 17/01/2010 18:57, Wayne Watson wrote: I was just looking at the (Win) Python documentation via the Help on IDLE, and a Global Module Index. Does anything like that exist for numpy, matplotlib, scipy? Wayne, folks,  

Re: [Numpy-discussion] indexing, searchsorting, ...

2010-01-25 Thread josef . pktd
On Mon, Jan 25, 2010 at 5:16 PM, Keith Goodman kwgood...@gmail.com wrote: On Mon, Jan 25, 2010 at 1:38 PM, Jan Strube curious...@gmail.com wrote: Dear List, I'm trying to speed up a piece of code that selects a subsample based on some criteria: Setup: I have two samples, raw and cut. Cut

Re: [Numpy-discussion] fromfile() for reading text (one more time!)

2010-01-23 Thread josef . pktd
On Sat, Jan 23, 2010 at 1:53 PM, Christopher Barker chris.bar...@noaa.gov wrote: a...@ajackson.org wrote: it doesn't support the reasonable (IMHO) behavior of treating quote delimited strings in the input file as a single field. I'd use the csv module for that. Which makes me wonder if it

Re: [Numpy-discussion] Broadcasting and indexing

2010-01-22 Thread josef . pktd
On Thu, Jan 21, 2010 at 1:03 PM, Emmanuelle Gouillart emmanuelle.gouill...@normalesup.org wrote: Hi Thomas, broadcasting rules are only for ufuncs (and by extension, some numpy functions using ufuncs). Indexing obeys different rules and always starts by the first dimension. Just a

Re: [Numpy-discussion] how to do a proper 2 column sort on a 2 dimensional array ??

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 1:09 PM, Robert Kern robert.k...@gmail.com wrote: On Tue, Jan 19, 2010 at 11:53, robert somerville somer...@telus.net wrote: Hi;  i am having trouble trying to sort the rows of a 2 dimensional array by the values in the first column .. does anybody know how or have an

Re: [Numpy-discussion] Efficient orthogonalisation with scipy/numpy

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 3:47 PM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Tue, Jan 19, 2010 at 02:22:30PM -0600, Robert Kern wrote: y = y - np.dot(confounds.T, linalg.lstsq(confounds.T, y)[0]) with np = numpy and linalg = scipy.linalg where scipy calls ATLAS. For

Re: [Numpy-discussion] Efficient orthogonalisation with scipy/numpy

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 4:11 PM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: you can  combine x, and y for one call to leastsq, if it makes a difference linalg.lstsq(confounds.T, [x,y])  #format? columnstack? Indeed! Thank you Joseph. That's a gain of 10% in the total computation time

Re: [Numpy-discussion] Efficient orthogonalisation with scipy/numpy

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 4:29 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Jan 19, 2010 at 1:47 PM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Tue, Jan 19, 2010 at 02:22:30PM -0600, Robert Kern wrote: y = y - np.dot(confounds.T, linalg.lstsq(confounds.T, y)[0])

Re: [Numpy-discussion] Efficient orthogonalisation with scipy/numpy

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 6:48 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Jan 19, 2010 at 2:34 PM, josef.p...@gmail.com wrote: On Tue, Jan 19, 2010 at 4:29 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Jan 19, 2010 at 1:47 PM, Gael Varoquaux

Re: [Numpy-discussion] Efficient orthogonalisation with scipy/numpy

2010-01-19 Thread josef . pktd
On Tue, Jan 19, 2010 at 9:47 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Jan 19, 2010 at 6:08 PM, josef.p...@gmail.com wrote: On Tue, Jan 19, 2010 at 6:48 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Tue, Jan 19, 2010 at 2:34 PM, josef.p...@gmail.com

Re: [Numpy-discussion] logic problem

2010-01-18 Thread josef . pktd
2010/1/18 Ernest Adrogué eadro...@gmx.net: Hi, This is hard to explain. In this code: reduce(np.logical_or, [m1 m2, m1 m3, m2 m3]) where m1, m2 and m3 are boolean arrays, I'm trying to figure out an expression that works with an arbitrary number of arrays, not just 3. Any idea??

Re: [Numpy-discussion] logic problem

2010-01-18 Thread josef . pktd
On Mon, Jan 18, 2010 at 2:18 PM, Warren Weckesser warren.weckes...@enthought.com wrote: Ernest Adrogué wrote: Hi, This is hard to explain. In this code: reduce(np.logical_or, [m1 m2, m1 m3, m2 m3]) where m1, m2 and m3 are boolean arrays, I'm trying to figure out an expression that

Re: [Numpy-discussion] Module Index for numpy?

2010-01-17 Thread josef . pktd
On Sun, Jan 17, 2010 at 2:19 PM, David Goldsmith d.l.goldsm...@gmail.com wrote: Hi, Wayne. They're not nearly as structured, but for the time being (indefinitely? unless a volunteer steps forward to build something for us more closely resembling the GMI), you could use the numpy and scipy

Re: [Numpy-discussion] Waf or scons/numscons for a C/Fortran/Cython/Python project -- what's your recommendation?

2010-01-16 Thread josef . pktd
On Sat, Jan 16, 2010 at 4:43 PM, Dag Sverre Seljebotn da...@student.matnat.uio.no wrote: Kurt Smith wrote: On Sat, Jan 16, 2010 at 2:38 PM, Dag Sverre Seljebotn da...@student.matnat.uio.no wrote: Not that I really know anything about it, but note that one of the purposes of David's toydist

Re: [Numpy-discussion] svn log and blank entries

2010-01-15 Thread josef . pktd
On Fri, Jan 15, 2010 at 9:50 AM, Nils Wagner nwag...@iam.uni-stuttgart.de wrote: Hi all, An svn log CHANGELOG in svn/numpy yields some blank entries Is that intended ? r8055 | ariver | 2010-01-15 03:02:30 +0100 (Fr,

[Numpy-discussion] linalg.eig getting the original matrix back ?

2010-01-15 Thread josef . pktd
I had a problem because linal.eig doesn't rebuild the original matrix, linalg.eigh does, see script below Whats the trick with linalg.eig to get the original (or the inverse) back ? None of my variations on the formulas worked. Thanks, Josef import numpy as np import scipy as sp import

Re: [Numpy-discussion] linalg.eig getting the original matrix back ?

2010-01-15 Thread josef . pktd
On Fri, Jan 15, 2010 at 11:32 AM, Sebastian Walter sebastian.wal...@gmail.com wrote: numpy.linalg.eig guarantees to return right eigenvectors. evec is not necessarily an orthonormal matrix when there are eigenvalues with multiplicity 1. For symmetrical matrices you'll have mutually orthogonal

Re: [Numpy-discussion] linalg.eig getting the original matrix back ?

2010-01-15 Thread josef . pktd
On Fri, Jan 15, 2010 at 12:07 PM, josef.p...@gmail.com wrote: On Fri, Jan 15, 2010 at 11:32 AM, Sebastian Walter sebastian.wal...@gmail.com wrote: numpy.linalg.eig guarantees to return right eigenvectors. evec is not necessarily an orthonormal matrix when there are eigenvalues with

Re: [Numpy-discussion] comparison operators (e.g. ==) on array with dtype object do not work

2010-01-14 Thread josef . pktd
On Thu, Jan 14, 2010 at 5:49 PM, Warren Weckesser warren.weckes...@enthought.com wrote: Yaroslav Halchenko wrote: Dear NumPy People, First I want to apologize if I misbehaved on NumPy Trac by reopening the closed ticket http://projects.scipy.org/numpy/ticket/1362 but I still feel strongly

Re: [Numpy-discussion] numpy1.4 dtype issues: scipy.stats pytables

2010-01-12 Thread josef . pktd
On Tue, Jan 12, 2010 at 11:33 AM, denis denis-bz...@t-online.de wrote: On 11/01/2010 18:10, josef.p...@gmail.com wrote: For this problem, it's supposed to be only those packages that have or import cython generated code. Right; is this a known bug, is there a known fix  for mac dmgs ?

Re: [Numpy-discussion] wrong casting of augmented assignment statements

2010-01-12 Thread josef . pktd
On Tue, Jan 12, 2010 at 1:05 PM, Sebastian Walter sebastian.wal...@gmail.com wrote: Hello, I have a question about the augmented assignment statements *=, +=, etc. Apparently, the casting of types is not working correctly. Is this known resp. intended behavior of numpy? (I'm using

Re: [Numpy-discussion] numpy sum table by category

2010-01-12 Thread josef . pktd
On Tue, Jan 12, 2010 at 3:33 PM, Marc Schwarzschild m...@thebrookhavengroup.com wrote: I have a csv file like this:    Account, Symbol, Quantity, Price    One,SPY,5,119.00    One,SPY,3,120.00    One,SPY,-2,125.00    One,GE,...    One,GE,...    Two,SPY, ...    Three,GE, ...     ...

Re: [Numpy-discussion] numpy1.4 dtype issues: scipy.stats pytables

2010-01-11 Thread josef . pktd
On Mon, Jan 11, 2010 at 11:56 AM, denis denis-bz...@t-online.de wrote: Only 2 of the 21 top-level subpackages draw that warning with numpy-1.4.0-py2.6-python.org.dmg scipy-0.7.1-py2.6-python.org.dmg on my mac 10.4 ppc, python 2.6.4: try:     import scipy.cluster except ValueError, e:    

Re: [Numpy-discussion] TypeError: 'module' object is not callable

2010-01-11 Thread josef . pktd
On Mon, Jan 11, 2010 at 7:16 PM, Jankins andyjian430...@gmail.com wrote: I am sorry. My bad.   File C:\test.py, line 7, in module     print linalg.eigen(M) TypeError: 'module' object is not callable I installed pythonxy. pythonxy has already included the scipy package. On 1/11/2010 6:12

Re: [Numpy-discussion] TypeError: 'module' object is not callable

2010-01-11 Thread josef . pktd
On Mon, Jan 11, 2010 at 9:03 PM, Jankins andyjian430...@gmail.com wrote: Here is the command line python:   import scipy.sparse.linalg as linalg     linalg.eigen() Traceback (most recent call last):   File stdin, line 1, in module TypeError: 'module' object is not callable  

Re: [Numpy-discussion] TypeError: 'module' object is not callable

2010-01-11 Thread josef . pktd
On Mon, Jan 11, 2010 at 10:31 PM, Jankins andyjian430...@gmail.com wrote: linalg has no attribute eigen. You should post full tracebacks. I don't understand this error, because before eigen seemed to exist. You could run the test suite to see if the installation is ok and sparse is working

Re: [Numpy-discussion] fromfile() for reading text (one more time!)

2010-01-07 Thread josef . pktd
On Thu, Jan 7, 2010 at 3:08 PM, Christopher Barker chris.bar...@noaa.gov wrote: Pauli Virtanen wrote: ma, 2010-01-04 kello 17:05 -0800, Christopher Barker kirjoitti: it also does odd things with spaces embedded in the separator: , $ # matches all of:  ,$#   , $#  ,$ # That's a documented

Re: [Numpy-discussion] fromfile() for reading text (one more time!)

2010-01-07 Thread josef . pktd
On Thu, Jan 7, 2010 at 4:45 PM, Christopher Barker chris.bar...@noaa.gov wrote: Bruce Southey wrote: chris.bar...@noaa.gov wrote: Using the numpy NaN or similar (noting R's approach to missing values which in turn allows it to have the above functionality) is just a very bad idea for missing

Re: [Numpy-discussion] fromfile() for reading text (one more time!)

2010-01-07 Thread josef . pktd
On Thu, Jan 7, 2010 at 11:10 PM, Bruce Southey bsout...@gmail.com wrote: On Thu, Jan 7, 2010 at 3:45 PM, Christopher Barker chris.bar...@noaa.gov wrote: Bruce Southey wrote: chris.bar...@noaa.gov wrote: Using the numpy NaN or similar (noting R's approach to missing values which in turn

Re: [Numpy-discussion] fromfile() for reading text (one more time!)

2010-01-04 Thread josef . pktd
On Mon, Jan 4, 2010 at 10:39 PM, a...@ajackson.org wrote: Hi folks, I'm taking a look once again at fromfile() for reading text files. I often have the need to read a LOT of numbers form a text file, and it can actually be pretty darn slow do i the normal python way: for line in file:    data

Re: [Numpy-discussion] Complex slicing and take

2009-12-30 Thread josef . pktd
On Wed, Dec 30, 2009 at 12:19 PM, Keith Goodman kwgood...@gmail.com wrote: On Wed, Dec 30, 2009 at 12:08 AM, Eric Emsellem eemse...@eso.org wrote: Hi thanks for the tips. Unfortunately this is not what I am after. ? import numpy as num ? startarray = random((1000,100)) ? take_sample =

Re: [Numpy-discussion] Announcing toydist, improving distribution and packaging situation

2009-12-30 Thread josef . pktd
On Wed, Dec 30, 2009 at 1:47 PM, René Dudfield ren...@gmail.com wrote: On Wed, Dec 30, 2009 at 2:26 PM, Ravi lists_r...@lavabit.com wrote: On Wednesday 30 December 2009 06:15:45 René Dudfield wrote: I agree with many things in that post.  Except your conclusion on multiple versions of

Re: [Numpy-discussion] [matplotlib-devel] Announcing toydist, improving distribution and packaging situation

2009-12-29 Thread josef . pktd
On Tue, Dec 29, 2009 at 10:55 AM, Gael Varoquaux gael.varoqu...@normalesup.org wrote: On Tue, Dec 29, 2009 at 11:34:44PM +0900, David Cournapeau wrote: Buildout, virtualenv all work by sandboxing from the system python: each of them do not see each other, which may be useful for development,

Re: [Numpy-discussion] nicest way to apply an arbitrary sequence of row deltas to an array

2009-12-21 Thread josef . pktd
On Mon, Dec 21, 2009 at 6:35 AM, Pauli Virtanen pav...@iki.fi wrote: Mon, 21 Dec 2009 09:35:08 +, Neil wrote: [clip] I'm also interested to see if there are any answers to this; I came across a similar problem recently. It would have been convenient to do something like U[idx] += dU, but

Re: [Numpy-discussion] cos -- NameError: global name 'cos' is not defined

2009-12-21 Thread josef . pktd
On Mon, Dec 21, 2009 at 6:11 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Yes, thanks. That's the what I finally changed to. This originated up a thread or so when I displayed the highly populated code with math. Some said I didn't need it, so I thought I'd give it a go. I just

Re: [Numpy-discussion] cos -- NameError: global name 'cos' is not defined

2009-12-21 Thread josef . pktd
On Mon, Dec 21, 2009 at 8:25 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: I may have inadvertently made a slip between using script versus shell. What I'm getting at it that the namespace is the same for both the editor window and shell window. I find that a little bizarre. I would

Re: [Numpy-discussion] indexing question

2009-12-20 Thread josef . pktd
On Sun, Dec 20, 2009 at 8:58 PM, Alan G Isaac ais...@american.edu wrote: Why is s3 F_CONTIGUOUS, and perhaps equivalently, why is its C_CONTIGUOUS data in s3.base (below)? Thanks, Alan Isaac a3 array([[ 0,  1,  2,  3,  4,  5],        [ 6,  7,  8,  9, 10, 11]]) a3.flags   C_CONTIGUOUS :

Re: [Numpy-discussion] Segmentation fault with argsort

2009-12-18 Thread josef . pktd
On Fri, Dec 18, 2009 at 1:01 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Fri, Dec 18, 2009 at 10:57 AM, Skipper Seabold jsseab...@gmail.com wrote: On Fri, Dec 18, 2009 at 12:52 PM, Robert Kern robert.k...@gmail.com wrote: On Fri, Dec 18, 2009 at 11:46, Keith Goodman

Re: [Numpy-discussion] Proposal for matrix_rank function in numpy

2009-12-18 Thread josef . pktd
On Fri, Dec 18, 2009 at 10:21 PM, Fernando Perez fperez@gmail.com wrote: On Wed, Dec 16, 2009 at 10:56 AM, Skipper Seabold jsseab...@gmail.com wrote: Presumably the doctests should be turned into actual tests (noting Robert's comment) to make it more likely that it gets in Just curious:

Re: [Numpy-discussion] ValueError for numpy when importing KDTree from scipy.spatial

2009-12-18 Thread josef . pktd
On Sat, Dec 19, 2009 at 12:42 AM, Howard Chong hgch...@berkeley.edu wrote: I'm getting an odd behavior when I try to load KDTree. In the interactive interpreter: the first time I try to load it, it gives me an error. Second time works fine. When trying to run from command line, same error. To

Re: [Numpy-discussion] Proposal for matrix_rank function in numpy

2009-12-15 Thread josef . pktd
On Tue, Dec 15, 2009 at 12:01 PM, Matthew Brett matthew.br...@gmail.com wrote: Hi, Following on from the occasional discussion on the list, can I propose a small matrix_rank function for inclusion in numpy/linalg? I suggest it because it seems rather a basic need for linear algebra, and

Re: [Numpy-discussion] Slicing slower than matrix multiplication?

2009-12-14 Thread josef . pktd
On Mon, Dec 14, 2009 at 12:51 PM, Francesc Alted fal...@pytables.org wrote: A Monday 14 December 2009 18:20:32 Jasper van de Gronde escrigué: Francesc Alted wrote: A Monday 14 December 2009 17:09:13 Francesc Alted escrigué: The things seems to be worst than 1.6x times slower for numpy, as

Re: [Numpy-discussion] Question on timeseries, for financial application

2009-12-13 Thread josef . pktd
On Sun, Dec 13, 2009 at 3:31 AM, Pierre GM pgmdevl...@gmail.com wrote: On Dec 13, 2009, at 12:11 AM, Robert Ferrell wrote: Have you considered creating a TimeSeries for each data series, and then putting them all together in a dict, keyed by symbol? That's an idea As far as I understand,

Re: [Numpy-discussion] Repeated dot products

2009-12-12 Thread josef . pktd
On Sat, Dec 12, 2009 at 4:55 PM, T J tjhn...@gmail.com wrote: Hi, Suppose I have an array of shape:  (n, k, k).  In this case, I have n k-by-k matrices.  My goal is to compute the product of a (potentially large) user-specified selection (with replacement) of these matrices. For example,  

Re: [Numpy-discussion] Question on timeseries, for financial application

2009-12-12 Thread josef . pktd
On Sat, Dec 12, 2009 at 8:08 PM, THOMAS BROWNE totalb...@mac.com wrote: Hello all, Quite new to numpy / timeseries module, please forgive the elementary question. I wish to do quite to do a bunch of multivariate analysis on 1000 different financial markets series, each holding about 1800

Re: [Numpy-discussion] Assigning complex values to a real array

2009-12-09 Thread josef . pktd
On Wed, Dec 9, 2009 at 9:46 AM, Ryan May rma...@gmail.com wrote: On Wed, Dec 9, 2009 at 3:51 AM, David Warde-Farley d...@cs.toronto.edu wrote: On 9-Dec-09, at 1:26 AM, Dr. Phillip M. Feldman wrote: Unfortunately, NumPy seems to be a sort of step-child of Python, tolerated, but not fully

Re: [Numpy-discussion] Release blockers for 1.4.0 ?

2009-12-08 Thread josef . pktd
On Tue, Dec 8, 2009 at 10:21 AM, David Cournapeau courn...@gmail.com wrote: On Wed, Dec 9, 2009 at 12:10 AM, Bruce Southey bsout...@gmail.com wrote: I understand the reason for the masked arrays behavior but changing the seterr default will be a problem until ma is changed. With Python 2.6

Re: [Numpy-discussion] slices of structured arrays?

2009-12-08 Thread josef . pktd
2009/12/8 Ernest Adrogué eadro...@gmx.net: Hello, Here's a structured array with fields 'a','b' and 'c': s=[(i,int) for i in 'abc'] t=np.zeros(1,s) It has the form: array([(0, 0, 0)] I was wondering if such an array can be accessed like an ordinary array (e.g., with a slice) in order to

Re: [Numpy-discussion] Flattening an array

2009-12-08 Thread josef . pktd
On Tue, Dec 8, 2009 at 7:29 PM, Jake VanderPlas jake...@gmail.com wrote: Hello, I have a function -- call it f() -- which takes a length-N 1D numpy array as an argument, and returns a length-N 1D array. I want to pass it the data in an N-D array, and obtain the N-D array of the result. I've

Re: [Numpy-discussion] Release blockers for 1.4.0 ?

2009-12-07 Thread josef . pktd
On Mon, Dec 7, 2009 at 1:24 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Dec 7, 2009 at 11:16 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Dec 7, 2009 at 10:31 AM, David Cournapeau courn...@gmail.com wrote: On Tue, Dec 8, 2009 at 1:48 AM, Charles R

Re: [Numpy-discussion] Release blockers for 1.4.0 ?

2009-12-07 Thread josef . pktd
On Mon, Dec 7, 2009 at 11:22 PM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: josef.p...@gmail.com wrote: On Mon, Dec 7, 2009 at 1:24 PM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Dec 7, 2009 at 11:16 AM, Charles R Harris charlesr.har...@gmail.com wrote: On Mon, Dec

Re: [Numpy-discussion] Release blockers for 1.4.0 ?

2009-12-07 Thread josef . pktd
On Mon, Dec 7, 2009 at 11:54 PM, David Cournapeau da...@ar.media.kyoto-u.ac.jp wrote: josef.p...@gmail.com wrote: warning is no problem, but I haven't figured out what the pattern is for repeated warnings. By default, warnings are only 'raised' once. You need to use the standard warnings

Re: [Numpy-discussion] histogram for discrete data

2009-12-06 Thread josef . pktd
2009/12/6 Ernest Adrogué eadro...@gmx.net: Hi, A few weeks ago there was a discussion about a histogram_discrete() function --sorry for starting a new thread but I have lost the mails. Somebody pointed out that bincount() already can be used to histogram discrete data (except that it

[Numpy-discussion] np.equal

2009-12-06 Thread josef . pktd
what's the difference in the implementation between np.equal and == ? np.equal raises NotImplemented for strings, while == works. aa array(['a', 'b', 'a', 'aa', 'a'], dtype='|S2') aa == 'a' array([ True, False, True, False, True], dtype=bool) np.equal(aa,'a') NotImplemented

Re: [Numpy-discussion] non-standard standard deviation

2009-12-06 Thread josef . pktd
On Sun, Dec 6, 2009 at 11:01 AM, Colin J. Williams c...@ncf.ca wrote: On 04-Dec-09 10:54 AM, Bruce Southey wrote: On 12/04/2009 06:18 AM, yogesh karpate wrote: @ Pauli and @ Colin:                                   Sorry for the late reply. I was busy in some other assignments. # As far as

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

2009-11-28 Thread josef . pktd
On Sat, Nov 28, 2009 at 6:18 AM, Sebastian seb...@gmail.com wrote: On Sat, Nov 28, 2009 at 1:01 AM, josef.p...@gmail.com wrote: On Fri, Nov 27, 2009 at 9:44 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Joseph, That got it by the fig problem but there is yet another one. value is

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

2009-11-27 Thread josef . pktd
On Fri, Nov 27, 2009 at 8:43 AM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Thanks. That sounds like it should help a lot. Finding meaningful examples anywhere hasn't been easy. I thought I'd look through Amazon for books on Python and scientific uses. I found almost all were written by

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

2009-11-27 Thread josef . pktd
On Fri, Nov 27, 2009 at 12:57 PM, Skipper Seabold jsseab...@gmail.com wrote: On Fri, Nov 27, 2009 at 12:41 PM, Christopher Barker chris.bar...@noaa.gov wrote: Wayne Watson wrote: Yes, I'm just beginning to deal with the contents of NumPy, SciLab, and SciPy. They all have seemed part of one

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

2009-11-27 Thread josef . pktd
On Fri, Nov 27, 2009 at 9:05 PM, Sebastian seb...@gmail.com wrote: Hi Chris, yeah there should, try the following: import numpy import matplotlib.pyplot as pylab regards On Fri, Nov 27, 2009 at 8:47 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: I tried this and it put ranges on y

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

2009-11-27 Thread josef . pktd
On Fri, Nov 27, 2009 at 9:44 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Joseph, That got it by the fig problem but there is yet another one. value is not defined on the very long line: range = ...    Wayne (values is the data array, ... no idea about

Re: [Numpy-discussion] Computing Simple Statistics When Only they Frequency Distribution is Known

2009-11-27 Thread josef . pktd
On Fri, Nov 27, 2009 at 9:47 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: How do I compute avg, std dev, min, max and other simple stats if I only know the frequency distribution? If you are willing to assign to all observations in a bin the value at the bin midpoint, then you could do

[Numpy-discussion] matrix inverse

2009-11-26 Thread josef . pktd
why is http://docs.scipy.org/numpy/docs/numpy.matrixlib.defmatrix.matrix.I/ classified as unimportant? .T ? I was looking at http://projects.scipy.org/numpy/ticket/1093 and didn't find any docs for the matrix inverse. Josef ___ NumPy-Discussion

Re: [Numpy-discussion] matrix inverse

2009-11-26 Thread josef . pktd
On Thu, Nov 26, 2009 at 11:34 AM, josef.p...@gmail.com wrote: why is http://docs.scipy.org/numpy/docs/numpy.matrixlib.defmatrix.matrix.I/ classified as unimportant?   .T ? I was looking at http://projects.scipy.org/numpy/ticket/1093 and didn't find any docs for the matrix inverse. Ok

Re: [Numpy-discussion] matrix inverse

2009-11-26 Thread josef . pktd
On Thu, Nov 26, 2009 at 11:49 AM, josef.p...@gmail.com wrote: On Thu, Nov 26, 2009 at 11:34 AM,  josef.p...@gmail.com wrote: why is http://docs.scipy.org/numpy/docs/numpy.matrixlib.defmatrix.matrix.I/ classified as unimportant?   .T ? I was looking at

Re: [Numpy-discussion] NumPy Histogram for Tentative NumPy Tutorial Questions

2009-11-26 Thread josef . pktd
On Thu, Nov 26, 2009 at 2:44 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: I decided to try some example code from Subject. import numpy import pylab # Build a vector of 1 normal deviates with variance 0.5^2 and mean 2 mu, sigma = 2, 0.5 v = numpy.random.normal(mu,sigma,1) #

Re: [Numpy-discussion] NumPy Histogram for Tentative NumPy Tutorial Questions

2009-11-26 Thread josef . pktd
On Thu, Nov 26, 2009 at 8:18 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Yes, I'm just beginning to deal with the contents of NumPy, SciLab, and SciPy. They all have seemed part of one another, but I think I see how they've divided up the game. NumPy has no graphics. PyLab has them,

Re: [Numpy-discussion] NumPy Histogram for Tentative NumPy Tutorial Questions

2009-11-26 Thread josef . pktd
On Thu, Nov 26, 2009 at 9:48 PM, josef.p...@gmail.com wrote: On Thu, Nov 26, 2009 at 8:18 PM, Wayne Watson sierra_mtnv...@sbcglobal.net wrote: Yes, I'm just beginning to deal with the contents of NumPy, SciLab, and SciPy. They all have seemed part of one another, but I think I see how

<    4   5   6   7   8   9   10   11   12   >