X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=chimeracheckrdp.h;h=cc23c2e9e8b0b794339a0dfa8fa91bae50f28c00;hb=70491a12902e89b85cfa6b44a7b7fbe066ee2ac1;hp=fe4620614b8ae2d78d985599c594185313b40c4c;hpb=63e089e0b3aad1741bab60119ed7ccc784dce347;p=mothur.git diff --git a/chimeracheckrdp.h b/chimeracheckrdp.h index fe46206..cc23c2e 100644 --- a/chimeracheckrdp.h +++ b/chimeracheckrdp.h @@ -25,45 +25,35 @@ class ChimeraCheckRDP : public Chimera { public: - ChimeraCheckRDP(string, string); + ChimeraCheckRDP(string, string, string, bool, int, int, string); //fasta, template, name, svg, increment, ksize, outputDir ~ChimeraCheckRDP(); - void getChimeras(); - void print(ostream&); - - void setCons(string){}; - void setQuantiles(string q) {}; + int getChimeras(Sequence*); + Sequence print(ostream&, ostream&); + #ifdef USE_MPI + Sequence print(MPI_File&, MPI_File&); + #endif private: - vector lines; - vector querySeqs; + Sequence* querySeq; AlignmentDB* templateDB; Kmer* kmer; - - vector< vector > IS; //IS[0] is the vector of IS values for each window for querySeqs[0] - float chimeraCutoff; - - - //map > >:: iterator it; - //map::iterator it2; - - vector closest; //closest[0] is the closest overall seq to querySeqs[0]. - - string fastafile, templateFile; + Sequence closest; //closest is the closest overall seq to query + + vector IS; //IS is the vector of IS values for each window for query + string fastafile; map names; + string name; + bool svg; + int kmerSize, increment; - vector findIS(int); + vector findIS(); int calcKmers(map, map); - void getCutoff(); - void makeSVGpic(vector, int); + void makeSVGpic(vector); void readName(string); - - void createProcessesIS(); - }; - /***********************************************************/ #endif