1 /* Generate table of tests in tst-strtod-round.c from
2 tst-strtod-round-data.
3 Copyright (C) 2012-2022 Free Software Foundation, Inc.
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 <https://www.gnu.org/licenses/>. */
19
20 /* Compile this program as:
21
22 gcc -std=gnu11 -O2 -Wall -Wextra gen-tst-strtod-round.c -lmpfr \
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
33 #define _GNU_SOURCE
34 #include <assert.h>
35 #include <stdbool.h>
36 #include <stdio.h>
37 #include <stdlib.h>
38 #include <string.h>
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
string_to_fp(mpfr_t f,const char * s,mpfr_rnd_t rnd)46 string_to_fp (mpfr_t f, const char *s, mpfr_rnd_t rnd)
47 {
48 mpfr_clear_overflow ();
49 #ifdef WORKAROUND
50 mpfr_t f2;
51 mpfr_init2 (f2, 100000);
52 int r0 = mpfr_strtofr (f2, s, NULL, 0, rnd);
53 int r = mpfr_set (f, f2, rnd);
54 r |= mpfr_subnormalize (f, r, rnd);
55 mpfr_clear (f2);
56 return r0 | r;
57 #else
58 int r = mpfr_strtofr (f, s, NULL, 0, rnd);
59 r |= mpfr_subnormalize (f, r, rnd);
60 return r;
61 #endif
62 }
63
64 void
print_fp(FILE * fout,mpfr_t f,const char * suffix)65 print_fp (FILE *fout, mpfr_t f, const char *suffix)
66 {
67 if (mpfr_inf_p (f))
68 mpfr_fprintf (fout, "\t%sINF%s", mpfr_signbit (f) ? "-" : "", suffix);
69 else
70 mpfr_fprintf (fout, "\t%Ra%s", f, suffix);
71 }
72
73 static void
round_str(FILE * fout,const char * s,int prec,int emin,int emax,bool ibm_ld)74 round_str (FILE *fout, const char *s, int prec, int emin, int emax,
75 bool ibm_ld)
76 {
77 mpfr_t max_value;
78 mpfr_t f;
79 mpfr_set_default_prec (prec);
80 mpfr_set_emin (emin);
81 mpfr_set_emax (emax);
82 mpfr_init (f);
83 int r = string_to_fp (f, s, MPFR_RNDD);
84 bool overflow = mpfr_overflow_p () != 0;
85 if (ibm_ld)
86 {
87 assert (prec == 106 && emin == -1073 && emax == 1024);
88 /* The maximum value in IBM long double has discontiguous
89 mantissa bits. */
90 mpfr_init2 (max_value, 107);
91 mpfr_set_str (max_value, "0x1.fffffffffffff7ffffffffffffcp+1023", 0,
92 MPFR_RNDN);
93 if (mpfr_cmpabs (f, max_value) > 0)
94 {
95 r = 1;
96 overflow = true;
97 }
98 }
99 mpfr_fprintf (fout, "\t%s,\n", r ? "false" : "true");
100 print_fp (fout, f, overflow ? ", true,\n" : ", false,\n");
101 string_to_fp (f, s, MPFR_RNDN);
102 overflow = (mpfr_overflow_p () != 0
103 || (ibm_ld && mpfr_cmpabs (f, max_value) > 0));
104 print_fp (fout, f, overflow ? ", true,\n" : ", false,\n");
105 string_to_fp (f, s, MPFR_RNDZ);
106 overflow = (mpfr_overflow_p () != 0
107 || (ibm_ld && mpfr_cmpabs (f, max_value) > 0));
108 print_fp (fout, f, overflow ? ", true,\n" : ", false,\n");
109 string_to_fp (f, s, MPFR_RNDU);
110 overflow = (mpfr_overflow_p () != 0
111 || (ibm_ld && mpfr_cmpabs (f, max_value) > 0));
112 print_fp (fout, f, overflow ? ", true" : ", false");
113 mpfr_clear (f);
114 if (ibm_ld)
115 mpfr_clear (max_value);
116 }
117
118 static void
round_for_all(FILE * fout,const char * s)119 round_for_all (FILE *fout, const char *s)
120 {
121 static const struct fmt {
122 int prec;
123 int emin;
124 int emax;
125 bool ibm_ld;
126 } formats[] = {
127 { 24, -148, 128, false },
128 { 53, -1073, 1024, false },
129 /* This is the Intel extended float format. */
130 { 64, -16444, 16384, false },
131 /* This is the Motorola extended float format. */
132 { 64, -16445, 16384, false },
133 { 106, -1073, 1024, true },
134 { 113, -16493, 16384, false },
135 };
136 mpfr_fprintf (fout, " TEST (\"");
137 const char *p;
138 for (p = s; *p; p++)
139 {
140 fputc (*p, fout);
141 if ((p - s) % 60 == 59 && p[1])
142 mpfr_fprintf (fout, "\"\n\t\"");
143 }
144 mpfr_fprintf (fout, "\",\n");
145 int i;
146 int n_formats = sizeof (formats) / sizeof (formats[0]);
147 for (i = 0; i < n_formats; i++)
148 {
149 round_str (fout, s, formats[i].prec, formats[i].emin,
150 formats[i].emax, formats[i].ibm_ld);
151 if (i < n_formats - 1)
152 mpfr_fprintf (fout, ",\n");
153 }
154 mpfr_fprintf (fout, "),\n");
155 }
156
157 int
main(int argc,char ** argv)158 main (int argc, char **argv)
159 {
160 char *p = NULL;
161 size_t len;
162 ssize_t nbytes;
163 FILE *fin, *fout;
164 char *fin_name, *fout_name;
165
166 if (argc < 3)
167 {
168 fprintf (stderr, "Usage: %s <input> <output>\n", basename (argv[0]));
169 return EXIT_FAILURE;
170 }
171
172 fin_name = argv[1];
173 fout_name = argv[2];
174
175 fin = fopen (fin_name, "r");
176 if (fin == NULL)
177 {
178 perror ("Could not open input for reading");
179 return EXIT_FAILURE;
180 }
181
182 fout = fopen (fout_name, "w");
183 if (fout == NULL)
184 {
185 perror ("Could not open output for writing");
186 return EXIT_FAILURE;
187 }
188
189 fprintf (fout, "/* This file was generated by %s from %s. */\n",
190 __FILE__, fin_name);
191 fputs ("static const struct test tests[] = {\n", fout);
192 while ((nbytes = getline (&p, &len, fin)) != -1)
193 {
194 if (p[nbytes - 1] == '\n')
195 p[nbytes - 1] = 0;
196 round_for_all (fout, p);
197 free (p);
198 p = NULL;
199 }
200 fputs ("};\n", fout);
201
202 return EXIT_SUCCESS;
203 }
204