fixing mtx bug

COMBINE-lab · Oct 4, 2019 · 167b494 · 167b494
1 parent 179191a
commit 167b494
Showing 1 changed file with 7 additions and 7 deletions.
diff --git a/src/CollapsedCellOptimizer.cpp b/src/CollapsedCellOptimizer.cpp
@@ -437,12 +437,6 @@ void optimizeCell(std::vector<std::string>& trueBarcodes,
                        << "\t" << totalExpGenes
                        << "\t" << numGenesOverMean;
 
-        if (dumpUmiGraph) {
-          featuresStream << arboString.rdbuf();
-        } else {
-          featuresStream << "\t" << averageNumMolPerArbo;
-        }
-
         if (mRnaGenes.size() > 1) {
           featureCode += 1;
           featuresStream << "\t" << mitoCount / totalUmiCount;
@@ -453,6 +447,12 @@ void optimizeCell(std::vector<std::string>& trueBarcodes,
           featuresStream << "\t" << riboCount / totalUmiCount;
         }
 
+        if (dumpUmiGraph) {
+          featuresStream << arboString.rdbuf();
+        } else {
+          featuresStream << "\t" << averageNumMolPerArbo;
+        }
+
         features = featuresStream.str();
       } // end making features
 
@@ -788,7 +788,7 @@ bool CollapsedCellOptimizer::optimize(EqMapT& fullEqMap,
       };
 
       uint32_t zerod_cells {0};
-      size_t numFlags = (numGenes/8)+1;
+      size_t numFlags = std::ceil(numGenes/8);
       std::vector<uint8_t> alphasFlag (numFlags, 0);
       size_t flagSize = sizeof(decltype(alphasFlag)::value_type);