use if/elsif to replace switch in rsem-calculate-expression

[rsem.git] / simulation.cpp

diff --git a/simulation.cpp b/simulation.cpp
index 838a88fcccae5bf0887c504ac5fef770da8458a7..1288c65178afa360f89f4cff3ee5a2f78b7de5cf 100644 (file)
--- 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);