//store in temps to avoid multiple repetitive function calls
Atotal += (double)shared[0]->getAbundance(i);
Btotal += (double)shared[1]->getAbundance(i);
}
//calculate the theta denominator sums
//store in temps to avoid multiple repetitive function calls
Atotal += (double)shared[0]->getAbundance(i);
Btotal += (double)shared[1]->getAbundance(i);
}
//calculate the theta denominator sums