I know there are probably alternatives for this with the standard library, but I think that would kill the speed I get with numarray:
Say I have two 2-dimensional numarrays (x_mat and y_mat, say), and a function f(x,y) that I would like to evaluate at every index. Basically I want to be able to say f(x_mat,y_mat) and have it return a numarray with the same shape and element wise evaluation of f. I know, I want a ufunc, but they look scary to write my own. Are there any functions that do this? Or is writing ufuncs easier than it seems?
It depends. If you write f(x,y) such that it operates on arrays as a whole (meaning you can't use if statements, etc.), then it may Just Work.
If you need to do something where you can't express it like that, then take a look at Scipy's[1] vectorize[2] function.
[1] http://www.scipy.org [2] http://oliphant.ee.byu.edu:81/scipy_base/vectorize/
-- Robert Kern [EMAIL PROTECTED]
"In the fields of hell where the grass grows high Are the graves of dreams allowed to die." -- Richard Harter -- http://mail.python.org/mailman/listinfo/python-list
