[SciPy-Dev] special.factorial accepting floats

The docstring for special.factorial says that the function accepts ints, and int arrays. However, it also seems to accept float arrays (including those that contain np.nan). Is there a rationale for it accepting floats/float arrays? One can't calculate factorials of non-integers.