stdlib/gen-fpioconst.c

   1 /* Generate data for fpioconst.c.
   2    Copyright (C) 2012-2014 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4
   5    The GNU C Library is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU Lesser General Public
   7    License as published by the Free Software Foundation; either
   8    version 2.1 of the License, or (at your option) any later version.
   9
  10    The GNU C Library is distributed in the hope that it will be useful,
  11    but WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13    Lesser General Public License for more details.
  14
  15    You should have received a copy of the GNU Lesser General Public
  16    License along with the GNU C Library; if not, see
  17    <http://www.gnu.org/licenses/>.  */
  18
  19 #include <inttypes.h>
  20 #include <stdio.h>
  21 #include <stdlib.h>
  22 #include <string.h>
  23 #include <gmp.h>
  24 #include <stdint.h>
  25
  26 int
  27 main (void)
  28 {
  29   FILE *out32 = fopen ("fpioconst-32", "w");
  30   if (out32 == NULL)
  31     abort ();
  32   FILE *out64 = fopen ("fpioconst-64", "w");
  33   if (out64 == NULL)
  34     abort ();
  35   FILE *outtable = fopen ("fpioconst-table", "w");
  36   if (outtable == NULL)
  37     abort ();
  38   mpz_t p;
  39   mpz_init (p);
  40   for (int i = 0; i <= 14; i++)
  41     {
  42       int j = 1 << i;
  43       mpz_ui_pow_ui (p, 10, j - 1);
  44       int exp_m = mpz_sizeinbase (p, 2);
  45       mpz_ui_pow_ui (p, 10, j);
  46       int exp_p = mpz_sizeinbase (p, 2);
  47       int size32 = 2 + (exp_p + 31) / 32;
  48       int size64 = 1 + (exp_p + 63) / 64;
  49       uint32_t data32[size32];
  50       uint64_t data64[size64];
  51       memset (data32, 0, sizeof data32);
  52       memset (data64, 0, sizeof data64);
  53       mpz_export (data32 + 2, NULL, -1, 4, 0, 0, p);
  54       mpz_export (data64 + 1, NULL, -1, 8, 0, 0, p);
  55       if (i == 0)
  56         {
  57           fprintf (out32, "#define TENS_P%d_IDX\t0\n", i);
  58           fprintf (out64, "#define TENS_P%d_IDX\t0\n", i);
  59         }
  60       else
  61         {
  62           fprintf (out32, "#define TENS_P%d_IDX\t"
  63                    "(TENS_P%d_IDX + TENS_P%d_SIZE)\n",
  64                    i, i - 1, i - 1);
  65           fprintf (out64, "#define TENS_P%d_IDX\t"
  66                    "(TENS_P%d_IDX + TENS_P%d_SIZE)\n",
  67                    i, i - 1, i - 1);
  68         }
  69       fprintf (out32, "#define TENS_P%d_SIZE\t%d\n", i, size32);
  70       fprintf (out64, "#define TENS_P%d_SIZE\t%d\n", i, size64);
  71       for (int k = 0; k < size32; k++)
  72         {
  73           if (k == 0)
  74             fprintf (out32, "  [TENS_P%d_IDX] = ", i);
  75           else if (k % 6 == 5)
  76             fprintf (out32, "\n  ");
  77           else
  78             fprintf (out32, " ");
  79           fprintf (out32, "0x%08"PRIx32",", data32[k]);
  80         }
  81       for (int k = 0; k < size64; k++)
  82         {
  83           if (k == 0)
  84             fprintf (out64, "  [TENS_P%d_IDX] = ", i);
  85           else if (k % 3 == 2)
  86             fprintf (out64, "\n  ");
  87           else
  88             fprintf (out64, " ");
  89           fprintf (out64, "0x%016"PRIx64"ull,", data64[k]);
  90         }
  91       fprintf (out32, "\n\n");
  92       fprintf (out64, "\n\n");
  93       const char *t = (i >= 10 ? "\t" : "\t\t");
  94       if (i == 0)
  95         fprintf (outtable, "  { TENS_P%d_IDX, TENS_P%d_SIZE,%s%d,\t      },\n",
  96                  i, i, t, exp_p);
  97       else
  98         fprintf (outtable, "  { TENS_P%d_IDX, TENS_P%d_SIZE,%s%d,\t%5d },\n",
  99                  i, i, t, exp_p, exp_m);
 100     }
 101   if (fclose (out32) != 0)
 102     abort ();
 103   if (fclose (out64) != 0)
 104     abort ();
 105   if (fclose (outtable) != 0)
 106     abort ();
 107   return 0;
 108 }