[Numpy-discussion] kind of a matrix multiplication

2011-10-10 Thread Martin Raspaud
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Hi all, I have a stack of vectors: v1 = np.arange(3) v2 = np.arange(3) + 3 stack = np.vstack(v1, v2) (now stack is : array([[0, 1, 2], [3, 4, 5]])) and a 3d matrix: mat = np.dstack((np.eye(3), np.eye(3) * 2)) (mat is now array([[[ 1., 2.],

Re: [Numpy-discussion] np.where and broadcasting

2011-10-10 Thread Pauli Virtanen
10.10.2011 23:02, Jesper Larsen kirjoitti: > Hi numpy-users > > I have a 2d array of shape (ny, nx). I want to broadcast (copy) this > array to a target array of shape (nt, nz, ny, nz) or (nt, ny, nx) so > that the 2d array is repeated for each t and z index (corresponding to > nt and nz). I am not

[Numpy-discussion] np.where and broadcasting

2011-10-10 Thread Jesper Larsen
Hi numpy-users I have a 2d array of shape (ny, nx). I want to broadcast (copy) this array to a target array of shape (nt, nz, ny, nz) or (nt, ny, nx) so that the 2d array is repeated for each t and z index (corresponding to nt and nz). I am not sure how to do this (generic solution, for different

[Numpy-discussion] failure to build

2011-10-10 Thread Amos Anderson
hello -- i'm trying to build numpy-1.6.1. so far, it works on several machines that i've tried, but now it's failing. I put some details on the machine below. I tried to change the compiler, but setup.py objected to the things I tried (--compiler=g++4, --compiler=icc). thanks for any advice! Amo

Re: [Numpy-discussion] Lookup array

2011-10-10 Thread Olivier Delalleau
2011/10/10 Andrey N. Sobolev > В Mon, 10 Oct 2011 10:03:48 +0100 > Bob Dowling пишет: > > > > > On 10/10/11 09:53, Andrey N. Sobolev wrote: > > > > > I have 2 arrays - A with the dimensions of 1000x4 and B with the > > > dimensions of 5000x4. B doesn't (hopefully) contain any rows that > > > are

[Numpy-discussion] scikits.image 0.3 release

2011-10-10 Thread Stéfan van der Walt
Announcement: scikits.image 0.3 === After a brief (!) absence, we're back with a new and shiny version of scikits.image, the image processing toolbox for SciPy. This release runs under all major operating systems where Python (>=2.6 or 3.x), NumPy and SciPy can be inst

Re: [Numpy-discussion] Lookup array

2011-10-10 Thread Andrey N. Sobolev
В Mon, 10 Oct 2011 10:03:48 +0100 Bob Dowling пишет: > > On 10/10/11 09:53, Andrey N. Sobolev wrote: > > > I have 2 arrays - A with the dimensions of 1000x4 and B with the > > dimensions of 5000x4. B doesn't (hopefully) contain any rows that > > are not in A. I need to create a lookup array C,

Re: [Numpy-discussion] Lookup array

2011-10-10 Thread Bob Dowling
On 10/10/11 09:53, Andrey N. Sobolev wrote: > I have 2 arrays - A with the dimensions of 1000x4 and B with the > dimensions of 5000x4. B doesn't (hopefully) contain any rows that are > not in A. I need to create a lookup array C, the i-th value of which > will be the index of B[i] in A. In the (v

[Numpy-discussion] Lookup array

2011-10-10 Thread Andrey N. Sobolev
Hi all, I have 2 arrays - A with the dimensions of 1000x4 and B with the dimensions of 5000x4. B doesn't (hopefully) contain any rows that are not in A. I need to create a lookup array C, the i-th value of which will be the index of B[i] in A. In the (very rare) case when B[i] is not in A C[i] sho

Re: [Numpy-discussion] numpy.lib.npyio.load

2011-10-10 Thread David Cournapeau
Hi Jean-Louis, On Sun, Oct 9, 2011 at 2:37 PM, Jean-Louis Durrieu wrote: > Hi everyone, > > I was just wondering something: lately, I had to use the load function, to > load arrays stored in npz files. > > During one session, I need to read quite a few times several files (or even > the same fi