Sheldon wrote: > average(compress(ravel(equal(wk,z)),ravel(sattmp)),axis=None) > This is much more compact and elegant. Thanks for pointing this out. > I don't know why average() returned a divide by zero error and to avoid > this I inserted this if statement. Now it works much better !
Probably you had a case where the array length was zero, but it wouldn't happen in the present case unless your input arrays are zero by zero. Carl Banks -- http://mail.python.org/mailman/listinfo/python-list