Home | History | Annotate | Line # | Download | only in libm 
      1 /* $NetBSD: t_casinh.c,v 1.2 2016/09/20 17:19:28 christos Exp $ */
      2 
      3 /*
      4  * Written by Maya Rashish
      5  * Public domain.
      6  *
      7  * Testing special values of casinh
      8  * Values from ISO/IEC 9899:201x G.6.2.2
      9  */
     10 
     11 #include <atf-c.h>
     12 #include <complex.h>
     13 #include <math.h>
     14 
     15 #define RE(z) (((double *)(&z))[0])
     16 #define IM(z) (((double *)(&z))[1])
     17 
     18 static const struct {
     19 	double input_re;
     20 	double input_im;
     21 	double result_re;
     22 	double result_im;
     23 } values[] = {
     24 	{ +0,		+0,		+0,		+0},
     25 	{ +5.032E3,	+INFINITY,	+INFINITY,	+M_PI/2},
     26 	{ +INFINITY,	+5.023E3,	+INFINITY,	+0},
     27 	{ +INFINITY,	+INFINITY,	+INFINITY,	+M_PI/4},
     28 #ifdef __HAVE_NANF
     29 	{ +INFINITY,	+NAN,		+INFINITY,	+NAN},
     30 	{ +5.032E3,	+NAN,		+NAN,		+NAN}, /* + FE_INVALID optionally raised */
     31 	{ +NAN,		+0,		+NAN,		+0},
     32 	{ +NAN,		-5.023E3,	+NAN,		+NAN}, /* + FE_INVALID optionally raised */
     33 	{ +NAN,		+INFINITY,	+INFINITY,	+NAN}, /* sign of real part of result unspecified */
     34 	{ +NAN,		+NAN,		+NAN,		+NAN},
     35 #endif
     36 };
     37 
     38 #ifdef __HAVE_NANF
     39 #define both_nan(a,b) (isnan(a) && isnan(b))
     40 #else
     41 #define both_nan(a,b) 0
     42 #endif
     43 
     44 #define crude_equality(a,b) ((a == b) || both_nan(a,b))
     45 
     46 #define ATF_COMPLEX_EQUAL(a,b) do { \
     47 	complex double ci = casinh(a); \
     48 	ATF_CHECK_MSG(crude_equality(creal(ci),creal(b)) && \
     49 	    crude_equality(cimag(ci), cimag(b)), \
     50 	    "for casinh([%g,%g]) = [%g,%g] != [%g,%g]", \
     51 	    creal(a), cimag(a), creal(ci), cimag(ci), creal(b), cimag(b)); \
     52 } while (0/*CONSTCOND*/)
     53 
     54 
     55 ATF_TC(casinh);
     56 ATF_TC_HEAD(casinh, tc)
     57 {
     58 	atf_tc_set_md_var(tc, "descr","Check casinh family - special values");
     59 }
     60 
     61 ATF_TC_BODY(casinh, tc)
     62 {
     63 	complex double input;
     64 	complex double result;
     65 	unsigned int i;
     66 	for (i = 0; i < __arraycount(values); i++) {
     67 		RE(input) = values[i].input_re;
     68 		IM(input) = values[i].input_im;
     69 		RE(result) = values[i].result_re;
     70 		IM(result) = values[i].result_im;
     71 		ATF_COMPLEX_EQUAL(input, result);
     72 	}
     73 }
     74 
     75 ATF_TP_ADD_TCS(tp)
     76 {
     77 
     78 	ATF_TP_ADD_TC(tp, casinh);
     79 
     80 	return atf_no_error();
     81 }
     82