[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.1 and 1.1.1.3

version 1.1.1.1, 2000/01/10 15:35:22 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 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.
   
 The GNU MP Library is free software; you can redistribute it and/or modify  The GNU MP Library is free software; you can redistribute it and/or modify
 it under the terms of the GNU Library General Public License as published by  it under the terms of the GNU Lesser General Public License as published by
 the Free Software Foundation; either version 2 of the License, or (at your  the Free Software Foundation; either version 2.1 of the License, or (at your
 option) any later version.  option) any later version.
   
 The GNU MP Library is distributed in the hope that it will be useful, but  The GNU MP Library is distributed in the hope that it will be useful, but
 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY  WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
 or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Library General Public  or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
 License for more details.  License for more details.
   
 You should have received a copy of the GNU Library General Public License  You should have received a copy of the GNU Lesser General Public License
 along with the GNU MP Library; see the file COPYING.LIB.  If not, write to  along with the GNU MP Library; see the file COPYING.LIB.  If not, write to
 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,  the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 MA 02111-1307, USA. */  MA 02111-1307, USA. */
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 50  mpf_ui_div (r, u, v)
Line 43  mpf_ui_div (r, u, v)
   prec = r->_mp_prec;    prec = r->_mp_prec;
   
   if (vsize == 0)    if (vsize == 0)
     vsize = 1 / vsize;          /* divide by zero as directed */      DIVIDE_BY_ZERO;
   
   if (u == 0)    if (u == 0)
     {      {
       r->_mp_size = 0;        r->_mp_size = 0;
Line 74  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 91  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 113  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.1  
changed lines
  Added in v.1.1.1.3
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>