11ec3feb6SAlex Richardson /* $NetBSD: t_sinh.c,v 1.7 2018/11/07 03:59:36 riastradh Exp $ */
257718be8SEnji Cooper
357718be8SEnji Cooper /*-
457718be8SEnji Cooper * Copyright (c) 2011 The NetBSD Foundation, Inc.
557718be8SEnji Cooper * All rights reserved.
657718be8SEnji Cooper *
757718be8SEnji Cooper * This code is derived from software contributed to The NetBSD Foundation
857718be8SEnji Cooper * by Jukka Ruohonen.
957718be8SEnji Cooper *
1057718be8SEnji Cooper * Redistribution and use in source and binary forms, with or without
1157718be8SEnji Cooper * modification, are permitted provided that the following conditions
1257718be8SEnji Cooper * are met:
1357718be8SEnji Cooper * 1. Redistributions of source code must retain the above copyright
1457718be8SEnji Cooper * notice, this list of conditions and the following disclaimer.
1557718be8SEnji Cooper * 2. Redistributions in binary form must reproduce the above copyright
1657718be8SEnji Cooper * notice, this list of conditions and the following disclaimer in the
1757718be8SEnji Cooper * documentation and/or other materials provided with the distribution.
1857718be8SEnji Cooper *
1957718be8SEnji Cooper * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS
2057718be8SEnji Cooper * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
2157718be8SEnji Cooper * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
2257718be8SEnji Cooper * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS
2357718be8SEnji Cooper * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
2457718be8SEnji Cooper * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
2557718be8SEnji Cooper * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
2657718be8SEnji Cooper * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
2757718be8SEnji Cooper * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
2857718be8SEnji Cooper * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
2957718be8SEnji Cooper * POSSIBILITY OF SUCH DAMAGE.
3057718be8SEnji Cooper */
3157718be8SEnji Cooper #include <sys/cdefs.h>
321ec3feb6SAlex Richardson __RCSID("$NetBSD: t_sinh.c,v 1.7 2018/11/07 03:59:36 riastradh Exp $");
3357718be8SEnji Cooper
3457718be8SEnji Cooper #include <atf-c.h>
351ec3feb6SAlex Richardson #include <float.h>
3657718be8SEnji Cooper #include <math.h>
3757718be8SEnji Cooper #include <stdio.h>
3857718be8SEnji Cooper
3957718be8SEnji Cooper static const struct {
4057718be8SEnji Cooper double x;
4157718be8SEnji Cooper double y;
4257718be8SEnji Cooper } values[] = {
431ec3feb6SAlex Richardson { -10, -11013.232874703393, },
441ec3feb6SAlex Richardson { -2, -3.626860407847019, },
451ec3feb6SAlex Richardson { -1, -1.1752011936438014, },
461ec3feb6SAlex Richardson { -0.05, -0.050020835937655016, },
471ec3feb6SAlex Richardson { -0.001,-0.0010000001666666751, },
481ec3feb6SAlex Richardson { 0.001, 0.0010000001666666751, },
491ec3feb6SAlex Richardson { 0.05, 0.050020835937655016, },
501ec3feb6SAlex Richardson { 1, 1.1752011936438014, },
511ec3feb6SAlex Richardson { 2, 3.626860407847019, },
521ec3feb6SAlex Richardson { 10, 11013.232874703393, },
5357718be8SEnji Cooper };
5457718be8SEnji Cooper
5557718be8SEnji Cooper /*
5657718be8SEnji Cooper * sinh(3)
5757718be8SEnji Cooper */
5857718be8SEnji Cooper ATF_TC(sinh_inrange);
ATF_TC_HEAD(sinh_inrange,tc)5957718be8SEnji Cooper ATF_TC_HEAD(sinh_inrange, tc)
6057718be8SEnji Cooper {
6157718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "sinh(x) for some values");
6257718be8SEnji Cooper }
6357718be8SEnji Cooper
ATF_TC_BODY(sinh_inrange,tc)6457718be8SEnji Cooper ATF_TC_BODY(sinh_inrange, tc)
6557718be8SEnji Cooper {
661ec3feb6SAlex Richardson const double eps = DBL_EPSILON;
6757718be8SEnji Cooper size_t i;
6857718be8SEnji Cooper
6957718be8SEnji Cooper for (i = 0; i < __arraycount(values); i++) {
701ec3feb6SAlex Richardson double x = values[i].x;
711ec3feb6SAlex Richardson double sinh_x = values[i].y;
7257718be8SEnji Cooper
731ec3feb6SAlex Richardson if (!(fabs((sinh(x) - sinh_x)/sinh_x) <= eps)) {
741ec3feb6SAlex Richardson atf_tc_fail_nonfatal("sinh(%.17g) = %.17g != %.17g\n",
751ec3feb6SAlex Richardson x, sinh(x), sinh_x);
761ec3feb6SAlex Richardson }
7757718be8SEnji Cooper }
7857718be8SEnji Cooper }
7957718be8SEnji Cooper
8057718be8SEnji Cooper ATF_TC(sinh_nan);
ATF_TC_HEAD(sinh_nan,tc)8157718be8SEnji Cooper ATF_TC_HEAD(sinh_nan, tc)
8257718be8SEnji Cooper {
8357718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinh(NaN) == NaN");
8457718be8SEnji Cooper }
8557718be8SEnji Cooper
ATF_TC_BODY(sinh_nan,tc)8657718be8SEnji Cooper ATF_TC_BODY(sinh_nan, tc)
8757718be8SEnji Cooper {
8857718be8SEnji Cooper const double x = 0.0L / 0.0L;
8957718be8SEnji Cooper
9057718be8SEnji Cooper ATF_CHECK(isnan(x) != 0);
9157718be8SEnji Cooper ATF_CHECK(isnan(sinh(x)) != 0);
9257718be8SEnji Cooper }
9357718be8SEnji Cooper
9457718be8SEnji Cooper ATF_TC(sinh_inf_neg);
ATF_TC_HEAD(sinh_inf_neg,tc)9557718be8SEnji Cooper ATF_TC_HEAD(sinh_inf_neg, tc)
9657718be8SEnji Cooper {
9757718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinh(-Inf) == -Inf");
9857718be8SEnji Cooper }
9957718be8SEnji Cooper
ATF_TC_BODY(sinh_inf_neg,tc)10057718be8SEnji Cooper ATF_TC_BODY(sinh_inf_neg, tc)
10157718be8SEnji Cooper {
10257718be8SEnji Cooper const double x = -1.0L / 0.0L;
10357718be8SEnji Cooper double y = sinh(x);
10457718be8SEnji Cooper
10557718be8SEnji Cooper ATF_CHECK(isinf(y) != 0);
10657718be8SEnji Cooper ATF_CHECK(signbit(y) != 0);
10757718be8SEnji Cooper }
10857718be8SEnji Cooper
10957718be8SEnji Cooper ATF_TC(sinh_inf_pos);
ATF_TC_HEAD(sinh_inf_pos,tc)11057718be8SEnji Cooper ATF_TC_HEAD(sinh_inf_pos, tc)
11157718be8SEnji Cooper {
11257718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinh(+Inf) == +Inf");
11357718be8SEnji Cooper }
11457718be8SEnji Cooper
ATF_TC_BODY(sinh_inf_pos,tc)11557718be8SEnji Cooper ATF_TC_BODY(sinh_inf_pos, tc)
11657718be8SEnji Cooper {
11757718be8SEnji Cooper const double x = 1.0L / 0.0L;
11857718be8SEnji Cooper double y = sinh(x);
11957718be8SEnji Cooper
12057718be8SEnji Cooper ATF_CHECK(isinf(y) != 0);
12157718be8SEnji Cooper ATF_CHECK(signbit(y) == 0);
12257718be8SEnji Cooper }
12357718be8SEnji Cooper
12457718be8SEnji Cooper ATF_TC(sinh_zero_neg);
ATF_TC_HEAD(sinh_zero_neg,tc)12557718be8SEnji Cooper ATF_TC_HEAD(sinh_zero_neg, tc)
12657718be8SEnji Cooper {
12757718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinh(-0.0) == -0.0");
12857718be8SEnji Cooper }
12957718be8SEnji Cooper
ATF_TC_BODY(sinh_zero_neg,tc)13057718be8SEnji Cooper ATF_TC_BODY(sinh_zero_neg, tc)
13157718be8SEnji Cooper {
13257718be8SEnji Cooper const double x = -0.0L;
13357718be8SEnji Cooper double y = sinh(x);
13457718be8SEnji Cooper
13557718be8SEnji Cooper if (fabs(y) > 0.0 || signbit(y) == 0)
13657718be8SEnji Cooper atf_tc_fail_nonfatal("sinh(-0.0) != -0.0");
13757718be8SEnji Cooper }
13857718be8SEnji Cooper
13957718be8SEnji Cooper ATF_TC(sinh_zero_pos);
ATF_TC_HEAD(sinh_zero_pos,tc)14057718be8SEnji Cooper ATF_TC_HEAD(sinh_zero_pos, tc)
14157718be8SEnji Cooper {
14257718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinh(+0.0) == +0.0");
14357718be8SEnji Cooper }
14457718be8SEnji Cooper
ATF_TC_BODY(sinh_zero_pos,tc)14557718be8SEnji Cooper ATF_TC_BODY(sinh_zero_pos, tc)
14657718be8SEnji Cooper {
14757718be8SEnji Cooper const double x = 0.0L;
14857718be8SEnji Cooper double y = sinh(x);
14957718be8SEnji Cooper
15057718be8SEnji Cooper if (fabs(y) > 0.0 || signbit(y) != 0)
15157718be8SEnji Cooper atf_tc_fail_nonfatal("sinh(+0.0) != +0.0");
15257718be8SEnji Cooper }
15357718be8SEnji Cooper
15457718be8SEnji Cooper /*
15557718be8SEnji Cooper * sinhf(3)
15657718be8SEnji Cooper */
15757718be8SEnji Cooper ATF_TC(sinhf_inrange);
ATF_TC_HEAD(sinhf_inrange,tc)15857718be8SEnji Cooper ATF_TC_HEAD(sinhf_inrange, tc)
15957718be8SEnji Cooper {
16057718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "sinhf(x) for some values");
16157718be8SEnji Cooper }
16257718be8SEnji Cooper
ATF_TC_BODY(sinhf_inrange,tc)16357718be8SEnji Cooper ATF_TC_BODY(sinhf_inrange, tc)
16457718be8SEnji Cooper {
1651ec3feb6SAlex Richardson const float eps = FLT_EPSILON;
16657718be8SEnji Cooper size_t i;
16757718be8SEnji Cooper
16857718be8SEnji Cooper for (i = 0; i < __arraycount(values); i++) {
1691ec3feb6SAlex Richardson float x = values[i].x;
1701ec3feb6SAlex Richardson float sinh_x = values[i].y;
17157718be8SEnji Cooper
1721ec3feb6SAlex Richardson if (!(fabsf((sinhf(x) - sinh_x)/sinh_x) <= eps)) {
1731ec3feb6SAlex Richardson atf_tc_fail_nonfatal("sinhf(%.8g) = %.8g != %.8g\n",
1741ec3feb6SAlex Richardson (double)x, (double)sinhf(x), (double)sinh_x);
1751ec3feb6SAlex Richardson }
17657718be8SEnji Cooper }
17757718be8SEnji Cooper }
17857718be8SEnji Cooper
17957718be8SEnji Cooper ATF_TC(sinhf_nan);
ATF_TC_HEAD(sinhf_nan,tc)18057718be8SEnji Cooper ATF_TC_HEAD(sinhf_nan, tc)
18157718be8SEnji Cooper {
18257718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinhf(NaN) == NaN");
18357718be8SEnji Cooper }
18457718be8SEnji Cooper
ATF_TC_BODY(sinhf_nan,tc)18557718be8SEnji Cooper ATF_TC_BODY(sinhf_nan, tc)
18657718be8SEnji Cooper {
18757718be8SEnji Cooper const float x = 0.0L / 0.0L;
18857718be8SEnji Cooper
18957718be8SEnji Cooper ATF_CHECK(isnan(x) != 0);
19057718be8SEnji Cooper ATF_CHECK(isnan(sinhf(x)) != 0);
19157718be8SEnji Cooper }
19257718be8SEnji Cooper
19357718be8SEnji Cooper ATF_TC(sinhf_inf_neg);
ATF_TC_HEAD(sinhf_inf_neg,tc)19457718be8SEnji Cooper ATF_TC_HEAD(sinhf_inf_neg, tc)
19557718be8SEnji Cooper {
19657718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinhf(-Inf) == -Inf");
19757718be8SEnji Cooper }
19857718be8SEnji Cooper
ATF_TC_BODY(sinhf_inf_neg,tc)19957718be8SEnji Cooper ATF_TC_BODY(sinhf_inf_neg, tc)
20057718be8SEnji Cooper {
20157718be8SEnji Cooper const float x = -1.0L / 0.0L;
20257718be8SEnji Cooper float y = sinhf(x);
20357718be8SEnji Cooper
20457718be8SEnji Cooper ATF_CHECK(isinf(y) != 0);
20557718be8SEnji Cooper ATF_CHECK(signbit(y) != 0);
20657718be8SEnji Cooper }
20757718be8SEnji Cooper
20857718be8SEnji Cooper ATF_TC(sinhf_inf_pos);
ATF_TC_HEAD(sinhf_inf_pos,tc)20957718be8SEnji Cooper ATF_TC_HEAD(sinhf_inf_pos, tc)
21057718be8SEnji Cooper {
21157718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinhf(+Inf) == +Inf");
21257718be8SEnji Cooper }
21357718be8SEnji Cooper
ATF_TC_BODY(sinhf_inf_pos,tc)21457718be8SEnji Cooper ATF_TC_BODY(sinhf_inf_pos, tc)
21557718be8SEnji Cooper {
21657718be8SEnji Cooper const float x = 1.0L / 0.0L;
21757718be8SEnji Cooper float y = sinhf(x);
21857718be8SEnji Cooper
21957718be8SEnji Cooper ATF_CHECK(isinf(y) != 0);
22057718be8SEnji Cooper ATF_CHECK(signbit(y) == 0);
22157718be8SEnji Cooper }
22257718be8SEnji Cooper
22357718be8SEnji Cooper ATF_TC(sinhf_zero_neg);
ATF_TC_HEAD(sinhf_zero_neg,tc)22457718be8SEnji Cooper ATF_TC_HEAD(sinhf_zero_neg, tc)
22557718be8SEnji Cooper {
22657718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinhf(-0.0) == -0.0");
22757718be8SEnji Cooper }
22857718be8SEnji Cooper
ATF_TC_BODY(sinhf_zero_neg,tc)22957718be8SEnji Cooper ATF_TC_BODY(sinhf_zero_neg, tc)
23057718be8SEnji Cooper {
23157718be8SEnji Cooper const float x = -0.0L;
23257718be8SEnji Cooper float y = sinhf(x);
23357718be8SEnji Cooper
23457718be8SEnji Cooper if (fabsf(y) > 0.0 || signbit(y) == 0)
23557718be8SEnji Cooper atf_tc_fail_nonfatal("sinhf(-0.0) != -0.0");
23657718be8SEnji Cooper }
23757718be8SEnji Cooper
23857718be8SEnji Cooper ATF_TC(sinhf_zero_pos);
ATF_TC_HEAD(sinhf_zero_pos,tc)23957718be8SEnji Cooper ATF_TC_HEAD(sinhf_zero_pos, tc)
24057718be8SEnji Cooper {
24157718be8SEnji Cooper atf_tc_set_md_var(tc, "descr", "Test sinhf(+0.0) == +0.0");
24257718be8SEnji Cooper }
24357718be8SEnji Cooper
ATF_TC_BODY(sinhf_zero_pos,tc)24457718be8SEnji Cooper ATF_TC_BODY(sinhf_zero_pos, tc)
24557718be8SEnji Cooper {
24657718be8SEnji Cooper const float x = 0.0L;
24757718be8SEnji Cooper float y = sinhf(x);
24857718be8SEnji Cooper
24957718be8SEnji Cooper if (fabsf(y) > 0.0 || signbit(y) != 0)
25057718be8SEnji Cooper atf_tc_fail_nonfatal("sinhf(+0.0) != +0.0");
25157718be8SEnji Cooper }
25257718be8SEnji Cooper
ATF_TP_ADD_TCS(tp)25357718be8SEnji Cooper ATF_TP_ADD_TCS(tp)
25457718be8SEnji Cooper {
25557718be8SEnji Cooper
25657718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_inrange);
25757718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_nan);
25857718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_inf_neg);
25957718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_inf_pos);
26057718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_zero_neg);
26157718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinh_zero_pos);
26257718be8SEnji Cooper
26357718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_inrange);
26457718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_nan);
26557718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_inf_neg);
26657718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_inf_pos);
26757718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_zero_neg);
26857718be8SEnji Cooper ATF_TP_ADD_TC(tp, sinhf_zero_pos);
26957718be8SEnji Cooper
27057718be8SEnji Cooper return atf_no_error();
27157718be8SEnji Cooper }
272