2  * SPDX-License-Identifier: BSD-3-Clause
45  * fall-through case statement annotations
51 /* Non-standard but supported by at least gcc and clang */
74         e = scalbnl(1.0L, -ilogbl(expect));  in within_error()
76     difference = fabsl(expect - result) * e;  in within_error()
88         long double diff = fabsl(expect - result);  in check_long_double()
108         long long diff = expect - result;  in check_long_long()
162 #define DEFAULT_PREC 0x1p-55L
163 #define SQRTL_PREC 0x1.0p-63L
167 #define DEFAULT_PREC 0x1p-105L
168 #define SQRTL_PREC 0x1.0p-112L
170 #define DEFAULT_PREC 0x1p-97L
171 #define SQRTL_PREC 0x1.0p-105L
174 #define DEFAULT_PREC 0x1p-48L
175 #define SQRTL_PREC 0x1.0p-52L
225     long double v = 0.0L;  in naive_strtold()  local
250             return -(long double)INFINITY;  in naive_strtold()
258             return -(long double)INFINITY;  in naive_strtold()
259         case '-':  in naive_strtold()
261                 exp_sign = -1;  in naive_strtold()
264             return -(long double)INFINITY;  in naive_strtold()
270             return -(long double)INFINITY;  in naive_strtold()
273             digit = c - '0';  in naive_strtold()
281                 return -(long double)INFINITY;  in naive_strtold()
288                 return -(long double)INFINITY;  in naive_strtold()
294                 digit = c - 'A' + 10;  in naive_strtold()
297             return -(long double)INFINITY;  in naive_strtold()
304                 return -(long double)INFINITY;  in naive_strtold()
310                 digit = c - 'a' + 10;  in naive_strtold()
313             return -(long double)INFINITY;  in naive_strtold()
315             return -(long double)INFINITY;  in naive_strtold()
319             v = v * base + digit;  in naive_strtold()
322             v = v + digit * frac_mul;  in naive_strtold()
332         long ebot = exp - etop;  in naive_strtold()
335         long double vpow = v * epow_top;  in naive_strtold()
339         return ldexpl(v, exp * exp_sign);  in naive_strtold()
357     return fabsl((have - want) / want) <= max_error;  in close()
361     __LDBL_MIN_EXP__ - __LDBL_MANT_DIG__ - 1,
362     __LDBL_MIN_EXP__ - __LDBL_MANT_DIG__,
363     __LDBL_MIN_EXP__ - __LDBL_MANT_DIG__ + 1,
364     __LDBL_MIN_EXP__ - __LDBL_MANT_DIG__ + 2,
365     __LDBL_MIN_EXP__ - __LDBL_MANT_DIG__ + 3,
366     __LDBL_MIN_EXP__ - 3,
367     __LDBL_MIN_EXP__ - 2,
368     __LDBL_MIN_EXP__ - 1,
373     -3,
374     -2,
375     -1,
380     __LDBL_MAX_EXP__ - 3,
381     __LDBL_MAX_EXP__ - 2,
382     __LDBL_MAX_EXP__ - 1,
390  * For 64-bit values, we may have exact conversions. Otherwise, allow
397 #  define MAX_DECIMAL_ERROR     1e-10L
401 #  define MAX_DECIMAL_ERROR       1e-5L
403 #  define MAX_DECIMAL_ERROR       1e-10L
421             long double v, r;  in test_io()  local
422             v = ldexpl(vals[i], test_exp[e]);  in test_io()
431                     max_error_naive = 1e-6L;  in test_io()
434                 sprintf(buf, formats[j], v);  in test_io()
435                 if (isinf(v)) {  in test_io()
440                 } else if (isnan(v)) {  in test_io()
447                     if (!close(r, v, max_error_naive)) {  in test_io()
448 …io naive i %d val %La exp %d: \"%s\", is %La should be %La\n", i, vals[i], test_exp[e], buf, r, v);  in test_io()
453                 if (!close(r, v, max_error) && !(isnan(v) && isnan(r))) {  in test_io()
454 …io scanf i %d val %La exp %d: \"%s\", is %La should be %La\n", i, vals[i], test_exp[e], buf, r, v);  in test_io()
458 …             if ((!close(r, v, max_error) && !(isnan(v) && isnan(r)))|| end != buf + strlen(buf)) {  in test_io()
459 … strtold i %d val %La exp %d: \"%s\", is %La should be %La\n", i, vals[i], test_exp[e], buf, r, v);  in test_io()