Fix for all homo-dimer reactions

ahellander · Aug 29, 2014 · 2dddcad · 2dddcad
1 parent d4d89df
commit 2dddcad
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/pyurdme/model.py b/pyurdme/model.py
@@ -326,7 +326,7 @@ def create_mass_action(self):
             raise ReactionError("Reaction: " +self.name + "A mass-action reaction cannot involve more than two species.")
 
         # Case EmptySet -> Y
-        propensity_function = 'return ' + self.marate.name;
+        propensity_function = self.marate.name;
 
         # There are only three ways to get 'total_stoch==2':
         for r in self.reactants:
@@ -345,7 +345,7 @@ def create_mass_action(self):
             propensity_function += "*vol"
 
 
-        self.propensity_function = propensity_function + ';'
+        self.propensity_function = "return " + propensity_function + ';'
 
     def set_type(self,type):
         if type not in {'mass-action','customized'}: