diff options
| author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2002-02-06 02:57:55 +0000 |
|---|---|---|
| committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2002-02-06 02:57:55 +0000 |
| commit | 2d9c4ac4e03b1a5de88271453cae9b067a87dc44 (patch) | |
| tree | 2438f12c2659b310aeaba7f2087254746e619b0d /pow.c | |
| parent | 0a6ebfd2487ea512bebedfb443fa21215317fa6d (diff) | |
| download | mpfr-2d9c4ac4e03b1a5de88271453cae9b067a87dc44.tar.gz | |
Some optimizations.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@1682 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'pow.c')
| -rw-r--r-- | pow.c | 44 |
1 files changed, 23 insertions, 21 deletions
@@ -40,13 +40,9 @@ mpfr_pow (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mp_rnd_t rnd_mode) if (MPFR_IS_INF(y)) { - mpfr_t px; + mpfr_t one; int cmp; - /* This is inefficient and should be optimized later. */ - mpfr_init2(px, MPFR_PREC(x)); - mpfr_abs(px, x, GMP_RNDN); - if (MPFR_SIGN(y) > 0) { if (MPFR_IS_INF(x)) @@ -57,31 +53,34 @@ mpfr_pow (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mp_rnd_t rnd_mode) else MPFR_SET_ZERO(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } - cmp = mpfr_cmp_ui(px, 1); + MPFR_CLEAR_FLAGS(z); + if (MPFR_IS_ZERO(x)) + { + MPFR_SET_ZERO(z); + MPFR_SET_POS(z); + MPFR_RET(0); + } + mpfr_init2(one, BITS_PER_MP_LIMB); + mpfr_set_ui(one, 1, GMP_RNDN); + cmp = mpfr_cmp_abs(x, one); + mpfr_clear(one); if (cmp > 0) { - MPFR_CLEAR_FLAGS(z); MPFR_SET_INF(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } else if (cmp < 0) { - MPFR_CLEAR_FLAGS(z); MPFR_SET_ZERO(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } else { - MPFR_CLEAR_FLAGS(z); MPFR_SET_NAN(z); - mpfr_clear(px); MPFR_RET_NAN; } } @@ -95,31 +94,34 @@ mpfr_pow (mpfr_ptr z, mpfr_srcptr x, mpfr_srcptr y, mp_rnd_t rnd_mode) else MPFR_SET_INF(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } - cmp = mpfr_cmp_ui(px, 1); + if (MPFR_IS_ZERO(x)) + { + MPFR_SET_INF(z); + MPFR_SET_POS(z); + MPFR_RET(0); + } + mpfr_init2(one, BITS_PER_MP_LIMB); + mpfr_set_ui(one, 1, GMP_RNDN); + cmp = mpfr_cmp_abs(x, one); + mpfr_clear(one); + MPFR_CLEAR_FLAGS(z); if (cmp > 0) { - MPFR_CLEAR_FLAGS(z); MPFR_SET_ZERO(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } else if (cmp < 0) { - MPFR_CLEAR_FLAGS(z); MPFR_SET_INF(z); MPFR_SET_POS(z); - mpfr_clear(px); MPFR_RET(0); } else { - MPFR_CLEAR_FLAGS(z); MPFR_SET_NAN(z); - mpfr_clear(px); MPFR_RET_NAN; } } |