X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=simulation.cpp;h=1288c65178afa360f89f4cff3ee5a2f78b7de5cf;hb=1d786dca8b4c9c1e2176a28a669d2d16cd93e395;hp=838a88fcccae5bf0887c504ac5fef770da8458a7;hpb=9c2e46183a19d661f0a618a8eabe8ce1f6a8e2d6;p=rsem.git

diff --git a/simulation.cpp b/simulation.cpp
index 838a88f..1288c65 100644
--- a/simulation.cpp
+++ b/simulation.cpp
@@ -133,7 +133,7 @@ void simulate(char* modelF, char* resultsF) {
 	  size_t pos2 = line.find_first_of('\t', pos);
 	  if (pos2 == string::npos) pos2 = line.length();
 	  tpm = atof(line.substr(pos, pos2 - pos).c_str());
-	  theta[i] = tpm * eel[i];
+	  theta[i] = tpm * eel[i]; // during simulation, there is no check for effL < 0. The reason is for that case, eel[i] here = 0 and therefore no chance to sample from it
 	  denom += theta[i];
 	}
 	assert(denom > EPSILON);