]>
Commit | Line | Data |
---|---|---|
af92131a JM |
1 | /* Generate table of tests in tst-strtod-round.c from |
2 | tst-strtod-round-data. | |
bfff8b1b | 3 | Copyright (C) 2012-2017 Free Software Foundation, Inc. |
af92131a JM |
4 | This file is part of the GNU C Library. |
5 | ||
6 | The GNU C Library is free software; you can redistribute it and/or | |
7 | modify it under the terms of the GNU Lesser General Public | |
8 | License as published by the Free Software Foundation; either | |
9 | version 2.1 of the License, or (at your option) any later version. | |
10 | ||
11 | The GNU C Library is distributed in the hope that it will be useful, | |
12 | but WITHOUT ANY WARRANTY; without even the implied warranty of | |
13 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
14 | Lesser General Public License for more details. | |
15 | ||
16 | You should have received a copy of the GNU Lesser General Public | |
17 | License along with the GNU C Library; if not, see | |
18 | <http://www.gnu.org/licenses/>. */ | |
19 | ||
118fbf0e PM |
20 | /* Compile this program as: |
21 | ||
1ced34c0 | 22 | gcc -std=gnu11 -O2 -Wall -Wextra gen-tst-strtod-round.c -lmpfr \ |
118fbf0e PM |
23 | -o gen-tst-strtod-round |
24 | ||
25 | (use of current MPFR version recommended) and run it as: | |
26 | ||
27 | gen-tst-strtod-round tst-strtod-round-data tst-strtod-round-data.h | |
28 | ||
29 | The output file will be generated as tst-strtod-round-data.h | |
30 | */ | |
31 | ||
32 | ||
af92131a | 33 | #define _GNU_SOURCE |
ed8c2ecd | 34 | #include <assert.h> |
af92131a JM |
35 | #include <stdbool.h> |
36 | #include <stdio.h> | |
37 | #include <stdlib.h> | |
118fbf0e | 38 | #include <string.h> |
af92131a JM |
39 | #include <mpfr.h> |
40 | ||
41 | /* Work around incorrect ternary value from mpfr_strtofr | |
42 | <https://sympa.inria.fr/sympa/arc/mpfr/2012-08/msg00005.html>. */ | |
43 | #define WORKAROUND | |
44 | ||
45 | static int | |
46 | string_to_fp (mpfr_t f, const char *s, mpfr_rnd_t rnd) | |
47 | { | |
48 | #ifdef WORKAROUND | |
49 | mpfr_t f2; | |
50 | mpfr_init2 (f2, 100000); | |
51 | int r0 = mpfr_strtofr (f2, s, NULL, 0, rnd); | |
52 | int r = mpfr_set (f, f2, rnd); | |
4725d33e | 53 | r |= mpfr_subnormalize (f, r, rnd); |
af92131a JM |
54 | mpfr_clear (f2); |
55 | return r0 | r; | |
56 | #else | |
57 | int r = mpfr_strtofr (f, s, NULL, 0, rnd); | |
4725d33e | 58 | r |= mpfr_subnormalize (f, r, rnd); |
af92131a JM |
59 | return r; |
60 | #endif | |
61 | } | |
62 | ||
1ced34c0 PM |
63 | void |
64 | print_fp (FILE *fout, mpfr_t f, const char *suffix) | |
af92131a JM |
65 | { |
66 | if (mpfr_inf_p (f)) | |
1ced34c0 | 67 | mpfr_fprintf (fout, "\t%sINF%s", mpfr_signbit (f) ? "-" : "", suffix); |
af92131a | 68 | else |
1ced34c0 | 69 | mpfr_fprintf (fout, "\t%Ra%s", f, suffix); |
af92131a JM |
70 | } |
71 | ||
72 | static void | |
1ced34c0 PM |
73 | round_str (FILE *fout, const char *s, int prec, int emin, int emax, |
74 | bool ibm_ld) | |
af92131a JM |
75 | { |
76 | mpfr_t f; | |
77 | mpfr_set_default_prec (prec); | |
78 | mpfr_set_emin (emin); | |
79 | mpfr_set_emax (emax); | |
80 | mpfr_init (f); | |
81 | int r = string_to_fp (f, s, MPFR_RNDD); | |
1f24b9ad | 82 | if (ibm_ld) |
ed8c2ecd JM |
83 | { |
84 | assert (prec == 106 && emin == -1073 && emax == 1024); | |
85 | /* The maximum value in IBM long double has discontiguous | |
86 | mantissa bits. */ | |
87 | mpfr_t max_value; | |
88 | mpfr_init2 (max_value, 107); | |
89 | mpfr_set_str (max_value, "0x1.fffffffffffff7ffffffffffffcp+1023", 0, | |
90 | MPFR_RNDN); | |
91 | if (mpfr_cmpabs (f, max_value) > 0) | |
92 | r = 1; | |
ed8c2ecd JM |
93 | mpfr_clear (max_value); |
94 | } | |
118fbf0e | 95 | mpfr_fprintf (fout, "\t%s,\n", r ? "false" : "true"); |
1ced34c0 | 96 | print_fp (fout, f, ",\n"); |
af92131a | 97 | string_to_fp (f, s, MPFR_RNDN); |
1ced34c0 | 98 | print_fp (fout, f, ",\n"); |
af92131a | 99 | string_to_fp (f, s, MPFR_RNDZ); |
1ced34c0 | 100 | print_fp (fout, f, ",\n"); |
af92131a | 101 | string_to_fp (f, s, MPFR_RNDU); |
1ced34c0 | 102 | print_fp (fout, f, ""); |
af92131a JM |
103 | mpfr_clear (f); |
104 | } | |
105 | ||
106 | static void | |
118fbf0e | 107 | round_for_all (FILE *fout, const char *s) |
af92131a JM |
108 | { |
109 | static const struct fmt { | |
af92131a JM |
110 | int prec; |
111 | int emin; | |
112 | int emax; | |
1f24b9ad | 113 | bool ibm_ld; |
1ced34c0 PM |
114 | } formats[] = { |
115 | { 24, -148, 128, false }, | |
116 | { 53, -1073, 1024, false }, | |
e1343020 | 117 | /* This is the Intel extended float format. */ |
1ced34c0 | 118 | { 64, -16444, 16384, false }, |
e1343020 | 119 | /* This is the Motorola extended float format. */ |
1ced34c0 PM |
120 | { 64, -16445, 16384, false }, |
121 | { 106, -1073, 1024, true }, | |
122 | { 113, -16493, 16384, false }, | |
af92131a | 123 | }; |
118fbf0e | 124 | mpfr_fprintf (fout, " TEST (\""); |
af92131a JM |
125 | const char *p; |
126 | for (p = s; *p; p++) | |
127 | { | |
118fbf0e | 128 | fputc (*p, fout); |
af92131a | 129 | if ((p - s) % 60 == 59 && p[1]) |
118fbf0e | 130 | mpfr_fprintf (fout, "\"\n\t\""); |
af92131a | 131 | } |
118fbf0e | 132 | mpfr_fprintf (fout, "\",\n"); |
af92131a | 133 | int i; |
1ced34c0 PM |
134 | int n_formats = sizeof (formats) / sizeof (formats[0]); |
135 | for (i = 0; i < n_formats; i++) | |
af92131a | 136 | { |
1ced34c0 PM |
137 | round_str (fout, s, formats[i].prec, formats[i].emin, |
138 | formats[i].emax, formats[i].ibm_ld); | |
139 | if (i < n_formats - 1) | |
118fbf0e | 140 | mpfr_fprintf (fout, ",\n"); |
af92131a | 141 | } |
118fbf0e | 142 | mpfr_fprintf (fout, "),\n"); |
af92131a JM |
143 | } |
144 | ||
145 | int | |
118fbf0e | 146 | main (int argc, char **argv) |
af92131a JM |
147 | { |
148 | char *p = NULL; | |
149 | size_t len; | |
150 | ssize_t nbytes; | |
118fbf0e PM |
151 | FILE *fin, *fout; |
152 | char *fin_name, *fout_name; | |
153 | ||
154 | if (argc < 3) | |
155 | { | |
156 | fprintf (stderr, "Usage: %s <input> <output>\n", basename (argv[0])); | |
157 | return EXIT_FAILURE; | |
158 | } | |
159 | ||
160 | fin_name = argv[1]; | |
161 | fout_name = argv[2]; | |
162 | ||
163 | fin = fopen (fin_name, "r"); | |
164 | if (fin == NULL) | |
165 | { | |
166 | perror ("Could not open input for reading"); | |
167 | return EXIT_FAILURE; | |
168 | } | |
169 | ||
170 | fout = fopen (fout_name, "w"); | |
171 | if (fout == NULL) | |
172 | { | |
173 | perror ("Could not open output for writing"); | |
174 | return EXIT_FAILURE; | |
175 | } | |
176 | ||
177 | fprintf (fout, "/* This file was generated by %s from %s. */\n", | |
178 | __FILE__, fin_name); | |
179 | fputs ("static const struct test tests[] = {\n", fout); | |
180 | while ((nbytes = getline (&p, &len, fin)) != -1) | |
af92131a JM |
181 | { |
182 | if (p[nbytes - 1] == '\n') | |
183 | p[nbytes - 1] = 0; | |
118fbf0e | 184 | round_for_all (fout, p); |
af92131a JM |
185 | free (p); |
186 | p = NULL; | |
187 | } | |
118fbf0e PM |
188 | fputs ("};\n", fout); |
189 | ||
190 | return EXIT_SUCCESS; | |
af92131a | 191 | } |