1/* $NetBSD: t_fpsetmask.c,v 1.16 2016/03/12 11:55:14 martin Exp $ */ 2 3/*- 4 * Copyright (c) 1995 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 <sys/param.h> 30 31#include <atf-c.h> 32 33#include <stdio.h> 34#include <signal.h> 35#include <float.h> 36#include <setjmp.h> 37#include <stdlib.h> 38#include <string.h> 39 40#include "isqemu.h" 41 42#ifndef _FLOAT_IEEE754 43 44ATF_TC(no_test); 45ATF_TC_HEAD(no_test, tc) 46{ 47 48 atf_tc_set_md_var(tc, "descr", "Dummy test case"); 49} 50 51ATF_TC_BODY(no_test, tc) 52{ 53 54 atf_tc_skip("Test not available on this architecture."); 55} 56 57#else /* defined(_FLOAT_IEEE754) */ 58 59#include <ieeefp.h> 60 61#if __arm__ && !__SOFTFP__ 62 /* 63 * Some NEON fpus do not implement IEEE exception handling, 64 * skip these tests if running on them and compiled for 65 * hard float. 66 */ 67#define FPU_PREREQ() \ 68 if (0 == fpsetmask(fpsetmask(FP_X_INV))) \ 69 atf_tc_skip("FPU does not implement exception handling"); 70#endif 71 72#ifndef FPU_PREREQ 73#define FPU_PREREQ() /* nothing */ 74#endif 75 76void sigfpe(int, siginfo_t *, void *); 77 78volatile sig_atomic_t signal_caught; 79volatile int sicode; 80 81static volatile const float f_one = 1.0; 82static volatile const float f_zero = 0.0; 83static volatile const double d_one = 1.0; 84static volatile const double d_zero = 0.0; 85static volatile const long double ld_one = 1.0; 86static volatile const long double ld_zero = 0.0; 87 88static volatile const float f_huge = FLT_MAX; 89static volatile const float f_tiny = FLT_MIN; 90static volatile const double d_huge = DBL_MAX; 91static volatile const double d_tiny = DBL_MIN; 92static volatile const long double ld_huge = LDBL_MAX; 93static volatile const long double ld_tiny = LDBL_MIN; 94 95static volatile float f_x; 96static volatile double d_x; 97static volatile long double ld_x; 98 99/* trip divide by zero */ 100static void 101f_dz(void) 102{ 103 104 f_x = f_one / f_zero; 105} 106 107static void 108d_dz(void) 109{ 110 111 d_x = d_one / d_zero; 112} 113 114static void 115ld_dz(void) 116{ 117 118 ld_x = ld_one / ld_zero; 119} 120 121/* trip invalid operation */ 122static void 123d_inv(void) 124{ 125 126 d_x = d_zero / d_zero; 127} 128 129static void 130ld_inv(void) 131{ 132 133 ld_x = ld_zero / ld_zero; 134} 135 136static void 137f_inv(void) 138{ 139 140 f_x = f_zero / f_zero; 141} 142 143/* trip overflow */ 144static void 145f_ofl(void) 146{ 147 148 f_x = f_huge * f_huge; 149} 150 151static void 152d_ofl(void) 153{ 154 155 d_x = d_huge * d_huge; 156} 157 158static void 159ld_ofl(void) 160{ 161 162 ld_x = ld_huge * ld_huge; 163} 164 165/* trip underflow */ 166static void 167f_ufl(void) 168{ 169 170 f_x = f_tiny * f_tiny; 171} 172 173static void 174d_ufl(void) 175{ 176 177 d_x = d_tiny * d_tiny; 178} 179 180static void 181ld_ufl(void) 182{ 183 184 ld_x = ld_tiny * ld_tiny; 185} 186 187struct ops { 188 void (*op)(void); 189 fp_except mask; 190 int sicode; 191}; 192 193static const struct ops float_ops[] = { 194 { f_dz, FP_X_DZ, FPE_FLTDIV }, 195 { f_inv, FP_X_INV, FPE_FLTINV }, 196 { f_ofl, FP_X_OFL, FPE_FLTOVF }, 197 { f_ufl, FP_X_UFL, FPE_FLTUND }, 198 { NULL, 0, 0 } 199}; 200 201static const struct ops double_ops[] = { 202 { d_dz, FP_X_DZ, FPE_FLTDIV }, 203 { d_inv, FP_X_INV, FPE_FLTINV }, 204 { d_ofl, FP_X_OFL, FPE_FLTOVF }, 205 { d_ufl, FP_X_UFL, FPE_FLTUND }, 206 { NULL, 0, 0 } 207}; 208 209static const struct ops long_double_ops[] = { 210 { ld_dz, FP_X_DZ, FPE_FLTDIV }, 211 { ld_inv, FP_X_INV, FPE_FLTINV }, 212 { ld_ofl, FP_X_OFL, FPE_FLTOVF }, 213 { ld_ufl, FP_X_UFL, FPE_FLTUND }, 214 { NULL, 0, 0 } 215}; 216 217static sigjmp_buf b; 218 219static void 220fpsetmask_masked(const struct ops *test_ops) 221{ 222 struct sigaction sa; 223 fp_except ex1, ex2; 224 const struct ops *t; 225 226 /* mask all exceptions, clear history */ 227 fpsetmask(0); 228 fpsetsticky(0); 229 230 /* set up signal handler */ 231 sa.sa_sigaction = sigfpe; 232 sigemptyset(&sa.sa_mask); 233 sa.sa_flags = SA_SIGINFO; 234 sigaction(SIGFPE, &sa, 0); 235 signal_caught = 0; 236 237 /* 238 * exceptions masked, check whether "sticky" bits are set correctly 239 */ 240 for (t = test_ops; t->op != NULL; t++) { 241 (*t->op)(); 242 ex1 = fpgetsticky(); 243 ATF_CHECK_EQ(ex1 & t->mask, t->mask); 244 ATF_CHECK_EQ(signal_caught, 0); 245 246 /* check correct fpsetsticky() behaviour */ 247 ex2 = fpsetsticky(0); 248 ATF_CHECK_EQ(fpgetsticky(), 0); 249 ATF_CHECK_EQ(ex1, ex2); 250 } 251} 252 253/* force delayed exceptions to be delivered */ 254#define BARRIER() fpsetmask(0); f_x = f_one * f_one 255 256static void 257fpsetmask_unmasked(const struct ops *test_ops) 258{ 259 struct sigaction sa; 260 int r; 261 const struct ops *volatile t; 262 263 /* mask all exceptions, clear history */ 264 fpsetmask(0); 265 fpsetsticky(0); 266 267 /* set up signal handler */ 268 sa.sa_sigaction = sigfpe; 269 sigemptyset(&sa.sa_mask); 270 sa.sa_flags = SA_SIGINFO; 271 sigaction(SIGFPE, &sa, 0); 272 signal_caught = 0; 273 274 /* 275 * exception unmasked, check SIGFPE delivery and correct siginfo 276 */ 277 for (t = test_ops; t->op != NULL; t++) { 278 fpsetmask(t->mask); 279 r = sigsetjmp(b, 1); 280 if (!r) { 281 (*t->op)(); 282 BARRIER(); 283 } 284 ATF_CHECK_EQ(signal_caught, 1); 285 ATF_CHECK_EQ(sicode, t->sicode); 286 signal_caught = 0; 287 } 288} 289 290void 291sigfpe(int s, siginfo_t *si, void *c) 292{ 293 signal_caught = 1; 294 sicode = si->si_code; 295 siglongjmp(b, 1); 296} 297 298#define TEST(m, t) \ 299 ATF_TC(m##_##t); \ 300 \ 301 ATF_TC_HEAD(m##_##t, tc) \ 302 { \ 303 \ 304 atf_tc_set_md_var(tc, "descr", \ 305 "Test " ___STRING(m) " exceptions for " \ 306 ___STRING(t) "values"); \ 307 } \ 308 \ 309 ATF_TC_BODY(m##_##t, tc) \ 310 { \ 311 \ 312 FPU_PREREQ(); \ 313 \ 314 if (strcmp(MACHINE, "macppc") == 0) \ 315 atf_tc_expect_fail("PR port-macppc/46319"); \ 316 \ 317 if (isQEMU()) \ 318 atf_tc_expect_fail("PR misc/44767"); \ 319 \ 320 m(t##_ops); \ 321 } 322 323TEST(fpsetmask_masked, float) 324TEST(fpsetmask_masked, double) 325TEST(fpsetmask_masked, long_double) 326TEST(fpsetmask_unmasked, float) 327TEST(fpsetmask_unmasked, double) 328TEST(fpsetmask_unmasked, long_double) 329 330ATF_TC(fpsetmask_basic); 331ATF_TC_HEAD(fpsetmask_basic, tc) 332{ 333 atf_tc_set_md_var(tc, "descr", "A basic test of fpsetmask(3)"); 334} 335 336ATF_TC_BODY(fpsetmask_basic, tc) 337{ 338 size_t i; 339 fp_except_t msk, lst[] = { FP_X_INV, FP_X_DZ, FP_X_OFL, FP_X_UFL }; 340 341 FPU_PREREQ(); 342 343 msk = fpgetmask(); 344 for (i = 0; i < __arraycount(lst); i++) { 345 fpsetmask(msk | lst[i]); 346 ATF_CHECK((fpgetmask() & lst[i]) != 0); 347 fpsetmask(msk & ~lst[i]); 348 ATF_CHECK((fpgetmask() & lst[i]) == 0); 349 } 350 351} 352 353#endif /* defined(_FLOAT_IEEE754) */ 354 355ATF_TP_ADD_TCS(tp) 356{ 357 358#ifndef _FLOAT_IEEE754 359 ATF_TP_ADD_TC(tp, no_test); 360#else 361 ATF_TP_ADD_TC(tp, fpsetmask_basic); 362 ATF_TP_ADD_TC(tp, fpsetmask_masked_float); 363 ATF_TP_ADD_TC(tp, fpsetmask_masked_double); 364 ATF_TP_ADD_TC(tp, fpsetmask_masked_long_double); 365 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_float); 366 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_double); 367 ATF_TP_ADD_TC(tp, fpsetmask_unmasked_long_double); 368#endif 369 370 return atf_no_error(); 371} 372