code: def fastcc_consistent_subset(model, epsilon, solver):
      
      reaction_set = set(model.reactions)
      return reaction_set.difference(fastcc(model, epsilon, solver))
docs: Return consistent subset of model.
  
      The largest consistent subset is returned as
      a set of reaction names.
  
      Args:
          model: :class:`MetabolicModel` to solve.
          epsilon: Flux threshold value.
          solver: LP solver instance to use.
  
      Returns:
          Set of reaction IDs in the consistent reaction subset.