| Index: gcc/gmp/gen-bases.c
|
| diff --git a/gcc/gmp/gen-bases.c b/gcc/gmp/gen-bases.c
|
| deleted file mode 100644
|
| index 31895e8b240939c1bacc634c4fad19be941e7aae..0000000000000000000000000000000000000000
|
| --- a/gcc/gmp/gen-bases.c
|
| +++ /dev/null
|
| @@ -1,180 +0,0 @@
|
| -/* Generate mp_bases data.
|
| -
|
| -Copyright 1991, 1993, 1994, 1996, 2000, 2002, 2004 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 <math.h>
|
| -
|
| -#include "dumbmp.c"
|
| -
|
| -
|
| -int chars_per_limb;
|
| -double chars_per_bit_exactly;
|
| -mpz_t big_base;
|
| -int normalization_steps;
|
| -mpz_t big_base_inverted;
|
| -
|
| -mpz_t t;
|
| -
|
| -#define POW2_P(n) (((n) & ((n) - 1)) == 0)
|
| -
|
| -unsigned int
|
| -ulog2 (unsigned int x)
|
| -{
|
| - unsigned int i;
|
| - for (i = 0; x != 0; i++)
|
| - x >>= 1;
|
| - return i;
|
| -}
|
| -
|
| -void
|
| -generate (int limb_bits, int nail_bits, int base)
|
| -{
|
| - int numb_bits = limb_bits - nail_bits;
|
| -
|
| - mpz_set_ui (t, 1L);
|
| - mpz_mul_2exp (t, t, numb_bits);
|
| - mpz_set_ui (big_base, 1L);
|
| - chars_per_limb = 0;
|
| - for (;;)
|
| - {
|
| - mpz_mul_ui (big_base, big_base, (long) base);
|
| - if (mpz_cmp (big_base, t) > 0)
|
| - break;
|
| - chars_per_limb++;
|
| - }
|
| -
|
| - chars_per_bit_exactly = 0.69314718055994530942 / log ((double) base);
|
| -
|
| - mpz_ui_pow_ui (big_base, (long) base, (long) chars_per_limb);
|
| -
|
| - normalization_steps = limb_bits - mpz_sizeinbase (big_base, 2);
|
| -
|
| - mpz_set_ui (t, 1L);
|
| - mpz_mul_2exp (t, t, 2*limb_bits - normalization_steps);
|
| - mpz_tdiv_q (big_base_inverted, t, big_base);
|
| - mpz_set_ui (t, 1L);
|
| - mpz_mul_2exp (t, t, limb_bits);
|
| - mpz_sub (big_base_inverted, big_base_inverted, t);
|
| -}
|
| -
|
| -void
|
| -header (int limb_bits, int nail_bits)
|
| -{
|
| - int numb_bits = limb_bits - nail_bits;
|
| -
|
| - generate (limb_bits, nail_bits, 10);
|
| -
|
| - printf ("/* This file generated by gen-bases.c - DO NOT EDIT. */\n");
|
| - printf ("\n");
|
| - printf ("#if GMP_NUMB_BITS != %d\n", numb_bits);
|
| - printf ("Error, error, this data is for %d bits\n", numb_bits);
|
| - printf ("#endif\n");
|
| - printf ("\n");
|
| - printf ("/* mp_bases[10] data, as literal values */\n");
|
| - printf ("#define MP_BASES_CHARS_PER_LIMB_10 %d\n", chars_per_limb);
|
| - printf ("#define MP_BASES_BIG_BASE_10 CNST_LIMB(0x");
|
| - mpz_out_str (stdout, 16, big_base);
|
| - printf (")\n");
|
| - printf ("#define MP_BASES_BIG_BASE_INVERTED_10 CNST_LIMB(0x");
|
| - mpz_out_str (stdout, 16, big_base_inverted);
|
| - printf (")\n");
|
| - printf ("#define MP_BASES_NORMALIZATION_STEPS_10 %d\n", normalization_steps);
|
| -}
|
| -
|
| -void
|
| -table (int limb_bits, int nail_bits)
|
| -{
|
| - int numb_bits = limb_bits - nail_bits;
|
| - int base;
|
| -
|
| - printf ("/* This file generated by gen-bases.c - DO NOT EDIT. */\n");
|
| - printf ("\n");
|
| - printf ("#include \"gmp.h\"\n");
|
| - printf ("#include \"gmp-impl.h\"\n");
|
| - printf ("\n");
|
| - printf ("#if GMP_NUMB_BITS != %d\n", numb_bits);
|
| - printf ("Error, error, this data is for %d bits\n", numb_bits);
|
| - printf ("#endif\n");
|
| - printf ("\n");
|
| - puts ("const struct bases mp_bases[257] =\n{");
|
| - puts (" /* 0 */ { 0, 0.0, 0 },");
|
| - puts (" /* 1 */ { 0, 1e37, 0 },");
|
| - for (base = 2; base <= 256; base++)
|
| - {
|
| - generate (limb_bits, nail_bits, base);
|
| -
|
| - printf (" /* %3u */ { ", base);
|
| - if (POW2_P (base))
|
| - {
|
| - printf ("%u, %.16f, 0x%x },\n",
|
| - chars_per_limb, chars_per_bit_exactly, ulog2 (base) - 1);
|
| - }
|
| - else
|
| - {
|
| - printf ("%u, %.16f, CNST_LIMB(0x",
|
| - chars_per_limb, chars_per_bit_exactly);
|
| - mpz_out_str (stdout, 16, big_base);
|
| - printf ("), CNST_LIMB(0x");
|
| - mpz_out_str (stdout, 16, big_base_inverted);
|
| - printf (") },\n");
|
| - }
|
| - }
|
| -
|
| - puts ("};");
|
| -}
|
| -
|
| -int
|
| -main (int argc, char **argv)
|
| -{
|
| - int limb_bits, nail_bits;
|
| -
|
| - mpz_init (big_base);
|
| - mpz_init (big_base_inverted);
|
| - mpz_init (t);
|
| -
|
| - if (argc != 4)
|
| - {
|
| - fprintf (stderr, "Usage: gen-bases <header|table> <limbbits> <nailbits>\n");
|
| - exit (1);
|
| - }
|
| -
|
| - limb_bits = atoi (argv[2]);
|
| - nail_bits = atoi (argv[3]);
|
| -
|
| - if (limb_bits <= 0
|
| - || nail_bits < 0
|
| - || nail_bits >= limb_bits)
|
| - {
|
| - fprintf (stderr, "Invalid limb/nail bits: %d %d\n",
|
| - limb_bits, nail_bits);
|
| - exit (1);
|
| - }
|
| -
|
| - if (strcmp (argv[1], "header") == 0)
|
| - header (limb_bits, nail_bits);
|
| - else if (strcmp (argv[1], "table") == 0)
|
| - table (limb_bits, nail_bits);
|
| - else
|
| - {
|
| - fprintf (stderr, "Invalid header/table choice: %s\n", argv[1]);
|
| - exit (1);
|
| - }
|
| -
|
| - return 0;
|
| -}
|
|
|