rsem v1.1.13, speed up EM by only updating model parameters for first 10 iterations...

[rsem.git] / SingleRead.h

diff --git a/SingleRead.h b/SingleRead.h
index 2531176b6cb470941537dce24fb58e26eb0967d8..3c61242202ab3f4bfaad173787767417b04858f3 100644 (file)
--- a/SingleRead.h
+++ b/SingleRead.h
@@ -61,6 +61,8 @@ void SingleRead::calc_lq() {
   int numA = 0, numT = 0, numAO = 0, numTO = 0; // numAO : number of A in overlap seed region
   int threshold_1, threshold_2;
 
+  if (len < OLEN) { low_quality = true; return; }
+
   threshold_1 = int(0.9 * len - 1.5 * sqrt(len * 1.0) + 0.5);
   threshold_2 = (OLEN - 1) / 2 + 1;
   for (int i = 0; i < len; i++) {