Fix allostery free enzyme ratio bug

It sounds like it should be a per-enzyme vector, but it's a per-edge one. Oops!
biosustain · Aug 14, 2024 · 9ed5064 · 9ed5064
1 parent 094d460
commit 9ed5064
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/maud/stan/functions.stan b/maud/stan/functions.stan
@@ -242,7 +242,7 @@ functions {
       }
       out[f] = inv(1
                    + tc_edge
-                     * (free_enzyme_ratio[edge_to_enzyme[f]] * Q_num / Q_denom) ^ subunits[edge_to_enzyme[f]]);
+                     * (free_enzyme_ratio[f] * Q_num / Q_denom) ^ subunits[edge_to_enzyme[f]]);
     }
     return out;
   }
@@ -680,7 +680,7 @@ functions {
       out[f] = 1
                / (1
                   + tc_edge
-                    * (free_enzyme_ratio[edge_to_enzyme[f]] * Q_num / Q_denom) ^ subunits[edge_to_enzyme[f]]);
+                    * (free_enzyme_ratio[f] * Q_num / Q_denom) ^ subunits[edge_to_enzyme[f]]);
     }
     return out;
   }