[BACK]Return to ui_div.c CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / gmp / mpf

Diff for /OpenXM_contrib/gmp/mpf/Attic/ui_div.c between version 1.1.1.2 and 1.1.1.3

version 1.1.1.2, 2000/09/09 14:13:15 version 1.1.1.3, 2003/08/25 16:06:35
Line 1 
Line 1 
 /* mpf_ui_div -- Divide an unsigned integer with a float.  /* mpf_ui_div -- Divide an unsigned integer with a float.
   
 Copyright (C) 1993, 1994, 1995, 1996, 2000 Free Software Foundation, Inc.  Copyright 1993, 1994, 1995, 1996, 2000, 2001, 2002 Free Software Foundation,
   Inc.
   
 This file is part of the GNU MP Library.  This file is part of the GNU MP Library.
   
Line 24  MA 02111-1307, USA. */
Line 25  MA 02111-1307, USA. */
 #include "longlong.h"  #include "longlong.h"
   
 void  void
 #if __STDC__  
 mpf_ui_div (mpf_ptr r, unsigned long int u, mpf_srcptr v)  mpf_ui_div (mpf_ptr r, unsigned long int u, mpf_srcptr v)
 #else  
 mpf_ui_div (r, u, v)  
      mpf_ptr r;  
      unsigned long int u;  
      mpf_srcptr v;  
 #endif  
 {  {
   mp_srcptr vp;    mp_srcptr vp;
   mp_ptr rp, tp;    mp_ptr rp, tp;
Line 39  mpf_ui_div (r, u, v)
Line 33  mpf_ui_div (r, u, v)
   mp_size_t rsize, tsize;    mp_size_t rsize, tsize;
   mp_size_t sign_quotient;    mp_size_t sign_quotient;
   mp_size_t prec;    mp_size_t prec;
   unsigned normalization_steps;  
   mp_limb_t q_limb;    mp_limb_t q_limb;
   mp_exp_t rexp;    mp_exp_t rexp;
   TMP_DECL (marker);    TMP_DECL (marker);
Line 75  mpf_ui_div (r, u, v)
Line 68  mpf_ui_div (r, u, v)
   tp = (mp_ptr) TMP_ALLOC ((tsize + 1) * BYTES_PER_MP_LIMB);    tp = (mp_ptr) TMP_ALLOC ((tsize + 1) * BYTES_PER_MP_LIMB);
   MPN_ZERO (tp, tsize);    MPN_ZERO (tp, tsize);
   
   count_leading_zeros (normalization_steps, vp[vsize - 1]);  
   
   /* Normalize the divisor and the dividend.  */    /* Normalize the divisor and the dividend.  */
   if (normalization_steps != 0)    if ((vp[vsize - 1] & GMP_NUMB_HIGHBIT) == 0)
     {      {
       mp_ptr tmp;        mp_ptr tmp;
       mp_limb_t dividend_high, dividend_low;        mp_limb_t dividend_high, dividend_low;
         unsigned normalization_steps;
   
         count_leading_zeros (normalization_steps, vp[vsize - 1]);
         normalization_steps -= GMP_NAIL_BITS;
   
       /* Shift up the divisor setting the most significant bit of        /* Shift up the divisor setting the most significant bit of
          the most significant limb.  Use temporary storage not to clobber           the most significant limb.  Use temporary storage not to clobber
          the original contents of the divisor.  */           the original contents of the divisor.  */
Line 92  mpf_ui_div (r, u, v)
Line 87  mpf_ui_div (r, u, v)
   
       /* Shift up the dividend, possibly introducing a new most        /* Shift up the dividend, possibly introducing a new most
          significant word.  */           significant word.  */
       dividend_high = (mp_limb_t) u >> (BITS_PER_MP_LIMB - normalization_steps);        dividend_high = (mp_limb_t) u >> (GMP_NUMB_BITS - normalization_steps);
       dividend_low = (mp_limb_t) u << normalization_steps;        dividend_low = ((mp_limb_t) u << normalization_steps) & GMP_NUMB_MASK;
   
       tp[tsize - 1] = dividend_low;        tp[tsize - 1] = dividend_low;
       if (dividend_high != 0)        if (dividend_high != 0)
         {          {
           tp[tsize] = dividend_high;  #if GMP_NAIL_BITS != 0
           tsize++;            if (dividend_high > vp[vsize - 1])
           rexp++;              {
                 tp[tsize - 2] = dividend_low;
                 tp[tsize - 1] = dividend_high & GMP_NUMB_MASK;
                 tp[tsize] = dividend_high >> GMP_NUMB_BITS;
                 tsize++;
                 rexp += 2;
               }
             else
   #endif
               {
                 tp[tsize] = dividend_high;
                 tsize++;
                 rexp++;
               }
         }          }
     }      }
   else    else
Line 114  mpf_ui_div (r, u, v)
Line 122  mpf_ui_div (r, u, v)
           vp = (mp_srcptr) tmp;            vp = (mp_srcptr) tmp;
         }          }
   
       tp[tsize - 1] = u;        tp[tsize - 1] = u & GMP_NUMB_MASK;
   #if GMP_NAIL_BITS != 0
         if (u > GMP_NUMB_MAX)
           {
             tp[tsize] = u >> GMP_NUMB_BITS;
             tsize++;
             rexp++;
           }
   #endif
     }      }
   
   q_limb = mpn_divmod (rp, tp, tsize, vp, vsize);    q_limb = mpn_divmod (rp, tp, tsize, vp, vsize);
Legend:
Removed from v.1.1.1.2  
changed lines
  Added in v.1.1.1.3
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>