Cast integer to double. Solves issue #10.

gavryushkina · gavryushkina · commit d1a84bf65177 · 2018-06-27T14:22:09.000+02:00
diff --git a/src/beast/evolution/substitutionmodel/LewisMK.java b/src/beast/evolution/substitutionmodel/LewisMK.java
@@ -105,7 +105,7 @@ public void getTransitionProbabilities(Node node, double fStartTime, double fEnd
                 matrix[r + i] += e1;
             }
         } else {
-            double fDelta = (nrOfStates / (nrOfStates - 1)) * (fStartTime - fEndTime);
+            double fDelta = (nrOfStates / (nrOfStates - 1.0)) * (fStartTime - fEndTime);
             double fPStay = (1.0 + (nrOfStates - 1) * Math.exp(-fDelta * fRate)) / nrOfStates;
             double fPMove = (1.0 - Math.exp(-fDelta * fRate)) / nrOfStates;
             // fill the matrix with move probabilities

Original file line number	Diff line number	Diff line change
`@@ -105,7 +105,7 @@ public void getTransitionProbabilities(Node node, double fStartTime, double fEnd`
`105`	`105`	`matrix[r + i] += e1;`
`106`	`106`	`}`
`107`	`107`	`} else {`
`108`		`- double fDelta = (nrOfStates / (nrOfStates - 1)) * (fStartTime - fEndTime);`
	`108`	`+ double fDelta = (nrOfStates / (nrOfStates - 1.0)) * (fStartTime - fEndTime);`
`109`	`109`	`double fPStay = (1.0 + (nrOfStates - 1) * Math.exp(-fDelta * fRate)) / nrOfStates;`
`110`	`110`	`double fPMove = (1.0 - Math.exp(-fDelta * fRate)) / nrOfStates;`
`111`	`111`	`// fill the matrix with move probabilities`