1*57718be8SEnji Cooper /* $NetBSD: t_tanh.c,v 1.7 2014/03/03 10:39:08 martin Exp $ */
2*57718be8SEnji Cooper
3*57718be8SEnji Cooper /*-
4*57718be8SEnji Cooper * Copyright (c) 2011 The NetBSD Foundation, Inc.
5*57718be8SEnji Cooper * All rights reserved.
6*57718be8SEnji Cooper *
7*57718be8SEnji Cooper * This code is derived from software contributed to The NetBSD Foundation
8*57718be8SEnji Cooper * by Jukka Ruohonen.
9*57718be8SEnji Cooper *
10*57718be8SEnji Cooper * Redistribution and use in source and binary forms, with or without
11*57718be8SEnji Cooper * modification, are permitted provided that the following conditions
12*57718be8SEnji Cooper * are met:
13*57718be8SEnji Cooper * 1. Redistributions of source code must retain the above copyright
14*57718be8SEnji Cooper * notice, this list of conditions and the following disclaimer.
15*57718be8SEnji Cooper * 2. Redistributions in binary form must reproduce the above copyright
16*57718be8SEnji Cooper * notice, this list of conditions and the following disclaimer in the
17*57718be8SEnji Cooper * documentation and/or other materials provided with the distribution.
18*57718be8SEnji Cooper *
19*57718be8SEnji Cooper * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
20*57718be8SEnji Cooper * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
21*57718be8SEnji Cooper * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
22*57718be8SEnji Cooper * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
23*57718be8SEnji Cooper * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
24*57718be8SEnji Cooper * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
25*57718be8SEnji Cooper * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
26*57718be8SEnji Cooper * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
27*57718be8SEnji Cooper * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
28*57718be8SEnji Cooper * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
29*57718be8SEnji Cooper * POSSIBILITY OF SUCH DAMAGE.
30*57718be8SEnji Cooper */
31*57718be8SEnji Cooper #include <sys/cdefs.h>
32*57718be8SEnji Cooper __RCSID("$NetBSD: t_tanh.c,v 1.7 2014/03/03 10:39:08 martin Exp $");
33*57718be8SEnji Cooper
34*57718be8SEnji Cooper #include <atf-c.h>
35*57718be8SEnji Cooper #include <math.h>
36*57718be8SEnji Cooper
37*57718be8SEnji Cooper /*
38*57718be8SEnji Cooper * tanh(3)
39*57718be8SEnji Cooper */
40*57718be8SEnji Cooper ATF_TC(tanh_nan);
ATF_TC_HEAD(tanh_nan,tc)41*57718be8SEnji Cooper ATF_TC_HEAD(tanh_nan, tc)
42*57718be8SEnji Cooper {
43*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanh(NaN) == NaN");
44*57718be8SEnji Cooper }
45*57718be8SEnji Cooper
ATF_TC_BODY(tanh_nan,tc)46*57718be8SEnji Cooper ATF_TC_BODY(tanh_nan, tc)
47*57718be8SEnji Cooper {
48*57718be8SEnji Cooper const double x = 0.0L / 0.0L;
49*57718be8SEnji Cooper
50*57718be8SEnji Cooper ATF_CHECK(isnan(x) != 0);
51*57718be8SEnji Cooper ATF_CHECK(isnan(tanh(x)) != 0);
52*57718be8SEnji Cooper }
53*57718be8SEnji Cooper
54*57718be8SEnji Cooper ATF_TC(tanh_inf_neg);
ATF_TC_HEAD(tanh_inf_neg,tc)55*57718be8SEnji Cooper ATF_TC_HEAD(tanh_inf_neg, tc)
56*57718be8SEnji Cooper {
57*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanh(-Inf) == -1.0");
58*57718be8SEnji Cooper }
59*57718be8SEnji Cooper
ATF_TC_BODY(tanh_inf_neg,tc)60*57718be8SEnji Cooper ATF_TC_BODY(tanh_inf_neg, tc)
61*57718be8SEnji Cooper {
62*57718be8SEnji Cooper const double x = -1.0L / 0.0L;
63*57718be8SEnji Cooper
64*57718be8SEnji Cooper ATF_CHECK(tanh(x) == -1.0);
65*57718be8SEnji Cooper }
66*57718be8SEnji Cooper
67*57718be8SEnji Cooper ATF_TC(tanh_inf_pos);
ATF_TC_HEAD(tanh_inf_pos,tc)68*57718be8SEnji Cooper ATF_TC_HEAD(tanh_inf_pos, tc)
69*57718be8SEnji Cooper {
70*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanh(+Inf) == +1.0");
71*57718be8SEnji Cooper }
72*57718be8SEnji Cooper
ATF_TC_BODY(tanh_inf_pos,tc)73*57718be8SEnji Cooper ATF_TC_BODY(tanh_inf_pos, tc)
74*57718be8SEnji Cooper {
75*57718be8SEnji Cooper const double x = 1.0L / 0.0L;
76*57718be8SEnji Cooper
77*57718be8SEnji Cooper ATF_CHECK(tanh(x) == 1.0);
78*57718be8SEnji Cooper }
79*57718be8SEnji Cooper
80*57718be8SEnji Cooper ATF_TC(tanh_zero_neg);
ATF_TC_HEAD(tanh_zero_neg,tc)81*57718be8SEnji Cooper ATF_TC_HEAD(tanh_zero_neg, tc)
82*57718be8SEnji Cooper {
83*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanh(-0.0) == -0.0");
84*57718be8SEnji Cooper }
85*57718be8SEnji Cooper
ATF_TC_BODY(tanh_zero_neg,tc)86*57718be8SEnji Cooper ATF_TC_BODY(tanh_zero_neg, tc)
87*57718be8SEnji Cooper {
88*57718be8SEnji Cooper const double x = -0.0L;
89*57718be8SEnji Cooper double y = tanh(x);
90*57718be8SEnji Cooper
91*57718be8SEnji Cooper ATF_CHECK(x == y);
92*57718be8SEnji Cooper ATF_CHECK(signbit(x) != 0);
93*57718be8SEnji Cooper
94*57718be8SEnji Cooper ATF_REQUIRE_MSG(signbit(y) != 0,
95*57718be8SEnji Cooper "compiler bug, waiting for newer gcc import, see PR lib/44057");
96*57718be8SEnji Cooper }
97*57718be8SEnji Cooper
98*57718be8SEnji Cooper ATF_TC(tanh_zero_pos);
ATF_TC_HEAD(tanh_zero_pos,tc)99*57718be8SEnji Cooper ATF_TC_HEAD(tanh_zero_pos, tc)
100*57718be8SEnji Cooper {
101*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanh(+0.0) == +0.0");
102*57718be8SEnji Cooper }
103*57718be8SEnji Cooper
ATF_TC_BODY(tanh_zero_pos,tc)104*57718be8SEnji Cooper ATF_TC_BODY(tanh_zero_pos, tc)
105*57718be8SEnji Cooper {
106*57718be8SEnji Cooper const double x = 0.0L;
107*57718be8SEnji Cooper double y = tanh(x);
108*57718be8SEnji Cooper
109*57718be8SEnji Cooper ATF_CHECK(x == y);
110*57718be8SEnji Cooper ATF_CHECK(signbit(x) == 0);
111*57718be8SEnji Cooper ATF_CHECK(signbit(y) == 0);
112*57718be8SEnji Cooper }
113*57718be8SEnji Cooper
114*57718be8SEnji Cooper /*
115*57718be8SEnji Cooper * tanhf(3)
116*57718be8SEnji Cooper */
117*57718be8SEnji Cooper ATF_TC(tanhf_nan);
ATF_TC_HEAD(tanhf_nan,tc)118*57718be8SEnji Cooper ATF_TC_HEAD(tanhf_nan, tc)
119*57718be8SEnji Cooper {
120*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanhf(NaN) == NaN");
121*57718be8SEnji Cooper }
122*57718be8SEnji Cooper
ATF_TC_BODY(tanhf_nan,tc)123*57718be8SEnji Cooper ATF_TC_BODY(tanhf_nan, tc)
124*57718be8SEnji Cooper {
125*57718be8SEnji Cooper const float x = 0.0L / 0.0L;
126*57718be8SEnji Cooper
127*57718be8SEnji Cooper ATF_CHECK(isnan(x) != 0);
128*57718be8SEnji Cooper ATF_CHECK(isnan(tanhf(x)) != 0);
129*57718be8SEnji Cooper }
130*57718be8SEnji Cooper
131*57718be8SEnji Cooper ATF_TC(tanhf_inf_neg);
ATF_TC_HEAD(tanhf_inf_neg,tc)132*57718be8SEnji Cooper ATF_TC_HEAD(tanhf_inf_neg, tc)
133*57718be8SEnji Cooper {
134*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanhf(-Inf) == -1.0");
135*57718be8SEnji Cooper }
136*57718be8SEnji Cooper
ATF_TC_BODY(tanhf_inf_neg,tc)137*57718be8SEnji Cooper ATF_TC_BODY(tanhf_inf_neg, tc)
138*57718be8SEnji Cooper {
139*57718be8SEnji Cooper const float x = -1.0L / 0.0L;
140*57718be8SEnji Cooper
141*57718be8SEnji Cooper ATF_CHECK(tanhf(x) == -1.0);
142*57718be8SEnji Cooper }
143*57718be8SEnji Cooper
144*57718be8SEnji Cooper ATF_TC(tanhf_inf_pos);
ATF_TC_HEAD(tanhf_inf_pos,tc)145*57718be8SEnji Cooper ATF_TC_HEAD(tanhf_inf_pos, tc)
146*57718be8SEnji Cooper {
147*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanhf(+Inf) == +1.0");
148*57718be8SEnji Cooper }
149*57718be8SEnji Cooper
ATF_TC_BODY(tanhf_inf_pos,tc)150*57718be8SEnji Cooper ATF_TC_BODY(tanhf_inf_pos, tc)
151*57718be8SEnji Cooper {
152*57718be8SEnji Cooper const float x = 1.0L / 0.0L;
153*57718be8SEnji Cooper
154*57718be8SEnji Cooper ATF_CHECK(tanhf(x) == 1.0);
155*57718be8SEnji Cooper }
156*57718be8SEnji Cooper
157*57718be8SEnji Cooper ATF_TC(tanhf_zero_neg);
ATF_TC_HEAD(tanhf_zero_neg,tc)158*57718be8SEnji Cooper ATF_TC_HEAD(tanhf_zero_neg, tc)
159*57718be8SEnji Cooper {
160*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanhf(-0.0) == -0.0");
161*57718be8SEnji Cooper }
162*57718be8SEnji Cooper
ATF_TC_BODY(tanhf_zero_neg,tc)163*57718be8SEnji Cooper ATF_TC_BODY(tanhf_zero_neg, tc)
164*57718be8SEnji Cooper {
165*57718be8SEnji Cooper const float x = -0.0L;
166*57718be8SEnji Cooper float y = tanh(x);
167*57718be8SEnji Cooper
168*57718be8SEnji Cooper ATF_CHECK(x == y);
169*57718be8SEnji Cooper ATF_CHECK(signbit(x) != 0);
170*57718be8SEnji Cooper
171*57718be8SEnji Cooper ATF_REQUIRE_MSG(signbit(y) != 0,
172*57718be8SEnji Cooper "compiler bug, waiting for newer gcc import, see PR lib/44057");
173*57718be8SEnji Cooper }
174*57718be8SEnji Cooper
175*57718be8SEnji Cooper ATF_TC(tanhf_zero_pos);
ATF_TC_HEAD(tanhf_zero_pos,tc)176*57718be8SEnji Cooper ATF_TC_HEAD(tanhf_zero_pos, tc)
177*57718be8SEnji Cooper {
178*57718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test tanhf(+0.0) == +0.0");
179*57718be8SEnji Cooper }
180*57718be8SEnji Cooper
ATF_TC_BODY(tanhf_zero_pos,tc)181*57718be8SEnji Cooper ATF_TC_BODY(tanhf_zero_pos, tc)
182*57718be8SEnji Cooper {
183*57718be8SEnji Cooper const float x = 0.0L;
184*57718be8SEnji Cooper float y = tanhf(x);
185*57718be8SEnji Cooper
186*57718be8SEnji Cooper ATF_CHECK(x == y);
187*57718be8SEnji Cooper ATF_CHECK(signbit(x) == 0);
188*57718be8SEnji Cooper ATF_CHECK(signbit(y) == 0);
189*57718be8SEnji Cooper }
190*57718be8SEnji Cooper
ATF_TP_ADD_TCS(tp)191*57718be8SEnji Cooper ATF_TP_ADD_TCS(tp)
192*57718be8SEnji Cooper {
193*57718be8SEnji Cooper
194*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanh_nan);
195*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanh_inf_neg);
196*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanh_inf_pos);
197*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanh_zero_neg);
198*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanh_zero_pos);
199*57718be8SEnji Cooper
200*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanhf_nan);
201*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanhf_inf_neg);
202*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanhf_inf_pos);
203*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanhf_zero_neg);
204*57718be8SEnji Cooper ATF_TP_ADD_TC(tp, tanhf_zero_pos);
205*57718be8SEnji Cooper
206*57718be8SEnji Cooper return atf_no_error();
207*57718be8SEnji Cooper }
208