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.

def fastcc_consistent_subset(model, epsilon, solver):
    """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.
    """
    reaction_set = set(model.reactions)
    return reaction_set.difference(fastcc(model, epsilon, solver))