Re: numpy : efficient sum computations

2007-10-16 Thread Robert Kern
TG wrote: > Okay, another one which I don't have answer for. > > it is the reverse case, sort of : > >> phi.shape (x,y) >> d.shape (a,b) > > I want to return m : >> m.shape = (x,y,a,b) > with > m[x,y] = d * phi[x,y] > > currently, my code is : >> m = empty(phi.shape + d.shape) >> m[:,:] = d >

Re: numpy : efficient sum computations

2007-10-16 Thread Robert Kern
TG wrote: > Hi there. > > I want to do some intensive computations with numpy, and I'm > struggling a bit to find my wayy. The best place to ask numpy questions is on the numpy mailing list. http://www.scipy.org/Mailing_Lists > Here is the problem : > > m and d are two matrices : >

Re: numpy : efficient sum computations

2007-10-16 Thread TG
Okay, another one which I don't have answer for. it is the reverse case, sort of : > phi.shape (x,y) > d.shape (a,b) I want to return m : > m.shape = (x,y,a,b) with m[x,y] = d * phi[x,y] currently, my code is : > m = empty(phi.shape + d.shape) > m[:,:] = d this repeats the matrix d x*y times,