def remove_reaction(self, reaction):
if reaction not in self._reaction_set:<EOL><INDENT>return<EOL><DEDENT>self._reaction_set.remove(reaction)<EOL>self._limits_lower.pop(reaction, None)<EOL>self._limits_upper.pop(reaction, None)<EOL>for compound, value in self._database.get_reaction_values(reaction):<EOL><INDENT>reactions = frozenset(<EOL>self._database.get_compound_reactions(compound))<EOL>if all(other_reaction not in self._reaction_set<EOL>for other_reaction in reactions):<EOL><INDENT>self._compound_set.remove(compound)<EOL><DEDENT><DEDENT>
Remove reaction from model