t_erf.c revision 272343
1272343Sngie/* $NetBSD: t_erf.c,v 1.2 2014/03/03 10:39:08 martin Exp $ */ 2272343Sngie 3272343Sngie/*- 4272343Sngie * Copyright (c) 2011 The NetBSD Foundation, Inc. 5272343Sngie * All rights reserved. 6272343Sngie * 7272343Sngie * This code is derived from software contributed to The NetBSD Foundation 8272343Sngie * by Jukka Ruohonen. 9272343Sngie * 10272343Sngie * Redistribution and use in source and binary forms, with or without 11272343Sngie * modification, are permitted provided that the following conditions 12272343Sngie * are met: 13272343Sngie * 1. Redistributions of source code must retain the above copyright 14272343Sngie * notice, this list of conditions and the following disclaimer. 15272343Sngie * 2. Redistributions in binary form must reproduce the above copyright 16272343Sngie * notice, this list of conditions and the following disclaimer in the 17272343Sngie * documentation and/or other materials provided with the distribution. 18272343Sngie * 19272343Sngie * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS 20272343Sngie * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 21272343Sngie * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 22272343Sngie * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS 23272343Sngie * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 24272343Sngie * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 25272343Sngie * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 26272343Sngie * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 27272343Sngie * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 28272343Sngie * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 29272343Sngie * POSSIBILITY OF SUCH DAMAGE. 30272343Sngie */ 31272343Sngie#include <sys/cdefs.h> 32272343Sngie__RCSID("$NetBSD: t_erf.c,v 1.2 2014/03/03 10:39:08 martin Exp $"); 33272343Sngie 34272343Sngie#include <atf-c.h> 35272343Sngie#include <math.h> 36272343Sngie 37272343Sngie/* 38272343Sngie * erf(3) 39272343Sngie */ 40272343SngieATF_TC(erf_nan); 41272343SngieATF_TC_HEAD(erf_nan, tc) 42272343Sngie{ 43272343Sngie atf_tc_set_md_var(tc, "descr", "Test erf(NaN) == NaN"); 44272343Sngie} 45272343Sngie 46272343SngieATF_TC_BODY(erf_nan, tc) 47272343Sngie{ 48272343Sngie const double x = 0.0L / 0.0L; 49272343Sngie 50272343Sngie ATF_CHECK(isnan(erf(x)) != 0); 51272343Sngie} 52272343Sngie 53272343SngieATF_TC(erf_inf_neg); 54272343SngieATF_TC_HEAD(erf_inf_neg, tc) 55272343Sngie{ 56272343Sngie atf_tc_set_md_var(tc, "descr", "Test erf(-Inf) == -1.0"); 57272343Sngie} 58272343Sngie 59272343SngieATF_TC_BODY(erf_inf_neg, tc) 60272343Sngie{ 61272343Sngie const double x = -1.0L / 0.0L; 62272343Sngie 63272343Sngie if (erf(x) != -1.0) 64272343Sngie atf_tc_fail_nonfatal("erf(-Inf) != -1.0"); 65272343Sngie} 66272343Sngie 67272343SngieATF_TC(erf_inf_pos); 68272343SngieATF_TC_HEAD(erf_inf_pos, tc) 69272343Sngie{ 70272343Sngie atf_tc_set_md_var(tc, "descr", "Test erf(+Inf) == 1.0"); 71272343Sngie} 72272343Sngie 73272343SngieATF_TC_BODY(erf_inf_pos, tc) 74272343Sngie{ 75272343Sngie const double x = 1.0L / 0.0L; 76272343Sngie 77272343Sngie if (erf(x) != 1.0) 78272343Sngie atf_tc_fail_nonfatal("erf(+Inf) != 1.0"); 79272343Sngie} 80272343Sngie 81272343SngieATF_TC(erf_zero_neg); 82272343SngieATF_TC_HEAD(erf_zero_neg, tc) 83272343Sngie{ 84272343Sngie atf_tc_set_md_var(tc, "descr", "Test erf(-0.0) == -0.0"); 85272343Sngie} 86272343Sngie 87272343SngieATF_TC_BODY(erf_zero_neg, tc) 88272343Sngie{ 89272343Sngie const double x = -0.0L; 90272343Sngie double y = erf(x); 91272343Sngie 92272343Sngie if (fabs(y) > 0.0 || signbit(y) == 0) 93272343Sngie atf_tc_fail_nonfatal("erf(-0.0) != -0.0"); 94272343Sngie} 95272343Sngie 96272343SngieATF_TC(erf_zero_pos); 97272343SngieATF_TC_HEAD(erf_zero_pos, tc) 98272343Sngie{ 99272343Sngie atf_tc_set_md_var(tc, "descr", "Test erf(+0.0) == +0.0"); 100272343Sngie} 101272343Sngie 102272343SngieATF_TC_BODY(erf_zero_pos, tc) 103272343Sngie{ 104272343Sngie const double x = 0.0L; 105272343Sngie double y = erf(x); 106272343Sngie 107272343Sngie if (fabs(y) > 0.0 || signbit(y) != 0) 108272343Sngie atf_tc_fail_nonfatal("erf(+0.0) != +0.0"); 109272343Sngie} 110272343Sngie 111272343Sngie/* 112272343Sngie * erff(3) 113272343Sngie */ 114272343SngieATF_TC(erff_nan); 115272343SngieATF_TC_HEAD(erff_nan, tc) 116272343Sngie{ 117272343Sngie atf_tc_set_md_var(tc, "descr", "Test erff(NaN) == NaN"); 118272343Sngie} 119272343Sngie 120272343SngieATF_TC_BODY(erff_nan, tc) 121272343Sngie{ 122272343Sngie const float x = 0.0L / 0.0L; 123272343Sngie 124272343Sngie ATF_CHECK(isnan(erff(x)) != 0); 125272343Sngie} 126272343Sngie 127272343SngieATF_TC(erff_inf_neg); 128272343SngieATF_TC_HEAD(erff_inf_neg, tc) 129272343Sngie{ 130272343Sngie atf_tc_set_md_var(tc, "descr", "Test erff(-Inf) == -1.0"); 131272343Sngie} 132272343Sngie 133272343SngieATF_TC_BODY(erff_inf_neg, tc) 134272343Sngie{ 135272343Sngie const float x = -1.0L / 0.0L; 136272343Sngie 137272343Sngie if (erff(x) != -1.0) 138272343Sngie atf_tc_fail_nonfatal("erff(-Inf) != -1.0"); 139272343Sngie} 140272343Sngie 141272343SngieATF_TC(erff_inf_pos); 142272343SngieATF_TC_HEAD(erff_inf_pos, tc) 143272343Sngie{ 144272343Sngie atf_tc_set_md_var(tc, "descr", "Test erff(+Inf) == 1.0"); 145272343Sngie} 146272343Sngie 147272343SngieATF_TC_BODY(erff_inf_pos, tc) 148272343Sngie{ 149272343Sngie const float x = 1.0L / 0.0L; 150272343Sngie 151272343Sngie if (erff(x) != 1.0) 152272343Sngie atf_tc_fail_nonfatal("erff(+Inf) != 1.0"); 153272343Sngie} 154272343Sngie 155272343SngieATF_TC(erff_zero_neg); 156272343SngieATF_TC_HEAD(erff_zero_neg, tc) 157272343Sngie{ 158272343Sngie atf_tc_set_md_var(tc, "descr", "Test erff(-0.0) == -0.0"); 159272343Sngie} 160272343Sngie 161272343SngieATF_TC_BODY(erff_zero_neg, tc) 162272343Sngie{ 163272343Sngie const float x = -0.0L; 164272343Sngie float y = erff(x); 165272343Sngie 166272343Sngie if (fabsf(y) > 0.0 || signbit(y) == 0) 167272343Sngie atf_tc_fail_nonfatal("erff(-0.0) != -0.0"); 168272343Sngie} 169272343Sngie 170272343SngieATF_TC(erff_zero_pos); 171272343SngieATF_TC_HEAD(erff_zero_pos, tc) 172272343Sngie{ 173272343Sngie atf_tc_set_md_var(tc, "descr", "Test erff(+0.0) == +0.0"); 174272343Sngie} 175272343Sngie 176272343SngieATF_TC_BODY(erff_zero_pos, tc) 177272343Sngie{ 178272343Sngie const float x = 0.0L; 179272343Sngie float y = erff(x); 180272343Sngie 181272343Sngie if (fabsf(y) > 0.0 || signbit(y) != 0) 182272343Sngie atf_tc_fail_nonfatal("erff(+0.0) != +0.0"); 183272343Sngie} 184272343Sngie 185272343Sngie/* 186272343Sngie * erfc(3) 187272343Sngie */ 188272343SngieATF_TC(erfc_nan); 189272343SngieATF_TC_HEAD(erfc_nan, tc) 190272343Sngie{ 191272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfc(NaN) == NaN"); 192272343Sngie} 193272343Sngie 194272343SngieATF_TC_BODY(erfc_nan, tc) 195272343Sngie{ 196272343Sngie const double x = 0.0L / 0.0L; 197272343Sngie 198272343Sngie ATF_CHECK(isnan(erfc(x)) != 0); 199272343Sngie} 200272343Sngie 201272343SngieATF_TC(erfc_inf_neg); 202272343SngieATF_TC_HEAD(erfc_inf_neg, tc) 203272343Sngie{ 204272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfc(-Inf) == 2.0"); 205272343Sngie} 206272343Sngie 207272343SngieATF_TC_BODY(erfc_inf_neg, tc) 208272343Sngie{ 209272343Sngie const double x = -1.0L / 0.0L; 210272343Sngie 211272343Sngie if (erfc(x) != 2.0) 212272343Sngie atf_tc_fail_nonfatal("erfc(-Inf) != 2.0"); 213272343Sngie} 214272343Sngie 215272343SngieATF_TC(erfc_inf_pos); 216272343SngieATF_TC_HEAD(erfc_inf_pos, tc) 217272343Sngie{ 218272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfc(+Inf) == +0.0"); 219272343Sngie} 220272343Sngie 221272343SngieATF_TC_BODY(erfc_inf_pos, tc) 222272343Sngie{ 223272343Sngie const double x = 1.0L / 0.0L; 224272343Sngie double y = erfc(x); 225272343Sngie 226272343Sngie if (fabs(y) > 0.0 || signbit(y) != 0) 227272343Sngie atf_tc_fail_nonfatal("erfc(+Inf) != +0.0"); 228272343Sngie} 229272343Sngie 230272343Sngie/* 231272343Sngie * erfcf(3) 232272343Sngie */ 233272343SngieATF_TC(erfcf_nan); 234272343SngieATF_TC_HEAD(erfcf_nan, tc) 235272343Sngie{ 236272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfcf(NaN) == NaN"); 237272343Sngie} 238272343Sngie 239272343SngieATF_TC_BODY(erfcf_nan, tc) 240272343Sngie{ 241272343Sngie const float x = 0.0L / 0.0L; 242272343Sngie 243272343Sngie ATF_CHECK(isnan(erfcf(x)) != 0); 244272343Sngie} 245272343Sngie 246272343SngieATF_TC(erfcf_inf_neg); 247272343SngieATF_TC_HEAD(erfcf_inf_neg, tc) 248272343Sngie{ 249272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfcf(-Inf) == 2.0"); 250272343Sngie} 251272343Sngie 252272343SngieATF_TC_BODY(erfcf_inf_neg, tc) 253272343Sngie{ 254272343Sngie const float x = -1.0L / 0.0L; 255272343Sngie 256272343Sngie if (erfcf(x) != 2.0) 257272343Sngie atf_tc_fail_nonfatal("erfcf(-Inf) != 2.0"); 258272343Sngie} 259272343Sngie 260272343SngieATF_TC(erfcf_inf_pos); 261272343SngieATF_TC_HEAD(erfcf_inf_pos, tc) 262272343Sngie{ 263272343Sngie atf_tc_set_md_var(tc, "descr", "Test erfcf(+Inf) == +0.0"); 264272343Sngie} 265272343Sngie 266272343SngieATF_TC_BODY(erfcf_inf_pos, tc) 267272343Sngie{ 268272343Sngie const float x = 1.0L / 0.0L; 269272343Sngie float y = erfcf(x); 270272343Sngie 271272343Sngie if (fabsf(y) > 0.0 || signbit(y) != 0) 272272343Sngie atf_tc_fail_nonfatal("erfcf(+Inf) != +0.0"); 273272343Sngie} 274272343Sngie 275272343SngieATF_TP_ADD_TCS(tp) 276272343Sngie{ 277272343Sngie 278272343Sngie ATF_TP_ADD_TC(tp, erf_nan); 279272343Sngie ATF_TP_ADD_TC(tp, erf_inf_neg); 280272343Sngie ATF_TP_ADD_TC(tp, erf_inf_pos); 281272343Sngie ATF_TP_ADD_TC(tp, erf_zero_neg); 282272343Sngie ATF_TP_ADD_TC(tp, erf_zero_pos); 283272343Sngie 284272343Sngie ATF_TP_ADD_TC(tp, erff_nan); 285272343Sngie ATF_TP_ADD_TC(tp, erff_inf_neg); 286272343Sngie ATF_TP_ADD_TC(tp, erff_inf_pos); 287272343Sngie ATF_TP_ADD_TC(tp, erff_zero_neg); 288272343Sngie ATF_TP_ADD_TC(tp, erff_zero_pos); 289272343Sngie 290272343Sngie ATF_TP_ADD_TC(tp, erfc_nan); 291272343Sngie ATF_TP_ADD_TC(tp, erfc_inf_neg); 292272343Sngie ATF_TP_ADD_TC(tp, erfc_inf_pos); 293272343Sngie 294272343Sngie ATF_TP_ADD_TC(tp, erfcf_nan); 295272343Sngie ATF_TP_ADD_TC(tp, erfcf_inf_neg); 296272343Sngie ATF_TP_ADD_TC(tp, erfcf_inf_pos); 297272343Sngie 298272343Sngie return atf_no_error(); 299272343Sngie} 300