1/* $NetBSD: t_asin.c,v 1.3 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 32#include <atf-c.h> 33#include <math.h> 34 35static const struct { 36 double x; 37 double y; 38} values[] = { 39 { -1.0, -M_PI / 2, }, 40 { -0.9, -1.119769514998634, }, 41 { -0.5, -M_PI / 6, }, 42 { -0.1, -0.1001674211615598, }, 43 { 0.1, 0.1001674211615598, }, 44 { 0.5, M_PI / 6, }, 45 { 0.9, 1.119769514998634, }, 46 { 1.0, M_PI / 2, }, 47}; 48 49/* 50 * asin(3) 51 */ 52ATF_TC(asin_nan); 53ATF_TC_HEAD(asin_nan, tc) 54{ 55 atf_tc_set_md_var(tc, "descr", "Test asin(NaN) == NaN"); 56} 57 58ATF_TC_BODY(asin_nan, tc) 59{ 60 const double x = 0.0L / 0.0L; 61 62 if (isnan(asin(x)) == 0) 63 atf_tc_fail_nonfatal("asin(NaN) != NaN"); 64} 65 66ATF_TC(asin_inf_neg); 67ATF_TC_HEAD(asin_inf_neg, tc) 68{ 69 atf_tc_set_md_var(tc, "descr", "Test asin(-Inf) == NaN"); 70} 71 72ATF_TC_BODY(asin_inf_neg, tc) 73{ 74 const double x = -1.0L / 0.0L; 75 76 if (isnan(asin(x)) == 0) 77 atf_tc_fail_nonfatal("asin(-Inf) != NaN"); 78} 79 80ATF_TC(asin_inf_pos); 81ATF_TC_HEAD(asin_inf_pos, tc) 82{ 83 atf_tc_set_md_var(tc, "descr", "Test asin(+Inf) == NaN"); 84} 85 86ATF_TC_BODY(asin_inf_pos, tc) 87{ 88 const double x = 1.0L / 0.0L; 89 90 if (isnan(asin(x)) == 0) 91 atf_tc_fail_nonfatal("asin(+Inf) != NaN"); 92} 93 94ATF_TC(asin_range); 95ATF_TC_HEAD(asin_range, tc) 96{ 97 atf_tc_set_md_var(tc, "descr", "Test asin(x) == NaN, x < -1, x > 1"); 98} 99 100ATF_TC_BODY(asin_range, tc) 101{ 102 const double x[] = { -1.1, -1.000000001, 1.1, 1.000000001 }; 103 size_t i; 104 105 for (i = 0; i < __arraycount(x); i++) { 106 107 if (isnan(asin(x[i])) == 0) 108 atf_tc_fail_nonfatal("asin(%f) != NaN", x[i]); 109 } 110} 111 112ATF_TC(asin_inrange); 113ATF_TC_HEAD(asin_inrange, tc) 114{ 115 atf_tc_set_md_var(tc, "descr", "Test asin(x) for some values"); 116} 117 118ATF_TC_BODY(asin_inrange, tc) 119{ 120 const double eps = 1.0e-15; 121 double y; 122 size_t i; 123 124 for (i = 0; i < __arraycount(values); i++) { 125 y = asin(values[i].x); 126 if (fabs(y - values[i].y) > eps) 127 atf_tc_fail_nonfatal("asin(%g) != %g", 128 values[i].x, values[i].y); 129 } 130} 131 132ATF_TC(asin_zero_neg); 133ATF_TC_HEAD(asin_zero_neg, tc) 134{ 135 atf_tc_set_md_var(tc, "descr", "Test asin(-0.0) == -0.0"); 136} 137 138ATF_TC_BODY(asin_zero_neg, tc) 139{ 140 const double x = -0.0L; 141 double y = asin(x); 142 143 if (fabs(y) > 0.0 || signbit(y) == 0) 144 atf_tc_fail_nonfatal("asin(-0.0) != -0.0"); 145} 146 147ATF_TC(asin_zero_pos); 148ATF_TC_HEAD(asin_zero_pos, tc) 149{ 150 atf_tc_set_md_var(tc, "descr", "Test asin(+0.0) == +0.0"); 151} 152 153ATF_TC_BODY(asin_zero_pos, tc) 154{ 155 const double x = 0.0L; 156 double y = asin(x); 157 158 if (fabs(y) > 0.0 || signbit(y) != 0) 159 atf_tc_fail_nonfatal("asin(+0.0) != +0.0"); 160} 161 162/* 163 * asinf(3) 164 */ 165ATF_TC(asinf_nan); 166ATF_TC_HEAD(asinf_nan, tc) 167{ 168 atf_tc_set_md_var(tc, "descr", "Test asinf(NaN) == NaN"); 169} 170 171ATF_TC_BODY(asinf_nan, tc) 172{ 173 const float x = 0.0L / 0.0L; 174 175 if (isnan(asinf(x)) == 0) 176 atf_tc_fail_nonfatal("asinf(NaN) != NaN"); 177} 178 179ATF_TC(asinf_inf_neg); 180ATF_TC_HEAD(asinf_inf_neg, tc) 181{ 182 atf_tc_set_md_var(tc, "descr", "Test asinf(-Inf) == NaN"); 183} 184 185ATF_TC_BODY(asinf_inf_neg, tc) 186{ 187 const float x = -1.0L / 0.0L; 188 189 if (isnan(asinf(x)) == 0) 190 atf_tc_fail_nonfatal("asinf(-Inf) != NaN"); 191} 192 193ATF_TC(asinf_inf_pos); 194ATF_TC_HEAD(asinf_inf_pos, tc) 195{ 196 atf_tc_set_md_var(tc, "descr", "Test asinf(+Inf) == NaN"); 197} 198 199ATF_TC_BODY(asinf_inf_pos, tc) 200{ 201 const float x = 1.0L / 0.0L; 202 203 if (isnan(asinf(x)) == 0) 204 atf_tc_fail_nonfatal("asinf(+Inf) != NaN"); 205} 206 207ATF_TC(asinf_range); 208ATF_TC_HEAD(asinf_range, tc) 209{ 210 atf_tc_set_md_var(tc, "descr", "Test asinf(x) == NaN, x < -1, x > 1"); 211} 212 213ATF_TC_BODY(asinf_range, tc) 214{ 215 const float x[] = { -1.1, -1.0000001, 1.1, 1.0000001 }; 216 size_t i; 217 218 for (i = 0; i < __arraycount(x); i++) { 219 220 if (isnan(asinf(x[i])) == 0) 221 atf_tc_fail_nonfatal("asinf(%f) != NaN", x[i]); 222 } 223} 224 225ATF_TC(asinf_inrange); 226ATF_TC_HEAD(asinf_inrange, tc) 227{ 228 atf_tc_set_md_var(tc, "descr", "Test asinf(x) for some values"); 229} 230 231ATF_TC_BODY(asinf_inrange, tc) 232{ 233 const float eps = 1.0e-6; 234 float x; 235 float y; 236 size_t i; 237 238 for (i = 0; i < __arraycount(values); i++) { 239 x = values[i].x; 240 y = values[i].y; 241 if (fabs(asinf(x) - y) > eps) 242 atf_tc_fail_nonfatal("asinf(%g) != %g", x, y); 243 } 244} 245 246ATF_TC(asinf_zero_neg); 247ATF_TC_HEAD(asinf_zero_neg, tc) 248{ 249 atf_tc_set_md_var(tc, "descr", "Test asinf(-0.0) == -0.0"); 250} 251 252ATF_TC_BODY(asinf_zero_neg, tc) 253{ 254 const float x = -0.0L; 255 float y = asinf(x); 256 257 if (fabsf(y) > 0.0 || signbit(y) == 0) 258 atf_tc_fail_nonfatal("asinf(-0.0) != -0.0"); 259} 260 261ATF_TC(asinf_zero_pos); 262ATF_TC_HEAD(asinf_zero_pos, tc) 263{ 264 atf_tc_set_md_var(tc, "descr", "Test asinf(+0.0) == +0.0"); 265} 266 267ATF_TC_BODY(asinf_zero_pos, tc) 268{ 269 const float x = 0.0L; 270 float y = asinf(x); 271 272 if (fabsf(y) > 0.0 || signbit(y) != 0) 273 atf_tc_fail_nonfatal("asinf(+0.0) != +0.0"); 274} 275 276ATF_TP_ADD_TCS(tp) 277{ 278 279 ATF_TP_ADD_TC(tp, asin_nan); 280 ATF_TP_ADD_TC(tp, asin_inf_neg); 281 ATF_TP_ADD_TC(tp, asin_inf_pos); 282 ATF_TP_ADD_TC(tp, asin_range); 283 ATF_TP_ADD_TC(tp, asin_inrange); 284 ATF_TP_ADD_TC(tp, asin_zero_neg); 285 ATF_TP_ADD_TC(tp, asin_zero_pos); 286 287 ATF_TP_ADD_TC(tp, asinf_nan); 288 ATF_TP_ADD_TC(tp, asinf_inf_neg); 289 ATF_TP_ADD_TC(tp, asinf_inf_pos); 290 ATF_TP_ADD_TC(tp, asinf_range); 291 ATF_TP_ADD_TC(tp, asinf_inrange); 292 ATF_TP_ADD_TC(tp, asinf_zero_neg); 293 ATF_TP_ADD_TC(tp, asinf_zero_pos); 294 295 return atf_no_error(); 296} 297