Home | History | Annotate | Line # | Download | only in libm 
      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