added abundance and sample parameters to get.coremicrobiome command.

[mothur.git] / classify.cpp

diff --git a/classify.cpp b/classify.cpp
index c92d9fdc6526788df19d9e6f487ffb641bbaca44..2d01183479a504c798a18071df045f723643f850 100644 (file)
--- a/classify.cpp
+++ b/classify.cpp
@@ -61,7 +61,7 @@ void Classify::generateDatabaseAndNames(string tfile, string tempFile, string me
                                        names.push_back(temp.getName());
                                        database->addSequence(temp);    
                                }
-//                             database->generateDB();
+                               database->generateDB();
                        }else if ((method == "kmer") && (!needToGenerate)) {    
                                ifstream kmerFileTest(kmerDBName.c_str());
                                database->readKmerDB(kmerFileTest);     
@@ -150,7 +150,7 @@ void Classify::generateDatabaseAndNames(string tfile, string tempFile, string me
                                        }
                                }
                                
-//                             database->generateDB();
+                               database->generateDB();
                                MPI_File_close(&inMPI);
                                MPI_Barrier(MPI_COMM_WORLD); //make everyone wait - just in case
                #else
@@ -200,7 +200,7 @@ void Classify::generateDatabaseAndNames(string tfile, string tempFile, string me
                                }
                                fastaFile.close();
 
-//                             database->generateDB();
+                               database->generateDB();
                                
                        }else if ((method == "kmer") && (!needToGenerate)) {    
                                ifstream kmerFileTest(kmerDBName.c_str());
@@ -223,9 +223,9 @@ void Classify::generateDatabaseAndNames(string tfile, string tempFile, string me
                        database->setNumSeqs(names.size());
                        
                        //sanity check
-                       //bool okay = phyloTree->ErrorCheck(names);
+                       bool okay = phyloTree->ErrorCheck(names);
                        
-                       //if (!okay) { m->control_pressed = true; }
+                       if (!okay) { m->control_pressed = true; }
                        
                        m->mothurOut("DONE."); m->mothurOutEndLine();
                        m->mothurOut("It took " + toString(time(NULL) - start) + " seconds generate search database. "); m->mothurOutEndLine();
@@ -238,7 +238,7 @@ void Classify::generateDatabaseAndNames(string tfile, string tempFile, string me
        }
 }
 /**************************************************************************************************/
-Classify::Classify() {         m = MothurOut::getInstance();   database = NULL;        }
+Classify::Classify() {         m = MothurOut::getInstance();   database = NULL;        flipped=false; }
 /**************************************************************************************************/
 
 int Classify::readTaxonomy(string file) {
@@ -296,7 +296,8 @@ int Classify::readTaxonomy(string file) {
                        delete buf4;
 
                        istringstream iss (tempBuf,istringstream::in);
-                       iss >> name >> taxInfo;
+                       iss >> name; m->gobble(iss);
+            iss >> taxInfo;
                        taxonomy[name] = taxInfo;
                        phyloTree->addSeqToTree(name, taxInfo);
                }
@@ -309,8 +310,9 @@ int Classify::readTaxonomy(string file) {
        
                //read template seqs and save
                while (!inTax.eof()) {
-                       inTax >> name >> taxInfo;
-                       
+                       inTax >> name; m->gobble(inTax);
+            inTax >> taxInfo;
+               
                        taxonomy[name] = taxInfo;
                        
                        phyloTree->addSeqToTree(name, taxInfo);