[Numpy-discussion] Any interest in a 'heaviside' ufunc?

I have an implementation of the Heaviside function as numpy ufunc. Is there any interest in adding this to numpy? The function is simply: 0 if x < 0 heaviside(x) = 0.5 if x == 0 1 if x > 0 Warren