code: def rescale_around1(self, times):
          
          if self._unit == self._UNIT_STEP:
              return times, 'step' # nothing to do
  
          m = np.mean(times)
          mult = 1.0
          cur_unit = self._unit
  
          # numbers are too small. Making them larger and reducing the unit:
          if (m < 0.001):
              while mult*m < 0.001 and cur_unit >= 0:
                  mult *= 1000
                  cur_unit -= 1
              return mult*times, self._unit_names[cur_unit]
  
          # numbers are too large. Making them smaller and increasing the unit:
          if (m > 1000):
              while mult*m > 1000 and cur_unit <= 5:
                  mult /= 1000
                  cur_unit += 1
              return mult*times, self._unit_names[cur_unit]
  
          # nothing to do
          return times, self._unit
docs: Suggests a rescaling factor and new physical time unit to balance the given time multiples around 1.
  
          Parameters
          ----------
          times : float array
              array of times in multiple of the present elementary unit