1 | /* Copyright (C) 2012-2022 Free Software Foundation, Inc. |
2 | This file is part of the GNU C Library. |
3 | |
4 | The GNU C Library is free software; you can redistribute it and/or |
5 | modify it under the terms of the GNU Lesser General Public |
6 | License as published by the Free Software Foundation; either |
7 | version 2.1 of the License, or (at your option) any later version. |
8 | |
9 | The GNU C Library is distributed in the hope that it will be useful, |
10 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
11 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
12 | Lesser General Public License for more details. |
13 | |
14 | You should have received a copy of the GNU Lesser General Public |
15 | License along with the GNU C Library; if not, see |
16 | <https://www.gnu.org/licenses/>. */ |
17 | |
18 | /* Test bug #13941. */ |
19 | |
20 | #include <float.h> |
21 | #include <math.h> |
22 | #include <stdio.h> |
23 | #include <string.h> |
24 | |
25 | static int |
26 | do_test (void) |
27 | { |
28 | #if LDBL_MANT_DIG >= 106 |
29 | volatile union { long double l; long long x[2]; } u, v; |
30 | char buf[64]; |
31 | #endif |
32 | int result = 0; |
33 | |
34 | #if LDBL_MANT_DIG == 106 || LDBL_MANT_DIG == 113 |
35 | # define COMPARE_LDBL(u, v) \ |
36 | ((u).l == (v).l && (u).x[0] == (v).x[0] && (u).x[1] == (v).x[1]) |
37 | #else |
38 | # define COMPARE_LDBL(u, v) ((u).l == (v).l) |
39 | #endif |
40 | |
41 | #define TEST_N(val, n) \ |
42 | do \ |
43 | { \ |
44 | u.l = (val); \ |
45 | snprintf (buf, sizeof buf, "%." #n "LgL", u.l); \ |
46 | if (strcmp (buf, #val) != 0) \ |
47 | { \ |
48 | printf ("Error on line %d: %s != %s\n", __LINE__, buf, #val); \ |
49 | result = 1; \ |
50 | } \ |
51 | if (sscanf (#val, "%Lg", &v.l) != 1 || !COMPARE_LDBL (u, v)) \ |
52 | { \ |
53 | printf ("Error sscanf on line %d: %." #n "Lg != %." #n "Lg\n", \ |
54 | __LINE__, u.l, v.l); \ |
55 | result = 1; \ |
56 | } \ |
57 | /* printf ("%s %Lg %016Lx %016Lx\n", #val, u.l, u.x[0], u.x[1]); */ \ |
58 | } \ |
59 | while (0) |
60 | |
61 | #define TEST(val) TEST_N (val,30) |
62 | |
63 | #if LDBL_MANT_DIG >= 106 |
64 | # if LDBL_MANT_DIG == 106 |
65 | TEST (2.22507385850719347803989925739e-308L); |
66 | TEST (2.22507385850719397210554509863e-308L); |
67 | TEST (2.22507385850720088902458687609e-308L); |
68 | |
69 | /* Verify precision is not lost for long doubles |
70 | of the form +1.pN,-1.pM. */ |
71 | TEST_N (3.32306998946228968225951765070082e+35L, 34); |
72 | # endif |
73 | TEST (2.22507385850720138309023271733e-308L); |
74 | TEST (2.22507385850720187715587855858e-308L); |
75 | TEST (2.2250738585074419930597574044e-308L); |
76 | TEST (4.45014771701440227211481959342e-308L); |
77 | TEST (4.45014771701440276618046543466e-308L); |
78 | TEST (4.45014771701440375431175711716e-308L); |
79 | TEST (4.45014771701440474244304879965e-308L); |
80 | TEST (7.12023634722304600689881138745e-307L); |
81 | TEST (1.13923781555569064960474854133e-305L); |
82 | TEST (1.13777777777777776389998996996L); |
83 | TEST (1.13777777777777765287768750745L); |
84 | TEST (20988295479420645138.2044444444L); |
85 | TEST (20988295479420643090.2044444444L); |
86 | TEST (2.14668699894294423266045294316e-292L); |
87 | # if LDBL_MANT_DIG == 106 |
88 | TEST (-2.35993711055432139266626434123e-292L); |
89 | TEST (6.26323524637968345414769634658e-302L); |
90 | TEST (1.49327164802066885331814201989e-308L); |
91 | TEST (3.71834550652787023640837473722e-308L); |
92 | TEST (9.51896449671134907001349268087e-306L); |
93 | # endif |
94 | #endif |
95 | return result; |
96 | } |
97 | |
98 | #define TEST_FUNCTION do_test () |
99 | #include "../test-skeleton.c" |
100 | |