数学计算

开发平台：
Unix_Linux

gen-bases.c：源码内容
							/* 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 != %dn", numb_bits);
  printf ("Error, error, this data is for %d bitsn", numb_bits);
  printf ("#endifn");
  printf ("n");
  printf ("/* mp_bases[10] data, as literal values */n");
  printf ("#define MP_BASES_CHARS_PER_LIMB_10      %dn", 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 %dn", 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 != %dn", numb_bits);
  printf ("Error, error, this data is for %d bitsn", numb_bits);
  printf ("#endifn");
  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 %dn",
               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: %sn", argv[1]);
      exit (1);
    }
  return 0;
}

						