diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2000-10-26 11:59:59 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2000-10-26 11:59:59 +0000 |
commit | da4b19c01ec0e27d42071d294e29fc56e1a1b16a (patch) | |
tree | 393496f9c6cc54d9bcde0fc0d77691fdcfdf255a /log.c | |
parent | 33336e9a85f1ae242191e1b2286cded95f11362b (diff) | |
download | mpfr-da4b19c01ec0e27d42071d294e29fc56e1a1b16a.tar.gz |
protected all macros: xxx -> MPFR_xxx
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@786 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'log.c')
-rw-r--r-- | log.c | 14 |
1 files changed, 7 insertions, 7 deletions
@@ -63,16 +63,16 @@ mpfr_log(r, a, rnd_mode) TMP_DECL(marker); /* If a is NaN or a is negative or null, the result is NaN */ - if (FLAG_NAN(a) || (NOTZERO(a)==0) || (MPFR_SIGN(a)<0)) - { SET_NAN(r); return 1; } + if (MPFR_IS_NAN(a) || (MPFR_NOTZERO(a)==0) || (MPFR_SIGN(a)<0)) + { MPFR_SET_NAN(r); return 1; } /* If a is 1, the result is 0 */ if (mpfr_cmp_ui_2exp(a,1,0)==0){ - SET_ZERO(r); + MPFR_SET_ZERO(r); return 0; /* only case where the result is exact */ } - q=PREC(r); + q=MPFR_PREC(r); ref=mpfr_get_d(a)-1.0; if (ref<0) @@ -92,7 +92,7 @@ mpfr_log(r, a, rnd_mode) printf("p=%d\n", p); #endif /* Calculus of m (depends on p) */ - m=(int) ceil(((double) p)/2.0) -EXP(a)+1; + m=(int) ceil(((double) p)/2.0) -MPFR_EXP(a)+1; /* All the mpfr_t needed have a precision of p */ TMP_MARK(marker); @@ -116,9 +116,9 @@ mpfr_log(r, a, rnd_mode) mpfr_div(tmp2,cst,tmp1,GMP_RNDN); /* pi/2*AG(1,4/s), err<=5ulps */ mpfr_const_log2(cst,GMP_RNDN); /* compute log(2), err<=1ulp */ mpfr_mul(tmp1,cst,mm,GMP_RNDN); /* I compute m*log(2), err<=2ulps */ - cancel = EXP(tmp2); + cancel = MPFR_EXP(tmp2); mpfr_sub(cst,tmp2,tmp1,GMP_RNDN); /* log(a), err<=7ulps+cancel */ - cancel -= EXP(cst); + cancel -= MPFR_EXP(cst); #ifdef DEBUG printf("cancelled bits=%d\n", cancel); printf("approx="); mpfr_print_raw(cst); putchar('\n'); |