//store in temps to avoid multiple repetitive function calls
pi = shared[0]->getAbundance(j) / Atotal;
qi = shared[1]->getAbundance(j) / Btotal;
//store in temps to avoid multiple repetitive function calls
pi = shared[0]->getAbundance(j) / Atotal;
qi = shared[1]->getAbundance(j) / Btotal;
if (isnan(data[1]) || isinf(data[1])) { data[1] = 0; }
if (isnan(data[2]) || isinf(data[2])) { data[2] = 0; }
if (isnan(data[1]) || isinf(data[1])) { data[1] = 0; }
if (isnan(data[2]) || isinf(data[2])) { data[2] = 0; }