1 /* Support code for testing libm functions (driver).
2 Copyright (C) 1997-2025 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>
23 /* Flags set by the including file. */
24 const int flag_test_errno
= TEST_ERRNO
;
25 const int flag_test_exceptions
= TEST_EXCEPTIONS
;
26 const int flag_test_mathvec
= TEST_MATHVEC
;
29 const int snan_tests_arg
= SNAN_TESTS (ARG_FLOAT
);
31 const int snan_tests_arg
= SNAN_TESTS (FLOAT
);
35 #define STR(x) STRX (x)
36 #define STR_FLOAT STR (FLOAT)
37 #define STR_ARG_FLOAT STR (ARG_FLOAT)
38 #define STR_VEC_LEN STR (VEC_LEN)
40 /* Informal description of the functions being tested. */
42 # define TEST_MSG "testing " STR_FLOAT " (vector length " STR_VEC_LEN ")\n"
44 # define TEST_MSG "testing " STR_FLOAT " (argument " STR_ARG_FLOAT ")\n"
46 # define TEST_MSG "testing " STR_FLOAT " (without inline functions)\n"
48 const char test_msg
[] = TEST_MSG
;
50 /* Allow platforms without all rounding modes to test properly,
51 assuming they provide an __FE_UNDEFINED in <bits/fenv.h> which
52 causes fesetround() to return failure. */
54 # define FE_TONEAREST __FE_UNDEFINED
57 # define FE_TOWARDZERO __FE_UNDEFINED
60 # define FE_UPWARD __FE_UNDEFINED
63 # define FE_DOWNWARD __FE_UNDEFINED
66 #define TEST_NAN_PAYLOAD_CANONICALIZE (SNAN_TESTS_PRESERVE_PAYLOAD \
70 const char qtype_str
[] = TYPE_STR
;
72 /* Various constants derived from pi. We must supply them precalculated for
73 accuracy. They are written as a series of postfix operations to keep
74 them concise yet somewhat readable. */
77 #define lit_pi_3_m_4_d LIT (2.356194490192344928846982537459627163)
78 /* pi * 3 / (4 * ln(10)) */
79 #define lit_pi_3_m_4_ln10_m_d LIT (1.023282265381381010614337719073516828)
80 /* pi / (2 * ln(10)) */
81 #define lit_pi_2_ln10_m_d LIT (0.682188176920920673742891812715677885)
82 /* pi / (4 * ln(10)) */
83 #define lit_pi_4_ln10_m_d LIT (0.341094088460460336871445906357838943)
85 #define lit_pi_ln10_d LIT (1.364376353841841347485783625431355770)
87 #define lit_pi_2_d LITM (M_PI_2)
89 #define lit_pi_4_d LITM (M_PI_4)
91 #define lit_pi LITM (M_PI)
93 /* Other useful constants. */
96 #define lit_e LITM (M_E)
98 #define plus_zero LIT (0.0)
99 #define minus_zero LIT (-0.0)
100 #define plus_infty FUNC (__builtin_inf) ()
101 #define minus_infty -(FUNC (__builtin_inf) ())
102 #define qnan_value_pl(S) FUNC (__builtin_nan) (S)
103 #define qnan_value qnan_value_pl ("")
104 #define snan_value_pl(S) FUNC (__builtin_nans) (S)
105 #define snan_value snan_value_pl ("")
106 #define max_value TYPE_MAX
107 #define min_value TYPE_MIN
108 #define min_subnorm_value TYPE_TRUE_MIN
110 #define arg_plus_zero ARG_LIT (0.0)
111 #define arg_minus_zero ARG_LIT (-0.0)
112 #define arg_plus_infty ARG_FUNC (__builtin_inf) ()
113 #define arg_minus_infty -(ARG_FUNC (__builtin_inf) ())
114 #define arg_qnan_value_pl(S) ARG_FUNC (__builtin_nan) (S)
115 #define arg_qnan_value arg_qnan_value_pl ("")
116 #define arg_snan_value_pl(S) ARG_FUNC (__builtin_nans) (S)
117 #define arg_snan_value arg_snan_value_pl ("")
118 #define arg_max_value ARG_TYPE_MAX
119 #define arg_min_value ARG_TYPE_MIN
120 #define arg_min_subnorm_value ARG_TYPE_TRUE_MIN
122 /* For nexttoward tests. */
123 #define snan_value_ld __builtin_nansl ("")
125 /* For pseudo-normal number tests. */
126 #if TEST_COND_intel96
127 # include <math_ldbl.h>
128 #define pseudo_inf { .parts = { 0x00000000, 0x00000000, 0x7fff }}
129 #define pseudo_zero { .parts = { 0x00000000, 0x00000000, 0x0100 }}
130 #define pseudo_qnan { .parts = { 0x00000001, 0x00000000, 0x7fff }}
131 #define pseudo_snan { .parts = { 0x00000001, 0x40000000, 0x7fff }}
132 #define pseudo_unnormal { .parts = { 0x00000001, 0x40000000, 0x0100 }}
135 /* Structures for each kind of test. */
136 /* Used for both RUN_TEST_LOOP_f_f and RUN_TEST_LOOP_fp_f. */
147 struct test_ff_f_data
157 /* Strictly speaking, a j type argument is one gen-libm-test.py will not
158 attempt to muck with. For now, it is only used to prevent it from
159 mucking up an explicitly long double argument. */
160 struct test_fj_f_data
182 struct test_aa_f_data
185 ARG_FLOAT arg1
, arg2
;
192 struct test_aaa_f_data
195 ARG_FLOAT arg1
, arg2
, arg3
;
203 struct test_fi_f_data
214 struct test_fl_f_data
225 struct test_fL_f_data
236 struct test_if_f_data
247 struct test_fff_f_data
250 FLOAT arg1
, arg2
, arg3
;
257 struct test_fiu_M_data
269 struct test_fiu_U_data
291 /* Used for both RUN_TEST_LOOP_f_f1 and RUN_TEST_LOOP_fI_f1. */
292 struct test_f_f1_data
304 struct test_fF_f1_data
313 FLOAT extra_expected
;
316 struct test_ffI_f1_data
338 struct test_cc_c_data
341 FLOAT arg1r
, arg1c
, arg2r
, arg2c
;
348 /* Used for all of RUN_TEST_LOOP_f_i, RUN_TEST_LOOP_f_i_tg,
349 RUN_TEST_LOOP_f_b and RUN_TEST_LOOP_f_b_tg. */
360 /* Used for RUN_TEST_LOOP_f_i_tg_u and RUN_TEST_LOOP_f_b_tg_u. */
361 #if TEST_COND_intel96
362 struct test_j_i_data_u
365 ieee_long_double_shape_type arg
;
373 /* Used for RUN_TEST_LOOP_ff_b, RUN_TEST_LOOP_fpfp_b and
374 RUN_TEST_LOOP_ff_i_tg. */
375 struct test_ff_i_data
401 long long int expected
;
405 struct test_fFF_11_data
413 FLOAT extra1_expected
;
415 FLOAT extra2_expected
;
418 /* Used for both RUN_TEST_LOOP_Ff_b1 and RUN_TEST_LOOP_Ffp_b1. */
419 struct test_Ff_b1_data
428 FLOAT extra_expected
;
432 /* Set the rounding mode, or restore the saved value. */
433 #define IF_ROUND_INIT_ /* Empty. */
434 #define IF_ROUND_INIT_FE_DOWNWARD \
435 int save_round_mode = fegetround (); \
436 if (ROUNDING_TESTS (FLOAT, FE_DOWNWARD) \
438 && fesetround (FE_DOWNWARD) == 0)
439 #define IF_ROUND_INIT_FE_TONEAREST \
440 int save_round_mode = fegetround (); \
441 if (ROUNDING_TESTS (FLOAT, FE_TONEAREST) \
442 && fesetround (FE_TONEAREST) == 0)
443 #define IF_ROUND_INIT_FE_TOWARDZERO \
444 int save_round_mode = fegetround (); \
445 if (ROUNDING_TESTS (FLOAT, FE_TOWARDZERO) \
447 && fesetround (FE_TOWARDZERO) == 0)
448 #define IF_ROUND_INIT_FE_UPWARD \
449 int save_round_mode = fegetround (); \
450 if (ROUNDING_TESTS (FLOAT, FE_UPWARD) \
452 && fesetround (FE_UPWARD) == 0)
453 #define ROUND_RESTORE_ /* Empty. */
454 #define ROUND_RESTORE_FE_DOWNWARD \
455 fesetround (save_round_mode)
456 #define ROUND_RESTORE_FE_TONEAREST \
457 fesetround (save_round_mode)
458 #define ROUND_RESTORE_FE_TOWARDZERO \
459 fesetround (save_round_mode)
460 #define ROUND_RESTORE_FE_UPWARD \
461 fesetround (save_round_mode)
463 /* Field name to use for a given rounding mode. */
465 #define RM_FE_DOWNWARD rd
466 #define RM_FE_TONEAREST rn
467 #define RM_FE_TOWARDZERO rz
468 #define RM_FE_UPWARD ru
470 /* Common setup for an individual test. */
471 #define COMMON_TEST_SETUP(ARG_STR) \
473 if (asprintf (&test_name, "%s (%s)", this_func, (ARG_STR)) == -1) \
476 /* Setup for a test with an extra output. */
477 #define EXTRA_OUTPUT_TEST_SETUP(ARG_STR, N) \
478 char *extra##N##_name; \
479 if (asprintf (&extra##N##_name, "%s (%s) extra output " #N, \
480 this_func, (ARG_STR)) == -1) \
483 /* Common cleanup after an individual test. */
484 #define COMMON_TEST_CLEANUP \
487 /* Cleanup for a test with an extra output. */
488 #define EXTRA_OUTPUT_TEST_CLEANUP(N) \
489 free (extra##N##_name)
491 /* Run an individual test, including any required setup and checking
492 of results, or loop over all tests in an array. */
493 #define RUN_TEST_1_f(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
496 if (enable_test (EXCEPTIONS)) \
498 COMMON_TEST_SETUP (ARG_STR); \
499 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG), \
500 EXPECTED, EXCEPTIONS); \
501 COMMON_TEST_CLEANUP; \
504 #define RUN_TEST_LOOP_1_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
505 IF_ROUND_INIT_ ## ROUNDING_MODE \
506 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
507 RUN_TEST_1_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
508 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
509 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
510 ROUND_RESTORE_ ## ROUNDING_MODE
511 #define RUN_TEST_LOOP_f_f RUN_TEST_LOOP_1_f
512 #define RUN_TEST_LOOP_a_f RUN_TEST_LOOP_1_f
513 #define RUN_TEST_fp_f(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
516 if (enable_test (EXCEPTIONS)) \
518 COMMON_TEST_SETUP (ARG_STR); \
519 check_float (test_name, FUNC_TEST (FUNC_NAME) (&(ARG)), \
520 EXPECTED, EXCEPTIONS); \
521 COMMON_TEST_CLEANUP; \
524 #define RUN_TEST_LOOP_fp_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
525 IF_ROUND_INIT_ ## ROUNDING_MODE \
526 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
527 RUN_TEST_fp_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
528 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
529 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
530 ROUND_RESTORE_ ## ROUNDING_MODE
531 #define RUN_TEST_2_f(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
534 if (enable_test (EXCEPTIONS)) \
536 COMMON_TEST_SETUP (ARG_STR); \
537 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2), \
538 EXPECTED, EXCEPTIONS); \
539 COMMON_TEST_CLEANUP; \
542 #define RUN_TEST_LOOP_2_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
543 IF_ROUND_INIT_ ## ROUNDING_MODE \
544 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
545 RUN_TEST_2_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
547 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
548 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
549 ROUND_RESTORE_ ## ROUNDING_MODE
550 #define RUN_TEST_ff_f RUN_TEST_2_f
551 #define RUN_TEST_LOOP_ff_f RUN_TEST_LOOP_2_f
552 #define RUN_TEST_LOOP_fj_f RUN_TEST_LOOP_2_f
553 #define RUN_TEST_LOOP_aa_f RUN_TEST_LOOP_2_f
554 #define RUN_TEST_fi_f RUN_TEST_2_f
555 #define RUN_TEST_LOOP_fi_f RUN_TEST_LOOP_2_f
556 #define RUN_TEST_fl_f RUN_TEST_2_f
557 #define RUN_TEST_LOOP_fl_f RUN_TEST_LOOP_2_f
558 #define RUN_TEST_fL_f RUN_TEST_2_f
559 #define RUN_TEST_LOOP_fL_f RUN_TEST_LOOP_2_f
560 #define RUN_TEST_if_f RUN_TEST_2_f
561 #define RUN_TEST_LOOP_if_f RUN_TEST_LOOP_2_f
562 #define RUN_TEST_3_f(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
563 EXPECTED, EXCEPTIONS) \
565 if (enable_test (EXCEPTIONS)) \
567 COMMON_TEST_SETUP (ARG_STR); \
568 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
569 EXPECTED, EXCEPTIONS); \
570 COMMON_TEST_CLEANUP; \
573 #define RUN_TEST_LOOP_3_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
574 IF_ROUND_INIT_ ## ROUNDING_MODE \
575 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
576 RUN_TEST_3_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
577 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
578 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
579 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
580 ROUND_RESTORE_ ## ROUNDING_MODE
581 #define RUN_TEST_LOOP_fff_f RUN_TEST_LOOP_3_f
582 #define RUN_TEST_LOOP_aaa_f RUN_TEST_LOOP_3_f
583 #define RUN_TEST_fiu_M(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
584 EXPECTED, EXCEPTIONS) \
586 if (enable_test (EXCEPTIONS)) \
588 COMMON_TEST_SETUP (ARG_STR); \
589 check_intmax_t (test_name, \
590 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
591 EXPECTED, EXCEPTIONS); \
592 COMMON_TEST_CLEANUP; \
595 #define RUN_TEST_LOOP_fiu_M(FUNC_NAME, ARRAY, ROUNDING_MODE) \
596 IF_ROUND_INIT_ ## ROUNDING_MODE \
597 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
598 RUN_TEST_fiu_M ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
599 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
600 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
601 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
602 ROUND_RESTORE_ ## ROUNDING_MODE
603 #define RUN_TEST_fiu_U(ARG_STR, FUNC_NAME, ARG1, ARG2, ARG3, \
604 EXPECTED, EXCEPTIONS) \
606 if (enable_test (EXCEPTIONS)) \
608 COMMON_TEST_SETUP (ARG_STR); \
609 check_uintmax_t (test_name, \
610 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, ARG3), \
611 EXPECTED, EXCEPTIONS); \
612 COMMON_TEST_CLEANUP; \
615 #define RUN_TEST_LOOP_fiu_U(FUNC_NAME, ARRAY, ROUNDING_MODE) \
616 IF_ROUND_INIT_ ## ROUNDING_MODE \
617 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
618 RUN_TEST_fiu_U ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1, \
619 (ARRAY)[i].arg2, (ARRAY)[i].arg3, \
620 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
621 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
622 ROUND_RESTORE_ ## ROUNDING_MODE
623 #define RUN_TEST_c_f(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
626 if (enable_test (EXCEPTIONS)) \
628 COMMON_TEST_SETUP (ARG_STR); \
629 check_float (test_name, \
630 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARG1, ARG2)),\
631 EXPECTED, EXCEPTIONS); \
632 COMMON_TEST_CLEANUP; \
635 #define RUN_TEST_LOOP_c_f(FUNC_NAME, ARRAY, ROUNDING_MODE) \
636 IF_ROUND_INIT_ ## ROUNDING_MODE \
637 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
638 RUN_TEST_c_f ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].argr, \
640 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
641 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
642 ROUND_RESTORE_ ## ROUNDING_MODE
643 #define RUN_TEST_f_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
644 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
647 if (enable_test (EXCEPTIONS)) \
649 COMMON_TEST_SETUP (ARG_STR); \
650 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
651 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
653 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
655 check_int (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
656 EXTRA_OUTPUT_TEST_CLEANUP (1); \
657 COMMON_TEST_CLEANUP; \
660 #define RUN_TEST_LOOP_f_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
661 IF_ROUND_INIT_ ## ROUNDING_MODE \
662 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
663 RUN_TEST_f_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
664 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
665 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
667 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
668 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
669 ROUND_RESTORE_ ## ROUNDING_MODE
670 #define RUN_TEST_fF_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
671 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
674 if (enable_test (EXCEPTIONS)) \
676 COMMON_TEST_SETUP (ARG_STR); \
677 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
678 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA_VAR)), \
679 EXPECTED, EXCEPTIONS); \
680 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
682 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
683 EXTRA_OUTPUT_TEST_CLEANUP (1); \
684 COMMON_TEST_CLEANUP; \
687 #define RUN_TEST_LOOP_fF_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
688 IF_ROUND_INIT_ ## ROUNDING_MODE \
689 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
690 RUN_TEST_fF_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
691 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
692 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
694 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
695 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
696 ROUND_RESTORE_ ## ROUNDING_MODE
697 #define RUN_TEST_fI_f1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
698 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
701 if (enable_test (EXCEPTIONS)) \
703 COMMON_TEST_SETUP (ARG_STR); \
704 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
705 check_float (test_name, FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA_VAR)), \
706 EXPECTED, EXCEPTIONS); \
707 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
709 check_int (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, 0); \
710 EXTRA_OUTPUT_TEST_CLEANUP (1); \
711 COMMON_TEST_CLEANUP; \
714 #define RUN_TEST_LOOP_fI_f1(FUNC_NAME, ARRAY, ROUNDING_MODE, EXTRA_VAR) \
715 IF_ROUND_INIT_ ## ROUNDING_MODE \
716 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
717 RUN_TEST_fI_f1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
718 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
719 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
721 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
722 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
723 ROUND_RESTORE_ ## ROUNDING_MODE
724 #define RUN_TEST_ffI_f1_mod8(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
725 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
728 if (enable_test (EXCEPTIONS)) \
730 COMMON_TEST_SETUP (ARG_STR); \
731 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
732 check_float (test_name, \
733 FUNC_TEST (FUNC_NAME) (ARG1, ARG2, &(EXTRA_VAR)), \
734 EXPECTED, EXCEPTIONS); \
735 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
737 check_int (extra1_name, (EXTRA_VAR) % 8, EXTRA_EXPECTED, 0); \
738 EXTRA_OUTPUT_TEST_CLEANUP (1); \
739 COMMON_TEST_CLEANUP; \
742 #define RUN_TEST_LOOP_ffI_f1_mod8(FUNC_NAME, ARRAY, ROUNDING_MODE, \
744 IF_ROUND_INIT_ ## ROUNDING_MODE \
745 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
746 RUN_TEST_ffI_f1_mod8 ((ARRAY)[i].arg_str, FUNC_NAME, \
747 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
748 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
749 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
751 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
752 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
753 ROUND_RESTORE_ ## ROUNDING_MODE
754 #define RUN_TEST_Ff_b1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
755 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
758 if (enable_test (EXCEPTIONS)) \
760 COMMON_TEST_SETUP (ARG_STR); \
761 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
762 /* Clear any exceptions from comparison involving sNaN \
764 feclearexcept (FE_ALL_EXCEPT); \
765 check_bool (test_name, FUNC_TEST (FUNC_NAME) (&(EXTRA_VAR), \
767 EXPECTED, EXCEPTIONS); \
768 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
770 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, \
771 (EXCEPTIONS) & TEST_NAN_PAYLOAD); \
772 EXTRA_OUTPUT_TEST_CLEANUP (1); \
773 COMMON_TEST_CLEANUP; \
776 #define RUN_TEST_LOOP_Ff_b1(FUNC_NAME, ARRAY, ROUNDING_MODE, \
778 IF_ROUND_INIT_ ## ROUNDING_MODE \
779 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
780 RUN_TEST_Ff_b1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
781 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
782 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
784 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
785 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
786 ROUND_RESTORE_ ## ROUNDING_MODE
787 #define RUN_TEST_Ffp_b1(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
788 EXCEPTIONS, EXTRA_VAR, EXTRA_TEST, \
791 if (enable_test (EXCEPTIONS)) \
793 COMMON_TEST_SETUP (ARG_STR); \
794 (EXTRA_VAR) = (EXTRA_EXPECTED) == 0 ? 1 : 0; \
795 check_bool (test_name, FUNC_TEST (FUNC_NAME) (&(EXTRA_VAR), \
797 EXPECTED, EXCEPTIONS); \
798 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
800 check_float (extra1_name, EXTRA_VAR, EXTRA_EXPECTED, \
801 (EXCEPTIONS) & TEST_NAN_PAYLOAD); \
802 EXTRA_OUTPUT_TEST_CLEANUP (1); \
803 COMMON_TEST_CLEANUP; \
806 #define RUN_TEST_LOOP_Ffp_b1(FUNC_NAME, ARRAY, ROUNDING_MODE, \
808 IF_ROUND_INIT_ ## ROUNDING_MODE \
809 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
810 RUN_TEST_Ffp_b1 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
811 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
812 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
814 (ARRAY)[i].RM_##ROUNDING_MODE.extra_test, \
815 (ARRAY)[i].RM_##ROUNDING_MODE.extra_expected); \
816 ROUND_RESTORE_ ## ROUNDING_MODE
817 #define RUN_TEST_c_c(ARG_STR, FUNC_NAME, ARGR, ARGC, EXPR, EXPC, \
820 if (enable_test (EXCEPTIONS)) \
822 COMMON_TEST_SETUP (ARG_STR); \
823 check_complex (test_name, \
824 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARGR, ARGC)), \
825 BUILD_COMPLEX (EXPR, EXPC), EXCEPTIONS); \
826 COMMON_TEST_CLEANUP; \
829 #define RUN_TEST_LOOP_c_c(FUNC_NAME, ARRAY, ROUNDING_MODE) \
830 IF_ROUND_INIT_ ## ROUNDING_MODE \
831 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
832 RUN_TEST_c_c ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].argr, \
834 (ARRAY)[i].RM_##ROUNDING_MODE.expr, \
835 (ARRAY)[i].RM_##ROUNDING_MODE.expc, \
836 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
837 ROUND_RESTORE_ ## ROUNDING_MODE
838 #define RUN_TEST_cc_c(ARG_STR, FUNC_NAME, ARG1R, ARG1C, ARG2R, ARG2C, \
839 EXPR, EXPC, EXCEPTIONS) \
841 if (enable_test (EXCEPTIONS)) \
843 COMMON_TEST_SETUP (ARG_STR); \
844 check_complex (test_name, \
845 FUNC_TEST (FUNC_NAME) (BUILD_COMPLEX (ARG1R, ARG1C), \
846 BUILD_COMPLEX (ARG2R, ARG2C)), \
847 BUILD_COMPLEX (EXPR, EXPC), EXCEPTIONS); \
848 COMMON_TEST_CLEANUP; \
851 #define RUN_TEST_LOOP_cc_c(FUNC_NAME, ARRAY, ROUNDING_MODE) \
852 IF_ROUND_INIT_ ## ROUNDING_MODE \
853 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
854 RUN_TEST_cc_c ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg1r, \
855 (ARRAY)[i].arg1c, (ARRAY)[i].arg2r, \
857 (ARRAY)[i].RM_##ROUNDING_MODE.expr, \
858 (ARRAY)[i].RM_##ROUNDING_MODE.expc, \
859 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
860 ROUND_RESTORE_ ## ROUNDING_MODE
861 #define RUN_TEST_f_i(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
863 if (enable_test (EXCEPTIONS)) \
865 COMMON_TEST_SETUP (ARG_STR); \
866 check_int (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
868 COMMON_TEST_CLEANUP; \
871 #define RUN_TEST_LOOP_f_i(FUNC_NAME, ARRAY, ROUNDING_MODE) \
872 IF_ROUND_INIT_ ## ROUNDING_MODE \
873 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
874 RUN_TEST_f_i ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
875 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
876 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
877 ROUND_RESTORE_ ## ROUNDING_MODE
878 #define RUN_TEST_f_i_tg(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
881 if (enable_test (EXCEPTIONS)) \
883 COMMON_TEST_SETUP (ARG_STR); \
884 check_int (test_name, FUNC_NAME (ARG), EXPECTED, EXCEPTIONS); \
885 COMMON_TEST_CLEANUP; \
888 #define RUN_TEST_LOOP_f_i_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
889 IF_ROUND_INIT_ ## ROUNDING_MODE \
890 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
891 RUN_TEST_f_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
892 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
893 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
894 ROUND_RESTORE_ ## ROUNDING_MODE
895 #define RUN_TEST_LOOP_j_b_tg_u(FUNC_NAME, ARRAY, ROUNDING_MODE) \
896 IF_ROUND_INIT_ ## ROUNDING_MODE \
897 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
898 RUN_TEST_f_b_tg ((ARRAY)[i].arg_str, FUNC_NAME, \
899 (FLOAT)(ARRAY)[i].arg.value, \
900 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
901 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
902 ROUND_RESTORE_ ## ROUNDING_MODE
903 #define RUN_TEST_LOOP_j_i_tg_u(FUNC_NAME, ARRAY, ROUNDING_MODE) \
904 IF_ROUND_INIT_ ## ROUNDING_MODE \
905 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
906 RUN_TEST_f_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, \
907 (FLOAT)(ARRAY)[i].arg.value, \
908 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
909 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
910 ROUND_RESTORE_ ## ROUNDING_MODE
911 #define RUN_TEST_ff_b(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
914 if (enable_test (EXCEPTIONS)) \
916 COMMON_TEST_SETUP (ARG_STR); \
917 check_bool (test_name, FUNC_TEST (FUNC_NAME) (ARG1, ARG2), \
918 EXPECTED, EXCEPTIONS); \
919 COMMON_TEST_CLEANUP; \
922 #define RUN_TEST_LOOP_ff_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
923 IF_ROUND_INIT_ ## ROUNDING_MODE \
924 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
925 RUN_TEST_ff_b ((ARRAY)[i].arg_str, FUNC_NAME, \
926 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
927 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
928 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
929 ROUND_RESTORE_ ## ROUNDING_MODE
930 #define RUN_TEST_fpfp_b(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
933 if (enable_test (EXCEPTIONS)) \
935 COMMON_TEST_SETUP (ARG_STR); \
936 check_bool (test_name, \
937 FUNC_TEST (FUNC_NAME) (&(ARG1), &(ARG2)), \
938 EXPECTED, EXCEPTIONS); \
939 COMMON_TEST_CLEANUP; \
942 #define RUN_TEST_LOOP_fpfp_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
943 IF_ROUND_INIT_ ## ROUNDING_MODE \
944 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
945 RUN_TEST_fpfp_b ((ARRAY)[i].arg_str, FUNC_NAME, \
946 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
947 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
948 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
949 ROUND_RESTORE_ ## ROUNDING_MODE
950 #define RUN_TEST_ff_i_tg(ARG_STR, FUNC_NAME, ARG1, ARG2, EXPECTED, \
953 if (enable_test (EXCEPTIONS)) \
955 COMMON_TEST_SETUP (ARG_STR); \
956 check_int (test_name, FUNC_NAME (ARG1, ARG2), EXPECTED, \
958 COMMON_TEST_CLEANUP; \
961 #define RUN_TEST_LOOP_ff_i_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
962 IF_ROUND_INIT_ ## ROUNDING_MODE \
963 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
964 RUN_TEST_ff_i_tg ((ARRAY)[i].arg_str, FUNC_NAME, \
965 (ARRAY)[i].arg1, (ARRAY)[i].arg2, \
966 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
967 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
968 ROUND_RESTORE_ ## ROUNDING_MODE
969 #define RUN_TEST_f_b(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
971 if (enable_test (EXCEPTIONS)) \
973 COMMON_TEST_SETUP (ARG_STR); \
974 check_bool (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
976 COMMON_TEST_CLEANUP; \
979 #define RUN_TEST_LOOP_f_b(FUNC_NAME, ARRAY, ROUNDING_MODE) \
980 IF_ROUND_INIT_ ## ROUNDING_MODE \
981 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
982 RUN_TEST_f_b ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
983 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
984 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
985 ROUND_RESTORE_ ## ROUNDING_MODE
986 #define RUN_TEST_f_b_tg(ARG_STR, FUNC_NAME, ARG, EXPECTED, \
989 if (enable_test (EXCEPTIONS)) \
991 COMMON_TEST_SETUP (ARG_STR); \
992 check_bool (test_name, FUNC_NAME (ARG), EXPECTED, EXCEPTIONS); \
993 COMMON_TEST_CLEANUP; \
996 #define RUN_TEST_LOOP_f_b_tg(FUNC_NAME, ARRAY, ROUNDING_MODE) \
997 IF_ROUND_INIT_ ## ROUNDING_MODE \
998 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
999 RUN_TEST_f_b_tg ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
1000 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
1001 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
1002 ROUND_RESTORE_ ## ROUNDING_MODE
1003 #define RUN_TEST_f_l(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
1005 if (enable_test (EXCEPTIONS)) \
1007 COMMON_TEST_SETUP (ARG_STR); \
1008 check_long (test_name, FUNC_TEST (FUNC_NAME) (ARG), EXPECTED, \
1010 COMMON_TEST_CLEANUP; \
1013 #define RUN_TEST_LOOP_f_l(FUNC_NAME, ARRAY, ROUNDING_MODE) \
1014 IF_ROUND_INIT_ ## ROUNDING_MODE \
1015 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
1016 RUN_TEST_f_l ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
1017 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
1018 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
1019 ROUND_RESTORE_ ## ROUNDING_MODE
1020 #define RUN_TEST_f_L(ARG_STR, FUNC_NAME, ARG, EXPECTED, EXCEPTIONS) \
1022 if (enable_test (EXCEPTIONS)) \
1024 COMMON_TEST_SETUP (ARG_STR); \
1025 check_longlong (test_name, FUNC_TEST (FUNC_NAME) (ARG), \
1026 EXPECTED, EXCEPTIONS); \
1027 COMMON_TEST_CLEANUP; \
1030 #define RUN_TEST_LOOP_f_L(FUNC_NAME, ARRAY, ROUNDING_MODE) \
1031 IF_ROUND_INIT_ ## ROUNDING_MODE \
1032 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
1033 RUN_TEST_f_L ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
1034 (ARRAY)[i].RM_##ROUNDING_MODE.expected, \
1035 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions); \
1036 ROUND_RESTORE_ ## ROUNDING_MODE
1037 #define RUN_TEST_fFF_11(ARG_STR, FUNC_NAME, ARG, EXCEPTIONS, \
1038 EXTRA1_VAR, EXTRA1_TEST, \
1039 EXTRA1_EXPECTED, EXTRA2_VAR, \
1040 EXTRA2_TEST, EXTRA2_EXPECTED) \
1042 if (enable_test (EXCEPTIONS)) \
1044 COMMON_TEST_SETUP (ARG_STR); \
1045 FUNC_TEST (FUNC_NAME) (ARG, &(EXTRA1_VAR), &(EXTRA2_VAR)); \
1046 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 1); \
1048 check_float (extra1_name, EXTRA1_VAR, EXTRA1_EXPECTED, \
1050 EXTRA_OUTPUT_TEST_CLEANUP (1); \
1051 EXTRA_OUTPUT_TEST_SETUP (ARG_STR, 2); \
1053 check_float (extra2_name, EXTRA2_VAR, EXTRA2_EXPECTED, 0); \
1054 EXTRA_OUTPUT_TEST_CLEANUP (2); \
1055 COMMON_TEST_CLEANUP; \
1058 #define RUN_TEST_LOOP_fFF_11(FUNC_NAME, ARRAY, ROUNDING_MODE, \
1059 EXTRA1_VAR, EXTRA2_VAR) \
1060 IF_ROUND_INIT_ ## ROUNDING_MODE \
1061 for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
1062 RUN_TEST_fFF_11 ((ARRAY)[i].arg_str, FUNC_NAME, (ARRAY)[i].arg, \
1063 (ARRAY)[i].RM_##ROUNDING_MODE.exceptions, \
1065 (ARRAY)[i].RM_##ROUNDING_MODE.extra1_test, \
1066 (ARRAY)[i].RM_##ROUNDING_MODE.extra1_expected, \
1068 (ARRAY)[i].RM_##ROUNDING_MODE.extra2_test, \
1069 (ARRAY)[i].RM_##ROUNDING_MODE.extra2_expected); \
1070 ROUND_RESTORE_ ## ROUNDING_MODE
1073 # define TEST_SUFF VEC_SUFF
1074 # define TEST_SUFF_STR
1077 # define TEST_SUFF_STR "_" ARG_TYPE_STR
1080 # define TEST_SUFF_STR
1083 #define STR_CONCAT(a, b, c) __STRING (a##b##c)
1084 #define STR_CON3(a, b, c) STR_CONCAT (a, b, c)
1087 # define TEST_COND_any_ibm128 (TEST_COND_ibm128 || TEST_COND_arg_ibm128)
1089 # define TEST_COND_any_ibm128 TEST_COND_ibm128
1092 /* Start and end the tests for a given function. */
1093 #define START(FUN, SUFF, EXACT) \
1095 const char *this_func \
1096 = STR_CON3 (FUN, SUFF, TEST_SUFF) TEST_SUFF_STR; \
1097 init_max_error (this_func, EXACT, TEST_COND_any_ibm128)
1099 check_max_error (this_func)
1100 #define END_COMPLEX \
1101 check_complex_max_error (this_func)
1103 /* Run tests for a given function in all rounding modes. */
1104 #define ALL_RM_TEST(FUNC, EXACT, ARRAY, LOOP_MACRO, END_MACRO, ...) \
1109 START (FUNC,, EXACT); \
1110 LOOP_MACRO (FUNC, ARRAY, , ## __VA_ARGS__); \
1116 START (FUNC, _downward, EXACT); \
1117 LOOP_MACRO (FUNC, ARRAY, FE_DOWNWARD, ## __VA_ARGS__); \
1123 START (FUNC, _towardzero, EXACT); \
1124 LOOP_MACRO (FUNC, ARRAY, FE_TOWARDZERO, ## __VA_ARGS__); \
1130 START (FUNC, _upward, EXACT); \
1131 LOOP_MACRO (FUNC, ARRAY, FE_UPWARD, ## __VA_ARGS__); \
1138 /* Short description of program. */
1139 const char doc
[] = "Math test suite: " TEST_MSG
;
1141 static void do_test (void);
1144 main (int argc
, char **argv
)
1146 libm_test_init (argc
, argv
);
1149 return libm_test_finish ();