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