]> git.donarmstrong.com Git - mothur.git/blobdiff - mothurout.cpp
modified trim.seqs to speed up the append with multiple processors
[mothur.git] / mothurout.cpp
index bd1a098c8d5662044552071c7c44a273214018ce..eb452b19a681e5ebe83d71ca7975a4140f703a55 100644 (file)
@@ -9,6 +9,7 @@
 
 #include "mothurout.h"
 
+
 /******************************************************/
 MothurOut* MothurOut::getInstance() {
        if( _uniqueInstance == 0) {
@@ -359,18 +360,21 @@ string MothurOut::getline(ifstream& fileHandle) {
 }
 /***********************************************************************/
 
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
 #ifdef USE_COMPRESSION
 inline bool endsWith(string s, const char * suffix){
   size_t suffixLength = strlen(suffix);
   return s.size() >= suffixLength && s.substr(s.size() - suffixLength, suffixLength).compare(suffix) == 0;
 }
 #endif
+#endif
 
 string MothurOut::getRootName(string longName){
        try {
        
                string rootName = longName;
 
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
 #ifdef USE_COMPRESSION
     if (endsWith(rootName, ".gz") || endsWith(rootName, ".bz2")) {
       int pos = rootName.find_last_of('.');
@@ -378,7 +382,7 @@ string MothurOut::getRootName(string longName){
       cerr << "shortening " << longName << " to " << rootName << "\n";
     }
 #endif
-
+#endif
                if(rootName.find_last_of(".") != rootName.npos){
                        int pos = rootName.find_last_of('.')+1;
                        rootName = rootName.substr(0, pos);
@@ -632,7 +636,7 @@ int MothurOut::openInputFile(string fileName, ifstream& fileHandle, string m){
        try {
                        //get full path name
                        string completeFileName = getFullPathName(fileName);
-
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
 #ifdef USE_COMPRESSION
       // check for gzipped or bzipped file
       if (endsWith(completeFileName, ".gz") || endsWith(completeFileName, ".bz2")) {
@@ -655,7 +659,7 @@ int MothurOut::openInputFile(string fileName, ifstream& fileHandle, string m){
         }
       }
 #endif
-
+#endif
                        fileHandle.open(completeFileName.c_str());
                        if(!fileHandle) {
                                //mothurOut("[ERROR]: Could not open " + completeFileName); mothurOutEndLine();
@@ -678,7 +682,7 @@ int MothurOut::openInputFile(string fileName, ifstream& fileHandle){
 
                //get full path name
                string completeFileName = getFullPathName(fileName);
-
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
 #ifdef USE_COMPRESSION
   // check for gzipped or bzipped file
   if (endsWith(completeFileName, ".gz") || endsWith(completeFileName, ".bz2")) {
@@ -701,7 +705,7 @@ int MothurOut::openInputFile(string fileName, ifstream& fileHandle){
     }
   }
 #endif
-
+#endif
 
                fileHandle.open(completeFileName.c_str());
                if(!fileHandle) {
@@ -756,7 +760,7 @@ int MothurOut::openOutputFile(string fileName, ofstream& fileHandle){
        try { 
        
                string completeFileName = getFullPathName(fileName);
-
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
 #ifdef USE_COMPRESSION
     // check for gzipped file
     if (endsWith(completeFileName, ".gz") || endsWith(completeFileName, ".bz2")) {
@@ -776,7 +780,7 @@ int MothurOut::openOutputFile(string fileName, ofstream& fileHandle){
       }
     }
 #endif
-
+#endif
                fileHandle.open(completeFileName.c_str(), ios::trunc);
                if(!fileHandle) {
                        mothurOut("[ERROR]: Could not open " + completeFileName); mothurOutEndLine();
@@ -1151,7 +1155,90 @@ vector<vector<double> > MothurOut::binomial(int maxOrder){
                exit(1);
        }
 }
-
+/**************************************************************************************************/
+unsigned int MothurOut::fromBase36(string base36){
+       try {
+               unsigned int num = 0;
+               
+               map<char, int> converts;
+               converts['A'] = 0;
+               converts['a'] = 0;
+               converts['B'] = 1;
+               converts['b'] = 1;
+               converts['C'] = 2;
+               converts['c'] = 2;
+               converts['D'] = 3;
+               converts['d'] = 3;
+               converts['E'] = 4;
+               converts['e'] = 4;
+               converts['F'] = 5;
+               converts['f'] = 5;
+               converts['G'] = 6;
+               converts['g'] = 6;
+               converts['H'] = 7;
+               converts['h'] = 7;
+               converts['I'] = 8;
+               converts['i'] = 8;
+               converts['J'] = 9;
+               converts['j'] = 9;
+               converts['K'] = 10;
+               converts['k'] = 10;
+               converts['L'] = 11;
+               converts['l'] = 11;
+               converts['M'] = 12;
+               converts['m'] = 12;
+               converts['N'] = 13;
+               converts['n'] = 13;
+               converts['O'] = 14;
+               converts['o'] = 14;
+               converts['P'] = 15;
+               converts['p'] = 15;
+               converts['Q'] = 16;
+               converts['q'] = 16;
+               converts['R'] = 17;
+               converts['r'] = 17;
+               converts['S'] = 18;
+               converts['s'] = 18;
+               converts['T'] = 19;
+               converts['t'] = 19;
+               converts['U'] = 20;
+               converts['u'] = 20;
+               converts['V'] = 21;
+               converts['v'] = 21;
+               converts['W'] = 22;
+               converts['w'] = 22;
+               converts['X'] = 23;
+               converts['x'] = 23;
+               converts['Y'] = 24;
+               converts['y'] = 24;
+               converts['Z'] = 25;
+               converts['z'] = 25;
+               converts['0'] = 26;
+               converts['1'] = 27;
+               converts['2'] = 28;
+               converts['3'] = 29;
+               converts['4'] = 30;
+               converts['5'] = 31;
+               converts['6'] = 32;
+               converts['7'] = 33;
+               converts['8'] = 34;
+               converts['9'] = 35;             
+               
+               int i = 0;
+               while (i < base36.length()) {
+                       char c = base36[i];
+                       num = 36 * num + converts[c];
+                       i++;
+               }
+               
+               return num;
+               
+       }
+       catch(exception& e) {
+               errorOut(e, "MothurOut", "fromBase36");
+               exit(1);
+       }
+}
 /***********************************************************************/
 
 int MothurOut::factorial(int num){
@@ -1217,17 +1304,6 @@ void MothurOut::splitAtChar(string& estim, vector<string>& container, char symbo
                }
                container.push_back(individual);
 
-               /*
-               
-               while (estim.find_first_of(symbol) != -1) {
-                       individual = estim.substr(0,estim.find_first_of(symbol));
-                       if ((estim.find_first_of(symbol)+1) <= estim.length()) { //checks to make sure you don't have dash at end of string
-                               estim = estim.substr(estim.find_first_of(symbol)+1, estim.length());
-                               container.push_back(individual);
-                       }
-               }
-               //get last one
-               container.push_back(estim); */
        }
        catch(exception& e) {
                errorOut(e, "MothurOut", "splitAtChar");
@@ -1462,20 +1538,21 @@ bool MothurOut::anyLabelsToProcess(string label, set<string>& userLabels, string
                }
                
                //go through users set and make them floats
-               for(it = userLabels.begin(); it != userLabels.end(); ++it) {
+               for(it = userLabels.begin(); it != userLabels.end();) {
                        
                        float temp;
                        if ((*it != "unique") && (convertTestFloat(*it, temp) == true)){
                                convert(*it, temp);
                                orderFloat.push_back(temp);
                                userMap[*it] = temp;
+                               it++;
                        }else if (*it == "unique") { 
                                orderFloat.push_back(-1.0);
                                userMap["unique"] = -1.0;
+                               it++;
                        }else {
-                               if (errorOff == "") {  cout << *it << " is not a valid label." << endl;  }
-                               userLabels.erase(*it); 
-                               it--;
+                               if (errorOff == "") {  mothurOut(*it + " is not a valid label."); mothurOutEndLine();  }
+                               userLabels.erase(it++); 
                        }
                }
                
@@ -1491,11 +1568,11 @@ bool MothurOut::anyLabelsToProcess(string label, set<string>& userLabels, string
                        if (orderFloat[i] < labelFloat) {
                                smaller = true;
                                if (orderFloat[i] == -1) { 
-                                       if (errorOff == "") { cout << "Your file does not include the label unique." << endl; }
+                                       if (errorOff == "") { mothurOut("Your file does not include the label unique."); mothurOutEndLine(); }
                                        userLabels.erase("unique");
                                }
                                else {  
-                                       if (errorOff == "") { cout << "Your file does not include the label " << endl; }
+                                       if (errorOff == "") { mothurOut("Your file does not include the label "); mothurOutEndLine(); }
                                        string s = "";
                                        for (it2 = userMap.begin(); it2!= userMap.end(); it2++) {  
                                                if (it2->second == orderFloat[i]) {  
@@ -1505,7 +1582,7 @@ bool MothurOut::anyLabelsToProcess(string label, set<string>& userLabels, string
                                                        break;
                                                }
                                        }
-                                       if (errorOff == "") {cout << s <<  ". I will use the next smallest distance. " << endl; }
+                                       if (errorOff == "") {mothurOut( s +  ". I will use the next smallest distance. "); mothurOutEndLine(); }
                                }
                        //since they are sorted once you find a bigger one stop looking
                        }else { break; }