| private static GammaRationalFunc ( double v0 ) : double | ||
| v0 | double | |
| return | double | |
private static double GammaRationalFunc(double v0) {
double numer = 0.0;
double denom = 0.0;
if (v0 < 1e15) {
numer = EvalPolynomial(v0, GammaNumerCoeffs);
denom = EvalPolynomial(v0, GammaDenomCoeffs);
} else {
double vRecip = 1.0 / v0;
numer = EvalPolynomial(vRecip, GammaNumerCoeffs, true);
denom = EvalPolynomial(vRecip, GammaDenomCoeffs, true);
}
return numer / denom;
}