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

[rsem.git] / SingleReadQ.h

diff --git a/SingleReadQ.h b/SingleReadQ.h
index a0099d9c41054b68da45880f94a8857e8f232f47..976637fc269d79fa484cb0033421e7f7a0769c49 100644 (file)
--- a/SingleReadQ.h
+++ b/SingleReadQ.h
@@ -67,6 +67,8 @@ void SingleReadQ::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++) {