- Added poly(A) tails to 'reference_name.transcripts.fa' so that the RSEM generated...

[rsem.git] / Gibbs.cpp

diff --git a/Gibbs.cpp b/Gibbs.cpp
index bf8563e16ea673b467e22b7d51343b7bd63b27e3..b0287a334cef604ea4385e250c7e22a8511c0424 100644 (file)
--- a/Gibbs.cpp
+++ b/Gibbs.cpp
@@ -254,7 +254,7 @@ void calcExpressionValues(const vector<double>& theta, const vector<double>& eel
            frac[i] = theta[i];
            denom += frac[i];
          }
-       general_assert(denom > 0, "No alignable reads?!");
+       general_assert(denom >= EPSILON, "No alignable reads?!");
        for (int i = 1; i <= M; i++) frac[i] /= denom;
   
        //calculate FPKM