1/* $NetBSD: t_sqrt.c,v 1.8 2018/11/07 03:59:36 riastradh 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_sqrt.c,v 1.8 2018/11/07 03:59:36 riastradh Exp $"); 33 34#include <atf-c.h> 35#include <math.h> 36#include <float.h> 37#include <stdio.h> 38 39/* 40 * sqrt(3) 41 */ 42ATF_TC(sqrt_nan); 43ATF_TC_HEAD(sqrt_nan, tc) 44{ 45 atf_tc_set_md_var(tc, "descr", "Test sqrt(NaN) == NaN"); 46} 47 48ATF_TC_BODY(sqrt_nan, tc) 49{ 50 const double x = 0.0L / 0.0L; 51 52 ATF_CHECK(isnan(x) != 0); 53 ATF_CHECK(isnan(sqrt(x)) != 0); 54} 55 56ATF_TC(sqrt_pow); 57ATF_TC_HEAD(sqrt_pow, tc) 58{ 59 atf_tc_set_md_var(tc, "descr", "Test sqrt(3) vs. pow(3)"); 60} 61 62ATF_TC_BODY(sqrt_pow, tc) 63{ 64 const double x[] = { 0.0, 0.005, 1.0, 99.0, 123.123, 9999.9999 }; 65 const double eps = DBL_EPSILON; 66 size_t i; 67 68 for (i = 0; i < __arraycount(x); i++) { 69 double x_sqrt = sqrt(x[i]); 70 double x_pow12 = pow(x[i], 1.0 / 2.0); 71 bool ok; 72 73 if (x[i] == 0) { 74 ok = (x_sqrt == x_pow12); 75 } else { 76 ok = (fabs((x_sqrt - x_pow12)/x_sqrt) <= eps); 77 } 78 79 if (!ok) { 80 atf_tc_fail_nonfatal("sqrt(%.17g) = %.17g != " 81 "pow(%.17g, 1/2) = %.17g\n", 82 x[i], x_sqrt, x[i], x_pow12); 83 } 84 } 85} 86 87ATF_TC(sqrt_inf_neg); 88ATF_TC_HEAD(sqrt_inf_neg, tc) 89{ 90 atf_tc_set_md_var(tc, "descr", "Test sqrt(-Inf) == NaN"); 91} 92 93ATF_TC_BODY(sqrt_inf_neg, tc) 94{ 95 const double x = -1.0L / 0.0L; 96 double y = sqrt(x); 97 98 ATF_CHECK(isnan(y) != 0); 99} 100 101ATF_TC(sqrt_inf_pos); 102ATF_TC_HEAD(sqrt_inf_pos, tc) 103{ 104 atf_tc_set_md_var(tc, "descr", "Test sqrt(+Inf) == +Inf"); 105} 106 107ATF_TC_BODY(sqrt_inf_pos, tc) 108{ 109 const double x = 1.0L / 0.0L; 110 double y = sqrt(x); 111 112 ATF_CHECK(isinf(y) != 0); 113 ATF_CHECK(signbit(y) == 0); 114} 115 116ATF_TC(sqrt_zero_neg); 117ATF_TC_HEAD(sqrt_zero_neg, tc) 118{ 119 atf_tc_set_md_var(tc, "descr", "Test sqrt(-0.0) == -0.0"); 120} 121 122ATF_TC_BODY(sqrt_zero_neg, tc) 123{ 124 const double x = -0.0L; 125 double y = sqrt(x); 126 127 if (fabs(y) > 0.0 || signbit(y) == 0) 128 atf_tc_fail_nonfatal("sqrt(-0.0) != -0.0"); 129} 130 131ATF_TC(sqrt_zero_pos); 132ATF_TC_HEAD(sqrt_zero_pos, tc) 133{ 134 atf_tc_set_md_var(tc, "descr", "Test sqrt(+0.0) == +0.0"); 135} 136 137ATF_TC_BODY(sqrt_zero_pos, tc) 138{ 139 const double x = 0.0L; 140 double y = sqrt(x); 141 142 if (fabs(y) > 0.0 || signbit(y) != 0) 143 atf_tc_fail_nonfatal("sqrt(+0.0) != +0.0"); 144} 145 146/* 147 * sqrtf(3) 148 */ 149ATF_TC(sqrtf_nan); 150ATF_TC_HEAD(sqrtf_nan, tc) 151{ 152 atf_tc_set_md_var(tc, "descr", "Test sqrtf(NaN) == NaN"); 153} 154 155ATF_TC_BODY(sqrtf_nan, tc) 156{ 157 const float x = 0.0L / 0.0L; 158 159 ATF_CHECK(isnan(x) != 0); 160 ATF_CHECK(isnan(sqrtf(x)) != 0); 161} 162 163ATF_TC(sqrtf_powf); 164ATF_TC_HEAD(sqrtf_powf, tc) 165{ 166 atf_tc_set_md_var(tc, "descr", "Test sqrtf(3) vs. powf(3)"); 167} 168 169ATF_TC_BODY(sqrtf_powf, tc) 170{ 171 const float x[] = { 0.0, 0.005, 1.0, 99.0, 123.123, 9999.9999 }; 172 const float eps = FLT_EPSILON; 173 size_t i; 174 175 for (i = 0; i < __arraycount(x); i++) { 176 float x_sqrt = sqrtf(x[i]); 177 float x_pow12 = powf(x[i], 1.0 / 2.0); 178 bool ok; 179 180 if (x[i] == 0) { 181 ok = (x_sqrt == x_pow12); 182 } else { 183 ok = (fabsf((x_sqrt - x_pow12)/x_sqrt) <= eps); 184 } 185 186 if (!ok) { 187 atf_tc_fail_nonfatal("sqrtf(%.8g) = %.8g != " 188 "powf(%.8g, 1/2) = %.8g\n", 189 (double)x[i], (double)x_sqrt, 190 (double)x[i], (double)x_pow12); 191 } 192 } 193} 194 195ATF_TC(sqrtf_inf_neg); 196ATF_TC_HEAD(sqrtf_inf_neg, tc) 197{ 198 atf_tc_set_md_var(tc, "descr", "Test sqrtf(-Inf) == NaN"); 199} 200 201ATF_TC_BODY(sqrtf_inf_neg, tc) 202{ 203 const float x = -1.0L / 0.0L; 204 float y = sqrtf(x); 205 206 ATF_CHECK(isnan(y) != 0); 207} 208 209ATF_TC(sqrtf_inf_pos); 210ATF_TC_HEAD(sqrtf_inf_pos, tc) 211{ 212 atf_tc_set_md_var(tc, "descr", "Test sqrtf(+Inf) == +Inf"); 213} 214 215ATF_TC_BODY(sqrtf_inf_pos, tc) 216{ 217 const float x = 1.0L / 0.0L; 218 float y = sqrtf(x); 219 220 ATF_CHECK(isinf(y) != 0); 221 ATF_CHECK(signbit(y) == 0); 222} 223 224ATF_TC(sqrtf_zero_neg); 225ATF_TC_HEAD(sqrtf_zero_neg, tc) 226{ 227 atf_tc_set_md_var(tc, "descr", "Test sqrtf(-0.0) == -0.0"); 228} 229 230ATF_TC_BODY(sqrtf_zero_neg, tc) 231{ 232 const float x = -0.0L; 233 float y = sqrtf(x); 234 235 if (fabsf(y) > 0.0 || signbit(y) == 0) 236 atf_tc_fail_nonfatal("sqrtf(-0.0) != -0.0"); 237} 238 239ATF_TC(sqrtf_zero_pos); 240ATF_TC_HEAD(sqrtf_zero_pos, tc) 241{ 242 atf_tc_set_md_var(tc, "descr", "Test sqrtf(+0.0) == +0.0"); 243} 244 245ATF_TC_BODY(sqrtf_zero_pos, tc) 246{ 247 const float x = 0.0L; 248 float y = sqrtf(x); 249 250 if (fabsf(y) > 0.0 || signbit(y) != 0) 251 atf_tc_fail_nonfatal("sqrtf(+0.0) != +0.0"); 252} 253 254/* 255 * sqrtl(3) 256 */ 257ATF_TC(sqrtl_nan); 258ATF_TC_HEAD(sqrtl_nan, tc) 259{ 260 atf_tc_set_md_var(tc, "descr", "Test sqrtl(NaN) == NaN"); 261} 262 263ATF_TC_BODY(sqrtl_nan, tc) 264{ 265 const long double x = 0.0L / 0.0L; 266 267 ATF_CHECK(isnan(x) != 0); 268 ATF_CHECK(isnan(sqrtl(x)) != 0); 269} 270 271ATF_TC(sqrtl_powl); 272ATF_TC_HEAD(sqrtl_powl, tc) 273{ 274 atf_tc_set_md_var(tc, "descr", "Test sqrtl(3) vs. powl(3)"); 275} 276 277ATF_TC_BODY(sqrtl_powl, tc) 278{ 279 const long double x[] = { 0.0, 0.005, 1.0, 99.0, 123.123, 9999.9999 }; 280 const long double eps = DBL_EPSILON; /* XXX powl == pow for now */ 281 size_t i; 282 283 for (i = 0; i < __arraycount(x); i++) { 284 long double x_sqrt = sqrtl(x[i]); 285 long double x_pow12 = powl(x[i], 1.0 / 2.0); 286 bool ok; 287 288 if (x[i] == 0) { 289 ok = (x_sqrt == x_pow12); 290 } else { 291 ok = (fabsl((x_sqrt - x_pow12)/x_sqrt) <= eps); 292 } 293 294 if (!ok) { 295 atf_tc_fail_nonfatal("sqrtl(%.35Lg) = %.35Lg != " 296 "powl(%.35Lg, 1/2) = %.35Lg\n", 297 x[i], x_sqrt, x[i], x_pow12); 298 } 299 } 300} 301 302ATF_TC(sqrtl_inf_neg); 303ATF_TC_HEAD(sqrtl_inf_neg, tc) 304{ 305 atf_tc_set_md_var(tc, "descr", "Test sqrtl(-Inf) == NaN"); 306} 307 308ATF_TC_BODY(sqrtl_inf_neg, tc) 309{ 310 const long double x = -1.0L / 0.0L; 311 long double y = sqrtl(x); 312 313 ATF_CHECK(isnan(y) != 0); 314} 315 316ATF_TC(sqrtl_inf_pos); 317ATF_TC_HEAD(sqrtl_inf_pos, tc) 318{ 319 atf_tc_set_md_var(tc, "descr", "Test sqrtl(+Inf) == +Inf"); 320} 321 322ATF_TC_BODY(sqrtl_inf_pos, tc) 323{ 324 const long double x = 1.0L / 0.0L; 325 long double y = sqrtl(x); 326 327 ATF_CHECK(isinf(y) != 0); 328 ATF_CHECK(signbit(y) == 0); 329} 330 331ATF_TC(sqrtl_zero_neg); 332ATF_TC_HEAD(sqrtl_zero_neg, tc) 333{ 334 atf_tc_set_md_var(tc, "descr", "Test sqrtl(-0.0) == -0.0"); 335} 336 337ATF_TC_BODY(sqrtl_zero_neg, tc) 338{ 339 const long double x = -0.0L; 340 long double y = sqrtl(x); 341 342 if (fabsl(y) > 0.0 || signbit(y) == 0) 343 atf_tc_fail_nonfatal("sqrtl(-0.0) != -0.0"); 344} 345 346ATF_TC(sqrtl_zero_pos); 347ATF_TC_HEAD(sqrtl_zero_pos, tc) 348{ 349 atf_tc_set_md_var(tc, "descr", "Test sqrtl(+0.0) == +0.0"); 350} 351 352ATF_TC_BODY(sqrtl_zero_pos, tc) 353{ 354 const long double x = 0.0L; 355 long double y = sqrtl(x); 356 357 if (fabsl(y) > 0.0 || signbit(y) != 0) 358 atf_tc_fail_nonfatal("sqrtl(+0.0) != +0.0"); 359} 360 361ATF_TP_ADD_TCS(tp) 362{ 363 364 ATF_TP_ADD_TC(tp, sqrt_nan); 365 ATF_TP_ADD_TC(tp, sqrt_pow); 366 ATF_TP_ADD_TC(tp, sqrt_inf_neg); 367 ATF_TP_ADD_TC(tp, sqrt_inf_pos); 368 ATF_TP_ADD_TC(tp, sqrt_zero_neg); 369 ATF_TP_ADD_TC(tp, sqrt_zero_pos); 370 371 ATF_TP_ADD_TC(tp, sqrtf_nan); 372 ATF_TP_ADD_TC(tp, sqrtf_powf); 373 ATF_TP_ADD_TC(tp, sqrtf_inf_neg); 374 ATF_TP_ADD_TC(tp, sqrtf_inf_pos); 375 ATF_TP_ADD_TC(tp, sqrtf_zero_neg); 376 ATF_TP_ADD_TC(tp, sqrtf_zero_pos); 377 378 ATF_TP_ADD_TC(tp, sqrtl_nan); 379 ATF_TP_ADD_TC(tp, sqrtl_powl); 380 ATF_TP_ADD_TC(tp, sqrtl_inf_neg); 381 ATF_TP_ADD_TC(tp, sqrtl_inf_pos); 382 ATF_TP_ADD_TC(tp, sqrtl_zero_neg); 383 ATF_TP_ADD_TC(tp, sqrtl_zero_pos); 384 385 return atf_no_error(); 386} 387