Fixed logP calcs in approximate steppers.

tgvaughan · Jun 18, 2014 · 237949d · 237949d
1 parent 0117b26
commit 237949d
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/master/steppers/SALStepper.java b/src/master/steppers/SALStepper.java
@@ -88,7 +88,7 @@ public void leap(Reaction reaction, PopulationState state,
 
         if (calcLogP) {
             if (rho>0)
-                stepLogP += -rho + q*Math.log(rho) - Gamma.logGamma(q+1);
+                stepLogP += -rho + q*Math.log(rho/thisdt) - Gamma.logGamma(q+1);
         }
 
         // Implement reactions:

diff --git a/src/master/steppers/TauLeapingStepper.java b/src/master/steppers/TauLeapingStepper.java
@@ -75,7 +75,7 @@ public void leap(Reaction reaction, PopulationState state, Model model, boolean
 
         if (calcLogP) {
             if (rho>0)
-                stepLogP += -rho + q*Math.log(rho) - Gamma.logGamma(q+1);
+                stepLogP += -rho + q*Math.log(rho/thisdt) - Gamma.logGamma(q+1);
         }
 
         // Implement reactions: