1 /* Support code for testing libm functions (driver).
2 Copyright (C) 1997-2020 Free Software Foundation, Inc.
3 This file is part of the GNU C Library.
5 The GNU C Library is free software; you can redistribute it and/or
6 modify it under the terms of the GNU Lesser General Public
7 License as published by the Free Software Foundation; either
8 version 2.1 of the License, or (at your option) any later version.
10 The GNU C Library is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 Lesser General Public License for more details.
15 You should have received a copy of the GNU Lesser General Public
16 License along with the GNU C Library; if not, see
17 <https://www.gnu.org/licenses/>. */
19 #include "libm-test-support.h"
21 #include <math-tests-arch.h>
22 #include <nan-pseudo-number.h>
24 /* Flags set by the including file. */
25 const int flag_test_errno
= TEST_ERRNO
;
26 const int flag_test_exceptions
= TEST_EXCEPTIONS
;
27 const int flag_test_mathvec
= TEST_MATHVEC
;
30 const int snan_tests_arg
= SNAN_TESTS (ARG_FLOAT
);
32 const int snan_tests_arg
= SNAN_TESTS (FLOAT
);
36 #define STR(x) STRX (x)
37 #define STR_FLOAT STR (FLOAT)
38 #define STR_ARG_FLOAT STR (ARG_FLOAT)
39 #define STR_VEC_LEN STR (VEC_LEN)
41 /* Informal description of the functions being tested. */
43 # define TEST_MSG "testing " STR_FLOAT " (vector length " STR_VEC_LEN ")\n"
45 # define TEST_MSG "testing " STR_FLOAT " (argument " STR_ARG_FLOAT ")\n"
47 # define TEST_MSG "testing " STR_FLOAT " (without inline functions)\n"
49 const char test_msg
[] = TEST_MSG
;
51 /* Allow platforms without all rounding modes to test properly,
52 assuming they provide an __FE_UNDEFINED in <bits/fenv.h> which
53 causes fesetround() to return failure. */
55 # define FE_TONEAREST __FE_UNDEFINED
58 # define FE_TOWARDZERO __FE_UNDEFINED
61 # define FE_UPWARD __FE_UNDEFINED
64 # define FE_DOWNWARD __FE_UNDEFINED
67 #define TEST_NAN_PAYLOAD_CANONICALIZE (SNAN_TESTS_PRESERVE_PAYLOAD \
71 const char qtype_str
[] = TYPE_STR
;
73 /* Various constants derived from pi. We must supply them precalculated for
74 accuracy. They are written as a series of postfix operations to keep
75 them concise yet somewhat readable. */
78 #define lit_pi_3_m_4_d LIT (2.356194490192344928846982537459627163)
79 /* pi * 3 / (4 * ln(10)) */
80 #define lit_pi_3_m_4_ln10_m_d LIT (1.023282265381381010614337719073516828)
81 /* pi / (2 * ln(10)) */
82 #define lit_pi_2_ln10_m_d LIT (0.682188176920920673742891812715677885)
83 /* pi / (4 * ln(10)) */
84 #define lit_pi_4_ln10_m_d LIT (0.341094088460460336871445906357838943)
86 #define lit_pi_ln10_d LIT (1.364376353841841347485783625431355770)
88 #define lit_pi_2_d LITM (M_PI_2)
90 #define lit_pi_4_d LITM (M_PI_4)
92 #define lit_pi LITM (M_PI)
94 /* Other useful constants. */
97 #define lit_e LITM (M_E)
99 #define plus_zero LIT (0.0)
100 #define minus_zero LIT (-0.0)
101 #define plus_infty FUNC (__builtin_inf) ()
102 #define minus_infty -(FUNC (__builtin_inf) ())
103 #define qnan_value_pl(S) FUNC (__builtin_nan) (S)
104 #define qnan_value qnan_value_pl ("")
105 #define snan_value_pl(S) FUNC (__builtin_nans) (S)
106 #define snan_value snan_value_pl ("")
107 #define max_value TYPE_MAX
108 #define min_value TYPE_MIN
109 #define min_subnorm_value TYPE_TRUE_MIN
111 #define arg_plus_zero ARG_LIT (0.0)
112 #define arg_minus_zero ARG_LIT (-0.0)
113 #define arg_plus_infty ARG_FUNC (__builtin_inf) ()
114 #define arg_minus_infty -(ARG_FUNC (__builtin_inf) ())
115 #define arg_qnan_value_pl(S) ARG_FUNC (__builtin_nan) (S)
116 #define arg_qnan_value arg_qnan_value_pl ("")
117 #define arg_snan_value_pl(S) ARG_FUNC (__builtin_nans) (S)
118 #define arg_snan_value arg_snan_value_pl ("")
119 #define arg_max_value ARG_TYPE_MAX
120 #define arg_min_value ARG_TYPE_MIN
121 #define arg_min_subnorm_value ARG_TYPE_TRUE_MIN
123 /* For nexttoward tests. */
124 #define snan_value_ld __builtin_nansl ("")
126 /* For pseudo-normal number tests. */
127 #if HANDLE_PSEUDO_NUMBERS
128 # include <math_ldbl.h>
129 #define pseudo_inf { .parts = { 0x00000000, 0x00000000, 0x7fff }}
130 #define pseudo_zero { .parts = { 0x00000000, 0x00000000, 0x0100 }}
131 #define pseudo_qnan { .parts = { 0x00000001, 0x00000000, 0x7fff }}
132 #define pseudo_snan { .parts = { 0x00000001, 0x40000000, 0x7fff }}
133 #define pseudo_unnormal { .parts = { 0x00000001, 0x40000000, 0x0100 }}
136 /* Structures for each kind of test. */
137 /* Used for both RUN_TEST_LOOP_f_f and RUN_TEST_LOOP_fp_f. */
148 struct test_ff_f_data
158 /* Strictly speaking, a j type argument is one gen-libm-test.py will not
159 attempt to muck with. For now, it is only used to prevent it from
160 mucking up an explicitly long double argument. */
161 struct test_fj_f_data
173 struct test_aa_f_data
176 ARG_FLOAT arg1
, arg2
;
184 struct test_fi_f_data
195 struct test_fl_f_data
206 struct test_if_f_data
217 struct test_fff_f_data
220 FLOAT arg1
, arg2
, arg3
;
227 struct test_fiu_M_data
239 struct test_fiu_U_data
261 /* Used for both RUN_TEST_LOOP_f_f1 and RUN_TEST_LOOP_fI_f1. */
262 struct test_f_f1_data
274 struct test_fF_f1_data
283 FLOAT extra_expected
;
286 struct test_ffI_f1_data
308 struct test_cc_c_data
311 FLOAT arg1r
, arg1c
, arg2r
, arg2c
;
318 /* Used for all of RUN_TEST_LOOP_f_i, RUN_TEST_LOOP_f_i_tg,
319 RUN_TEST_LOOP_f_b and RUN_TEST_LOOP_f_b_tg. */
330 /* Used for RUN_TEST_LOOP_f_i_tg_u and RUN_TEST_LOOP_f_b_tg_u. */
331 #if HANDLE_PSEUDO_NUMBERS
332 struct test_f_i_data_u
335 ieee_long_double_shape_type arg
;
343 /* Used for RUN_TEST_LOOP_ff_b, RUN_TEST_LOOP_fpfp_b and
344 RUN_TEST_LOOP_ff_i_tg. */
345 struct test_ff_i_data
371 long long int expected
;
375 struct test_fFF_11_data
383 FLOAT extra1_expected
;
385 FLOAT extra2_expected
;
388 /* Used for both RUN_TEST_LOOP_Ff_b1 and RUN_TEST_LOOP_Ffp_b1. */
389 struct test_Ff_b1_data
398 FLOAT extra_expected
;
402 /* Set the rounding mode, or restore the saved value. */
403 #define IF_ROUND_INIT_ /* Empty. */
404 #define IF_ROUND_INIT_FE_DOWNWARD \
405 int save_round_mode = fegetround (); \
406 if (ROUNDING_TESTS (FLOAT, FE_DOWNWARD) \
408 && fesetround (FE_DOWNWARD) == 0)
409 #define IF_ROUND_INIT_FE_TONEAREST \
410 int save_round_mode = fegetround (); \
411 if (ROUNDING_TESTS (FLOAT, FE_TONEAREST) \
412 && fesetround (FE_TONEAREST) == 0)
413 #define IF_ROUND_INIT_FE_TOWARDZERO \
414 int save_round_mode = fegetround (); \
415 if (ROUNDING_TESTS (FLOAT, FE_TOWARDZERO) \
417 && fesetround (FE_TOWARDZERO) == 0)
418 #define IF_ROUND_INIT_FE_UPWARD \
419 int save_round_mode = fegetround (); \
420 if (ROUNDING_TESTS (FLOAT, FE_UPWARD) \
422 && fesetround (FE_UPWARD) == 0)
423 #define ROUND_RESTORE_ /* Empty. */
424 #define ROUND_RESTORE_FE_DOWNWARD \
425 fesetround (save_round_mode)
426 #define ROUND_RESTORE_FE_TONEAREST \
427 fesetround (save_round_mode)
428 #define ROUND_RESTORE_FE_TOWARDZERO \
429 fesetround (save_round_mode)
430 #define ROUND_RESTORE_FE_UPWARD \
431 fesetround (save_round_mode)
433 /* Field name to use for a given rounding mode. */
435 #define RM_FE_DOWNWARD rd
436 #define RM_FE_TONEAREST rn
437 #define RM_FE_TOWARDZERO rz
438 #define RM_FE_UPWARD ru
440 /* Common setup for an individual test. */
441 #define COMMON_TEST_SETUP(ARG_STR) \
443 if (asprintf (&test_name, "%s (%s)", this_func, (ARG_STR)) == -1) \
446 /* Setup for a test with an extra output. */
447 #define EXTRA_OUTPUT_TEST_SETUP(ARG_STR, N) \
448 char *extra##N##_name; \
449 if (asprintf (&extra##N##_name, "%s (%s) extra output " #N, \
450 this_func, (ARG_STR)) == -1) \
453 /* Common cleanup after an individual test. */
454 #define COMMON_TEST_CLEANUP \
457 /* Cleanup for a test with an extra output. */
458 #define EXTRA_OUTPUT_TEST_CLEANUP(N) \
459 free (extra##N##_name)
461 /* Run an individual test, including any required setup and checking
462 of results, or loop over all tests in an array. */
463 #define RUN_TEST_f_f(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
466 if (enable_test (EXCEPTIONS)) \
468 COMMON_TEST_SETUP (ARG_STR); \
469 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG), \
470 EXPECTED, EXCEPTIONS); \
471 COMMON_TEST_CLEANUP; \
474 #define RUN_TEST_LOOP_f_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
475 IF_ROUND_INIT_ ## ROUNDING_MODE \
476 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
477 RUN_TEST_f_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
478 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
479 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
480 ROUND_RESTORE_ ## ROUNDING_MODE
481 #define RUN_TEST_fp_f(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
484 if (enable_test (EXCEPTIONS)) \
486 COMMON_TEST_SETUP (ARG_STR); \
487 check_float (test_name, FUNC_TEST (FUNC_NAME) (&(ARG)), \
488 EXPECTED, EXCEPTIONS); \
489 COMMON_TEST_CLEANUP; \
492 #define RUN_TEST_LOOP_fp_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
493 IF_ROUND_INIT_ ## ROUNDING_MODE \
494 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
495 RUN_TEST_fp_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
496 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
497 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
498 ROUND_RESTORE_ ## ROUNDING_MODE
499 #define RUN_TEST_2_f(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
502 if (enable_test (EXCEPTIONS)) \
504 COMMON_TEST_SETUP (ARG_STR); \
505 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2), \
506 EXPECTED, EXCEPTIONS); \
507 COMMON_TEST_CLEANUP; \
510 #define RUN_TEST_LOOP_2_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
511 IF_ROUND_INIT_ ## ROUNDING_MODE \
512 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
513 RUN_TEST_2_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
515 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
516 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
517 ROUND_RESTORE_ ## ROUNDING_MODE
518 #define RUN_TEST_ff_f RUN_TEST_2_f
519 #define RUN_TEST_LOOP_ff_f RUN_TEST_LOOP_2_f
520 #define RUN_TEST_LOOP_fj_f RUN_TEST_LOOP_2_f
521 #define RUN_TEST_LOOP_aa_f RUN_TEST_LOOP_2_f
522 #define RUN_TEST_fi_f RUN_TEST_2_f
523 #define RUN_TEST_LOOP_fi_f RUN_TEST_LOOP_2_f
524 #define RUN_TEST_fl_f RUN_TEST_2_f
525 #define RUN_TEST_LOOP_fl_f RUN_TEST_LOOP_2_f
526 #define RUN_TEST_if_f RUN_TEST_2_f
527 #define RUN_TEST_LOOP_if_f RUN_TEST_LOOP_2_f
528 #define RUN_TEST_fff_f(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
529 EXPECTED, EXCEPTIONS) \
531 if (enable_test (EXCEPTIONS)) \
533 COMMON_TEST_SETUP (ARG_STR); \
534 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
535 EXPECTED, EXCEPTIONS); \
536 COMMON_TEST_CLEANUP; \
539 #define RUN_TEST_LOOP_fff_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
540 IF_ROUND_INIT_ ## ROUNDING_MODE \
541 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
542 RUN_TEST_fff_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
543 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
544 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
545 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
546 ROUND_RESTORE_ ## ROUNDING_MODE
547 #define RUN_TEST_fiu_M(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
548 EXPECTED, EXCEPTIONS) \
550 if (enable_test (EXCEPTIONS)) \
552 COMMON_TEST_SETUP (ARG_STR); \
553 check_intmax_t (test_name, \
554 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
555 EXPECTED, EXCEPTIONS); \
556 COMMON_TEST_CLEANUP; \
559 #define RUN_TEST_LOOP_fiu_M(FUNC_NAME, ARRAY, ROUNDING_MODE) \
560 IF_ROUND_INIT_ ## ROUNDING_MODE \
561 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
562 RUN_TEST_fiu_M ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
563 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
564 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
565 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
566 ROUND_RESTORE_ ## ROUNDING_MODE
567 #define RUN_TEST_fiu_U(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
568 EXPECTED, EXCEPTIONS) \
570 if (enable_test (EXCEPTIONS)) \
572 COMMON_TEST_SETUP (ARG_STR); \
573 check_uintmax_t (test_name, \
574 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
575 EXPECTED, EXCEPTIONS); \
576 COMMON_TEST_CLEANUP; \
579 #define RUN_TEST_LOOP_fiu_U(FUNC_NAME, ARRAY, ROUNDING_MODE) \
580 IF_ROUND_INIT_ ## ROUNDING_MODE \
581 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
582 RUN_TEST_fiu_U ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
583 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
584 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
585 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
586 ROUND_RESTORE_ ## ROUNDING_MODE
587 #define RUN_TEST_c_f(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
590 if (enable_test (EXCEPTIONS)) \
592 COMMON_TEST_SETUP (ARG_STR); \
593 check_float (test_name, \
594 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARG1, ARG2)),\
595 EXPECTED, EXCEPTIONS); \
596 COMMON_TEST_CLEANUP; \
599 #define RUN_TEST_LOOP_c_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
600 IF_ROUND_INIT_ ## ROUNDING_MODE \
601 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
602 RUN_TEST_c_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].argr, \
604 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
605 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
606 ROUND_RESTORE_ ## ROUNDING_MODE
607 #define RUN_TEST_f_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
608 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
611 if (enable_test (EXCEPTIONS)) \
613 COMMON_TEST_SETUP (ARG_STR); \
614 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
615 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
617 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
619 check_int (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
620 EXTRA_OUTPUT_TEST_CLEANUP (1); \
621 COMMON_TEST_CLEANUP; \
624 #define RUN_TEST_LOOP_f_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
625 IF_ROUND_INIT_ ## ROUNDING_MODE \
626 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
627 RUN_TEST_f_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
628 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
629 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
631 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
632 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
633 ROUND_RESTORE_ ## ROUNDING_MODE
634 #define RUN_TEST_fF_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
635 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
638 if (enable_test (EXCEPTIONS)) \
640 COMMON_TEST_SETUP (ARG_STR); \
641 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
642 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA_VAR)), \
643 EXPECTED, EXCEPTIONS); \
644 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
646 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
647 EXTRA_OUTPUT_TEST_CLEANUP (1); \
648 COMMON_TEST_CLEANUP; \
651 #define RUN_TEST_LOOP_fF_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
652 IF_ROUND_INIT_ ## ROUNDING_MODE \
653 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
654 RUN_TEST_fF_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
655 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
656 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
658 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
659 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
660 ROUND_RESTORE_ ## ROUNDING_MODE
661 #define RUN_TEST_fI_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
662 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
665 if (enable_test (EXCEPTIONS)) \
667 COMMON_TEST_SETUP (ARG_STR); \
668 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
669 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA_VAR)), \
670 EXPECTED, EXCEPTIONS); \
671 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
673 check_int (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
674 EXTRA_OUTPUT_TEST_CLEANUP (1); \
675 COMMON_TEST_CLEANUP; \
678 #define RUN_TEST_LOOP_fI_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
679 IF_ROUND_INIT_ ## ROUNDING_MODE \
680 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
681 RUN_TEST_fI_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
682 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
683 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
685 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
686 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
687 ROUND_RESTORE_ ## ROUNDING_MODE
688 #define RUN_TEST_ffI_f1_mod8(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
689 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
692 if (enable_test (EXCEPTIONS)) \
694 COMMON_TEST_SETUP (ARG_STR); \
695 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
696 check_float (test_name, \
697 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, &(EXTRA_VAR)), \
698 EXPECTED, EXCEPTIONS); \
699 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
701 check_int (extra1_name, (EXTRA_VAR) % 8, EXTRA_EXPECTED, 0); \
702 EXTRA_OUTPUT_TEST_CLEANUP (1); \
703 COMMON_TEST_CLEANUP; \
706 #define RUN_TEST_LOOP_ffI_f1_mod8(FUNC_NAME, ARRAY, ROUNDING_MODE, \
708 IF_ROUND_INIT_ ## ROUNDING_MODE \
709 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
710 RUN_TEST_ffI_f1_mod8 ((ARRAY)[i].arg_str, FUNC_NAME, \
711 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
712 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
713 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
715 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
716 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
717 ROUND_RESTORE_ ## ROUNDING_MODE
718 #define RUN_TEST_Ff_b1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
719 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
722 if (enable_test (EXCEPTIONS)) \
724 COMMON_TEST_SETUP (ARG_STR); \
725 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
726 /* Clear any exceptions from comparison involving sNaN \
728 feclearexcept (FE_ALL_EXCEPT); \
729 check_bool (test_name, FUNC_TEST (FUNC_NAME) (&(EXTRA_VAR), \
731 EXPECTED, EXCEPTIONS); \
732 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
734 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, \
735 (EXCEPTIONS) & TEST_NAN_PAYLOAD); \
736 EXTRA_OUTPUT_TEST_CLEANUP (1); \
737 COMMON_TEST_CLEANUP; \
740 #define RUN_TEST_LOOP_Ff_b1(FUNC_NAME, ARRAY, ROUNDING_MODE, \
742 IF_ROUND_INIT_ ## ROUNDING_MODE \
743 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
744 RUN_TEST_Ff_b1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
745 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
746 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
748 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
749 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
750 ROUND_RESTORE_ ## ROUNDING_MODE
751 #define RUN_TEST_Ffp_b1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
752 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
755 if (enable_test (EXCEPTIONS)) \
757 COMMON_TEST_SETUP (ARG_STR); \
758 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
759 check_bool (test_name, FUNC_TEST (FUNC_NAME) (&(EXTRA_VAR), \
761 EXPECTED, EXCEPTIONS); \
762 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
764 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, \
765 (EXCEPTIONS) & TEST_NAN_PAYLOAD); \
766 EXTRA_OUTPUT_TEST_CLEANUP (1); \
767 COMMON_TEST_CLEANUP; \
770 #define RUN_TEST_LOOP_Ffp_b1(FUNC_NAME, ARRAY, ROUNDING_MODE, \
772 IF_ROUND_INIT_ ## ROUNDING_MODE \
773 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
774 RUN_TEST_Ffp_b1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
775 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
776 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
778 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
779 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
780 ROUND_RESTORE_ ## ROUNDING_MODE
781 #define RUN_TEST_c_c(ARG_STR, FUNC_NAME, ARGR, ARGC, EXPR, EXPC, \
784 if (enable_test (EXCEPTIONS)) \
786 COMMON_TEST_SETUP (ARG_STR); \
787 check_complex (test_name, \
788 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARGR, ARGC)), \
789 BUILD_COMPLEX (EXPR, EXPC), EXCEPTIONS); \
790 COMMON_TEST_CLEANUP; \
793 #define RUN_TEST_LOOP_c_c(FUNC_NAME, ARRAY, ROUNDING_MODE) \
794 IF_ROUND_INIT_ ## ROUNDING_MODE \
795 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
796 RUN_TEST_c_c ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].argr, \
798 (ARRAY)[i].RM_##ROUNDING_MODE.expr, \
799 (ARRAY)[i].RM_##ROUNDING_MODE.expc, \
800 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
801 ROUND_RESTORE_ ## ROUNDING_MODE
802 #define RUN_TEST_cc_c(ARG_STR, FUNC_NAME, ARG1R, ARG1C, ARG2R, ARG2C, \
803 EXPR, EXPC, EXCEPTIONS) \
805 if (enable_test (EXCEPTIONS)) \
807 COMMON_TEST_SETUP (ARG_STR); \
808 check_complex (test_name, \
809 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARG1R, ARG1C), \
810 BUILD_COMPLEX (ARG2R, ARG2C)), \
811 BUILD_COMPLEX (EXPR, EXPC), EXCEPTIONS); \
812 COMMON_TEST_CLEANUP; \
815 #define RUN_TEST_LOOP_cc_c(FUNC_NAME, ARRAY, ROUNDING_MODE) \
816 IF_ROUND_INIT_ ## ROUNDING_MODE \
817 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
818 RUN_TEST_cc_c ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1r, \
819 (ARRAY)[i].arg1c, (ARRAY)[i].arg2r, \
821 (ARRAY)[i].RM_##ROUNDING_MODE.expr, \
822 (ARRAY)[i].RM_##ROUNDING_MODE.expc, \
823 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
824 ROUND_RESTORE_ ## ROUNDING_MODE
825 #define RUN_TEST_f_i(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
827 if (enable_test (EXCEPTIONS)) \
829 COMMON_TEST_SETUP (ARG_STR); \
830 check_int (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
832 COMMON_TEST_CLEANUP; \
835 #define RUN_TEST_LOOP_f_i(FUNC_NAME, ARRAY, ROUNDING_MODE) \
836 IF_ROUND_INIT_ ## ROUNDING_MODE \
837 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
838 RUN_TEST_f_i ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
839 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
840 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
841 ROUND_RESTORE_ ## ROUNDING_MODE
842 #define RUN_TEST_f_i_tg(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
845 if (enable_test (EXCEPTIONS)) \
847 COMMON_TEST_SETUP (ARG_STR); \
848 check_int (test_name, FUNC_NAME (ARG), EXPECTED, EXCEPTIONS); \
849 COMMON_TEST_CLEANUP; \
852 #define RUN_TEST_LOOP_f_i_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
853 IF_ROUND_INIT_ ## ROUNDING_MODE \
854 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
855 RUN_TEST_f_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
856 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
857 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
858 ROUND_RESTORE_ ## ROUNDING_MODE
859 #define RUN_TEST_LOOP_f_b_tg_u(FUNC_NAME, ARRAY, ROUNDING_MODE) \
860 IF_ROUND_INIT_ ## ROUNDING_MODE \
861 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
862 RUN_TEST_f_b_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg.value, \
863 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
864 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
865 ROUND_RESTORE_ ## ROUNDING_MODE
866 #define RUN_TEST_LOOP_f_i_tg_u(FUNC_NAME, ARRAY, ROUNDING_MODE) \
867 IF_ROUND_INIT_ ## ROUNDING_MODE \
868 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
869 RUN_TEST_f_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg.value, \
870 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
871 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
872 ROUND_RESTORE_ ## ROUNDING_MODE
873 #define RUN_TEST_ff_b(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
876 if (enable_test (EXCEPTIONS)) \
878 COMMON_TEST_SETUP (ARG_STR); \
879 check_bool (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2), \
880 EXPECTED, EXCEPTIONS); \
881 COMMON_TEST_CLEANUP; \
884 #define RUN_TEST_LOOP_ff_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
885 IF_ROUND_INIT_ ## ROUNDING_MODE \
886 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
887 RUN_TEST_ff_b ((ARRAY)[i].arg_str, FUNC_NAME, \
888 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
889 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
890 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
891 ROUND_RESTORE_ ## ROUNDING_MODE
892 #define RUN_TEST_fpfp_b(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
895 if (enable_test (EXCEPTIONS)) \
897 COMMON_TEST_SETUP (ARG_STR); \
898 check_bool (test_name, \
899 FUNC_TEST (FUNC_NAME) (&(ARG1), &(ARG2)), \
900 EXPECTED, EXCEPTIONS); \
901 COMMON_TEST_CLEANUP; \
904 #define RUN_TEST_LOOP_fpfp_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
905 IF_ROUND_INIT_ ## ROUNDING_MODE \
906 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
907 RUN_TEST_fpfp_b ((ARRAY)[i].arg_str, FUNC_NAME, \
908 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
909 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
910 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
911 ROUND_RESTORE_ ## ROUNDING_MODE
912 #define RUN_TEST_ff_i_tg(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
915 if (enable_test (EXCEPTIONS)) \
917 COMMON_TEST_SETUP (ARG_STR); \
918 check_int (test_name, FUNC_NAME (ARG1, ARG2), EXPECTED, \
920 COMMON_TEST_CLEANUP; \
923 #define RUN_TEST_LOOP_ff_i_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
924 IF_ROUND_INIT_ ## ROUNDING_MODE \
925 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
926 RUN_TEST_ff_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, \
927 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
928 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
929 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
930 ROUND_RESTORE_ ## ROUNDING_MODE
931 #define RUN_TEST_f_b(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
933 if (enable_test (EXCEPTIONS)) \
935 COMMON_TEST_SETUP (ARG_STR); \
936 check_bool (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
938 COMMON_TEST_CLEANUP; \
941 #define RUN_TEST_LOOP_f_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
942 IF_ROUND_INIT_ ## ROUNDING_MODE \
943 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
944 RUN_TEST_f_b ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
945 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
946 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
947 ROUND_RESTORE_ ## ROUNDING_MODE
948 #define RUN_TEST_f_b_tg(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
951 if (enable_test (EXCEPTIONS)) \
953 COMMON_TEST_SETUP (ARG_STR); \
954 check_bool (test_name, FUNC_NAME (ARG), EXPECTED, EXCEPTIONS); \
955 COMMON_TEST_CLEANUP; \
958 #define RUN_TEST_LOOP_f_b_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
959 IF_ROUND_INIT_ ## ROUNDING_MODE \
960 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
961 RUN_TEST_f_b_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
962 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
963 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
964 ROUND_RESTORE_ ## ROUNDING_MODE
965 #define RUN_TEST_f_l(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
967 if (enable_test (EXCEPTIONS)) \
969 COMMON_TEST_SETUP (ARG_STR); \
970 check_long (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
972 COMMON_TEST_CLEANUP; \
975 #define RUN_TEST_LOOP_f_l(FUNC_NAME, ARRAY, ROUNDING_MODE) \
976 IF_ROUND_INIT_ ## ROUNDING_MODE \
977 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
978 RUN_TEST_f_l ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
979 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
980 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
981 ROUND_RESTORE_ ## ROUNDING_MODE
982 #define RUN_TEST_f_L(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
984 if (enable_test (EXCEPTIONS)) \
986 COMMON_TEST_SETUP (ARG_STR); \
987 check_longlong (test_name, FUNC_TEST (FUNC_NAME) (ARG), \
988 EXPECTED, EXCEPTIONS); \
989 COMMON_TEST_CLEANUP; \
992 #define RUN_TEST_LOOP_f_L(FUNC_NAME, ARRAY, ROUNDING_MODE) \
993 IF_ROUND_INIT_ ## ROUNDING_MODE \
994 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
995 RUN_TEST_f_L ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
996 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
997 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
998 ROUND_RESTORE_ ## ROUNDING_MODE
999 #define RUN_TEST_fFF_11(ARG_STR, FUNC_NAME, ARG, EXCEPTIONS, \
1000 EXTRA1_VAR, EXTRA1_TEST, \
1001 EXTRA1_EXPECTED, EXTRA2_VAR, \
1002 EXTRA2_TEST, EXTRA2_EXPECTED) \
1004 if (enable_test (EXCEPTIONS)) \
1006 COMMON_TEST_SETUP (ARG_STR); \
1007 FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA1_VAR), &(EXTRA2_VAR)); \
1008 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
1010 check_float (extra1_name, EXTRA1_VAR, EXTRA1_EXPECTED, \
1012 EXTRA_OUTPUT_TEST_CLEANUP (1); \
1013 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 2); \
1015 check_float (extra2_name, EXTRA2_VAR, EXTRA2_EXPECTED, 0); \
1016 EXTRA_OUTPUT_TEST_CLEANUP (2); \
1017 COMMON_TEST_CLEANUP; \
1020 #define RUN_TEST_LOOP_fFF_11(FUNC_NAME, ARRAY, ROUNDING_MODE, \
1021 EXTRA1_VAR, EXTRA2_VAR) \
1022 IF_ROUND_INIT_ ## ROUNDING_MODE \
1023 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
1024 RUN_TEST_fFF_11 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
1025 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
1027 (ARRAY)[i].RM_##ROUNDING_MODE.extra1_test, \
1028 (ARRAY)[i].RM_##ROUNDING_MODE.extra1_expected, \
1030 (ARRAY)[i].RM_##ROUNDING_MODE.extra2_test, \
1031 (ARRAY)[i].RM_##ROUNDING_MODE.extra2_expected); \
1032 ROUND_RESTORE_ ## ROUNDING_MODE
1035 # define TEST_SUFF VEC_SUFF
1036 # define TEST_SUFF_STR
1039 # define TEST_SUFF_STR "_" ARG_TYPE_STR
1042 # define TEST_SUFF_STR
1045 #define STR_CONCAT(a, b, c) __STRING (a##b##c)
1046 #define STR_CON3(a, b, c) STR_CONCAT (a, b, c)
1049 # define TEST_COND_any_ibm128 (TEST_COND_ibm128 || TEST_COND_arg_ibm128)
1051 # define TEST_COND_any_ibm128 TEST_COND_ibm128
1054 /* Start and end the tests for a given function. */
1055 #define START(FUN, SUFF, EXACT) \
1057 const char *this_func \
1058 = STR_CON3 (FUN, SUFF, TEST_SUFF) TEST_SUFF_STR; \
1059 init_max_error (this_func, EXACT, TEST_COND_any_ibm128)
1061 print_max_error (this_func)
1062 #define END_COMPLEX \
1063 print_complex_max_error (this_func)
1065 /* Run tests for a given function in all rounding modes. */
1066 #define ALL_RM_TEST(FUNC, EXACT, ARRAY, LOOP_MACRO, END_MACRO, ...) \
1071 START (FUNC,, EXACT); \
1072 LOOP_MACRO (FUNC, ARRAY, , ## __VA_ARGS__); \
1078 START (FUNC, _downward, EXACT); \
1079 LOOP_MACRO (FUNC, ARRAY, FE_DOWNWARD, ## __VA_ARGS__); \
1085 START (FUNC, _towardzero, EXACT); \
1086 LOOP_MACRO (FUNC, ARRAY, FE_TOWARDZERO, ## __VA_ARGS__); \
1092 START (FUNC, _upward, EXACT); \
1093 LOOP_MACRO (FUNC, ARRAY, FE_UPWARD, ## __VA_ARGS__); \
1100 /* Short description of program. */
1101 const char doc
[] = "Math test suite: " TEST_MSG
;
1103 static void do_test (void);
1106 main (int argc
, char **argv
)
1108 libm_test_init (argc
, argv
);
1111 return libm_test_finish ();