5 * Created by Pat Schloss on 10/25/08.
6 * Copyright 2008 Patrick D. Schloss. All rights reserved.
10 #include "command.hpp"
11 #include "readdistcommand.h"
12 #include "readtreecommand.h"
13 #include "readotucommand.h"
14 #include "clustercommand.h"
15 #include "collectcommand.h"
16 #include "collectsharedcommand.h"
17 #include "getgroupcommand.h"
18 #include "getlabelcommand.h"
19 #include "rarefactcommand.h"
20 #include "summarycommand.h"
21 #include "summarysharedcommand.h"
22 #include "rarefactsharedcommand.h"
23 #include "quitcommand.h"
24 #include "helpcommand.h"
25 #include "commandfactory.hpp"
26 #include "deconvolutecommand.h"
27 #include "parsimonycommand.h"
28 #include "unifracunweightedcommand.h"
29 #include "unifracweightedcommand.h"
30 #include "libshuffcommand.h"
31 #include "heatmapcommand.h"
32 #include "heatmapsimcommand.h"
33 #include "filterseqscommand.h"
34 #include "venncommand.h"
35 #include "nocommands.h"
36 #include "binsequencecommand.h"
37 #include "getoturepcommand.h"
38 #include "treegroupscommand.h"
39 #include "bootstrapsharedcommand.h"
40 //#include "consensuscommand.h"
41 #include "distancecommand.h"
42 #include "aligncommand.h"
43 #include "matrixoutputcommand.h"
44 #include "getsabundcommand.h"
45 #include "getrabundcommand.h"
46 #include "seqsummarycommand.h"
47 #include "screenseqscommand.h"
48 #include "reversecommand.h"
49 #include "trimseqscommand.h"
50 #include "mergefilecommand.h"
51 #include "chimeraseqscommand.h"
52 #include "listseqscommand.h"
53 #include "getseqscommand.h"
54 #include "removeseqscommand.h"
55 #include "systemcommand.h"
56 #include "secondarystructurecommand.h"
57 #include "getsharedotucommand.h"
58 #include "getlistcountcommand.h"
59 #include "hclustercommand.h"
60 #include "classifyseqscommand.h"
61 #include "phylotypecommand.h"
62 #include "mgclustercommand.h"
63 #include "preclustercommand.h"
64 #include "pcacommand.h"
65 #include "otuhierarchycommand.h"
66 #include "setdircommand.h"
67 #include "parselistscommand.h"
68 #include "parsesffcommand.h"
69 #include "chimeraccodecommand.h"
70 #include "chimeracheckcommand.h"
71 #include "chimeraslayercommand.h"
72 #include "chimerapintailcommand.h"
73 #include "chimerabellerophoncommand.h"
74 #include "setlogfilecommand.h"
75 #include "phylodiversitycommand.h"
77 /*******************************************************/
79 /******************************************************/
80 CommandFactory* CommandFactory::getInstance() {
81 if( _uniqueInstance == 0) {
82 _uniqueInstance = new CommandFactory();
84 return _uniqueInstance;
86 /***********************************************************/
88 /***********************************************************/
89 //note: This class is resposible for knowing which commands are mpiEnabled,
90 //If a command is not enabled only process 0 will execute the command.
91 //This avoids redundant outputs on pieces of code we have not paralellized.
92 //If you add mpi code to a existing command you need to modify the list below or the code will hang on MPI blocking commands like FIle_open.
93 //example: commands["dist.seqs"] = "MPIEnabled";
95 CommandFactory::CommandFactory(){
97 m = MothurOut::getInstance();
99 command = new NoCommand(s);
101 outputDir = ""; inputDir = "";
105 //initialize list of valid commands
106 commands["read.dist"] = "read.dist";
107 commands["read.otu"] = "read.otu";
108 commands["read.tree"] = "read.tree";
109 commands["bin.seqs"] = "bin.seqs";
110 commands["get.oturep"] = "get.oturep";
111 commands["cluster"] = "cluster";
112 commands["unique.seqs"] = "unique.seqs";
113 commands["dist.shared"] = "dist.shared";
114 commands["collect.single"] = "collect.single";
115 commands["collect.shared"] = "collect.shared";
116 commands["rarefaction.single"] = "rarefaction.single";
117 commands["rarefaction.shared"] = "rarefaction.shared";
118 commands["summary.single"] = "summary.single";
119 commands["summary.shared"] = "summary.shared";
120 commands["parsimony"] = "parsimony";
121 commands["unifrac.weighted"] = "unifrac.weighted";
122 commands["unifrac.unweighted"] = "unifrac.unweighted";
123 commands["libshuff"] = "libshuff";
124 commands["tree.shared"] = "tree.shared";
125 commands["heatmap.bin"] = "heatmap.bin";
126 commands["heatmap.sim"] = "heatmap.sim";
127 commands["venn"] = "venn";
128 commands["get.group"] = "get.group";
129 commands["get.label"] = "get.label";
130 commands["get.sabund"] = "get.sabund";
131 commands["get.rabund"] = "get.rabund";
132 commands["bootstrap.shared"] = "bootstrap.shared";
133 //commands["consensus"] = "consensus";
134 commands["help"] = "help";
135 commands["reverse.seqs"] = "reverse.seqs";
136 commands["trim.seqs"] = "trim.seqs";
137 commands["list.seqs"] = "list.seqs";
138 commands["get.seqs"] = "get.seqs";
139 commands["remove.seqs"] = "get.seqs";
140 commands["system"] = "system";
141 commands["align.check"] = "align.check";
142 commands["get.sharedseqs"] = "get.sharedseqs";
143 commands["get.otulist"] = "get.otulist";
144 commands["hcluster"] = "hcluster";
145 commands["phylotype"] = "phylotype";
146 commands["mgcluster"] = "mgcluster";
147 commands["pre.cluster"] = "pre.cluster";
148 commands["pcoa"] = "pcoa";
149 commands["otu.hierarchy"] = "otu.hierarchy";
150 commands["set.dir"] = "set.dir";
151 commands["merge.files"] = "merge.files";
152 commands["parse.list"] = "parse.list";
153 commands["parse.sff"] = "parse.sff";
154 commands["set.logfile"] = "set.logfile";
155 commands["phylo.diversity"] = "phylo.diversity";
156 commands["classify.seqs"] = "MPIEnabled";
157 commands["dist.seqs"] = "MPIEnabled";
158 commands["filter.seqs"] = "MPIEnabled";
159 commands["align.seqs"] = "MPIEnabled";
160 commands["chimera.seqs"] = "chimera.seqs";
161 commands["chimera.ccode"] = "MPIEnabled";
162 commands["chimera.check"] = "MPIEnabled";
163 commands["chimera.slayer"] = "MPIEnabled";
164 commands["chimera.pintail"] = "MPIEnabled";
165 commands["chimera.bellerophon"] = "MPIEnabled";
166 commands["screen.seqs"] = "MPIEnabled";
167 commands["summary.seqs"] = "MPIEnabled";
168 commands["quit"] = "MPIEnabled";
171 /***********************************************************/
173 /***********************************************************/
174 bool CommandFactory::MPIEnabled(string commandName) {
176 it = commands.find(commandName);
177 if (it != commands.end()) {
178 if (it->second == "MPIEnabled") { return true; }
182 /***********************************************************/
184 /***********************************************************/
185 CommandFactory::~CommandFactory(){
190 /***********************************************************/
192 /***********************************************************/
193 //This function calls the appropriate command fucntions based on user input.
194 Command* CommandFactory::getCommand(string commandName, string optionString){
196 delete command; //delete the old command
198 //user has opted to redirect output from dir where input files are located to some other place
199 if (outputDir != "") {
200 if (optionString != "") { optionString += ", outputdir=" + outputDir; }
201 else { optionString += "outputdir=" + outputDir; }
204 //user has opted to redirect input from dir where mothur.exe is located to some other place
205 if (inputDir != "") {
206 if (optionString != "") { optionString += ", inputdir=" + inputDir; }
207 else { optionString += "inputdir=" + inputDir; }
210 if(commandName == "read.dist") { command = new ReadDistCommand(optionString); }
211 else if(commandName == "read.otu") { command = new ReadOtuCommand(optionString); }
212 else if(commandName == "read.tree") { command = new ReadTreeCommand(optionString); }
213 else if(commandName == "cluster") { command = new ClusterCommand(optionString); }
214 else if(commandName == "unique.seqs") { command = new DeconvoluteCommand(optionString); }
215 else if(commandName == "parsimony") { command = new ParsimonyCommand(optionString); }
216 else if(commandName == "help") { command = new HelpCommand(optionString); }
217 else if(commandName == "quit") { command = new QuitCommand(optionString); }
218 else if(commandName == "collect.single") { command = new CollectCommand(optionString); }
219 else if(commandName == "collect.shared") { command = new CollectSharedCommand(optionString); }
220 else if(commandName == "rarefaction.single") { command = new RareFactCommand(optionString); }
221 else if(commandName == "rarefaction.shared") { command = new RareFactSharedCommand(optionString); }
222 else if(commandName == "summary.single") { command = new SummaryCommand(optionString); }
223 else if(commandName == "summary.shared") { command = new SummarySharedCommand(optionString); }
224 else if(commandName == "unifrac.weighted") { command = new UnifracWeightedCommand(optionString); }
225 else if(commandName == "unifrac.unweighted") { command = new UnifracUnweightedCommand(optionString); }
226 else if(commandName == "get.group") { command = new GetgroupCommand(optionString); }
227 else if(commandName == "get.label") { command = new GetlabelCommand(optionString); }
228 else if(commandName == "get.sabund") { command = new GetSAbundCommand(optionString); }
229 else if(commandName == "get.rabund") { command = new GetRAbundCommand(optionString); }
230 else if(commandName == "libshuff") { command = new LibShuffCommand(optionString); }
231 else if(commandName == "heatmap.bin") { command = new HeatMapCommand(optionString); }
232 else if(commandName == "heatmap.sim") { command = new HeatMapSimCommand(optionString); }
233 else if(commandName == "filter.seqs") { command = new FilterSeqsCommand(optionString); }
234 else if(commandName == "venn") { command = new VennCommand(optionString); }
235 else if(commandName == "bin.seqs") { command = new BinSeqCommand(optionString); }
236 else if(commandName == "get.oturep") { command = new GetOTURepCommand(optionString); }
237 else if(commandName == "tree.shared") { command = new TreeGroupCommand(optionString); }
238 else if(commandName == "dist.shared") { command = new MatrixOutputCommand(optionString); }
239 else if(commandName == "bootstrap.shared") { command = new BootSharedCommand(optionString); }
240 //else if(commandName == "consensus") { command = new ConcensusCommand(optionString); }
241 else if(commandName == "dist.seqs") { command = new DistanceCommand(optionString); }
242 else if(commandName == "align.seqs") { command = new AlignCommand(optionString); }
243 else if(commandName == "summary.seqs") { command = new SeqSummaryCommand(optionString); }
244 else if(commandName == "screen.seqs") { command = new ScreenSeqsCommand(optionString); }
245 else if(commandName == "reverse.seqs") { command = new ReverseSeqsCommand(optionString); }
246 else if(commandName == "trim.seqs") { command = new TrimSeqsCommand(optionString); }
247 else if(commandName == "chimera.seqs") { command = new ChimeraSeqsCommand(optionString); }
248 else if(commandName == "list.seqs") { command = new ListSeqsCommand(optionString); }
249 else if(commandName == "get.seqs") { command = new GetSeqsCommand(optionString); }
250 else if(commandName == "remove.seqs") { command = new RemoveSeqsCommand(optionString); }
251 else if(commandName == "merge.files") { command = new MergeFileCommand(optionString); }
252 else if(commandName == "system") { command = new SystemCommand(optionString); }
253 else if(commandName == "align.check") { command = new AlignCheckCommand(optionString); }
254 else if(commandName == "get.sharedseqs") { command = new GetSharedOTUCommand(optionString); }
255 else if(commandName == "get.otulist") { command = new GetListCountCommand(optionString); }
256 else if(commandName == "hcluster") { command = new HClusterCommand(optionString); }
257 else if(commandName == "classify.seqs") { command = new ClassifySeqsCommand(optionString); }
258 else if(commandName == "chimera.ccode") { command = new ChimeraCcodeCommand(optionString); }
259 else if(commandName == "chimera.check") { command = new ChimeraCheckCommand(optionString); }
260 else if(commandName == "chimera.slayer") { command = new ChimeraSlayerCommand(optionString); }
261 else if(commandName == "chimera.pintail") { command = new ChimeraPintailCommand(optionString); }
262 else if(commandName == "chimera.bellerophon") { command = new ChimeraBellerophonCommand(optionString); }
263 else if(commandName == "phylotype") { command = new PhylotypeCommand(optionString); }
264 else if(commandName == "mgcluster") { command = new MGClusterCommand(optionString); }
265 else if(commandName == "pre.cluster") { command = new PreClusterCommand(optionString); }
266 else if(commandName == "pcoa") { command = new PCACommand(optionString); }
267 else if(commandName == "otu.hierarchy") { command = new OtuHierarchyCommand(optionString); }
268 else if(commandName == "set.dir") { command = new SetDirectoryCommand(optionString); }
269 else if(commandName == "set.logfile") { command = new SetLogFileCommand(optionString); }
270 else if(commandName == "parse.list") { command = new ParseListCommand(optionString); }
271 else if(commandName == "parse.sff") { command = new ParseSFFCommand(optionString); }
272 else if(commandName == "phylo.diversity") { command = new PhyloDiversityCommand(optionString); }
273 else { command = new NoCommand(optionString); }
277 catch(exception& e) {
278 m->errorOut(e, "CommandFactory", "getCommand");
282 /***********************************************************/
283 //This function is used to interrupt a command
284 Command* CommandFactory::getCommand(){
286 delete command; //delete the old command
289 command = new NoCommand(s);
293 catch(exception& e) {
294 m->errorOut(e, "CommandFactory", "getCommand");
298 /***********************************************************************/
299 bool CommandFactory::isValidCommand(string command) {
302 //is the command in the map
303 if ((commands.find(command)) != (commands.end())) {
306 m->mothurOut(command + " is not a valid command in Mothur. Valid commands are ");
307 for (it = commands.begin(); it != commands.end(); it++) {
308 m->mothurOut(it->first + ", ");
310 m->mothurOutEndLine();
315 catch(exception& e) {
316 m->errorOut(e, "CommandFactory", "isValidCommand");
321 /***********************************************************************/
322 void CommandFactory::printCommands(ostream& out) {
324 out << "Valid commands are ";
325 for (it = commands.begin(); it != commands.end(); it++) {
326 out << it->first << ", ";
330 catch(exception& e) {
331 m->errorOut(e, "CommandFactory", "printCommands");
335 /***********************************************************************/