]> git.donarmstrong.com Git - mothur.git/blobdiff - screenseqscommand.h
removed read.dist, read.otu, read.tree and globaldata. added current to defaults...
[mothur.git] / screenseqscommand.h
index 425b473d4e1e7204ff460a00220bacc023dd0333..b8ab242e187b2f1b0fccc8b087e7c17f2cc0c903 100644 (file)
  */
 #include "mothur.h"
 #include "command.hpp"
-#include "globaldata.hpp"
-#include "readfasta.h"
-#include "readnexus.h"
-#include "readclustal.h"
-#include "readseqsphylip.h"
-
 
 class ScreenSeqsCommand : public Command {
        
 public:
+       ScreenSeqsCommand(string);
        ScreenSeqsCommand();
-       ~ScreenSeqsCommand();
-       int execute();
+       ~ScreenSeqsCommand() {}
+       
+       vector<string> setParameters();
+       string getCommandName()                 { return "screen.seqs";                         }
+       string getCommandCategory()             { return "Sequence Processing";         }
+       string getHelpString(); 
+       
+       int execute(); 
+       void help() { m->mothurOut(getHelpString()); }  
+       
+       
 private:
-       void screenNameGroupFile(set<string>);
-       void screenGroupFile(set<string>);
 
-       GlobalData* globaldata; 
+       struct linePair {
+               unsigned long int start;
+               unsigned long int end;
+               linePair(unsigned long int i, unsigned long int j) : start(i), end(j) {}
+       };
+
+       vector<int> processIDS;   //processid
+       vector<linePair*> lines;
+
+       int screenNameGroupFile(set<string>);
+       int screenGroupFile(set<string>);
+       int screenAlignReport(set<string>);
+       
+       int driver(linePair*, string, string, string, set<string>&);
+       int createProcesses(string, string, string, set<string>&);
+       
+       #ifdef USE_MPI
+       int driverMPI(int, int, MPI_File&, MPI_File&, MPI_File&, vector<unsigned long int>&, set<string>&);
+       #endif
+
+       bool abort;
+       string fastafile, namefile, groupfile, alignreport, outputDir;
+       int startPos, endPos, maxAmbig, maxHomoP, minLength, maxLength, processors, criteria;
+       vector<string> outputNames;
+       vector<string> optimize;
+       map<string, int> nameMap;
+       int readNames();
+       
+       int getSummary(vector<unsigned long int>&);
+       int createProcessesCreateSummary(vector<int>&, vector<int>&, vector<int>&, vector<int>&, vector<int>&, string);
+       int driverCreateSummary(vector<int>&, vector<int>&, vector<int>&, vector<int>&, vector<int>&, string, linePair*);       
 };
 
 #endif