Limits maxIters to numOtus sparcc

mothur · 8497ed352be4d8f309f55c5bff615325f17aff2d · mothur-westcott committed Sep 27, 2016 · Sep 27, 2016 · 8497ed3
1 parent 7bc9807
commit 8497ed352be4d8f309f55c5bff615325f17aff2d
Unified Split

Showing with 3 additions and 0 deletions.

+3 −0 source/commands/sparcccommand.cpp
diff --git a/source/commands/sparcccommand.cpp b/source/commands/sparcccommand.cpp
@@ -199,6 +199,9 @@ int SparccCommand::execute(){
         vector<SharedRAbundVector*> lookup = input.getSharedRAbundVectors();
         string lastLabel = lookup[0]->getLabel();
         
+        int numOTUs = lookup[0]->getNumBins();
+        if (numOTUs >= maxIterations) {  m->mothurOut("[WARNING]: The number of iterations is set higher than the number of OTUs, reducing to " + toString(numOTUs-1) + "\n"); maxIterations = numOTUs-1; }
+        
         //if the users enters label "0.06" and there is no "0.06" in their file use the next lowest label.
         set<string> processedLabels;
         set<string> userLabels = labels;