KineticLaw.py 文件源码

def getRawVelocities(self, subs={}, include_fast_reaction=True, include_slow_reaction=True):

        definition = self.__definition.getDeveloppedInternalMathFormula()
        if not KineticLawIdentifier.isReversible(self, definition):
            res = zeros(1)
            if (self.reaction.fast and include_fast_reaction) or (not self.reaction.fast and include_slow_reaction):
                res[0] += definition.subs(subs)
            return res

        else:
            res_front = zeros(1)
            res_backward = zeros(1)
            if (self.reaction.fast and include_fast_reaction) or (not self.reaction.fast and include_slow_reaction):
                (definition_front, definition_back) = KineticLawIdentifier.getReversibleRates(self, definition)
                res_front[0] += definition_front.subs(subs)
                res_backward[0] += definition_back.subs(subs)

            return res_front.col_join(res_backward)