aboutsummaryrefslogtreecommitdiff
path: root/math/libm-test.inc
diff options
context:
space:
mode:
Diffstat (limited to 'math/libm-test.inc')
-rw-r--r--math/libm-test.inc346
1 files changed, 179 insertions, 167 deletions
diff --git a/math/libm-test.inc b/math/libm-test.inc
index af6d28c..0843482 100644
--- a/math/libm-test.inc
+++ b/math/libm-test.inc
@@ -1001,6 +1001,14 @@ struct test_ff_i_data
int max_ulp;
int exceptions;
};
+struct test_f_l_data
+{
+ const char *test_name;
+ FLOAT arg;
+ long int expected;
+ long int max_ulp;
+ int exceptions;
+};
/* Set the rounding mode, or restore the saved value. */
#define IF_ROUND_INIT_ /* Empty. */
@@ -1262,6 +1270,13 @@ struct test_ff_i_data
MAX_ULP, EXCEPTIONS) \
check_long (TEST_NAME, FUNC (FUNC_NAME) (ARG), EXPECTED, \
MAX_ULP, EXCEPTIONS)
+#define RUN_TEST_LOOP_f_l(FUNC_NAME, ARRAY, ROUNDING_MODE) \
+ IF_ROUND_INIT_ ## ROUNDING_MODE \
+ for (size_t i = 0; i < sizeof (ARRAY) / sizeof (ARRAY)[0]; i++) \
+ RUN_TEST_f_l ((ARRAY)[i].test_name, FUNC_NAME, (ARRAY)[i].arg, \
+ (ARRAY)[i].expected, (ARRAY)[i].max_ulp, \
+ (ARRAY)[i].exceptions); \
+ ROUND_RESTORE_ ## ROUNDING_MODE
#define RUN_TEST_f_L(TEST_NAME, FUNC_NAME, ARG, EXPECTED, \
MAX_ULP, EXCEPTIONS) \
check_longlong (TEST_NAME, FUNC (FUNC_NAME) (ARG), EXPECTED, \
@@ -10346,195 +10361,187 @@ lgamma_test (void)
}
+static const struct test_f_l_data lrint_test_data[] =
+ {
+ START_DATA (lrint),
+ /* XXX this test is incomplete. We need to have a way to specifiy
+ the rounding method and test the critical cases. So far, only
+ unproblematic numbers are tested. */
+ /* TODO: missing +/-Inf as well as qNaN tests. */
+ TEST_f_l (lrint, 0.0, 0),
+ TEST_f_l (lrint, minus_zero, 0),
+ TEST_f_l (lrint, 0.2L, 0),
+ TEST_f_l (lrint, -0.2L, 0),
+
+ TEST_f_l (lrint, 1.4L, 1),
+ TEST_f_l (lrint, -1.4L, -1),
+
+ TEST_f_l (lrint, 8388600.3L, 8388600),
+ TEST_f_l (lrint, -8388600.3L, -8388600),
+
+ TEST_f_l (lrint, 1071930.0008, 1071930),
+#ifndef TEST_FLOAT
+ TEST_f_l (lrint, 1073741824.01, 1073741824),
+# if LONG_MAX > 281474976710656
+ TEST_f_l (lrint, 281474976710656.025, 281474976710656),
+# endif
+#endif
+ END_DATA (lrint)
+ };
+
static void
lrint_test (void)
{
- /* XXX this test is incomplete. We need to have a way to specifiy
- the rounding method and test the critical cases. So far, only
- unproblematic numbers are tested. */
- /* TODO: missing +/-Inf as well as qNaN tests. */
-
START (lrint);
+ RUN_TEST_LOOP_f_l (lrint, lrint_test_data, );
+ END (lrint);
+}
+
- TEST_f_l (lrint, 0.0, 0);
- TEST_f_l (lrint, minus_zero, 0);
- TEST_f_l (lrint, 0.2L, 0);
- TEST_f_l (lrint, -0.2L, 0);
+static const struct test_f_l_data lrint_tonearest_test_data[] =
+ {
+ START_DATA (lrint_tonearest),
+ TEST_f_l (lrint, 0.0, 0),
+ TEST_f_l (lrint, minus_zero, 0),
+ TEST_f_l (lrint, 0.2L, 0),
+ TEST_f_l (lrint, -0.2L, 0),
+ TEST_f_l (lrint, 0.5L, 0),
+ TEST_f_l (lrint, -0.5L, 0),
+ TEST_f_l (lrint, 0.8L, 1),
+ TEST_f_l (lrint, -0.8L, -1),
- TEST_f_l (lrint, 1.4L, 1);
- TEST_f_l (lrint, -1.4L, -1);
+ TEST_f_l (lrint, 1.4L, 1),
+ TEST_f_l (lrint, -1.4L, -1),
- TEST_f_l (lrint, 8388600.3L, 8388600);
- TEST_f_l (lrint, -8388600.3L, -8388600);
+ TEST_f_l (lrint, 8388600.3L, 8388600),
+ TEST_f_l (lrint, -8388600.3L, -8388600),
- TEST_f_l (lrint, 1071930.0008, 1071930);
+ TEST_f_l (lrint, 1071930.0008, 1071930),
#ifndef TEST_FLOAT
- TEST_f_l (lrint, 1073741824.01, 1073741824);
+ TEST_f_l (lrint, 1073741824.01, 1073741824),
# if LONG_MAX > 281474976710656
- TEST_f_l (lrint, 281474976710656.025, 281474976710656);
+ TEST_f_l (lrint, 281474976710656.025, 281474976710656),
# endif
#endif
-
- END (lrint);
-}
-
+ END_DATA (lrint_tonearest)
+ };
static void
lrint_test_tonearest (void)
{
- int save_round_mode;
START (lrint_tonearest);
+ RUN_TEST_LOOP_f_l (lrint, lrint_tonearest_test_data, FE_TONEAREST);
+ END (lrint_tonearest);
+}
- save_round_mode = fegetround ();
- if (!fesetround (FE_TONEAREST))
- {
- TEST_f_l (lrint, 0.0, 0);
- TEST_f_l (lrint, minus_zero, 0);
- TEST_f_l (lrint, 0.2L, 0);
- TEST_f_l (lrint, -0.2L, 0);
- TEST_f_l (lrint, 0.5L, 0);
- TEST_f_l (lrint, -0.5L, 0);
- TEST_f_l (lrint, 0.8L, 1);
- TEST_f_l (lrint, -0.8L, -1);
+static const struct test_f_l_data lrint_towardzero_test_data[] =
+ {
+ START_DATA (lrint_towardzero),
+ TEST_f_l (lrint, 0.0, 0),
+ TEST_f_l (lrint, minus_zero, 0),
+ TEST_f_l (lrint, 0.2L, 0),
+ TEST_f_l (lrint, -0.2L, 0),
+ TEST_f_l (lrint, 0.5L, 0),
+ TEST_f_l (lrint, -0.5L, 0),
+ TEST_f_l (lrint, 0.8L, 0),
+ TEST_f_l (lrint, -0.8L, 0),
- TEST_f_l (lrint, 1.4L, 1);
- TEST_f_l (lrint, -1.4L, -1);
+ TEST_f_l (lrint, 1.4L, 1),
+ TEST_f_l (lrint, -1.4L, -1),
- TEST_f_l (lrint, 8388600.3L, 8388600);
- TEST_f_l (lrint, -8388600.3L, -8388600);
+ TEST_f_l (lrint, 8388600.3L, 8388600),
+ TEST_f_l (lrint, -8388600.3L, -8388600),
- TEST_f_l (lrint, 1071930.0008, 1071930);
+ TEST_f_l (lrint, 1071930.0008, 1071930),
#ifndef TEST_FLOAT
- TEST_f_l (lrint, 1073741824.01, 1073741824);
+ TEST_f_l (lrint, 1073741824.01, 1073741824),
# if LONG_MAX > 281474976710656
- TEST_f_l (lrint, 281474976710656.025, 281474976710656);
+ TEST_f_l (lrint, 281474976710656.025, 281474976710656),
# endif
#endif
- }
-
- fesetround (save_round_mode);
-
- END (lrint_tonearest);
-}
-
+ END_DATA (lrint_towardzero)
+ };
static void
lrint_test_towardzero (void)
{
- int save_round_mode;
START (lrint_towardzero);
+ RUN_TEST_LOOP_f_l (lrint, lrint_towardzero_test_data, FE_TOWARDZERO);
+ END (lrint_towardzero);
+}
- save_round_mode = fegetround ();
- if (!fesetround (FE_TOWARDZERO))
- {
- TEST_f_l (lrint, 0.0, 0);
- TEST_f_l (lrint, minus_zero, 0);
- TEST_f_l (lrint, 0.2L, 0);
- TEST_f_l (lrint, -0.2L, 0);
- TEST_f_l (lrint, 0.5L, 0);
- TEST_f_l (lrint, -0.5L, 0);
- TEST_f_l (lrint, 0.8L, 0);
- TEST_f_l (lrint, -0.8L, 0);
+static const struct test_f_l_data lrint_downward_test_data[] =
+ {
+ START_DATA (lrint_downward),
+ TEST_f_l (lrint, 0.0, 0),
+ TEST_f_l (lrint, minus_zero, 0),
+ TEST_f_l (lrint, 0.2L, 0),
+ TEST_f_l (lrint, -0.2L, -1),
+ TEST_f_l (lrint, 0.5L, 0),
+ TEST_f_l (lrint, -0.5L, -1),
+ TEST_f_l (lrint, 0.8L, 0),
+ TEST_f_l (lrint, -0.8L, -1),
- TEST_f_l (lrint, 1.4L, 1);
- TEST_f_l (lrint, -1.4L, -1);
+ TEST_f_l (lrint, 1.4L, 1),
+ TEST_f_l (lrint, -1.4L, -2),
- TEST_f_l (lrint, 8388600.3L, 8388600);
- TEST_f_l (lrint, -8388600.3L, -8388600);
+ TEST_f_l (lrint, 8388600.3L, 8388600),
+ TEST_f_l (lrint, -8388600.3L, -8388601),
- TEST_f_l (lrint, 1071930.0008, 1071930);
+ TEST_f_l (lrint, 1071930.0008, 1071930),
#ifndef TEST_FLOAT
- TEST_f_l (lrint, 1073741824.01, 1073741824);
+ TEST_f_l (lrint, 1073741824.01, 1073741824),
# if LONG_MAX > 281474976710656
- TEST_f_l (lrint, 281474976710656.025, 281474976710656);
+ TEST_f_l (lrint, 281474976710656.025, 281474976710656),
# endif
#endif
- }
-
- fesetround (save_round_mode);
-
- END (lrint_towardzero);
-}
-
+ END_DATA (lrint_downward)
+ };
static void
lrint_test_downward (void)
{
- int save_round_mode;
START (lrint_downward);
+ RUN_TEST_LOOP_f_l (lrint, lrint_downward_test_data, FE_DOWNWARD);
+ END (lrint_downward);
+}
- save_round_mode = fegetround ();
- if (!fesetround (FE_DOWNWARD))
- {
- TEST_f_l (lrint, 0.0, 0);
- TEST_f_l (lrint, minus_zero, 0);
- TEST_f_l (lrint, 0.2L, 0);
- TEST_f_l (lrint, -0.2L, -1);
- TEST_f_l (lrint, 0.5L, 0);
- TEST_f_l (lrint, -0.5L, -1);
- TEST_f_l (lrint, 0.8L, 0);
- TEST_f_l (lrint, -0.8L, -1);
+static const struct test_f_l_data lrint_upward_test_data[] =
+ {
+ START_DATA (lrint_upward),
+ TEST_f_l (lrint, 0.0, 0),
+ TEST_f_l (lrint, minus_zero, 0),
+ TEST_f_l (lrint, 0.2L, 1),
+ TEST_f_l (lrint, -0.2L, 0),
+ TEST_f_l (lrint, 0.5L, 1),
+ TEST_f_l (lrint, -0.5L, 0),
+ TEST_f_l (lrint, 0.8L, 1),
+ TEST_f_l (lrint, -0.8L, 0),
- TEST_f_l (lrint, 1.4L, 1);
- TEST_f_l (lrint, -1.4L, -2);
+ TEST_f_l (lrint, 1.4L, 2),
+ TEST_f_l (lrint, -1.4L, -1),
- TEST_f_l (lrint, 8388600.3L, 8388600);
- TEST_f_l (lrint, -8388600.3L, -8388601);
+ TEST_f_l (lrint, 8388600.3L, 8388601),
+ TEST_f_l (lrint, -8388600.3L, -8388600),
- TEST_f_l (lrint, 1071930.0008, 1071930);
#ifndef TEST_FLOAT
- TEST_f_l (lrint, 1073741824.01, 1073741824);
-# if LONG_MAX > 281474976710656
- TEST_f_l (lrint, 281474976710656.025, 281474976710656);
+ TEST_f_l (lrint, 1071930.0008, 1071931),
+ TEST_f_l (lrint, 1073741824.01, 1073741825),
+# if LONG_MAX > 281474976710656 && defined (TEST_LDOUBLE)
+ TEST_f_l (lrint, 281474976710656.025, 281474976710656),
# endif
#endif
- }
-
- fesetround (save_round_mode);
-
- END (lrint_downward);
-}
-
+ END_DATA (lrint_upward)
+ };
static void
lrint_test_upward (void)
{
- int save_round_mode;
START (lrint_upward);
-
- save_round_mode = fegetround ();
-
- if (!fesetround (FE_UPWARD))
- {
- TEST_f_l (lrint, 0.0, 0);
- TEST_f_l (lrint, minus_zero, 0);
- TEST_f_l (lrint, 0.2L, 1);
- TEST_f_l (lrint, -0.2L, 0);
- TEST_f_l (lrint, 0.5L, 1);
- TEST_f_l (lrint, -0.5L, 0);
- TEST_f_l (lrint, 0.8L, 1);
- TEST_f_l (lrint, -0.8L, 0);
-
- TEST_f_l (lrint, 1.4L, 2);
- TEST_f_l (lrint, -1.4L, -1);
-
- TEST_f_l (lrint, 8388600.3L, 8388601);
- TEST_f_l (lrint, -8388600.3L, -8388600);
-
-#ifndef TEST_FLOAT
- TEST_f_l (lrint, 1071930.0008, 1071931);
- TEST_f_l (lrint, 1073741824.01, 1073741825);
-# if LONG_MAX > 281474976710656 && defined (TEST_LDOUBLE)
- TEST_f_l (lrint, 281474976710656.025, 281474976710656);
-# endif
-#endif
- }
-
- fesetround (save_round_mode);
-
+ RUN_TEST_LOOP_f_l (lrint, lrint_upward_test_data, FE_UPWARD);
END (lrint_upward);
}
@@ -11442,46 +11449,51 @@ logb_test_downward (void)
END (logb_downward);
}
-static void
-lround_test (void)
-{
- /* TODO: missing +/-Inf as well as qNaN tests. */
-
- START (lround);
-
- TEST_f_l (lround, 0, 0);
- TEST_f_l (lround, minus_zero, 0);
- TEST_f_l (lround, 0.2L, 0.0);
- TEST_f_l (lround, -0.2L, 0);
- TEST_f_l (lround, 0.5, 1);
- TEST_f_l (lround, -0.5, -1);
- TEST_f_l (lround, 0.8L, 1);
- TEST_f_l (lround, -0.8L, -1);
- TEST_f_l (lround, 1.5, 2);
- TEST_f_l (lround, -1.5, -2);
- TEST_f_l (lround, 22514.5, 22515);
- TEST_f_l (lround, -22514.5, -22515);
- TEST_f_l (lround, 1071930.0008, 1071930);
-#ifndef TEST_FLOAT
- TEST_f_l (lround, 1073741824.01, 1073741824);
+static const struct test_f_l_data lround_test_data[] =
+ {
+ START_DATA (lround),
+ /* TODO: missing +/-Inf as well as qNaN tests. */
+ TEST_f_l (lround, 0, 0),
+ TEST_f_l (lround, minus_zero, 0),
+ TEST_f_l (lround, 0.2L, 0.0),
+ TEST_f_l (lround, -0.2L, 0),
+ TEST_f_l (lround, 0.5, 1),
+ TEST_f_l (lround, -0.5, -1),
+ TEST_f_l (lround, 0.8L, 1),
+ TEST_f_l (lround, -0.8L, -1),
+ TEST_f_l (lround, 1.5, 2),
+ TEST_f_l (lround, -1.5, -2),
+ TEST_f_l (lround, 22514.5, 22515),
+ TEST_f_l (lround, -22514.5, -22515),
+ TEST_f_l (lround, 1071930.0008, 1071930),
+#ifndef TEST_FLOAT
+ TEST_f_l (lround, 1073741824.01, 1073741824),
# if LONG_MAX > 281474976710656
- TEST_f_l (lround, 281474976710656.025, 281474976710656);
- TEST_f_l (lround, 18014398509481974, 18014398509481974);
+ TEST_f_l (lround, 281474976710656.025, 281474976710656),
+ TEST_f_l (lround, 18014398509481974, 18014398509481974),
# endif
- TEST_f_l (lround, 2097152.5, 2097153);
- TEST_f_l (lround, -2097152.5, -2097153);
- /* nextafter(0.5,-1) */
- TEST_f_l (lround, 0x1.fffffffffffffp-2, 0);
- /* nextafter(-0.5,1) */
- TEST_f_l (lround, -0x1.fffffffffffffp-2, 0);
+ TEST_f_l (lround, 2097152.5, 2097153),
+ TEST_f_l (lround, -2097152.5, -2097153),
+ /* nextafter(0.5,-1) */
+ TEST_f_l (lround, 0x1.fffffffffffffp-2, 0),
+ /* nextafter(-0.5,1) */
+ TEST_f_l (lround, -0x1.fffffffffffffp-2, 0),
#else
- /* nextafter(0.5,-1) */
- TEST_f_l (lround, 0x1.fffffp-2, 0);
- /* nextafter(-0.5,1) */
- TEST_f_l (lround, -0x1.fffffp-2, 0);
- TEST_f_l (lround, 0x1.fffffep+23, 16777215);
- TEST_f_l (lround, -0x1.fffffep+23, -16777215);
+ /* nextafter(0.5,-1) */
+ TEST_f_l (lround, 0x1.fffffp-2, 0),
+ /* nextafter(-0.5,1) */
+ TEST_f_l (lround, -0x1.fffffp-2, 0),
+ TEST_f_l (lround, 0x1.fffffep+23, 16777215),
+ TEST_f_l (lround, -0x1.fffffep+23, -16777215),
#endif
+ END_DATA (lround)
+ };
+
+static void
+lround_test (void)
+{
+ START (lround);
+ RUN_TEST_LOOP_f_l (lround, lround_test_data, );
END (lround);
}