numerator = 0.0; denominator = 0.0; thetaN = 0.0; sumSharedA = 0.0; sumSharedB = 0.0; a = 0.0; b = 0.0; d = 0.0;
//get the total values we need to calculate the theta denominator sums
numerator = 0.0; denominator = 0.0; thetaN = 0.0; sumSharedA = 0.0; sumSharedB = 0.0; a = 0.0; b = 0.0; d = 0.0;
//get the total values we need to calculate the theta denominator sums
//store in temps to avoid multiple repetitive function calls
Atotal += shared[0]->getAbundance(i);
Btotal += shared[1]->getAbundance(i);
}
//calculate the theta denominator sums
//store in temps to avoid multiple repetitive function calls
Atotal += shared[0]->getAbundance(i);
Btotal += shared[1]->getAbundance(i);
}
//calculate the theta denominator sums
//store in temps to avoid multiple repetitive function calls
tempA = shared[0]->getAbundance(j);
tempB = shared[1]->getAbundance(j);
//store in temps to avoid multiple repetitive function calls
tempA = shared[0]->getAbundance(j);
tempB = shared[1]->getAbundance(j);