Peak to peak (maximum - minimum) along a given axis.

This is similar to ptp(), but adapted to ensure that the full precision given by the two doubles jd1 and jd2 is used.

Note that the out argument is present only for compatibility with ~numpy.ptp; since Time instances are immutable, it is not possible to have an actual out to store the result in.