docs: def _lemke_howson_capping(payoff_matrices, tableaux, bases, init_pivot,
                            max_iter, capping):
      """
      Execute the Lemke-Howson algorithm with the heuristics proposed by
      Codenotti et al.
  
      Parameters
      ----------
      payoff_matrices : tuple(ndarray(ndim=2))
          Tuple of two arrays representing payoff matrices, of shape
          (m, n) and (n, m), respectively.
  
      tableaux : tuple(ndarray(float, ndim=2))
          Tuple of two arrays to be used to store the tableaux, of shape
          (n, m+n+1) and (m, m+n+1), respectively. Modified in place.
  
      bases : tuple(ndarray(int, ndim=1))
          Tuple of two arrays to be used to store the bases, of shape (n,)
          and (m,), respectively. Modified in place.
  
      init_pivot : scalar(int)
          Integer k such that 0 <= k < m+n.
  
      max_iter : scalar(int)
          Maximum number of pivoting steps.
  
      capping : scalar(int)
          Value for capping. If set equal to `max_iter`, then the routine
          is equivalent to the standard Lemke-Howson algorithm.
  
      """
      m, n = tableaux[1].shape[0], tableaux[0].shape[0]
      init_pivot_curr = init_pivot
      max_iter_curr = max_iter
      total_num_iter = 0
  
      for k in range(m+n-1):
          capping_curr = min(max_iter_curr, capping)
  
          _initialize_tableaux(payoff_matrices, tableaux, bases)
          converged, num_iter = \
              _lemke_howson_tbl(tableaux, bases, init_pivot_curr, capping_curr)
  
          total_num_iter += num_iter
  
          if converged or total_num_iter >= max_iter:
              return converged, total_num_iter, init_pivot_curr
  
          init_pivot_curr += 1
          if init_pivot_curr >= m + n:
              init_pivot_curr -= m + n
          max_iter_curr -= num_iter
  
      _initialize_tableaux(payoff_matrices, tableaux, bases)
      converged, num_iter = \
          _lemke_howson_tbl(tableaux, bases, init_pivot_curr, max_iter_curr)
      total_num_iter += num_iter
  
      return converged, total_num_iter, init_pivot_curr