1 | /* Generate data for fpioconst.c. |
2 | Copyright (C) 2012-2022 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 | <https://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 | } |
109 | |