1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
/* $NetBSD: t_libm.h,v 1.6 2014/03/25 17:30:14 joerg Exp $ */
/*
* Check result of fn(arg) is correct within the bounds.
* Should be ok to do the checks using 'double' for 'float' functions.
* On i386 float and double values are returned on the x87 stack and might
* be out of range for the function - so save and print as 'long double'.
* (otherwise you can get 'inf != inf' reported!)
*/
#define T_LIBM_CHECK(subtest, fn, arg, expect_, epsilon_) do { \
long double epsilon = epsilon_; \
long double expect = expect_; \
long double r = fn(arg); \
long double e = fabsl(r - expect); \
if (r != expect && e > epsilon) \
atf_tc_fail_nonfatal( \
"subtest %u: " #fn "(%g) is %Lg (%.14La) " \
"not %Lg (%.13La), error %Lg (%.6La) > %Lg", \
subtest, arg, r, r, expect, expect, e, e, epsilon); \
} while (0)
/* Check that the result of fn(arg) is NaN */
#ifndef __vax__
#define T_LIBM_CHECK_NAN(subtest, fn, arg) do { \
double r = fn(arg); \
if (!isnan(r)) \
atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not NaN", \
subtest, arg, r); \
} while (0)
#else
/* vax doesn't support NaN */
#define T_LIBM_CHECK_NAN(subtest, fn, arg) (void)(arg)
#endif
/* Check that the result of fn(arg) is +0.0 */
#define T_LIBM_CHECK_PLUS_ZERO(subtest, fn, arg) do { \
double r = fn(arg); \
if (fabs(r) > 0.0 || signbit(r) != 0) \
atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not +0.0", \
subtest, arg, r); \
} while (0)
/* Check that the result of fn(arg) is -0.0 */
#define T_LIBM_CHECK_MINUS_ZERO(subtest, fn, arg) do { \
double r = fn(arg); \
if (fabs(r) > 0.0 || signbit(r) == 0) \
atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not -0.0", \
subtest, arg, r); \
} while (0)
/* Some useful constants (for test vectors) */
#ifndef __vax__ /* no NAN nor +/- INF on vax */
#define T_LIBM_NAN (0.0 / 0.0)
#define T_LIBM_PLUS_INF (+1.0 / 0.0)
#define T_LIBM_MINUS_INF (-1.0 / 0.0)
#endif
/* One line definition of a simple test */
#define ATF_LIBM_TEST(name, description) \
ATF_TC(name); \
ATF_TC_HEAD(name, tc) { atf_tc_set_md_var(tc, "descr", description); } \
ATF_TC_BODY(name, tc)
|