Array support for np.datetime64 and np.timedelta64

This includes support for the basic arithmetic ufuncs.
It also augments ufunc testing with more values (including interesting ones such as 0 or NaN), and fixes a couple buglets found when doing so.
For the record, I'm not sure about the handling of -0.0 here...