X-Git-Url: http://git.megacz.com/?a=blobdiff_plain;f=rts%2Fgmp%2Fmpn%2Fpa64%2Fudiv_qrnnd.c;fp=rts%2Fgmp%2Fmpn%2Fpa64%2Fudiv_qrnnd.c;h=0000000000000000000000000000000000000000;hb=fdf1cd0399158308769fcb2ab7e46e215a68d865;hp=1c9fe084db7c4dba8eff28ace1c411bdeb1a57d9;hpb=e552cfc427d2734b9a9629f2ab1d22f493e775f6;p=ghc-hetmet.git

diff --git a/rts/gmp/mpn/pa64/udiv_qrnnd.c b/rts/gmp/mpn/pa64/udiv_qrnnd.c
deleted file mode 100644
index 1c9fe08..0000000
--- a/rts/gmp/mpn/pa64/udiv_qrnnd.c
+++ /dev/null
@@ -1,111 +0,0 @@
-/*
-Copyright (C) 1999, 2000 Free Software Foundation, Inc.
-
-This file is part of the GNU MP Library.
-
-The GNU MP Library is free software; you can redistribute it and/or modify
-it under the terms of the GNU Lesser General Public License as published by
-the Free Software Foundation; either version 2.1 of the License, or (at your
-option) any later version.
-
-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
-or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
-License for more details.
-
-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
-the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
-MA 02111-1307, USA.
-*/
-
-#include "gmp.h"
-#include "gmp-impl.h"
-#include "longlong.h"
-
-#define TWO64 18446744073709551616.0
-
-mp_limb_t
-#if __STDC__
-__MPN(udiv_qrnnd) (mp_limb_t n1, mp_limb_t n0, mp_limb_t d, mp_limb_t *r)
-#else
-__MPN(udiv_qrnnd) (n1, n0, d, r)
-     mp_limb_t n1;
-     mp_limb_t n0;
-     mp_limb_t d;
-     mp_limb_t *r;
-#endif
-{
-  mp_limb_t q1, q2, q;
-  mp_limb_t p1, p0;
-  double di, dq;
-
-  di = 1.0 / d;
-
-  /* Generate upper 53 bits of quotient.  Be careful here; the `double'
-     quotient may be rounded to 2^64 which we cannot safely convert back
-     to a 64-bit integer.  */
-  dq = (TWO64 * (double) n1 + (double) n0) * di;
-  if (dq >= TWO64)
-    q1 = 0xfffffffffffff800LL;
-  else
-    q1 = (mp_limb_t) dq;
-
-  /* Multiply back in order to compare the product to the dividend.  */
-  umul_ppmm (p1, p0, q1, d);
-
-  /* Was the 53-bit quotient greater that our sought quotient?  Test the
-     sign of the partial remainder to find out.  */
-  if (n1 < p1 || (n1 == p1 && n0 < p0))
-    {
-      /* 53-bit quotient too large.  Partial remainder is negative.
-	 Compute the absolute value of the remainder in n1,,n0.  */
-      n1 = p1 - (n1 + (p0 < n0));
-      n0 = p0 - n0;
-
-      /* Now use the partial remainder as new dividend to compute more bits of
-	 quotient.  This is an adjustment for the one we got previously.  */
-      q2 = (mp_limb_t) ((TWO64 * (double) n1 + (double) n0) * di);
-      umul_ppmm (p1, p0, q2, d);
-
-      q = q1 - q2;
-      if (n1 < p1 || (n1 == p1 && n0 <= p0))
-	{
-	  n0 = p0 - n0;
-	}
-      else
-	{
-	  n0 = p0 - n0;
-	  n0 += d;
-	  q--;
-	}
-    }
-  else
-    {
-      n1 = n1 - (p1 + (n0 < p0));
-      n0 = n0 - p0;
-
-      q2 = (mp_limb_t) ((TWO64 * (double) n1 + (double) n0) * di);
-      umul_ppmm (p1, p0, q2, d);
-
-      q = q1 + q2;
-      if (n1 < p1 || (n1 == p1 && n0 < p0))
-	{
-	  n0 = n0 - p0;
-	  n0 += d;
-	  q--;
-	}
-      else
-	{
-	  n0 = n0 - p0;
-	  if (n0 >= d)
-	    {
-	      n0 -= d;
-	      q++;
-	    }
-	}
-    }
-
-  *r = n0;
-  return q;
-}