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);
}
}
-// database->generateDB();
+ database->generateDB();
MPI_File_close(&inMPI);
MPI_Barrier(MPI_COMM_WORLD); //make everyone wait - just in case
#else
}
fastaFile.close();
-// database->generateDB();
+ database->generateDB();
}else if ((method == "kmer") && (!needToGenerate)) {
ifstream kmerFileTest(kmerDBName.c_str());
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();
}
}
/**************************************************************************************************/
-Classify::Classify() { m = MothurOut::getInstance(); database = NULL; }
+Classify::Classify() { m = MothurOut::getInstance(); database = NULL; flipped=false; }
/**************************************************************************************************/
int Classify::readTaxonomy(string file) {
m->mothurOutEndLine();
m->mothurOut("Reading in the " + file + " taxonomy...\t"); cout.flush();
-
+ if (m->debug) { m->mothurOut("[DEBUG]: Taxonomies read in...\n"); }
+
#ifdef USE_MPI
int pid, num, processors;
vector<unsigned long long> positions;
delete buf4;
istringstream iss (tempBuf,istringstream::in);
- iss >> name >> taxInfo;
+ iss >> name; m->gobble(iss);
+ iss >> taxInfo;
+ if (m->debug) { m->mothurOut("[DEBUG]: name = " + name + " tax = " + taxInfo + "\n"); }
taxonomy[name] = taxInfo;
phyloTree->addSeqToTree(name, taxInfo);
}
//read template seqs and save
while (!inTax.eof()) {
- inTax >> name >> taxInfo;
-
+ inTax >> name; m->gobble(inTax);
+ inTax >> taxInfo;
+
+ if (m->debug) { m->mothurOut("[DEBUG]: name = '" + name + "' tax = '" + taxInfo + "'\n"); }
+
taxonomy[name] = taxInfo;
phyloTree->addSeqToTree(name, taxInfo);
}
inTax.close();
#endif
+
+
phyloTree->assignHeirarchyIDs(0);