- //store in temps to avoid multiple repetitive function calls
- double tempA = shared[0]->getAbundance(i);
- double tempB = shared[1]->getAbundance(i);
+ totalA += countA[i];
+ totalB += countB[i];
+ }
+
+ for (int i = 0; i < shared[0]->getNumBins(); i++) {
+ double tempA = countA[i] / totalA;
+ double tempB = countB[i] / totalB;