1311970Sngie/* $NetBSD: t_casinh.c,v 1.2 2016/09/20 17:19:28 christos Exp $ */ 2311970Sngie 3311970Sngie/* 4311970Sngie * Written by Maya Rashish 5311970Sngie * Public domain. 6311970Sngie * 7311970Sngie * Testing special values of casinh 8311970Sngie * Values from ISO/IEC 9899:201x G.6.2.2 9311970Sngie */ 10311970Sngie 11311970Sngie#include <atf-c.h> 12311970Sngie#include <complex.h> 13311970Sngie#include <math.h> 14311970Sngie 15311970Sngie#define RE(z) (((double *)(&z))[0]) 16311970Sngie#define IM(z) (((double *)(&z))[1]) 17311970Sngie 18311970Sngiestatic const struct { 19311970Sngie double input_re; 20311970Sngie double input_im; 21311970Sngie double result_re; 22311970Sngie double result_im; 23311970Sngie} values[] = { 24311970Sngie { +0, +0, +0, +0}, 25311970Sngie { +5.032E3, +INFINITY, +INFINITY, +M_PI/2}, 26311970Sngie { +INFINITY, +5.023E3, +INFINITY, +0}, 27311970Sngie { +INFINITY, +INFINITY, +INFINITY, +M_PI/4}, 28311970Sngie#ifdef __HAVE_NANF 29311970Sngie { +INFINITY, +NAN, +INFINITY, +NAN}, 30311970Sngie { +5.032E3, +NAN, +NAN, +NAN}, /* + FE_INVALID optionally raised */ 31311970Sngie { +NAN, +0, +NAN, +0}, 32311970Sngie { +NAN, -5.023E3, +NAN, +NAN}, /* + FE_INVALID optionally raised */ 33311970Sngie { +NAN, +INFINITY, +INFINITY, +NAN}, /* sign of real part of result unspecified */ 34311970Sngie { +NAN, +NAN, +NAN, +NAN}, 35311970Sngie#endif 36311970Sngie}; 37311970Sngie 38311970Sngie#ifdef __HAVE_NANF 39311970Sngie#define both_nan(a,b) (isnan(a) && isnan(b)) 40311970Sngie#else 41311970Sngie#define both_nan(a,b) 0 42311970Sngie#endif 43311970Sngie 44311970Sngie#define crude_equality(a,b) ((a == b) || both_nan(a,b)) 45311970Sngie 46311970Sngie#define ATF_COMPLEX_EQUAL(a,b) do { \ 47311970Sngie complex double ci = casinh(a); \ 48311970Sngie ATF_CHECK_MSG(crude_equality(creal(ci),creal(b)) && \ 49311970Sngie crude_equality(cimag(ci), cimag(b)), \ 50311970Sngie "for casinh([%g,%g]) = [%g,%g] != [%g,%g]", \ 51311970Sngie creal(a), cimag(a), creal(ci), cimag(ci), creal(b), cimag(b)); \ 52311970Sngie} while (0/*CONSTCOND*/) 53311970Sngie 54311970Sngie 55311970SngieATF_TC(casinh); 56311970SngieATF_TC_HEAD(casinh, tc) 57311970Sngie{ 58311970Sngie atf_tc_set_md_var(tc, "descr","Check casinh family - special values"); 59311970Sngie} 60311970Sngie 61311970SngieATF_TC_BODY(casinh, tc) 62311970Sngie{ 63311970Sngie complex double input; 64311970Sngie complex double result; 65311970Sngie unsigned int i; 66311970Sngie for (i = 0; i < __arraycount(values); i++) { 67311970Sngie RE(input) = values[i].input_re; 68311970Sngie IM(input) = values[i].input_im; 69311970Sngie RE(result) = values[i].result_re; 70311970Sngie IM(result) = values[i].result_im; 71311970Sngie ATF_COMPLEX_EQUAL(input, result); 72311970Sngie } 73311970Sngie} 74311970Sngie 75311970SngieATF_TP_ADD_TCS(tp) 76311970Sngie{ 77311970Sngie 78311970Sngie ATF_TP_ADD_TC(tp, casinh); 79311970Sngie 80311970Sngie return atf_no_error(); 81311970Sngie} 82