diff options
author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2002-02-01 21:35:34 +0000 |
---|---|---|
committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2002-02-01 21:35:34 +0000 |
commit | 6e879c6f3361346ede165d0733a26e9145cccf69 (patch) | |
tree | ba27f0eb5a934b1be7da19ed5937cb8f7b076608 /pow_si.c | |
parent | d2a8d03f413280358567c13baa9b2e1f6e3349bf (diff) | |
download | mpfr-6e879c6f3361346ede165d0733a26e9145cccf69.tar.gz |
Misc bug fixes and code clean-up.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@1678 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'pow_si.c')
-rw-r--r-- | pow_si.c | 78 |
1 files changed, 24 insertions, 54 deletions
@@ -1,6 +1,6 @@ /* mpfr_pow_si -- power function x^y with y an unsigned int -Copyright (C) 2001 Free Software Foundation, Inc. +Copyright (C) 2001-2002 Free Software Foundation, Inc. This file is part of the MPFR Library. @@ -34,77 +34,47 @@ MA 02111-1307, USA. */ int mpfr_pow_si (mpfr_ptr y, mpfr_srcptr x, long int n, mp_rnd_t rnd_mode) { - - if (n>0) - return mpfr_pow_ui(y,x,(unsigned long int)n,rnd_mode); + if (n > 0) + return mpfr_pow_ui(y, x, n, rnd_mode); else { - int inexact = 0; - n=-n; - - /* x is NaN*/ - if (MPFR_IS_NAN(x)) + if (MPFR_IS_NAN(x)) { - MPFR_SET_NAN(y); - return 1; + MPFR_SET_NAN(y); + MPFR_RET_NAN; } + MPFR_CLEAR_NAN(y); - /* n=0 */ - if(n==0) - return mpfr_set_ui(y,1,GMP_RNDN);; + if (n == 0) + return mpfr_set_ui(y, 1, GMP_RNDN); - /* case x is INF */ - if(MPFR_IS_INF(x)) + if (MPFR_IS_INF(x)) { - if(MPFR_SIGN(x)>0) /* +Inf */ - { - MPFR_SET_ZERO(y); - if(MPFR_SIGN(y) < 0) - MPFR_CHANGE_SIGN(y); - return 0; - } + MPFR_SET_ZERO(y); + if (MPFR_SIGN(x) > 0 || ((unsigned) n & 1) == 0) + MPFR_SET_POS(y); else - { - MPFR_SET_ZERO(y); /* -Inf */ - if(!(n%2)) /* n is odd */ - { - if(MPFR_SIGN(y) > 0) - MPFR_CHANGE_SIGN(y); - } - else /* n is not odd */ - { - if(MPFR_SIGN(y) < 0) - MPFR_CHANGE_SIGN(y); - } - return 0; - } + MPFR_SET_NEG(y); + MPFR_RET(0); } - /* case x=0 */ - if(mpfr_cmp_ui(x,0) == 0) + if (MPFR_IS_ZERO(x)) { - if(!(n%2)) /* n is odd */ - { - MPFR_SET_INF(y); - MPFR_SET_SAME_SIGN(y,x); - DIVIDE_BY_ZERO; /* Execption GMP*/ - return 0; - } - else /* n is not odd */ - { - MPFR_SET_INF(y); - if(MPFR_SIGN(y) < 0) - MPFR_CHANGE_SIGN(y); - DIVIDE_BY_ZERO; /* Execption GMP*/ - return 0; - } + MPFR_SET_INF(y); + if (MPFR_SIGN(x) > 0 || ((unsigned) n & 1) == 0) + MPFR_SET_POS(y); + else + MPFR_SET_NEG(y); + MPFR_RET(0); } MPFR_CLEAR_INF(y); + n = -n; + /* General case */ { /* Declaration of the intermediary variable */ |