Re: [Numpy-discussion] Forbidden charcter in the "names" argument of genfromtxt?

2012-02-20 Thread Skipper Seabold
On Mon, Feb 20, 2012 at 1:35 PM, Brett Olsen wrote: > On Sat, Feb 18, 2012 at 8:12 PM, Adam Hughes wrote: >> Hey everyone, >> >> I have timeseries data in which the column label is simply a filename from >> which the original data was taken.  Here's some sample data: >> >> name1.txt  name2.txt  n

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

2012-02-25 Thread Skipper Seabold
On Sat, Feb 25, 2012 at 5:13 PM, 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

[Numpy-discussion] speed of array creation from tuples

2012-02-27 Thread Skipper Seabold
I am surprised by this (though maybe I shouldn't be?) It's always faster to use list comprehension to unpack lists of tuples than np.array/asarray? [~/] [1]: X = [tuple(np.random.randint(10,size=2)) for _ in range(100)] [~/] [2]: timeit np.array([x1 for _,x1 in X]) 1 loops, best of 3: 26.4 us

Re: [Numpy-discussion] Missing data again

2012-03-03 Thread Skipper Seabold
On Sat, Mar 3, 2012 at 4:46 PM, Mark Wiebe wrote: > On Sat, Mar 3, 2012 at 12:30 PM, Travis Oliphant >> >>        * the reduction operations need to default to "skipna" --- this is >> the most common use case which has been re-inforced again to me today by a >> new user to Python who is using ma

[Numpy-discussion] checking for c compiler during build

2012-03-07 Thread Skipper Seabold
Is there a way to use numpy.distuils to programmatically check for a C compiler at build time in a platform independent way? ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] checking for c compiler during build

2012-03-07 Thread Skipper Seabold
On Wed, Mar 7, 2012 at 12:35 PM, Skipper Seabold wrote: > Is there a way to use numpy.distuils to programmatically check for a C > compiler at build time in a platform independent way? Wading through the numpy/distutils code some more. Would something as simple as this work all the time?

Re: [Numpy-discussion] numpy videos

2012-03-12 Thread Skipper Seabold
On Mon, Mar 12, 2012 at 6:04 PM, Abhishek Pratap wrote: > > Hey Guys > > Few days with folks at my first pycon has made me wonder how much of > cool things I was missing .. > > I am looking to do some quick catch up on numpy and wondering if there > are any set of videos that I can refer to. I lea

Re: [Numpy-discussion] numpy doc for percentile function

2012-04-04 Thread Skipper Seabold
On Wed, Apr 4, 2012 at 6:10 AM, Pierre Haessig wrote: > Hi, > > I'm looking for the entry point in Numpy doc for the percentile function. > I'm assuming it should sit in routines.statistics but do not see it : > http://docs.scipy.org/doc/numpy/reference/routines.statistics.html > > I don't see it

[Numpy-discussion] adding a cut function to numpy

2012-04-16 Thread Skipper Seabold
Hi, I have a pull request here [1] to add a cut function similar to R's [2]. It seems there are often requests for similar functionality. It's something I'm making use of for my own work and would like to use in statstmodels and in generating instances of pandas' Factor class, but is this generall

Re: [Numpy-discussion] adding a cut function to numpy

2012-04-16 Thread Skipper Seabold
On Mon, Apr 16, 2012 at 5:51 PM, Tony Yu wrote: > > > On Mon, Apr 16, 2012 at 5:27 PM, Skipper Seabold > wrote: >> >> Hi, >> >> I have a pull request here [1] to add a cut function similar to R's >> [2]. It seems there are often requests for simi

Re: [Numpy-discussion] adding a cut function to numpy

2012-04-16 Thread Skipper Seabold
On Mon, Apr 16, 2012 at 8:08 PM, Tony Yu wrote: > > > On Mon, Apr 16, 2012 at 6:01 PM, Skipper Seabold > wrote: >> >> On Mon, Apr 16, 2012 at 5:51 PM, Tony Yu wrote: >> > >> > >> > On Mon, Apr 16, 2012 at 5:27 PM, Skipper Seabold >> &

Re: [Numpy-discussion] linalg.lstsq

2012-04-25 Thread Skipper Seabold
On Wed, Apr 25, 2012 at 12:13 PM, Till Stensitzki wrote: > Hello, > is there weighted version of linalg.lstsq available? > In my case, b is a (N,K) matrix, so i can't use manual scaling of x and b. > What shape are the weights in this case? I'm not that familiar with problems with an N x K b matr

Re: [Numpy-discussion] Scipy build can't find BLAS when using numpy master? (Was: Should arr.diagonal() return a copy or a view? (1.7 compatibility issue))

2012-05-16 Thread Skipper Seabold
On Wed, May 16, 2012 at 11:50 AM, Robert Kern wrote: > On Wed, May 16, 2012 at 4:35 PM, Nathaniel Smith wrote: >> On Wed, May 16, 2012 at 4:24 PM, Robert Kern wrote: >>> On Wed, May 16, 2012 at 4:21 PM, Nathaniel Smith wrote: >>> I built some pristine python 2.7 installs from scratch (no v

Re: [Numpy-discussion] Fwd: Named dtype array: Difference between a[0]['name'] and a['name'][0]?

2012-05-21 Thread Skipper Seabold
On Mon, May 21, 2012 at 4:37 PM, Travis Oliphant wrote: > This is the right place to ask, it's just that it can take time to get an > answer because people who might know the answer may not have the time to > respond immediately. > > The short answer is that this is not really a "normal" bug, but

Re: [Numpy-discussion] NumPy 1.7 release delays

2012-06-26 Thread Skipper Seabold
On Tue, Jun 26, 2012 at 7:59 PM, Fernando Perez wrote: > On Tue, Jun 26, 2012 at 1:10 PM, Travis Oliphant wrote: >> One issues is the one that Sage identified about the array interface >> regression as noted by Jason.    Any other regressions from 1.5.x need to be >> addressed as well.    We'll h

Re: [Numpy-discussion] Multivariate hypergeometric distribution?

2012-07-02 Thread Skipper Seabold
On Mon, Jul 2, 2012 at 9:35 PM, wrote: > On Mon, Jul 2, 2012 at 8:08 PM, wrote: > > On Mon, Jul 2, 2012 at 4:16 PM, Fernando Perez > wrote: > >> Hi all, > >> > >> in recent work with a colleague, the need came up for a multivariate > >> hypergeometric sampler; I had a look in the numpy code an

Re: [Numpy-discussion] view of recarray issue

2012-07-25 Thread Skipper Seabold
On Sun, Jul 22, 2012 at 2:15 PM, Ralf Gommers wrote: > Hi, > > Just a heads up that right now views of recarrays seem to be problematic, > this doesn't work anymore: > import statsmodels.api as sm dta = sm.datasets.macrodata.load() # returns a record array with 14 fields dta.d

Re: [Numpy-discussion] ANN: NumPy 1.7.0b1 release

2012-08-21 Thread Skipper Seabold
On Tue, Aug 21, 2012 at 3:59 PM, Christoph Gohlke wrote: > On 8/21/2012 9:24 AM, Ondřej Čertík wrote: > > Hi, > > > > I'm pleased to announce the availability of the first beta release of > > NumPy 1.7.0b1. > > > > Sources and binary installers can be found at > > https://sourceforge.net/projects

Re: [Numpy-discussion] the mean, var, std of empty arrays

2012-11-21 Thread Skipper Seabold
On Wed, Nov 21, 2012 at 9:22 PM, Olivier Delalleau wrote: > Current behavior looks sensible to me. I personally would prefer no > warning but I think it makes sense to have one as it can be helpful to > detect issues faster. > > -=- Olivier > It's configurable. [~/] [1]: np.seterr(all='ignore')

[Numpy-discussion] result shape from dot for 0d, 1d, 2d scalar

2012-11-26 Thread Skipper Seabold
I discovered this because scipy.optimize.fmin_powell appears to squeeze 1d argmin to 0d unlike the other optimizers, but that's a different story. I would expect the 0d array to behave like the 1d array not the 2d as it does below. Thoughts? Maybe too big of a pain to change this behavior if indee

Re: [Numpy-discussion] result shape from dot for 0d, 1d, 2d scalar

2012-11-28 Thread Skipper Seabold
On Tue, Nov 27, 2012 at 11:16 AM, Sebastian Berg wrote: > On Mon, 2012-11-26 at 13:54 -0500, Skipper Seabold wrote: > > I discovered this because scipy.optimize.fmin_powell appears to > > squeeze 1d argmin to 0d unlike the other optimizers, but that's a > > different

Re: [Numpy-discussion] result shape from dot for 0d, 1d, 2d scalar

2012-11-28 Thread Skipper Seabold
On Wed, Nov 28, 2012 at 12:31 PM, Sebastian Berg wrote: > Maybe a strict matrix product would make sense too, but the dot function > behavior cannot be changed in any case, so its pointless to argue about > it. Just make sure your arrays are 2-d (or matrices) if you want a > matrix product, which

Re: [Numpy-discussion] Proposal to drop python 2.4 support in numpy 1.8

2012-12-13 Thread Skipper Seabold
On Thu, Dec 13, 2012 at 12:00 PM, David Cournapeau wrote: > I would even go as far as dropping 2.5 as well then (RHEL 6 > uses python 2.6). +1 Skipper ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/

Re: [Numpy-discussion] New numpy functions: filled, filled_like

2013-01-13 Thread Skipper Seabold
On Sun, Jan 13, 2013 at 6:39 PM, Nathaniel Smith wrote: > On Sun, Jan 13, 2013 at 11:24 PM, Robert Kern > wrote: > > On Sun, Jan 13, 2013 at 6:27 PM, Nathaniel Smith wrote: > >> Hi all, > >> > >> PR 2875 adds two new functions, that generalize zeros(), ones(), > >> zeros_like(), ones_like(), by

Re: [Numpy-discussion] Adding .abs() method to the array object

2013-02-25 Thread Skipper Seabold
On Mon, Feb 25, 2013 at 10:43 AM, Till Stensitzki wrote: > > First, sorry that i didnt search for an old thread, but because i disagree with > conclusion i would at least address my reason: > >> I don't like >> np.abs(arr).max() >> because I have to concentrate to much on the braces, especially if

Re: [Numpy-discussion] should get rid of the annoying numpy STDERR output

2011-03-24 Thread Skipper Seabold
On Thu, Mar 24, 2011 at 9:45 AM, Ralf Gommers wrote: > 2011/3/24 Dmitrey : >> >> Hi >> >> 2011/3/24 Dmitrey >>> >>> >>> from numpy import inf, array >>> >>> inf*0 >>> nan >>> >>> (ok) >>> >>> >>> array(inf) * 0.0 >>> StdErr: Warning: invalid value encountered in multiply >>> nan >>> >>> My cycled

Re: [Numpy-discussion] should get rid of the annoying numpy STDERR output

2011-03-24 Thread Skipper Seabold
On Thu, Mar 24, 2011 at 10:52 AM, eat wrote: > Hi > > On Thu, Mar 24, 2011 at 4:17 PM, Skipper Seabold >> Also, as Robert pointed out to me before np.errstate is a >> context-manager for ignoring these warnings. I often wrap optimization >> code with it. > >

[Numpy-discussion] "Best" dot product along an axis?

2011-04-04 Thread Skipper Seabold
Going through some of the recent threads on similar problems, I'm trying to discern which is "best." I have X is T x i, theta is T x i x j. I want a T by j array that contains X[t].T.dot(theta[t]) along axis 0. Say, T,i,j = 166,7,3 X = np.random.random((T,i)) theta = np.random.random((T,i,j))

Re: [Numpy-discussion] Remove deprecated skiprows

2011-04-05 Thread Skipper Seabold
On Sun, Apr 3, 2011 at 8:20 PM, Charles R Harris wrote: > Should skiprows be removed? > >     if skiprows: >     warnings.warn(\ >     "The use of `skiprows` is deprecated, it will be removed in > numpy 2.0.\n" \ >     "Please use `skip_header` instead.", >     Deprecat

Re: [Numpy-discussion] Remove deprecated skiprows

2011-04-05 Thread Skipper Seabold
On Tue, Apr 5, 2011 at 6:05 PM, Pierre GM wrote: > > On Apr 5, 2011, at 11:52 PM, Ralf Gommers wrote: > >> On Tue, Apr 5, 2011 at 11:45 PM, Skipper Seabold wrote: >>> On Sun, Apr 3, 2011 at 8:20 PM, Charles R Harris >>> wrote: >>>> Should s

[Numpy-discussion] Odd numerical difference between Numpy 1.5.1 and Numpy > 1.5.1

2011-04-11 Thread Skipper Seabold
All, We noticed some failing tests for statsmodels between numpy 1.5.1 and numpy >= 1.6.0. These are the versions where I noticed the change. It seems that when you divide a float array and multiply by a boolean array the answers are different (unless the others are also off by some floating point

Re: [Numpy-discussion] Bug: np.linalg.qr crash the python interpreter if getting a shape (0, x) array.

2011-05-08 Thread Skipper Seabold
On Fri, May 6, 2011 at 11:30 AM, Till Stensitzki wrote: > Hi, > discovered another small bug. Windows 7 32 bit, Python 2.6. > > In [1]: np.__version__ > Out[1]: '1.5.1' > > In [2]: a=np.zeros((0,2)) > > In [3]: np.linalg.qr(a) >  ** On entry to DGEQRF parameter number  4 had an illegal value > In

Re: [Numpy-discussion] Bug: np.linalg.qr crash the python interpreter if getting a shape (0, x) array.

2011-05-08 Thread Skipper Seabold
On Sun, May 8, 2011 at 11:42 AM, Skipper Seabold wrote: > On Fri, May 6, 2011 at 11:30 AM, Till Stensitzki wrote: >> Hi, >> discovered another small bug. Windows 7 32 bit, Python 2.6. >> >> In [1]: np.__version__ >> Out[1]: '1.5.1' >> >>

Re: [Numpy-discussion] finding elements that match any in a set

2011-05-27 Thread Skipper Seabold
On Fri, May 27, 2011 at 12:48 PM, Michael Katz wrote: > I have a numpy array, records, with named fields including a field named > "integer_field". I have an array (or list) of values of interest, and I want > to get the indexes where integer_field has any of those values. > > Because I can do > >

Re: [Numpy-discussion] New functions.

2011-05-31 Thread Skipper Seabold
On Tue, May 31, 2011 at 9:31 PM, Benjamin Root wrote: > > > On Tue, May 31, 2011 at 8:18 PM, Warren Weckesser > wrote: >> >> >> On Tue, May 31, 2011 at 8:08 PM, Charles R Harris >> wrote: >>> >>> Hi All, >>> >>> I've been contemplating new functions that could be added to numpy and >>> thought I

Re: [Numpy-discussion] New functions.

2011-05-31 Thread Skipper Seabold
On Tue, May 31, 2011 at 9:53 PM, Warren Weckesser wrote: > > > On Tue, May 31, 2011 at 8:36 PM, Skipper Seabold > wrote: >> I don't know if it's one pass off the top of my head, but I've used >> percentile for interpercentile ranges. >> >> [d

Re: [Numpy-discussion] New functions.

2011-06-01 Thread Skipper Seabold
On Wed, Jun 1, 2011 at 11:31 AM, Mark Miller wrote: > I'd love to see something like a "count_unique" function included. The > numpy.unique function is handy, but it can be a little awkward to > efficiently go back and get counts of each unique value after the > fact. > Does bincount do what you'

Re: [Numpy-discussion] bincount limitations

2011-06-02 Thread Skipper Seabold
On Wed, Jun 1, 2011 at 10:10 PM, Alan G Isaac wrote: >> On Thu, Jun 2, 2011 at 1:49 AM, Mark Miller   >> wrote: >>> Not quite. Bincount is fine if you have a set of approximately >>> sequential numbers. But if you don't > > > On 6/1/2011 9:35 PM, David Cournapeau wrote: >> Even worse, it fails

Re: [Numpy-discussion] bincount limitations

2011-06-02 Thread Skipper Seabold
On Thu, Jun 2, 2011 at 1:11 PM, Robert Kern wrote: > On Thu, Jun 2, 2011 at 12:08, Skipper Seabold wrote: >> On Wed, Jun 1, 2011 at 10:10 PM, Alan G Isaac wrote: >>>> On Thu, Jun 2, 2011 at 1:49 AM, Mark Miller   >>>> wrote: >>>>> Not quite. Bin

[Numpy-discussion] Numpy/Scipy Testing Guidelines URL?

2011-06-20 Thread Skipper Seabold
Are the testing guidelines included in the HTML docs anywhere? If I recall, they used to be, and I couldn't find them with a brief look/google. I'd like to link to them. Maybe the rendered rst page is considered their new home? https://github.com/numpy/numpy/blob/master/doc/TESTS.rst.txt Skipper

Re: [Numpy-discussion] Numpy/Scipy Testing Guidelines URL?

2011-06-20 Thread Skipper Seabold
On Mon, Jun 20, 2011 at 5:41 PM, Ralf Gommers wrote: > > > On Mon, Jun 20, 2011 at 11:36 PM, Skipper Seabold > wrote: >> >> Are the testing guidelines included in the HTML docs anywhere? If I >> recall, they used to be, and I couldn't find them with a brief

[Numpy-discussion] numpy submodules missing in intersphinx objects.inv?

2011-06-20 Thread Skipper Seabold
I was just trying to link to the numpy.testing module but can't using intersphinx. Does anyone know why certain submodules aren't included in objects.inv? It looks as though it has something to do either with having a reference at the top of the rst file (so you can :ref: link to it) or having that

Re: [Numpy-discussion] Multiply along axis

2011-06-29 Thread Skipper Seabold
On Wed, Jun 29, 2011 at 10:32 AM, Robert Elsner wrote: > > Hello everyone, > > I would like to solve the following problem (preferably without > reshaping / flipping the array a). > > Assume I have a vector v of length x and an n-dimensional array a where > one dimension has length x as well. Now

[Numpy-discussion] Two bugs in recfunctions.join_by with patch

2011-06-29 Thread Skipper Seabold
These two cases failed in recfunctions.join_by 1) the case for having either r1postfix or r2postfix as an empty string was not handled. 2) If there is more than one key and more than variable with a name collision. Patch and tests in a pull request here: https://github.com/numpy/numpy/pull/100

[Numpy-discussion] Moving lib.recfunctions?

2011-07-01 Thread Skipper Seabold
lib.recfunctions has never been fully advertised. The two bugs I just discovered lead me to believe that it's not that well vetted, but it is useful. I can't be the only one using these? What do people think of either deprecating lib.recfunctions or at least importing them into the numpy.rec names

Re: [Numpy-discussion] Moving lib.recfunctions?

2011-07-01 Thread Skipper Seabold
On Fri, Jul 1, 2011 at 2:22 PM, wrote: > On Fri, Jul 1, 2011 at 1:59 PM, Skipper Seabold wrote: >> lib.recfunctions has never been fully advertised. The two bugs I just >> discovered lead me to believe that it's not that well vetted, but it >> is useful. I can

Re: [Numpy-discussion] Moving lib.recfunctions?

2011-07-01 Thread Skipper Seabold
On Fri, Jul 1, 2011 at 3:05 PM, Benjamin Root wrote: > > > On Fri, Jul 1, 2011 at 12:59 PM, Skipper Seabold > wrote: >> >> lib.recfunctions has never been fully advertised. The two bugs I just >> discovered lead me to believe that it's not that well vetted,

Re: [Numpy-discussion] An NA compromise idea -- many-NA

2011-07-01 Thread Skipper Seabold
On Fri, Jul 1, 2011 at 3:46 PM, Dag Sverre Seljebotn wrote: > I propose a simple idea *for the long term* for generalizing Mark's > proposal, that I hope may perhaps put some people behind Mark's concrete > proposal in the short term. > > If key feature missing in Mark's proposal is the ability to

Re: [Numpy-discussion] Moving lib.recfunctions?

2011-07-05 Thread Skipper Seabold
On Fri, Jul 1, 2011 at 2:32 PM, Skipper Seabold wrote: > On Fri, Jul 1, 2011 at 2:22 PM,   wrote: >> On Fri, Jul 1, 2011 at 1:59 PM, Skipper Seabold wrote: >>> lib.recfunctions has never been fully advertised. The two bugs I just >>> discovered lead me to believe tha

Re: [Numpy-discussion] Moving lib.recfunctions?

2011-07-05 Thread Skipper Seabold
On Tue, Jul 5, 2011 at 2:46 PM, Pierre GM wrote: > > On Jul 5, 2011, at 8:33 PM, Skipper Seabold wrote: > >> On Fri, Jul 1, 2011 at 2:32 PM, Skipper Seabold wrote: >>> On Fri, Jul 1, 2011 at 2:22 PM,   wrote: >>>> On Fri, Jul 1, 2011 at 1

Re: [Numpy-discussion] NA/Missing Data Conference Call Summary

2011-07-06 Thread Skipper Seabold
On Wed, Jul 6, 2011 at 7:14 PM, Christopher Jordan-Squire wrote: > On Wed, Jul 6, 2011 at 3:47 PM, wrote: >> On Wed, Jul 6, 2011 at 4:38 PM,   wrote: >> > On Wed, Jul 6, 2011 at 4:22 PM, Christopher Jordan-Squire >> >> Mean value replacement, or more generally single scalar value >> >> replaceme

Re: [Numpy-discussion] code review/build & test for datetime business day API

2011-07-25 Thread Skipper Seabold
2011/7/25 Stéfan van der Walt : > On Mon, Jul 25, 2011 at 2:11 PM, Charles R Harris > wrote: >> It's just asking for import problems and general confusion to shadow a >> Python module, that's why we renamed io to npyio. > > Why?  Users can simply do > > import numpy.io as npyio ? > IIRC this was

Re: [Numpy-discussion] limit to number of fields in recarray

2011-08-02 Thread Skipper Seabold
On Tue, Aug 2, 2011 at 3:19 PM, Craig Yoshioka wrote: > duplicate column in dtype? > "Duplicate field names given."? Can you post code to replicate? > I just consolidated some of the columns and the error went away... none had > duplicate field names... hence the question. > I don't think this

Re: [Numpy-discussion] limit to number of fields in recarray

2011-08-02 Thread Skipper Seabold
On Tue, Aug 2, 2011 at 4:09 PM, Craig Yoshioka wrote: > yup, duplicate field names given.  I didn't commit the non-working version > and I didn't want to mess up my working code so I tried duplicating the dtype > in a new file and couldn't recreate the error.   I suppose the answer to my > ques

Re: [Numpy-discussion] numpydoc - latex longtables error

2011-08-10 Thread Skipper Seabold
On Wed, Aug 10, 2011 at 3:28 PM, Matthew Brett wrote: > Hi, > > I think this one might be for Pauli. > > I've run into an odd problem that seems to be an interaction of > numpydoc and autosummary and large classes. > > In summary, large classes and numpydoc lead to large tables of class > methods,

Re: [Numpy-discussion] summing an array

2011-08-18 Thread Skipper Seabold
On Thu, Aug 18, 2011 at 10:19 AM, Chris Withers wrote: > Hi All, > > Hopefully a simple newbie question, if I have an array such as : > > array([0, 1, 2, 3, 4]) > > ...what's the best way to cummulatively sum it so that I end up with: > > array([0, 1, 3, 6, 10]) > > How would I do this both in-pla

Re: [Numpy-discussion] Identifying Colinear Columns of a Matrix

2011-08-26 Thread Skipper Seabold
On Fri, Aug 26, 2011 at 1:10 PM, Mark Janikas wrote: > Hello All, > > > > I am trying to identify columns of a matrix that are perfectly collinear. > It is not that difficult to identify when two columns are identical are have > zero variance, but I do not know how to ID when the culprit is of a h

[Numpy-discussion] Bug in numpy std, etc. with other data structures?

2011-09-17 Thread Skipper Seabold
Just ran into this. Any objections for having numpy.std and other functions in core/fromnumeric.py call asanyarray before trying to use the array's method? Other data structures like pandas and larry define their own std method, for instance, and this doesn't allow them to pass through. I'm incline

[Numpy-discussion] why doesn't numpy.distutils find ATLAS?

2011-09-21 Thread Skipper Seabold
I have installed numpy with my own ATLAS, but trying to install PyMC, it can't find the ATLAS libs. I also have an older package that formerly installed but no longer does. I'm being a bit lazy, but am I missing something? Briefly checking, it looks like the conventions/assumptions for site.cfg in

Re: [Numpy-discussion] speeding up operations on small vectors

2011-10-11 Thread Skipper Seabold
On Tue, Oct 11, 2011 at 11:57 AM, Christoph Groth wrote: > Pauli Virtanen writes: > >>> Thank you for your suggestion.  It doesn't help me however, because >>> the algorithm I'm _really_ trying to speed up cannot be vectorized >>> with numpy in the way you vectorized my toy example. >>> >>> Any o

Re: [Numpy-discussion] speeding up operations on small vectors

2011-10-11 Thread Skipper Seabold
On Tue, Oct 11, 2011 at 12:41 PM, Christoph Groth wrote: > Skipper Seabold writes: > >> So it's the dot function being called repeatedly on smallish arrays >> that's the bottleneck? I've run into this as well. See this thread >> [1]. >> (...) >

Re: [Numpy-discussion] wanted: decent matplotlib alternative

2011-10-13 Thread Skipper Seabold
On Thu, Oct 13, 2011 at 6:22 PM, Gökhan Sever wrote: > > > On Thu, Oct 13, 2011 at 4:15 PM, Benjamin Root wrote: >> >> Myself and other developers would greatly appreciate help from the >> community to point out which examples are too confusing or out of date. We > > It would be nice to have a so

Re: [Numpy-discussion] simulate AR

2011-10-14 Thread Skipper Seabold
On Fri, Oct 14, 2011 at 2:18 PM, Alan G Isaac wrote: > > On 10/14/2011 1:42 PM, josef.p...@gmail.com wrote: > > If I remember correctly, signal.lfilter doesn't require stationarity, > > but handling of the starting values is a bit difficult. > > > Hmm.  Yes. > AR(1) is trivial, but how do you hand

[Numpy-discussion] adding unsigned int and int

2011-12-06 Thread Skipper Seabold
Hi, Is this intended? [~/] [1]: np.result_type(np.uint, np.int) [1]: dtype('float64') [~/] [2]: np.version.version [2]: '2.0.0.dev-aded70c' Skipper ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/nu

Re: [Numpy-discussion] adding unsigned int and int

2011-12-06 Thread Skipper Seabold
On Tue, Dec 6, 2011 at 7:53 AM, Matthew Brett wrote: > Hi, > > On Tue, Dec 6, 2011 at 4:45 AM, Skipper Seabold wrote: >> Hi, >> >> Is this intended? >> >> [~/] >> [1]: np.result_type(np.uint, np.int) >> [1]: dtype('float64') > &

Re: [Numpy-discussion] Numpy correlate

2013-03-18 Thread Skipper Seabold
On Mon, Mar 18, 2013 at 1:00 PM, Pierre Haessig wrote: > Hi Sudheer, > > Le 14/03/2013 10:18, Sudheer Joseph a écrit : > > Dear Numpy/Scipy experts, > Attached is a script which I > made to test the numpy.correlate ( which is called py plt.xcorr) to s

[Numpy-discussion] Unable to building numpy with openblas using bento or distutils

2013-03-23 Thread Skipper Seabold
Some help on this would be greatly appreciated. It's been recommended to use OpenBlas over ATLAS, so I've been trying to build numpy with openblas and have run into a few problems. 1) Build fails using bento master and waf 1.7.9, see below. 2) Distutils doesn't seem to be able to find lapack as pa

Re: [Numpy-discussion] Unable to building numpy with openblas using bento or distutils

2013-03-23 Thread Skipper Seabold
com > > Sent from my mobile device, please excuse my brevity. > > On 23.03.2013, at 19:19, Skipper Seabold wrote: > > Some help on this would be greatly appreciated. It's been recommended to > use OpenBlas over ATLAS, so I've been trying to build numpy with openblas &g

Re: [Numpy-discussion] Unable to building numpy with openblas using bento or distutils

2013-03-23 Thread Skipper Seabold
On Sat, Mar 23, 2013 at 7:26 PM, Ake Sandgren wrote: > On Sat, 2013-03-23 at 14:19 -0400, Skipper Seabold wrote: >> Some help on this would be greatly appreciated. It's been recommended >> to use OpenBlas over ATLAS, so I've been trying to build numpy with >> o

Re: [Numpy-discussion] Unable to building numpy with openblas using bento or distutils

2013-03-23 Thread Skipper Seabold
On Sat, Mar 23, 2013 at 8:44 PM, Skipper Seabold wrote: > On Sat, Mar 23, 2013 at 7:26 PM, Ake Sandgren > wrote: >> On Sat, 2013-03-23 at 14:19 -0400, Skipper Seabold wrote: >>> Some help on this would be greatly appreciated. It's been recommended >>> to us

Re: [Numpy-discussion] Unable to building numpy with openblas usingbento or distutils

2013-03-25 Thread Skipper Seabold
On Mon, Mar 25, 2013 at 2:40 PM, Dinesh B Vadhia wrote: > ** > Caveat: Not tested but it did look interesting: > http://osdf.github.com/blog/numpyscipy-with-openblas-for-ubuntu-1204-second-try.html > . > Would be interested to know if it worked out as want to try out OpenBlas > in the future. > >

Re: [Numpy-discussion] the mean, var, std of non-arrays

2013-07-18 Thread Skipper Seabold
On Thu, Jul 18, 2013 at 10:49 PM, Yaroslav Halchenko wrote: > Hi everyone, > > Some of my elderly code stopped working upon upgrades of numpy and > upcoming pandas: https://github.com/pydata/pandas/issues/4290 so I have > looked at the code of > > 2481 def mean(a, axis=None, dtype=None, out=No

Re: [Numpy-discussion] Deprecation of financial routines

2013-08-19 Thread Skipper Seabold
On Mon, Aug 19, 2013 at 9:38 AM, Cera, Tim wrote: > On Mon, Aug 19, 2013 at 2:37 AM, Juan Luis Cano > wrote: > > > > As now master is open for 1.9, following the discussion opened here > > > > https://github.com/numpy/numpy/issues/2880 > > > > it was suggested that we deprecate and eventually re

Re: [Numpy-discussion] Behavior of nan{max, min} and nanarg{max, min} for all-nan slices.

2013-10-03 Thread Skipper Seabold
On Thu, Oct 3, 2013 at 9:10 AM, Benjamin Root wrote: > On Thu, Oct 3, 2013 at 6:06 AM, Daniele Nicolodi wrote: >> >> Hello, >> >> sorry, I don't know where exactly jump in in the thread, it is getting >> quite long and articulated... >> >> On 02/10/2013 21:19, Charles R Harris wrote: >> > The mai

Re: [Numpy-discussion] deprecate numpy.matrix

2014-02-10 Thread Skipper Seabold
On Mon, Feb 10, 2014 at 2:49 PM, Matthew Brett wrote: > Hi, > > On Mon, Feb 10, 2014 at 11:44 AM, wrote: > > > > > > On Mon, Feb 10, 2014 at 2:12 PM, eat wrote: > >> > >> > >> > >> > >> On Mon, Feb 10, 2014 at 9:08 PM, alex wrote: > >>> > >>> On Mon, Feb 10, 2014 at 2:03 PM, eat wrote: > >>>

Re: [Numpy-discussion] svd error checking vs. speed

2014-02-15 Thread Skipper Seabold
On Sat, Feb 15, 2014 at 5:08 PM, wrote: > On Sat, Feb 15, 2014 at 4:56 PM, Sebastian Berg > wrote: > > On Sa, 2014-02-15 at 16:37 -0500, alex wrote: > >> Hello list, > >> > >> Here's another idea resurrection from numpy github comments that I've > >> been advised could be posted here for re-disc

Re: [Numpy-discussion] 1.8.1 release

2014-03-06 Thread Skipper Seabold
Hi, Should [1] be considered a release blocker for 1.8.1? Skipper [1] https://github.com/numpy/numpy/issues/4442 ___ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion

Re: [Numpy-discussion] Is there a pure numpy recipe for this?

2014-03-26 Thread Skipper Seabold
On Wed, Mar 26, 2014 at 3:48 PM, Slaunger wrote: > I am working on solving a recent recreational mathematical problem on > Project Euler . I have a solution, which works > fine for small N up to 10^5 but it takes too long to compute for the actual > problem, where N is o

Re: [Numpy-discussion] Is there a pure numpy recipe for this?

2014-03-26 Thread Skipper Seabold
On Wed, Mar 26, 2014 at 4:28 PM, Slaunger wrote: > jseabold wrote >> IIUC, >> >> [~/] >> [1]: np.logical_and([True, False, True], [False, False, True]) >> [1]: array([False, False, True], dtype=bool) >> >> You can avoid looping over k since they're all the same length >> >> [~/] >> [3]: np.logica

Re: [Numpy-discussion] Find indices of largest elements

2010-04-14 Thread Skipper Seabold
On Wed, Apr 14, 2010 at 11:16 AM, Nikolaus Rath wrote: > Hello, > > How do I best find out the indices of the largest x elements in an > array? > > Example: > > a = [ [1,8,2], [2,1,3] ] > magic_function(a, 2) == [ (0,1), (1,2) ] > > Since the largest 2 elements are at positions (0,1) and (1,2). S

Re: [Numpy-discussion] proposing a "beware of [as]matrix()" warning

2010-04-28 Thread Skipper Seabold
On Wed, Apr 28, 2010 at 2:12 PM, Alan G Isaac wrote: > On 4/28/2010 12:05 PM, Travis Oliphant wrote: >> A proposal was made to allow "calling a NumPy array" to infer dot >> product: >> >> a(b) is equivalent to dot(a,b) >> >> a(b)(c) would be equivalent to dot(dot(a,b),c) > > > Here is a related ti

Re: [Numpy-discussion] savetxt not working with python3.1

2010-05-17 Thread Skipper Seabold
On Thu, May 13, 2010 at 8:18 AM, Nadav Horesh wrote: > > in module npyio.py lines 794,796 "file" should be replaced by "_file" > What version of Numpy? Can you file a bug ticket with a failing example? I couldn't replicated with r8417 on Python 3. Skipper __

Re: [Numpy-discussion] Installing numpy on py 3.1.2 , osx

2010-06-07 Thread Skipper Seabold
On Mon, Jun 7, 2010 at 4:02 PM, Vincent Davis wrote: > Going down the prerequisite list I have. > > Python 3.1.2 (r312:79360M, Mar 24 2010, 01:33:18) > [GCC 4.0.1 (Apple Inc. build 5493)] on darwin > Type "help", "copyright", "credits" or "license" for more information. > nose.__version__ > '

Re: [Numpy-discussion] stacking record arrays

2010-06-21 Thread Skipper Seabold
On Mon, Jun 21, 2010 at 2:44 PM, Benjamin Root wrote: > Hello, > > I ran into a somewhat counter-intuitive situation that probably should be > documented somewhere with respect to record (structured?) arrays.  I wanted > to stack multiple arrays together that had the same names for the columns. >

Re: [Numpy-discussion] issue #1431 accessing multiple fields in a recarray

2010-06-27 Thread Skipper Seabold
On Sun, Jun 27, 2010 at 8:26 AM, sam tygier wrote: > Hello > > A while ago there was a discussion about field order when accessing > recarrays. eg: > x = np.array([(1.5,2.5,(1.0,2.0)),(3.,4.,(4.,5.)),(1.,3.,(2.,6.))], >         dtype=[('x','f4'),('y',np.float32),('value','f4',(2,2))]) x

Re: [Numpy-discussion] numpy.all docstring reality check

2010-06-29 Thread Skipper Seabold
On Tue, Jun 29, 2010 at 8:50 PM, David Goldsmith wrote: > Hi, folks.  Under Parameters, the docstring for numpy.core.fromnumeric.all > says: > > "out : ndarray, optionalAlternative output array in which to place the > result. It must have the same shape as the expected output and the type is > pre

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-06 Thread Skipper Seabold
On Tue, Jul 6, 2010 at 10:47 AM, Joshua Holbrook wrote: > I really really really want to work on this. I already forked datarray > on github and did some research on What Other People Have Done ( > http://jesusabdullah.github.com/2010/07/02/datarray.html ). With any > luck I'll contribute somethin

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-06 Thread Skipper Seabold
On Tue, Jul 6, 2010 at 11:55 AM, Keith Goodman wrote: > On Tue, Jul 6, 2010 at 7:47 AM, Joshua Holbrook > wrote: >> I really really really want to work on this. I already forked datarray >> on github and did some research on What Other People Have Done ( >> http://jesusabdullah.github.com/2010/0

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-06 Thread Skipper Seabold
probably not a deal-breaker, but I would imagine that using (a) > separate method(s) for label-based indexing may make allowing > integer-datatyped labels. > > Thoughts? Would you mind bottom-posting/ posting in-line to make the thread easier to follow? > > --Josh > > On T

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-08 Thread Skipper Seabold
On Thu, Jul 8, 2010 at 12:02 PM, Rob Speer wrote: >> While I haven't had a chance to really look in-depth at the changes >> myself (I'm a busy man! So many mailing lists!), I so far like the >> look and sound of them. That's just my opinion, though. > > If people are okay with the attribute magic,

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-08 Thread Skipper Seabold
On Thu, Jul 8, 2010 at 1:35 PM, Rob Speer wrote: >> Forgive me if this is has already been addressed, but my question is >> what happens when we have more than one "label" (not as in a labeled >> axis but an observation label -- but not a tick because they're not >> unique!) per say row axis and h

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-08 Thread Skipper Seabold
On Thu, Jul 8, 2010 at 1:38 PM, Lluís wrote: > Skipper Seabold writes: > >> On Thu, Jul 8, 2010 at 12:02 PM, Rob Speer wrote: > [...] >>> My proposal is that datarray.row should be equivalent to >>> datarray.axes[0], and datarray.column should be equivalent to

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-08 Thread Skipper Seabold
On Thu, Jul 8, 2010 at 2:41 PM, Rob Speer wrote: > On Thu, Jul 8, 2010 at 2:27 PM, Skipper Seabold wrote: >> On Thu, Jul 8, 2010 at 1:35 PM, Rob Speer wrote: >>> Your labels are unique if you look at them the right way. Here's how I >>> would represent that in

Re: [Numpy-discussion] BOF notes: Fernando's proposal: NumPy ndarray with named axes

2010-07-08 Thread Skipper Seabold
On Thu, Jul 8, 2010 at 4:05 PM, Lluís wrote: > Another reason to have multiple variables, is that the insertion of NaNs to > maintain shape homogeneity will make these "synthetic" NaNs undistinguishable > from other NaNs that might be on your original input data, unless you use a > masked array or

Re: [Numpy-discussion] Array concatenation performance

2010-07-15 Thread Skipper Seabold
On Thu, Jul 15, 2010 at 5:54 AM, John Porter wrote: > Has anyone got any advice about array creation. I've been using numpy > for a long time and have just noticed something unexpected about array > concatenation. > > It seems that using numpy.array([a,b,c]) is around 20 times slower > than creati

Re: [Numpy-discussion] Array concatenation performance

2010-07-15 Thread Skipper Seabold
. > > > On Thu, Jul 15, 2010 at 2:41 PM, Skipper Seabold wrote: >> On Thu, Jul 15, 2010 at 5:54 AM, John Porter >> wrote: >>> Has anyone got any advice about array creation. I've been using numpy >>> for a long time and have just noticed somethi

Re: [Numpy-discussion] Array concatenation performance

2010-07-15 Thread Skipper Seabold
y to use hstack and company on lists. In [1]: import numpy as np In [2]: a = np.arange(1000*1000) In [3]: b = a.tolist() In [4]: timeit b0 = np.array((a,a,a)) 1 loops, best of 3: 217 ms per loop In [5]: timeit b1 = np.vstack((b,b,b)) 1 loops, best of 3: 380 ms per loop Skipper > On Thu, J

[Numpy-discussion] build problem on windows: python26.dll not found

2010-07-16 Thread Skipper Seabold
I am way out of my element trying to build numpy on Windows, but we are having a replicable (though somewhat inconsistently on some machines) problem using fmin_cg on Windows that I am trying to debug. I'm on Windows 7 64-bit with Python 2.6.4 First, I try to build numpy. It complains about havi

Re: [Numpy-discussion] Stacking arrays along new dimension in front

2010-07-20 Thread Skipper Seabold
On Tue, Jul 20, 2010 at 5:11 AM, Gael Varoquaux wrote: > Is there in numpy a function that does: > >    np.concatenate([a_[np.newaxis] for a_ in a]) > > ? > > ie: add a dimension in front and stack along this dimension, just like > >    np.array(a) > > would do, but more efficient. > > This is som

Re: [Numpy-discussion] lstsq functionality

2010-07-20 Thread Skipper Seabold
On Tue, Jul 20, 2010 at 1:08 AM, Charles R Harris wrote: > > > On Mon, Jul 19, 2010 at 9:40 PM, Keith Goodman wrote: >> >> On Mon, Jul 19, 2010 at 8:27 PM, Charles R Harris >> wrote: >> > >> > >> > On Mon, Jul 19, 2010 at 9:02 PM, Keith Goodman >> > wrote: >> >> >> >> On Mon, Jul 19, 2010 at 6:

  1   2   3   >