[SciPy-User] Numpy/Scipy: Avoiding nested loops to operate on matrix-valued images

David Warde-Farley wardefar@iro.umontreal...
Wed Mar 21 22:34:48 CDT 2012


On 2012-03-20, at 11:33 AM, Nathaniel Smith wrote:

>  (Are there any generalized
> ufuncs actually defined in the standard packages? For instance, is
> np.dot implemented as a generalized ufunc? Should it be?)

Ideally, so long as it still made use of BLAS for the actual matrix products.

I tried my hand at implementing a gufunc for log(sum(exp(...))), with the sum being the "generalized" part. Did not have much luck...

David


More information about the SciPy-User mailing list