2008/12/4 Stéfan van der Walt <[EMAIL PROTECTED]>: > Hi Olivier > > 2008/12/4 Olivier Grisel <[EMAIL PROTECTED]>: >> To avoid the python level loop I then tried to use broadcasting as follows: >> >>>>> c = sum((a[:,newaxis,:] - b) ** 2, axis=2) >> >> But this build a useless and huge (d1, d2, d3) temporary array that >> does not fit in memory for large values of d1, d2 and d3... > > Does numpy.lib.broadcast_arrays do what you need?
That looks exactly what I am looking for. Apparently this is new in 1.2 since I cannot find it in the 1.1 version of my system. Thanks, -- Olivier _______________________________________________ Numpy-discussion mailing list [email protected] http://projects.scipy.org/mailman/listinfo/numpy-discussion
