fixed bug in cluster.split with classify method

[mothur.git] / decalc.cpp

diff --git a/decalc.cpp b/decalc.cpp
index 857bcf2e09c5e20f1f6a03cc95be674898411e6e..74cd8652909d17f0cd0e4e899b1472c6eea2ef93 100644 (file)
--- a/decalc.cpp
+++ b/decalc.cpp
@@ -428,7 +428,7 @@ vector< vector<quanMember> > DeCalculator::getQuantiles(vector<Sequence*> seqs,
                                float de = calcDE(obsi, exp);
                                                                
                                float dist = calcDist(query, subject, front, back); 
-       //o << i << '\t' <<  j << '\t' << dist << '\t' << de << endl;                   
+       //cout << i << '\t' <<  j << '\t' << dist << '\t' << de << endl;                        
                                dist = ceil(dist);
                                
                                quanMember newScore(de, i, j);
@@ -824,6 +824,8 @@ vector<Sequence*> DeCalculator::findClosest(Sequence* querySeq, vector<Sequence*
                        i++;
                }
 
+               if (numWanted > dists.size()) { m->mothurOut("numwanted is larger than the number of template sequences, adjusting numwanted."); m->mothurOutEndLine(); numWanted = dists.size(); }
+
 //cout << numWanted << endl;
                for (int i = 0; i < numWanted; i++) {
 //cout << dists[i].seq->getName() << '\t' << dists[i].dist << endl;