]> git.donarmstrong.com Git - mothur.git/blobdiff - engine.cpp
changed normalize.shared and sub.sample commands to fix bug in normalize.shared and...
[mothur.git] / engine.cpp
index 65349352d791f91aceeb850b114b697037984d41..af88e4f143924e5938740b52d4e76dbcd641f910 100644 (file)
@@ -25,14 +25,83 @@ Engine::Engine(){
                exit(1);
        }
 }
-
+/***********************************************************************/
+string Engine::findMothursPath(){
+       try { 
+               
+               string envPath = getenv("PATH");
+               string mothurPath = "";
+               
+               //delimiting path char
+               char delim;
+               #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
+                       delim = ':';
+               #else
+                       delim = ';';
+               #endif
+               
+               //break apart path variable by ':'
+               vector<string> dirs;
+               mout->splitAtChar(envPath, dirs, delim);
+               
+               //get path related to mothur
+               for (int i = 0; i < dirs.size(); i++) {
+                       //to lower so we can find it
+                       string tempLower = "";
+                       for (int j = 0; j < dirs[i].length(); j++) {  tempLower += tolower(dirs[i][j]);  }
+                       
+                       //is this mothurs path?
+                       if (tempLower.find("mothur") != -1) {  mothurPath = dirs[i]; break;  }
+               }
+               
+               if (mothurPath != "") {
+                       //add mothur so it looks like what argv would look like
+                       #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
+                               mothurPath += "/mothur";
+                       #else
+                               mothurPath += "\\mothur";
+                       #endif
+               }else {
+                       //okay mothur is not in the path, so the folder mothur is in must be in the path
+                       //lets find out which one
+                       
+                       //get path related to mothur
+                       for (int i = 0; i < dirs.size(); i++) {
+                                                               
+                               //is this mothurs path?
+                               ifstream in;
+                               string tempIn = dirs[i];
+                               #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
+                                       tempIn += "/mothur";
+                               #else
+                                       tempIn += "\\mothur";
+                               #endif
+                               mout->openInputFile(tempIn, in, "");
+                               
+                               //if this file exists
+                               if (in) { in.close(); mothurPath = tempIn;  break;  }
+                       }
+               }
+               
+               return mothurPath;
+               
+       }
+       catch(exception& e) {
+               mout->errorOut(e, "Engine", "findMothursPath");
+               exit(1);
+       }
+}
 /***********************************************************************/
 
 InteractEngine::InteractEngine(string path){
 
-       globaldata = GlobalData::getInstance();
-       globaldata->argv = path;
        
+       string temppath = path.substr(0, (path.find_last_of("othur")-5));
+       
+       //this will happen if you set the path variable to contain mothur's exe location
+       if (temppath == "") { path = findMothursPath(); }
+       
+       mout->argv = path;
 }
 
 /***********************************************************************/
@@ -51,6 +120,16 @@ bool InteractEngine::getInput(){
                
                while(quitCommandCalled != 1){
 
+                       #ifdef USE_MPI
+                               int pid, processors;
+                               MPI_Status status;
+                               MPI_Comm_rank(MPI_COMM_WORLD, &pid); 
+                               MPI_Comm_size(MPI_COMM_WORLD, &processors);
+                               
+                               if (pid == 0) {
+                               
+                       #endif
+                       
                        mout->mothurOutEndLine();
                        
                        input = getCommand();   
@@ -60,7 +139,31 @@ bool InteractEngine::getInput(){
                        
                        //allow user to omit the () on the quit command
                        if (input == "quit") { input = "quit()"; }
+
+                       
+                       #ifdef USE_MPI
+                               //send commandName
+                               for(int i = 1; i < processors; i++) { 
+                                               int length = input.length();
+                                               MPI_Send(&length, 1, MPI_INT, i, 2001, MPI_COMM_WORLD);
+                                               MPI_Send(&input[0], length, MPI_CHAR, i, 2001, MPI_COMM_WORLD);
+       
+                                       }
+                               }else {
+                                       int length;
+                                       MPI_Recv(&length, 1, MPI_INT, 0, 2001, MPI_COMM_WORLD, &status);
+                                       //recieve container
+                                       char* tempBuf = new char[length];
+                                       MPI_Recv(&tempBuf[0], length, MPI_CHAR, 0, 2001, MPI_COMM_WORLD, &status);
+                                       
+                                       input = tempBuf;
+                                       if (input.length() > length) { input = input.substr(0, length);  }
+                                       delete tempBuf; 
+                               }
+
                        
+                       #endif
+               
                        CommandOptionParser parser(input);
                        commandName = parser.getCommandString();
        
@@ -71,13 +174,27 @@ bool InteractEngine::getInput(){
                                        #ifdef USE_MPI
                                                int pid;
                                                MPI_Comm_rank(MPI_COMM_WORLD, &pid); 
-//cout << pid << " is here " << commandName << endl;
+                                               
                                                if ((cFactory->MPIEnabled(commandName)) || (pid == 0)) {
+                                       //cout << pid << " is in execute " << commandName << endl;
                                        #endif
                                        //executes valid command
+                                       mout->runParse = true;
+                                       mout->Groups.clear();
+                                       mout->namesOfGroups.clear();
+                                       mout->Treenames.clear();
+                                       mout->names.clear();
+                                       mout->saveNextLabel = "";
+                                       mout->printedHeaders = false;
+                                       mout->currentBinLabels.clear();
+                                       mout->binLabelsInFile.clear();
+                                                       
                                        Command* command = cFactory->getCommand(commandName, options);
                                        quitCommandCalled = command->execute();
+                                                       
+                                       //if we aborted command
+                                       if (quitCommandCalled == 2) {  mout->mothurOut("[ERROR]: did not complete " + commandName + "."); mout->mothurOutEndLine(); }
+
                                        mout->control_pressed = 0;
                                        mout->executing = false;
                                                                                
@@ -99,9 +216,6 @@ bool InteractEngine::getInput(){
 /***********************************************************************/
 string Engine::getCommand()  {
        try {
-               #ifdef USE_MPI
-                       MPI_Barrier(MPI_COMM_WORLD);
-               #endif
        
                #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
                        #ifdef USE_READLINE
@@ -110,7 +224,7 @@ string Engine::getCommand()  {
                                
                                if(nextCommand != NULL) {  add_history(nextCommand);  } 
                                else{ //^D causes null string and we want it to quit mothur
-                                       nextCommand = "quit"; 
+                                       nextCommand = strdup("quit");
                                        mout->mothurOut(nextCommand);
                                }       
                                
@@ -145,10 +259,15 @@ string Engine::getCommand()  {
 //This function opens the batchfile to be used by BatchEngine::getInput.
 BatchEngine::BatchEngine(string path, string batchFileName){
        try {
-               globaldata = GlobalData::getInstance();
        
-               openedBatch = openInputFile(batchFileName, inputBatchFile);
-               globaldata->argv = path;
+               openedBatch = mout->openInputFile(batchFileName, inputBatchFile);
+               
+               string temppath = path.substr(0, (path.find_last_of("othur")-5));
+       
+               //this will happen if you set the path variable to contain mothur's exe location
+               if (temppath == "") { path = findMothursPath(); }
+               
+               mout->argv = path;
                                
        }
        catch(exception& e) {
@@ -179,11 +298,46 @@ bool BatchEngine::getInput(){
                //CommandFactory cFactory;
                int quitCommandCalled = 0;
            int count = 0;
-               while(quitCommandCalled == 0){
-       
+               while(quitCommandCalled != 1){
+                       
+                       #ifdef USE_MPI
+                               int pid, processors;
+                               MPI_Status status;
+                               MPI_Comm_rank(MPI_COMM_WORLD, &pid); 
+                               MPI_Comm_size(MPI_COMM_WORLD, &processors);
+                               
+                               if (pid == 0) {
+                               
+                       #endif
+                       
                        input = getNextCommand(inputBatchFile);
                        count++;
                        
+                       #ifdef USE_MPI
+                               //send commandName
+                               for(int i = 1; i < processors; i++) { 
+                                               int length = input.length();
+                                               MPI_Send(&length, 1, MPI_INT, i, 2001, MPI_COMM_WORLD);
+                                               MPI_Send(&input[0], length, MPI_CHAR, i, 2001, MPI_COMM_WORLD);
+       
+                                       }
+                               }else {
+                                       int length;
+                                       MPI_Recv(&length, 1, MPI_INT, 0, 2001, MPI_COMM_WORLD, &status);
+                                       //recieve container
+                                       char* tempBuf = new char[length];
+                                       MPI_Recv(&tempBuf[0], length, MPI_CHAR, 0, 2001, MPI_COMM_WORLD, &status);
+                                       
+                                       input = tempBuf;
+                                       if (input.length() > length) { input = input.substr(0, length);  }
+                                       delete tempBuf; 
+                               }
+
+                       
+                       #endif
+
+                       
+                       
                        if (input[0] != '#') {
                                
                                mout->mothurOutEndLine();
@@ -209,8 +363,23 @@ bool BatchEngine::getInput(){
                                                if ((cFactory->MPIEnabled(commandName)) || (pid == 0)) {
                                        #endif
                                        //executes valid command
+                                       mout->runParse = true;
+                                       mout->Groups.clear();
+                                       mout->namesOfGroups.clear();
+                                       mout->Treenames.clear();
+                                       mout->names.clear();
+                                       mout->saveNextLabel = "";
+                                       mout->printedHeaders = false;
+                                       mout->currentBinLabels.clear();
+                                       mout->binLabelsInFile.clear();
+
+                                                       
                                        Command* command = cFactory->getCommand(commandName, options);
                                        quitCommandCalled = command->execute();
+                                                       
+                                       //if we aborted command
+                                       if (quitCommandCalled == 2) {  mout->mothurOut("[ERROR]: did not complete " + commandName + "."); mout->mothurOutEndLine(); }
+
                                        mout->control_pressed = 0;
                                        mout->executing = false;
                                                                                
@@ -223,7 +392,7 @@ bool BatchEngine::getInput(){
                                }
                                
                        }
-                       gobble(inputBatchFile);
+                       mout->gobble(inputBatchFile);
                }
                
                inputBatchFile.close();
@@ -238,15 +407,10 @@ bool BatchEngine::getInput(){
 string BatchEngine::getNextCommand(ifstream& inputBatchFile) {
        try {
                
-               #ifdef USE_MPI
-                       int err = MPI_Barrier(MPI_COMM_WORLD);
-//cout << "barrier = " << err << '\t' << MPI_SUCCESS << endl;
-               #endif
-               
                string nextcommand = "";
                
                if (inputBatchFile.eof()) { nextcommand = "quit()"; }
-               else { nextcommand = getline(inputBatchFile); }
+               else { nextcommand = mout->getline(inputBatchFile); }
                
                return nextcommand;
        }
@@ -261,12 +425,16 @@ string BatchEngine::getNextCommand(ifstream& inputBatchFile) {
 //This function opens the batchfile to be used by BatchEngine::getInput.
 ScriptEngine::ScriptEngine(string path, string commandString){
        try {
-               globaldata = GlobalData::getInstance();
                
                //remove quotes
                listOfCommands = commandString.substr(1, (commandString.length()-1));
+               
+               string temppath = path.substr(0, (path.find_last_of("othur")-5));
 
-               globaldata->argv = path;
+               //this will happen if you set the path variable to contain mothur's exe location
+               if (temppath == "") { path = findMothursPath(); }
+               
+               mout->argv = path;
                                
        }
        catch(exception& e) {
@@ -292,15 +460,53 @@ bool ScriptEngine::getInput(){
                //CommandFactory cFactory;
                int quitCommandCalled = 0;
        
-               while(quitCommandCalled == 0){
-               
+               while(quitCommandCalled != 1){
+                       
+                       #ifdef USE_MPI
+                               int pid, processors;
+                               MPI_Status status;
+                               MPI_Comm_rank(MPI_COMM_WORLD, &pid); 
+                               MPI_Comm_size(MPI_COMM_WORLD, &processors);
+                               
+                               if (pid == 0) {
+                               
+                       #endif
+                       
                        input = getNextCommand(listOfCommands); 
                        
                        if (input == "") { input = "quit()"; }
                        
-                       mout->mothurOutEndLine();
-                       mout->mothurOut("mothur > " + input);
-                       mout->mothurOutEndLine();
+                       if (mout->gui) {
+                               if ((input.find("quit") != string::npos) || (input.find("set.logfile") != string::npos)) {}
+                               else if ((input.find("get.current") != string::npos) && (!mout->hasCurrentFiles())) {}
+                               else { mout->mothurOutEndLine(); mout->mothurOut("mothur > " + input); mout->mothurOutEndLine(); }
+                       }else{
+                               mout->mothurOutEndLine(); mout->mothurOut("mothur > " + input); mout->mothurOutEndLine();
+                       }
+                       
+                       #ifdef USE_MPI
+                               //send commandName
+                               for(int i = 1; i < processors; i++) { 
+                                               int length = input.length();
+                                               MPI_Send(&length, 1, MPI_INT, i, 2001, MPI_COMM_WORLD);
+                                               MPI_Send(&input[0], length, MPI_CHAR, i, 2001, MPI_COMM_WORLD);
+       
+                                       }
+                               }else {
+                                       int length;
+                                       MPI_Recv(&length, 1, MPI_INT, 0, 2001, MPI_COMM_WORLD, &status);
+                                       //recieve container
+                                       char* tempBuf = new char[length];
+                                       MPI_Recv(&tempBuf[0], length, MPI_CHAR, 0, 2001, MPI_COMM_WORLD, &status);
+                                       
+                                       input = tempBuf;
+                                       if (input.length() > length) { input = input.substr(0, length);  }
+                                       delete tempBuf; 
+                               }
+
+                       
+                       #endif
+                       
                        
                        if (mout->control_pressed) { input = "quit()"; }
                                
@@ -315,7 +521,6 @@ bool ScriptEngine::getInput(){
                                        mout->executing = true;
                                        #ifdef USE_MPI
                                                int pid, numProcesses;
-                                               MPI_Status status; 
                                                
                                                MPI_Comm_rank(MPI_COMM_WORLD, &pid); 
                                                MPI_Comm_size(MPI_COMM_WORLD, &numProcesses); 
@@ -325,13 +530,28 @@ bool ScriptEngine::getInput(){
                                                        //cout << pid << " is in execute" << endl;      
                                        #endif
                                        //executes valid command
+                                       mout->runParse = true;
+                                       mout->Groups.clear();
+                                       mout->namesOfGroups.clear();
+                                       mout->Treenames.clear();
+                                       mout->names.clear();
+                                       mout->saveNextLabel = "";
+                                       mout->printedHeaders = false;
+                                       mout->currentBinLabels.clear();
+                                       mout->binLabelsInFile.clear();
+
+                                                       
                                        Command* command = cFactory->getCommand(commandName, options);
                                        quitCommandCalled = command->execute();
+                                       
+                                       //if we aborted command
+                                       if (quitCommandCalled == 2) {  mout->mothurOut("[ERROR]: did not complete " + commandName + "."); mout->mothurOutEndLine(); }
+                                                       
                                        mout->control_pressed = 0;
                                        mout->executing = false;
                                                                        
                                        #ifdef USE_MPI
-                                       cout << pid << " is done in execute" << endl;
+                                       //cout << pid << " is done in execute" << endl;
                                                }
                                        #endif
                                }else {         
@@ -353,17 +573,18 @@ bool ScriptEngine::getInput(){
 string ScriptEngine::getNextCommand(string& commandString) {
        try {
                
-               #ifdef USE_MPI
-                       MPI_Barrier(MPI_COMM_WORLD);
-               #endif
-               
                string nextcommand = "";
                int count = 0;
+               bool ignoreSemiColons = false;
                
                //go through string until you reach ; or end
                while (count < commandString.length()) { 
                        
-                       if (commandString[count] == ';') {  break;   }
+                        //you want to ignore any ; until you reach the next '
+                       if ((commandString[count] == '\'') && (!ignoreSemiColons)) {  ignoreSemiColons = true;  } 
+                       else if ((commandString[count] == '\'') && (ignoreSemiColons)) {  ignoreSemiColons = false;  } 
+                               
+                       if ((commandString[count] == ';') && (!ignoreSemiColons)) {  break;   }
                        else {          nextcommand += commandString[count];    }
                        
                        count++;