//get the total values we need to calculate the theta denominator sums
for (int i = 0; i < shared[0]->size(); i++) {
//store in temps to avoid multiple repetitive function calls
- Atotal += (float)shared[0]->getAbundance(i);
- Btotal += (float)shared[1]->getAbundance(i);
+ Atotal += (double)shared[0]->getAbundance(i);
+ Btotal += (double)shared[1]->getAbundance(i);
}
//calculate the theta denominator sums
if (isnan(data[1]) || isinf(data[1])) { data[1] = 0; }
if (isnan(data[2]) || isinf(data[2])) { data[2] = 0; }
+ data[0] = 1.0 - data[0];
return data;
}
catch(exception& e) {