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.