Home | History | Annotate | Line # | Download | only in libm
t_libm.h revision 1.2
      1  1.2  dsl /* $NetBSD: t_libm.h,v 1.2 2014/03/05 20:14:46 dsl Exp $ */
      2  1.1  dsl 
      3  1.1  dsl /*
      4  1.1  dsl  * Check result of fn(arg) is correct within the bounds.
      5  1.1  dsl  * Should be ok to do the checks using 'double' for 'float' functions.
      6  1.1  dsl  */
      7  1.1  dsl #define T_LIBM_CHECK(subtest, fn, arg, expect, epsilon) do { \
      8  1.1  dsl 	double r = fn(arg); \
      9  1.1  dsl 	double e = fabs(r - expect); \
     10  1.1  dsl 	if (e > epsilon) \
     11  1.1  dsl 		atf_tc_fail_nonfatal( \
     12  1.1  dsl 		    "subtest %u: " #fn "(%g) is %g not %g (error %g > %g)", \
     13  1.1  dsl 		    subtest, arg, r, expect, e, epsilon); \
     14  1.1  dsl     } while (0)
     15  1.1  dsl 
     16  1.1  dsl /* Check that the result of fn(arg) is NaN */
     17  1.1  dsl #ifndef __vax__
     18  1.1  dsl #define T_LIBM_CHECK_NAN(subtest, fn, arg) do { \
     19  1.1  dsl 	double r = fn(arg); \
     20  1.1  dsl 	if (!isnan(r)) \
     21  1.1  dsl 		atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not NaN", \
     22  1.1  dsl 		    subtest, arg, r); \
     23  1.1  dsl     } while (0)
     24  1.1  dsl #else
     25  1.1  dsl /* vax doesn't support NaN */
     26  1.1  dsl #define T_LIBM_CHECK_NAN(subtest, fn, arg) (void)(arg)
     27  1.1  dsl #endif
     28  1.1  dsl 
     29  1.1  dsl /* Check that the result of fn(arg) is +0.0 */
     30  1.1  dsl #define T_LIBM_CHECK_PLUS_ZERO(subtest, fn, arg) do { \
     31  1.1  dsl 	double r = fn(arg); \
     32  1.1  dsl 	if (fabs(r) > 0.0 || signbit(r) != 0) \
     33  1.1  dsl 		atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not +0.0", \
     34  1.1  dsl 		    subtest, arg, r); \
     35  1.1  dsl     } while (0)
     36  1.1  dsl 
     37  1.1  dsl /* Check that the result of fn(arg) is -0.0 */
     38  1.1  dsl #define T_LIBM_CHECK_MINUS_ZERO(subtest, fn, arg) do { \
     39  1.1  dsl 	double r = fn(arg); \
     40  1.1  dsl 	if (fabs(r) > 0.0 || signbit(r) == 0) \
     41  1.1  dsl 		atf_tc_fail_nonfatal("subtest %u: " #fn "(%g) is %g not -0.0", \
     42  1.1  dsl 		    subtest, arg, r); \
     43  1.1  dsl     } while (0)
     44  1.1  dsl 
     45  1.1  dsl /* Some useful constants (for test vectors) */
     46  1.2  dsl #define T_LIBM_NAN	(0.0 / 0.0)
     47  1.2  dsl #define T_LIBM_PLUS_INF	(+1.0 / 0.0)
     48  1.2  dsl #define T_LIBM_MINUS_INF (-1.0 / 0.0)
     49  1.1  dsl 
     50  1.1  dsl /* One line definition of a simple test */
     51  1.2  dsl #define ATF_LIBM_TEST(name, description) \
     52  1.1  dsl ATF_TC(name); \
     53  1.1  dsl ATF_TC_HEAD(name, tc) { atf_tc_set_md_var(tc, "descr", description); } \
     54  1.1  dsl ATF_TC_BODY(name, tc)
     55