change to phylo.diversity to make output = total branch length seen so far for each...

[mothur.git] / phylodiversity.cpp

diff --git a/phylodiversity.cpp b/phylodiversity.cpp
index 96d49d489e192cc3f9b1cf52f4091dc32cab9f20..37884c483619887bcc478b6fdf7e77be08424cd4 100644 (file)
--- a/phylodiversity.cpp
+++ b/phylodiversity.cpp
@@ -66,10 +66,11 @@ EstOutput PhyloDiversity::getValues(Tree* t, vector<int> treeNodes) {
                
        
                for (int i=0; i<globaldata->Groups.size(); i++) {   
-                       if (groupTotals[globaldata->Groups[i]] != 0.0) { //avoid divide by zero error
-                               float percent = DScore[globaldata->Groups[i]] / groupTotals[globaldata->Groups[i]];
-                               data.push_back(percent);  
-                       }else {  data.push_back(0.0);  }
+                       //if (groupTotals[globaldata->Groups[i]] != 0.0) { //avoid divide by zero error
+                               //float percent = DScore[globaldata->Groups[i]] / groupTotals[globaldata->Groups[i]];
+                       float percent = DScore[globaldata->Groups[i]]; 
+                       data.push_back(percent);  
+                       //}else {  data.push_back(0.0);  }
                }
                
                return data;