docs: def run_cm(cm, time_scale):
      """Iterate a connectivity matrix the specified number of steps.
  
      Args:
          cm (np.ndarray): A connectivity matrix.
          time_scale (int): The number of steps to run.
  
      Returns:
          np.ndarray: The connectivity matrix at the new timescale.
      """
      cm = np.linalg.matrix_power(cm, time_scale)
      # Round non-unitary values back to 1
      cm[cm > 1] = 1
      return cm