]> git.donarmstrong.com Git - mothur.git/blob - mothur.cpp
fixed sub.sample groups bug
[mothur.git] / mothur.cpp
1 /*
2  *  interface.cpp
3  *  
4  *
5  *  Created by Pat Schloss on 8/14/08.
6  *  Copyright 2008 Patrick D. Schloss. All rights reserved.
7  *
8  */
9  
10 #include "mothur.h"
11 #include "engine.hpp"
12 #include "globaldata.hpp"
13 #include "mothurout.h"
14
15 /**************************************************************************************************/
16
17 CommandFactory* CommandFactory::_uniqueInstance = 0;
18 MothurOut* MothurOut::_uniqueInstance = 0;
19
20 /***********************************************************************/
21 volatile int ctrlc_pressed = 0;
22 void ctrlc_handler ( int sig ) {
23         MothurOut* m = MothurOut::getInstance();
24     ctrlc_pressed = 1;
25         m->control_pressed = ctrlc_pressed;
26         
27         if (m->executing) { //if mid command quit execution, else quit mothur
28                 m->mothurOutEndLine(); m->mothurOut("quitting command...");  m->mothurOutEndLine();
29         }else{
30                 m->mothurOut("quitting mothur");  m->mothurOutEndLine();
31                 exit(1);
32         }
33 }
34 /***********************************************************************/
35 int main(int argc, char *argv[]){
36         MothurOut* m = MothurOut::getInstance();
37         try {   
38                 signal(SIGINT, ctrlc_handler );
39                                 
40                 time_t ltime = time(NULL); /* calendar time */  
41                 string logFileName = "mothur." + toString(ltime) + ".logfile";
42                 
43                 #ifdef USE_MPI
44                         MPI_Init(&argc, &argv); 
45                 #endif
46
47                 m->setFileName(logFileName);
48                 
49                 #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
50                         system("clear");
51                 #else
52                         system("CLS");
53                 #endif
54                 
55                 #ifdef MOTHUR_FILES
56                         string temp = MOTHUR_FILES; 
57                 
58                         //add / to name if needed
59                         string lastChar = temp.substr(temp.length()-1);
60                         #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
61                                 if (lastChar != "/") { temp += "/"; }
62                         #else
63                                 if (lastChar != "\\") { temp += "\\"; } 
64                         #endif
65                 
66                         temp = m->getFullPathName(temp);
67                         m->setDefaultPath(temp);
68                 #endif
69                 
70                 #ifdef USE_MPI
71                         int version, subversion;
72                         MPI_Get_version(&version, &subversion);
73                 #endif
74                 
75                 //get releaseDate from Make
76                 string releaseDate = RELEASE_DATE; 
77                 string mothurVersion = VERSION; 
78                 m->setReleaseDate(releaseDate);
79                 m->setVersion(mothurVersion);
80                 
81                 //will make the gui output "pretty"
82                 bool guiMode = false;
83                 if (argc>1) {
84                         string guiInput = argv[1];
85                         if (guiInput[0] == '+') { guiMode = true; }
86                 }
87                 
88                 if (!guiMode) {
89                         //version
90                         #if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
91                                 #if defined (__APPLE__) || (__MACH__)
92                                         m->mothurOutJustToLog("Mac version");
93                                         m->mothurOutEndLine(); m->mothurOutEndLine();
94                                 #else
95                                         m->mothurOutJustToLog("Linux version");
96                                         m->mothurOutEndLine(); m->mothurOutEndLine();
97                                 #endif
98
99                         #else
100                                 m->mothurOutJustToLog("Windows version");
101                                 m->mothurOutEndLine(); m->mothurOutEndLine();
102                         #endif          
103                         
104                         #ifdef USE_READLINE
105                                 m->mothurOutJustToLog("Using ReadLine");
106                                 m->mothurOutEndLine(); m->mothurOutEndLine();
107                         #endif
108                         
109                         #ifdef MOTHUR_FILES
110                                 m->mothurOutJustToLog("Using default file location " + temp);
111                                 m->mothurOutEndLine(); m->mothurOutEndLine();
112                         #endif
113                         
114                         #ifdef BIT_VERSION
115                                 m->mothurOutJustToLog("Running 64Bit Version");
116                                 m->mothurOutEndLine(); m->mothurOutEndLine();
117                         #else
118                                 m->mothurOutJustToLog("Running 32Bit Version");
119                                 m->mothurOutEndLine(); m->mothurOutEndLine();
120                         #endif
121                         
122                         //header
123                         m->mothurOut("mothur v." + mothurVersion);
124                         m->mothurOutEndLine();          
125                         m->mothurOut("Last updated: " + releaseDate);
126                         m->mothurOutEndLine();  
127                         m->mothurOutEndLine();          
128                         m->mothurOut("by");
129                         m->mothurOutEndLine();          
130                         m->mothurOut("Patrick D. Schloss");
131                         m->mothurOutEndLine();
132                         m->mothurOutEndLine();                  
133                         m->mothurOut("Department of Microbiology & Immunology");
134                         m->mothurOutEndLine();  
135                         m->mothurOut("University of Michigan");
136                         m->mothurOutEndLine();                  
137                         m->mothurOut("pschloss@umich.edu");
138                         m->mothurOutEndLine();          
139                         m->mothurOut("http://www.mothur.org");
140                         m->mothurOutEndLine();
141                         m->mothurOutEndLine();
142                         m->mothurOut("When using, please cite:");
143                         m->mothurOutEndLine();
144                         m->mothurOut("Schloss, P.D., et al., Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol, 2009. 75(23):7537-41.");
145                         m->mothurOutEndLine();  
146                         m->mothurOutEndLine();          
147                         m->mothurOut("Distributed under the GNU General Public License");
148                         m->mothurOutEndLine();
149                         m->mothurOutEndLine();                  
150                         m->mothurOut("Type 'help()' for information on the commands that are available");
151                         m->mothurOutEndLine();
152                         m->mothurOutEndLine();                  
153                         m->mothurOut("Type 'quit()' to exit program");
154                         m->mothurOutEndLine();  
155                         
156                         #ifdef USE_MPI
157                                 m->mothurOutJustToLog("Using MPI\tversion ");
158                                 m->mothurOutJustToLog(toString(version) + "." + toString(subversion) + "\n");
159                         #endif
160                 }
161                 
162                 //srand(54321);
163                 srand( (unsigned)time( NULL ) );
164                 
165                 Engine* mothur;
166                 bool bail = 0;
167                 string input;
168  
169                 if(argc>1){
170                         input = argv[1];
171                         //m->mothurOut("input = " + input); m->mothurOutEndLine();
172
173                         if (input[0] == '#') {
174                                 m->mothurOutJustToLog("Script Mode");
175                                 m->mothurOutEndLine(); m->mothurOutEndLine();
176
177                                 mothur = new ScriptEngine(argv[0], argv[1]);
178                         }else if (input[0] == '+') {
179                                         mothur = new ScriptEngine(argv[0], argv[1]);
180                                         m->gui = true;
181                         }else{
182                                 m->mothurOutJustToLog("Batch Mode");
183                                 m->mothurOutEndLine(); m->mothurOutEndLine();
184                                 
185                                 mothur = new BatchEngine(argv[0], argv[1]);
186                         }
187                 }
188                 else{
189                         m->mothurOutJustToLog("Interactive Mode");
190                         m->mothurOutEndLine(); m->mothurOutEndLine();
191                         
192                         mothur = new InteractEngine(argv[0]);   
193                 }
194                 
195                 while(bail == 0)        {       bail = mothur->getInput();      }
196                 
197                 //closes logfile so we can rename
198                 m->closeLog();
199                 
200                 string outputDir = mothur->getOutputDir();
201                 string tempLog = mothur->getLogFileName();
202                 bool append = mothur->getAppend();
203                 
204                 string newlogFileName;
205                 if (tempLog != "") {
206                         newlogFileName = outputDir + tempLog;
207                         
208                         if (!append) {  
209                                 //need this because m->mothurOut makes the logfile, but doesn't know where to put it
210                                 rename(logFileName.c_str(), newlogFileName.c_str()); //logfile with timestamp
211
212                         }else {
213                                 ofstream outNewLog;
214                                 m->openOutputFileAppend(newlogFileName, outNewLog);
215                                 
216                                 if (!m->gui) {
217                                         outNewLog << endl << endl << "*********************************************************************************" << endl << endl;
218                                 }
219                                 outNewLog.close();
220                                 
221                                 m->appendFiles(logFileName, newlogFileName);
222                                 remove(logFileName.c_str());
223                         }
224                 }else{  
225                         newlogFileName = outputDir + logFileName;
226                         //need this because m->mothurOut makes the logfile, but doesn't know where to put it
227                         rename(logFileName.c_str(), newlogFileName.c_str()); //logfile with timestamp
228                 }
229                 
230                                 
231                 delete mothur;
232                 
233                 #ifdef USE_MPI
234                         MPI_Finalize();
235                 #endif
236                 
237                 return 0;
238         }
239         catch(exception& e) {
240                 m->errorOut(e, "mothur", "main");
241                 exit(1);
242         }
243 }
244
245 /**************************************************************************************************/
246