]> git.donarmstrong.com Git - mothur.git/blobdiff - summarysharedcommand.h
fixes while testing 1.33.0
[mothur.git] / summarysharedcommand.h
index f486d520216bd8ec4858d50cd7f2d92b7fd23537..7a2a1909368fc9c71edbcbc90add7d3a966ac95d 100644 (file)
@@ -56,6 +56,8 @@
 #include "memchord.h"
 #include "memeuclidean.h"
 #include "mempearson.h"
+#include "sharedjsd.h"
+#include "sharedrjsd.h"
 
 class SummarySharedCommand : public Command {
 
@@ -112,6 +114,7 @@ struct summarySharedData {
        unsigned long long end;
        MothurOut* m;
        string sumFile;
+    int count;
        
        summarySharedData(){}
        summarySharedData(string sf, MothurOut* mout, unsigned long long st, unsigned long long en, vector<string> est, vector<SharedRAbundVector*> lu) {
@@ -121,6 +124,7 @@ struct summarySharedData {
                end = en;
         Estimators = est;
         thisLookup = lu;
+        count=0;
        }
 };
 /**************************************************************************************************/
@@ -214,6 +218,10 @@ static DWORD WINAPI MySummarySharedThreadFunction(LPVOID lpParam){
                     sumCalculators.push_back(new MemEuclidean());
                 }else if (pDataArray->Estimators[i] == "mempearson") { 
                     sumCalculators.push_back(new MemPearson());
+                }else if (pDataArray->Estimators[i] == "jsd") {
+                    sumCalculators.push_back(new JSD());
+                }else if (pDataArray->Estimators[i] == "rjsd") {
+                    sumCalculators.push_back(new RJSD());
                 }
             }
         }
@@ -225,7 +233,7 @@ static DWORD WINAPI MySummarySharedThreadFunction(LPVOID lpParam){
                
                vector<SharedRAbundVector*> subset;
                for (int k = pDataArray->start; k < pDataArray->end; k++) { // pass cdd each set of groups to compare
-            
+            pDataArray->count++;
                        for (int l = 0; l < k; l++) {
                                
                                outputFileHandle << pDataArray->thisLookup[0]->getLabel() << '\t';