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
18static 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
55ATF_TC(casinh);
56ATF_TC_HEAD(casinh, tc)
57{
58	atf_tc_set_md_var(tc, "descr","Check casinh family - special values");
59}
60
61ATF_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
75ATF_TP_ADD_TCS(tp)
76{
77
78	ATF_TP_ADD_TC(tp, casinh);
79
80	return atf_no_error();
81}
82