diff options
Diffstat (limited to 'src/digamma.c')
| -rw-r--r-- | src/digamma.c | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/digamma.c b/src/digamma.c index 001b47e3f..2613724fa 100644 --- a/src/digamma.c +++ b/src/digamma.c @@ -44,7 +44,7 @@ mpfr_digamma_approx (mpfr_ptr s, mpfr_srcptr x) mpfr_log (s, x, MPFR_RNDN); /* error <= 1/2 ulp */ mpfr_ui_div (t, 1, x, MPFR_RNDN); /* error <= 1/2 ulp */ - mpfr_div_2exp (t, t, 1, MPFR_RNDN); /* exact */ + mpfr_div_2ui (t, t, 1, MPFR_RNDN); /* exact */ mpfr_sub (s, s, t, MPFR_RNDN); /* error <= 1/2 + 1/2*2^(EXP(olds)-EXP(s)) + 1/2*2^(EXP(t)-EXP(s)). For x >= 2, log(x) >= 2*(1/(2x)), thus olds >= 2t, and olds - t >= olds/2, @@ -136,9 +136,9 @@ mpfr_digamma_reflection (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode) MPFR_ASSERTN(inex == 0); /* if x is half an integer, cot(Pi*x) = 0, thus Digamma(x) = Digamma(1-x) */ - mpfr_mul_2exp (u, u, 1, MPFR_RNDN); + mpfr_mul_2ui (u, u, 1, MPFR_RNDN); inex = mpfr_integer_p (u); - mpfr_div_2exp (u, u, 1, MPFR_RNDN); + mpfr_div_2ui (u, u, 1, MPFR_RNDN); if (inex) { inex = mpfr_digamma (y, u, rnd_mode); |