Skip to content
代码片段 群组 项目
提交 167b4945 编辑于 作者: k3yavi's avatar k3yavi
浏览文件

fixing mtx bug

上级 179191a0
分支
标签 v0.14.2
无相关合并请求
......@@ -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);
......
0% .
You are about to add 0 people to the discussion. Proceed with caution.
先完成此消息的编辑!
想要评论请 注册