Index: gcc/gmp/tune/modlinv.c |
diff --git a/gcc/gmp/tune/modlinv.c b/gcc/gmp/tune/modlinv.c |
deleted file mode 100644 |
index 3ed0c8a03f3f6153957f61acd3e61f845971dfaf..0000000000000000000000000000000000000000 |
--- a/gcc/gmp/tune/modlinv.c |
+++ /dev/null |
@@ -1,167 +0,0 @@ |
-/* Alternate implementations of binvert_limb to compare speeds. */ |
- |
-/* |
-Copyright 2000, 2002 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 3 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. If not, see http://www.gnu.org/licenses/. */ |
- |
-#include <stdio.h> |
-#include "gmp.h" |
-#include "gmp-impl.h" |
-#include "longlong.h" |
-#include "speed.h" |
- |
- |
-/* Like the standard version in gmp-impl.h, but with the expressions using a |
- "1-" form. This has the same number of steps, but "1-" is on the |
- dependent chain, whereas the "2*" in the standard version isn't. |
- Depending on the CPU this should be the same or a touch slower. */ |
- |
-#if BITS_PER_MP_LIMB <= 32 |
-#define binvert_limb_mul1(inv,n) \ |
- do { \ |
- mp_limb_t __n = (n); \ |
- mp_limb_t __inv; \ |
- ASSERT ((__n & 1) == 1); \ |
- __inv = binvert_limb_table[(__n&0xFF)/2]; /* 8 */ \ |
- __inv = (1 - __n * __inv) * __inv + __inv; /* 16 */ \ |
- __inv = (1 - __n * __inv) * __inv + __inv; /* 32 */ \ |
- ASSERT (__inv * __n == 1); \ |
- (inv) = __inv; \ |
- } while (0) |
-#endif |
- |
-#if BITS_PER_MP_LIMB > 32 && BITS_PER_MP_LIMB <= 64 |
-#define binvert_limb_mul1(inv,n) \ |
- do { \ |
- mp_limb_t __n = (n); \ |
- mp_limb_t __inv; \ |
- ASSERT ((__n & 1) == 1); \ |
- __inv = binvert_limb_table[(__n&0xFF)/2]; /* 8 */ \ |
- __inv = (1 - __n * __inv) * __inv + __inv; /* 16 */ \ |
- __inv = (1 - __n * __inv) * __inv + __inv; /* 32 */ \ |
- __inv = (1 - __n * __inv) * __inv + __inv; /* 64 */ \ |
- ASSERT (__inv * __n == 1); \ |
- (inv) = __inv; \ |
- } while (0) |
-#endif |
- |
- |
-/* The loop based version used in GMP 3.0 and earlier. Usually slower than |
- multiplying, due to the number of steps that must be performed. Much |
- slower when the processor has a good multiply. */ |
- |
-#define binvert_limb_loop(inv,n) \ |
- do { \ |
- mp_limb_t __v = (n); \ |
- mp_limb_t __v_orig = __v; \ |
- mp_limb_t __make_zero = 1; \ |
- mp_limb_t __two_i = 1; \ |
- mp_limb_t __v_inv = 0; \ |
- \ |
- ASSERT ((__v & 1) == 1); \ |
- \ |
- do \ |
- { \ |
- while ((__two_i & __make_zero) == 0) \ |
- __two_i <<= 1, __v <<= 1; \ |
- __v_inv += __two_i; \ |
- __make_zero -= __v; \ |
- } \ |
- while (__make_zero); \ |
- \ |
- ASSERT (__v_orig * __v_inv == 1); \ |
- (inv) = __v_inv; \ |
- } while (0) |
- |
- |
-/* Another loop based version with conditionals, but doing a fixed number of |
- steps. */ |
- |
-#define binvert_limb_cond(inv,n) \ |
- do { \ |
- mp_limb_t __n = (n); \ |
- mp_limb_t __rem = (1 - __n) >> 1; \ |
- mp_limb_t __inv = GMP_LIMB_HIGHBIT; \ |
- int __count; \ |
- \ |
- ASSERT ((__n & 1) == 1); \ |
- \ |
- __count = BITS_PER_MP_LIMB-1; \ |
- do \ |
- { \ |
- __inv >>= 1; \ |
- if (__rem & 1) \ |
- { \ |
- __inv |= GMP_LIMB_HIGHBIT; \ |
- __rem -= __n; \ |
- } \ |
- __rem >>= 1; \ |
- } \ |
- while (-- __count); \ |
- \ |
- ASSERT (__inv * __n == 1); \ |
- (inv) = __inv; \ |
- } while (0) |
- |
- |
-/* Another loop based bitwise version, but purely arithmetic, no |
- conditionals. */ |
- |
-#define binvert_limb_arith(inv,n) \ |
- do { \ |
- mp_limb_t __n = (n); \ |
- mp_limb_t __rem = (1 - __n) >> 1; \ |
- mp_limb_t __inv = GMP_LIMB_HIGHBIT; \ |
- mp_limb_t __lowbit; \ |
- int __count; \ |
- \ |
- ASSERT ((__n & 1) == 1); \ |
- \ |
- __count = BITS_PER_MP_LIMB-1; \ |
- do \ |
- { \ |
- __lowbit = __rem & 1; \ |
- __inv = (__inv >> 1) | (__lowbit << (BITS_PER_MP_LIMB-1)); \ |
- __rem = (__rem - (__n & -__lowbit)) >> 1; \ |
- } \ |
- while (-- __count); \ |
- \ |
- ASSERT (__inv * __n == 1); \ |
- (inv) = __inv; \ |
- } while (0) |
- |
- |
-double |
-speed_binvert_limb_mul1 (struct speed_params *s) |
-{ |
- SPEED_ROUTINE_MODLIMB_INVERT (binvert_limb_mul1); |
-} |
-double |
-speed_binvert_limb_loop (struct speed_params *s) |
-{ |
- SPEED_ROUTINE_MODLIMB_INVERT (binvert_limb_loop); |
-} |
-double |
-speed_binvert_limb_cond (struct speed_params *s) |
-{ |
- SPEED_ROUTINE_MODLIMB_INVERT (binvert_limb_cond); |
-} |
-double |
-speed_binvert_limb_arith (struct speed_params *s) |
-{ |
- SPEED_ROUTINE_MODLIMB_INVERT (binvert_limb_arith); |
-} |