//cout << numKmers << '\t' << genusNodes.size() << endl;
for (int j = 0; j < wordGenusProb.size(); j++) { wordGenusProb[j].resize(genusNodes.size()); }
//cout << numKmers << '\t' << genusNodes.size() << endl;
- ofstream out;
- ofstream out2;
+ //ofstream out;
+ //ofstream out2;
#ifdef USE_MPI
int pid;
#endif
- m->openOutputFile(probFileName, out);
+ //m->openOutputFile(probFileName, out);
//output mothur version
- out << "#" << m->getVersion() << endl;
+ //out << "#" << m->getVersion() << endl;
- out << numKmers << endl;
+ //out << numKmers << endl;
- m->openOutputFile(probFileName2, out2);
+ //m->openOutputFile(probFileName2, out2);
//output mothur version
- out2 << "#" << m->getVersion() << endl;
+ //out2 << "#" << m->getVersion() << endl;
#ifdef USE_MPI
}
if (pid == 0) {
#endif
- out << i << '\t';
+ //out << i << '\t';
#ifdef USE_MPI
}
if (pid == 0) {
#endif
- out << k << '\t' << wordGenusProb[i][k] << '\t';
+ //out << k << '\t' << wordGenusProb[i][k] << '\t';
#ifdef USE_MPI
}
if (pid == 0) {
#endif
- out << endl;
- out2 << probabilityInTemplate << '\t' << numNotZero << endl;
+ //out << endl;
+ //out2 << probabilityInTemplate << '\t' << numNotZero << endl;
#ifdef USE_MPI
}
if (pid == 0) {
#endif
- out.close();
- out2.close();
+ //out.close();
+ //out2.close();
#ifdef USE_MPI
}
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();
ifeq ($(strip $(64BIT_VERSION)),yes)
#if you are using centos uncomment the following lines
- #CXX = g++44
+ CXX = g++44
#if you are a mac user use the following line
- TARGET_ARCH += -arch x86_64
+# TARGET_ARCH += -arch x86_64
#if you using cygwin to build Windows the following line
#CXX = x86_64-w64-mingw32-g++
#TARGET_ARCH += -m64 -static
#if you are a linux user use the following line
- #CXXFLAGS += -mtune=native -march=native -m64
+ CXXFLAGS += -mtune=native -march=native -m64
CXXFLAGS += -DBIT_VERSION
endif
void ShhherCommand::initPyroCluster(){
try{
-
if (numOTUs < processors) { processors = 1; }
-
+
dist.assign(numSeqs * numOTUs, 0);
change.assign(numOTUs, 1);
centroids.assign(numOTUs, -1);