'Return indices of the minimum values along the given axis.
This is similar to :meth:`~numpy.ndarray.argmin`, but adapted to ensure
that the full precision given by the two doubles ``jd1`` and ``jd2``
is used.  See :func:`~numpy.argmin` for detailed documentation.'

def argmin(self, axis=None, out=None):

    jd = (self.jd1 + self.jd2)
    approx = jd.min(axis, keepdims=True)
    dt = ((self.jd1 - approx) + self.jd2)
    return dt.argmin(axis, out)
