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