summaryrefslogtreecommitdiff
path: root/numpy/random/src/legacy/legacy-distributions.c
diff options
context:
space:
mode:
authormattip <matti.picus@gmail.com>2019-10-02 22:02:03 +0300
committermattip <matti.picus@gmail.com>2019-10-11 15:08:46 +0300
commitaeac7d56e91577d54ee48a6ae928709414f176e5 (patch)
tree88dcd45fb67b5938c857f12fb351c1885ceb73db /numpy/random/src/legacy/legacy-distributions.c
parentfbda4849525316223d984b27f9f48bdd4ad477bb (diff)
downloadnumpy-aeac7d56e91577d54ee48a6ae928709414f176e5.tar.gz
API: remove unused functions from distributions.h
Diffstat (limited to 'numpy/random/src/legacy/legacy-distributions.c')
-rw-r--r--numpy/random/src/legacy/legacy-distributions.c8
1 files changed, 4 insertions, 4 deletions
diff --git a/numpy/random/src/legacy/legacy-distributions.c b/numpy/random/src/legacy/legacy-distributions.c
index 404c3e301..fd067fe8d 100644
--- a/numpy/random/src/legacy/legacy-distributions.c
+++ b/numpy/random/src/legacy/legacy-distributions.c
@@ -294,8 +294,8 @@ static RAND_INT_TYPE random_hypergeometric_hrua(bitgen_t *bitgen_state,
d7 = sqrt((double)(popsize - m) * sample * d4 * d5 / (popsize - 1) + 0.5);
d8 = D1 * d7 + D2;
d9 = (RAND_INT_TYPE)floor((double)(m + 1) * (mingoodbad + 1) / (popsize + 2));
- d10 = (loggam(d9 + 1) + loggam(mingoodbad - d9 + 1) + loggam(m - d9 + 1) +
- loggam(maxgoodbad - m + d9 + 1));
+ d10 = (random_loggam(d9 + 1) + random_loggam(mingoodbad - d9 + 1) +
+ random_loggam(m - d9 + 1) + random_loggam(maxgoodbad - m + d9 + 1));
d11 = MIN(MIN(m, mingoodbad) + 1.0, floor(d6 + 16 * d7));
/* 16 for 16-decimal-digit precision in D1 and D2 */
@@ -309,8 +309,8 @@ static RAND_INT_TYPE random_hypergeometric_hrua(bitgen_t *bitgen_state,
continue;
Z = (RAND_INT_TYPE)floor(W);
- T = d10 - (loggam(Z + 1) + loggam(mingoodbad - Z + 1) + loggam(m - Z + 1) +
- loggam(maxgoodbad - m + Z + 1));
+ T = d10 - (random_loggam(Z + 1) + random_loggam(mingoodbad - Z + 1) +
+ random_loggam(m - Z + 1) + random_loggam(maxgoodbad - m + Z + 1));
/* fast acceptance: */
if ((X * (4.0 - X) - 3.0) <= T)