t_fpclassify.c revision 1.8
1 /* $NetBSD: t_fpclassify.c,v 1.8 2025/04/08 00:29:09 riastradh Exp $ */ 2 3 /*- 4 * Copyright (c) 2011 The NetBSD Foundation, Inc. 5 * All rights reserved. 6 * 7 * Redistribution and use in source and binary forms, with or without 8 * modification, are permitted provided that the following conditions 9 * are met: 10 * 1. Redistributions of source code must retain the above copyright 11 * notice, this list of conditions and the following disclaimer. 12 * 2. Redistributions in binary form must reproduce the above copyright 13 * notice, this list of conditions and the following disclaimer in the 14 * documentation and/or other materials provided with the distribution. 15 * 16 * THIS SOFTWARE IS PROVIDED BY THE NETBSD FOUNDATION, INC. AND CONTRIBUTORS 17 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED 18 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR 19 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE FOUNDATION OR CONTRIBUTORS 20 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 21 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF 22 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS 23 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN 24 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) 25 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 26 * POSSIBILITY OF SUCH DAMAGE. 27 */ 28 29 #include <atf-c.h> 30 31 #include <float.h> 32 #include <math.h> 33 #include <stdio.h> 34 #include <string.h> 35 36 ATF_TC(fpclassify_float); 37 ATF_TC_HEAD(fpclassify_float, tc) 38 { 39 40 atf_tc_set_md_var(tc, "descr", "Test float operations"); 41 } 42 43 ATF_TC_BODY(fpclassify_float, tc) 44 { 45 float d0, d1, f; 46 int e; 47 48 d0 = FLT_MIN; 49 ATF_CHECK_EQ_MSG(fpclassify(d0), FP_NORMAL, 50 "fpclassify(%a)=%d FP_NORMAL=%d", 51 d0, fpclassify(d0), FP_NORMAL); 52 f = frexpf(d0, &e); 53 ATF_CHECK_EQ_MSG(e, FLT_MIN_EXP, 54 "frexpf(%a) returned normalized %a, exponent %d;" 55 " expected normalized %a, exponent %d", 56 d0, f, e, 0.5, FLT_MIN_EXP); 57 ATF_CHECK_EQ_MSG(f, 0.5, 58 "frexpf(%a) returned normalized %a, exponent %d;" 59 " expected normalized %a, exponent %d", 60 d0, f, e, 0.5, FLT_MIN_EXP); 61 d1 = d0; 62 63 #if __FLT_HAS_DENORM__ 64 /* shift a "1" bit through the mantissa (skip the implicit bit) */ 65 for (int i = 1; i < FLT_MANT_DIG; i++) { 66 float d2, ip; 67 68 d1 /= 2; 69 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_SUBNORMAL, 70 "[%d] fpclassify(%a)=%d FP_SUBNORMAL=%d", 71 i, d1, fpclassify(d1), FP_SUBNORMAL); 72 ATF_CHECK_MSG(d1 > 0 && d1 < d0, 73 "[%d] d1=%a d0=%a", i, d1, d0); 74 75 d2 = ldexpf(d0, -i); 76 ATF_CHECK_EQ_MSG(d2, d1, "[%d] ldexpf(%a, -%d)=%a != %a", 77 i, d0, i, d2, d1); 78 79 d2 = modff(d1, &ip); 80 ATF_CHECK_EQ_MSG(d2, d1, 81 "[%d] modff(%a) returned int %a, frac %a;" 82 " expected int %a, frac %a", 83 i, d1, ip, d2, 0., d1); 84 ATF_CHECK_EQ_MSG(ip, 0, 85 "[%d] modff(%a) returned int %a, frac %a;" 86 " expected int %a, frac %a", 87 i, d1, ip, d2, 0., d1); 88 89 f = frexpf(d1, &e); 90 ATF_CHECK_EQ_MSG(e, FLT_MIN_EXP - i, 91 "[%d] frexpf(%a) returned normalized %a, exponent %d;" 92 " expected normalized %a, exponent %d", 93 i, d1, f, e, 0.5, FLT_MIN_EXP - i); 94 ATF_CHECK_EQ_MSG(f, 0.5, 95 "[%d] frexpf(%a) returned normalized %a, exponent %d;" 96 " expected normalized %a, exponent %d", 97 i, d1, f, e, 0.5, FLT_MIN_EXP - i); 98 } 99 #endif 100 101 d1 /= 2; 102 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_ZERO, 103 "fpclassify(%a)=%d FP_ZERO=%d", 104 d1, fpclassify(d1), FP_ZERO); 105 f = frexpf(d1, &e); 106 ATF_CHECK_EQ_MSG(e, 0, 107 "frexpf(%a) returned normalized %a, exponent %d;" 108 " expected normalized %a, exponent %d", 109 d1, f, e, 0., 0); 110 ATF_CHECK_EQ_MSG(f, 0, 111 "frexpf(%a) returned normalized %a, exponent %d;" 112 " expected normalized %a, exponent %d", 113 d1, f, e, 0., 0); 114 } 115 116 ATF_TC(fpclassify_double); 117 ATF_TC_HEAD(fpclassify_double, tc) 118 { 119 120 atf_tc_set_md_var(tc, "descr", "Test double operations"); 121 } 122 123 ATF_TC_BODY(fpclassify_double, tc) 124 { 125 double d0, d1, f; 126 int e; 127 128 d0 = DBL_MIN; 129 ATF_CHECK_EQ_MSG(fpclassify(d0), FP_NORMAL, 130 "fpclassify(%a)=%d FP_NORMAL=%d", 131 d0, fpclassify(d0), FP_NORMAL); 132 f = frexp(d0, &e); 133 ATF_CHECK_EQ_MSG(e, DBL_MIN_EXP, 134 "frexp(%a) returned normalized %a, exponent %d;" 135 " expected normalized %a, exponent %d", 136 d0, f, e, 0.5, DBL_MIN_EXP); 137 ATF_CHECK_EQ_MSG(f, 0.5, 138 "frexp(%a) returned normalized %a, exponent %d;" 139 " expected normalized %a, exponent %d", 140 d0, f, e, 0.5, DBL_MIN_EXP); 141 d1 = d0; 142 143 #if __DBL_HAS_DENORM__ 144 /* shift a "1" bit through the mantissa (skip the implicit bit) */ 145 for (int i = 1; i < DBL_MANT_DIG; i++) { 146 double d2, ip; 147 148 d1 /= 2; 149 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_SUBNORMAL, 150 "[%d] fpclassify(%a)=%d FP_SUBNORMAL=%d", 151 i, d1, fpclassify(d1), FP_SUBNORMAL); 152 ATF_CHECK_MSG(d1 > 0 && d1 < d0, 153 "[%d] d1=%a d0=%a", i, d1, d0); 154 155 d2 = ldexp(d0, -i); 156 ATF_CHECK_EQ_MSG(d2, d1, "[%d] ldexp(%a, -%d)=%a != %a", 157 i, d0, i, d2, d1); 158 159 d2 = modf(d1, &ip); 160 ATF_CHECK_EQ_MSG(d2, d1, 161 "[%d] modf(%a) returned int %a, frac %a;" 162 " expected int %a, frac %a", 163 i, d1, ip, d2, 0., d1); 164 ATF_CHECK_EQ_MSG(ip, 0, 165 "[%d] modf(%a) returned int %a, frac %a;" 166 " expected int %a, frac %a", 167 i, d1, ip, d2, 0., d1); 168 169 f = frexp(d1, &e); 170 ATF_CHECK_EQ_MSG(e, DBL_MIN_EXP - i, 171 "[%d] frexp(%a) returned normalized %a, exponent %d;" 172 " expected normalized %a, exponent %d", 173 i, d1, f, e, 0.5, DBL_MIN_EXP - i); 174 ATF_CHECK_EQ_MSG(f, 0.5, 175 "[%d] frexp(%a) returned normalized %a, exponent %d;" 176 " expected normalized %a, exponent %d", 177 i, d1, f, e, 0.5, DBL_MIN_EXP - i); 178 } 179 #endif 180 181 d1 /= 2; 182 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_ZERO, 183 "fpclassify(%a)=%d FP_ZERO=%d", 184 d1, fpclassify(d1), FP_ZERO); 185 f = frexp(d1, &e); 186 ATF_CHECK_EQ_MSG(e, 0, 187 "frexp(%a) returned normalized %a, exponent %d;" 188 " expected normalized %a, exponent %d", 189 d1, f, e, 0., 0); 190 ATF_CHECK_EQ_MSG(f, 0, 191 "frexp(%a) returned normalized %a, exponent %d;" 192 " expected normalized %a, exponent %d", 193 d1, f, e, 0., 0); 194 } 195 196 ATF_TC(fpclassify_long_double); 197 ATF_TC_HEAD(fpclassify_long_double, tc) 198 { 199 200 atf_tc_set_md_var(tc, "descr", "Test long double operations"); 201 } 202 203 ATF_TC_BODY(fpclassify_long_double, tc) 204 { 205 long double d0, d1, f; 206 int e; 207 208 d0 = LDBL_MIN; 209 ATF_CHECK_EQ_MSG(fpclassify(d0), FP_NORMAL, 210 "fpclassify(%La)=%d FP_NORMAL=%d", 211 d0, fpclassify(d0), FP_NORMAL); 212 f = frexpl(d0, &e); 213 ATF_CHECK_EQ_MSG(e, LDBL_MIN_EXP, 214 "frexpl(%La) returned normalized %La, exponent %d;" 215 " expected normalized %La, exponent %d", 216 d0, f, e, 0.5L, LDBL_MIN_EXP); 217 ATF_CHECK_EQ_MSG(f, 0.5, 218 "frexpl(%La) returned normalized %La, exponent %d;" 219 " expected normalized %La, exponent %d", 220 d0, f, e, 0.5L, LDBL_MIN_EXP); 221 d1 = d0; 222 223 #if __LDBL_HAS_DENORM__ 224 /* shift a "1" bit through the mantissa (skip the implicit bit) */ 225 for (int i = 1; i < LDBL_MANT_DIG; i++) { 226 long double d2, ip; 227 228 d1 /= 2; 229 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_SUBNORMAL, 230 "[%d] fpclassify(%La)=%d FP_SUBNORMAL=%d", 231 i, d1, fpclassify(d1), FP_SUBNORMAL); 232 ATF_CHECK_MSG(d1 > 0 && d1 < d0, 233 "[%d] d1=%La d0=%La", i, d1, d0); 234 235 d2 = ldexpl(d0, -i); 236 ATF_CHECK_EQ_MSG(d2, d1, "[%d] ldexpl(%La, -%d)=%La != %La", 237 i, d0, i, d2, d1); 238 239 d2 = modfl(d1, &ip); 240 ATF_CHECK_EQ_MSG(d2, d1, 241 "[%d] modfl(%La) returned int %La, frac %La;" 242 " expected int %La, frac %La", 243 i, d1, ip, d2, 0.L, d1); 244 ATF_CHECK_EQ_MSG(ip, 0, 245 "[%d] modfl(%La) returned int %La, frac %La;" 246 " expected int %La, frac %La", 247 i, d1, ip, d2, 0.L, d1); 248 249 f = frexpl(d1, &e); 250 ATF_CHECK_EQ_MSG(e, LDBL_MIN_EXP - i, 251 "[%d] frexpl(%La) returned normalized %La, exponent %d;" 252 " expected normalized %La, exponent %d", 253 i, d1, f, e, 0.5L, LDBL_MIN_EXP - i); 254 ATF_CHECK_EQ_MSG(f, 0.5, 255 "[%d] frexpl(%La) returned normalized %La, exponent %d;" 256 " expected normalized %La, exponent %d", 257 i, d1, f, e, 0.5L, LDBL_MIN_EXP - i); 258 } 259 #endif 260 261 d1 /= 2; 262 ATF_CHECK_EQ_MSG(fpclassify(d1), FP_ZERO, 263 "fpclassify(%La)=%d FP_ZERO=%d", 264 d1, fpclassify(d1), FP_ZERO); 265 f = frexpl(d1, &e); 266 ATF_CHECK_EQ_MSG(e, 0, 267 "frexpl(%La) returned normalized %La, exponent %d;" 268 " expected normalized %La, exponent %d", 269 d1, f, e, 0.L, 0); 270 ATF_CHECK_EQ_MSG(f, 0, 271 "frexpl(%La) returned normalized %La, exponent %d;" 272 " expected normalized %La, exponent %d", 273 d1, f, e, 0.L, 0); 274 } 275 276 ATF_TP_ADD_TCS(tp) 277 { 278 279 ATF_TP_ADD_TC(tp, fpclassify_float); 280 ATF_TP_ADD_TC(tp, fpclassify_double); 281 ATF_TP_ADD_TC(tp, fpclassify_long_double); 282 283 return atf_no_error(); 284 } 285
Indexes created Sat Sep 20 22:09:52 GMT 2025