logGammaMinusLogGammaSum(ared, b)); } else { double prod1 = 1.0; return FastMath.log(prod1) + Gamma.logGamma(ared) + logGammaMinusLogGammaSum(ared, b); } else { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); if (b >= 10.0) { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); } else {
logGammaMinusLogGammaSum(ared, b)); } else { double prod1 = 1.0; return Math.log(prod1) + Gamma.logGamma(ared) + logGammaMinusLogGammaSum(ared, b); } else { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); if (b >= 10.0) { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); } else {
logGammaMinusLogGammaSum(ared, b)); } else { double prod1 = 1.0; return FastMath.log(prod1) + Gamma.logGamma(ared) + logGammaMinusLogGammaSum(ared, b); } else { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); if (b >= 10.0) { return Gamma.logGamma(a) + logGammaMinusLogGammaSum(a, b); } else {