code: def _return_path(acc_cost_mat):
      n_timestamps_1, n_timestamps_2 = acc_cost_mat.shape
      path = [(n_timestamps_1 - 1, n_timestamps_2 - 1)]
      while path[-1] != (0, 0):
          i, j = path[-1]
          if i == 0:
              path.append((0, j - 1))
          elif j == 0:
              path.append((i - 1, 0))
          else:
              arr = np.array(
                  [
                      acc_cost_mat[i - 1][j - 1],
                      acc_cost_mat[i - 1][j],
                      acc_cost_mat[i][j - 1],
                  ]
              )
              argmin = np.argmin(arr)
              if argmin == 0:
                  path.append((i - 1, j - 1))
              elif argmin == 1:
                  path.append((i - 1, j))
              else:
                  path.append((i, j - 1))
      return np.transpose(np.array(path)[::-1])
  
imports: import numpy
  from math import ceil, log2, sqrt
  from numba import njit, prange
  from sklearn.utils import check_array
label_desc: Use image node_id 13 for calling a global function with example usage: _return_path(acc_cost_mat) and returns: np