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