xref: /freebsd/contrib/netbsd-tests/lib/libm/t_tanh.c (revision 9268022b74279434ed6300244e3f977e56a8ceb5)
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