Re: [Numpy-discussion] numpy-izing a loop

2009-02-12 Thread Paul Rudin
Paul Rudin writes: > def compute_voxels2(depth_buffers): > dim = depth_buffers[0].shape[0] > znear, zfar, ynear, yfar, xnear, xfar = depth_buffers > z = numpy.arange(dim) > y = numpy.arange(dim)[:, None] > x = numpy.arange(dim)[:, None, None] > > return ((xnear[y,z] < xf

Re: [Numpy-discussion] numpy-izing a loop

2009-02-11 Thread Paul Rudin
Stéfan van der Walt writes: > 2009/2/10 Stéfan van der Walt : >> x = np.arange(dim) >> y = np.arange(dim)[:, None] >> z = np.arange(dim)[:, None, None] > > Do not operate heavy machinery or attempt broadcasting while tired or > under the influence. That order was incorrect: > >> z = np.arange(di

Re: [Numpy-discussion] numpy-izing a loop

2009-02-10 Thread Stéfan van der Walt
2009/2/10 Stéfan van der Walt : > x = np.arange(dim) > y = np.arange(dim)[:, None] > z = np.arange(dim)[:, None, None] Do not operate heavy machinery or attempt broadcasting while tired or under the influence. That order was incorrect: > z = np.arange(dim) > y = np.arange(dim)[:, None] > x = np.

Re: [Numpy-discussion] numpy-izing a loop

2009-02-10 Thread Stéfan van der Walt
2009/2/10 Robert Kern : > x = np.arange(dim)[:,np.newaxis,np.newaxis] > y = np.arange(dim)[np.newaxis,:,np.newaxis] > z = np.arange(dim)[np.newaxis,np.newaxis,:] Yes, sorry, I should have copied from my terminal. I think I had x = np.arange(dim) y = np.arange(dim)[:, None] z = np.arange(dim)[

Re: [Numpy-discussion] numpy-izing a loop

2009-02-10 Thread Robert Kern
On Tue, Feb 10, 2009 at 10:19, Stéfan van der Walt wrote: > Hi Paul > > 2009/2/10 Paul Rudin : >> >> I've just written this snippet of code: >> >> result = numpy.empty((dim, dim, dim), numpy.bool) >> for x in xrange(dim): >>for y in xrange(dim): >>for z in xrange(dim): >>re

Re: [Numpy-discussion] numpy-izing a loop

2009-02-10 Thread Stéfan van der Walt
Hi Paul 2009/2/10 Paul Rudin : > > I've just written this snippet of code: > > result = numpy.empty((dim, dim, dim), numpy.bool) > for x in xrange(dim): >for y in xrange(dim): >for z in xrange(dim): >result[x, y, z] = ((xnear[y, z] < xfar[y, z]) and >

[Numpy-discussion] numpy-izing a loop

2009-02-10 Thread Paul Rudin
I've just written this snippet of code: result = numpy.empty((dim, dim, dim), numpy.bool) for x in xrange(dim): for y in xrange(dim): for z in xrange(dim): result[x, y, z] = ((xnear[y, z] < xfar[y, z]) and (ynear[x, z] < yfar[x, z]) and