Re: [Numpy-discussion] OpenBLAS and dotblas

2014-08-07 Thread Charles R Harris
On Thu, Aug 7, 2014 at 8:32 PM, Charles R Harris wrote: > Hi All, > > It looks like numpy dot only uses BLAS if ATLAS is present, see > numpy/core/setup.py. Has anyone done the mods needed to use OpenBLAS? What > is the current status of using OpenBLAS with numpy? > NVM, the comments and "NO_ATL

[Numpy-discussion] OpenBLAS and dotblas

2014-08-07 Thread Charles R Harris
Hi All, It looks like numpy dot only uses BLAS if ATLAS is present, see numpy/core/setup.py. Has anyone done the mods needed to use OpenBLAS? What is the current status of using OpenBLAS with numpy? Also, I'm thinking of moving linalg/lapack_lite/blas_lite.c down into the core. Thoughts? Chuck _

Re: [Numpy-discussion] ANN: NumPy 1.8.2 release candidate

2014-08-07 Thread Frédéric Bastien
All Theano tests work. thanks! Fred On Tue, Aug 5, 2014 at 8:46 PM, Matthew Brett wrote: > Hi, > > On Tue, Aug 5, 2014 at 2:27 PM, Matthew Brett > wrote: > > Hi, > > > > On Tue, Aug 5, 2014 at 1:57 PM, Julian Taylor > > wrote: > >> On 05.08.2014 22:32, Christoph Gohlke wrote: > >>> On 8/5/20

Re: [Numpy-discussion] Preliminary thoughts on implementing __matmul__

2014-08-07 Thread cjw
On 07/08/2014 5:40 AM, Sebastian Berg wrote: > On Mi, 2014-08-06 at 14:05 -0700, Chris Barker wrote: >> On Wed, Aug 6, 2014 at 8:32 AM, Charles R Harris >> wrote: >> Should also mention that we don't have the ability to operate >> on stacked vectors because they can't be identif

Re: [Numpy-discussion] Inverted indices

2014-08-07 Thread Nicolas P. Rougier
Oh thanks, I would never have imagined a one-line solution... Here are the benchmarks: In [2]: %timeit stefan(S) 10 loops, best of 3: 10.8 µs per loop In [3]: %timeit gregor(S) 1 loops, best of 3: 48.1 µs per loop In [4]: %timeit gregor2(S) 10 loops, best of 3: 3.23 µs per loop

Re: [Numpy-discussion] Inverted indices

2014-08-07 Thread Gregor Thalhammer
Am 07.08.2014 um 13:59 schrieb Gregor Thalhammer : > > Am 07.08.2014 um 13:16 schrieb Nicolas P. Rougier : > >> >> Hi, >> >> I've a small problem for which I cannot find a solution and I'm quite sure >> there is an obvious one: >> >> I've an array Z (any dtype) with some data. >> I've a (so

Re: [Numpy-discussion] Inverted indices

2014-08-07 Thread Gregor Thalhammer
Am 07.08.2014 um 13:16 schrieb Nicolas P. Rougier : > > Hi, > > I've a small problem for which I cannot find a solution and I'm quite sure > there is an obvious one: > > I've an array Z (any dtype) with some data. > I've a (sorted) array I (of integer, same size as Z) that tells me the index

Re: [Numpy-discussion] Inverted indices

2014-08-07 Thread Nicolas P. Rougier
Nice ! Thanks Stéfan. I will add it to the numpy 100 problems. Nicolas On 07 Aug 2014, at 13:31, Stéfan van der Walt wrote: > Hi Nicolas > > On Thu, Aug 7, 2014 at 1:16 PM, Nicolas P. Rougier > wrote: >> Here is a small example: >> >> Z = [(0,0), (1,1), (2,2), (3,3), (4,4)) >> I = [0,

Re: [Numpy-discussion] Inverted indices

2014-08-07 Thread Stéfan van der Walt
Hi Nicolas On Thu, Aug 7, 2014 at 1:16 PM, Nicolas P. Rougier wrote: > Here is a small example: > > Z = [(0,0), (1,1), (2,2), (3,3), (4,4)) > I = [0, 20, 23, 24, 37] > > S = [ 20,20,0,24] > -> Result should be [(1,1), (1,1), (0,0),(3,3)] > > S = [15,15] > -> Wrong (15 not in I) but ideally, I wo

[Numpy-discussion] Inverted indices

2014-08-07 Thread Nicolas P. Rougier
Hi, I've a small problem for which I cannot find a solution and I'm quite sure there is an obvious one: I've an array Z (any dtype) with some data. I've a (sorted) array I (of integer, same size as Z) that tells me the index of Z[i] (if necessary, the index can be stored in Z). Now, I have

Re: [Numpy-discussion] Preliminary thoughts on implementing __matmul__

2014-08-07 Thread Sebastian Berg
On Mi, 2014-08-06 at 14:05 -0700, Chris Barker wrote: > On Wed, Aug 6, 2014 at 8:32 AM, Charles R Harris > wrote: > Should also mention that we don't have the ability to operate > on stacked vectors because they can't be identified by > dimension info. One workaround is to

Re: [Numpy-discussion] Preliminary thoughts on implementing __matmul__

2014-08-07 Thread Eelco Hoogendoorn
I don't expect stacked matrices/vectors to be used often, although there are some areas that might make heavy use of them, so I think we could live with the simple implementation, it's just a bit of a wart when there is broadcasting of arrays. Just to be clear, the '@' broadcasting differs from the