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