code: def _subtract_path_flux(net_flux, path):
      
  
      net_flux = copy.copy(net_flux)
  
      net_flux[path[:-1], path[1:]] -= net_flux[path[:-1], path[1:]].min()
  
      # The above *should* make the bottleneck have zero flux, but
      # numerically that may not be the case, so just set it to zero
      # to be sure.
      bottleneck_ind = net_flux[path[:-1], path[1:]].argmin()
      net_flux[path[bottleneck_ind], path[bottleneck_ind + 1]] = 0.0
  
      return net_flux
docs: Internal function for modifying the net flux matrix by subtracting
      a path's flux from every edge in the path.