From 6e879c6f3361346ede165d0733a26e9145cccf69 Mon Sep 17 00:00:00 2001
From: vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>
Date: Fri, 1 Feb 2002 21:35:34 +0000
Subject: Misc bug fixes and code clean-up.

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@1678 280ebfd0-de03-0410-8827-d642c229c3f4
---
 pow_si.c | 78 ++++++++++++++++++++--------------------------------------------
 1 file changed, 24 insertions(+), 54 deletions(-)

(limited to 'pow_si.c')

diff --git a/pow_si.c b/pow_si.c
index 8646e1c13..0380b2fa7 100644
--- a/pow_si.c
+++ b/pow_si.c
@@ -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 */
-- 
cgit v1.2.1