(M_PI_6l, M_E2l, M_E3l, M_2_SQRT_PIl, M_SQRT_PIl): Increase precision.
[kopensolaris-gnu/glibc.git] / math / libm-test.inc
index d2ceb2f..f4df889 100644 (file)
@@ -1,4 +1,4 @@
-/* Copyright (C) 1997, 1998, 1999 Free Software Foundation, Inc.
+/* Copyright (C) 1997, 1998, 1999, 2000 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
    Contributed by Andreas Jaeger <aj@arthur.rhein-neckar.de>, 1997.
 
@@ -77,7 +77,7 @@
    testing.  Handling of specific inputs (e.g. infinity, not-a-number)
    is also tested.  Correct handling of exceptions is checked
    against.  These implemented tests should check all cases that are
-   specified in ISO C 9X.
+   specified in ISO C99.
 
    Exception testing: At the moment only divide-by-zero and invalid
    exceptions are tested.  Overflow/underflow and inexact exceptions
 #include <errno.h>
 #include <stdlib.h>
 #include <stdio.h>
+#include <string.h>
 #include <argp.h>
 
 /* Possible exceptions */
 #define IGNORE_ZERO_INF_SIGN           0x10
 
 /* Various constants (we must supply them precalculated for accuracy).  */
-#define M_PI_6l                        .52359877559829887308L
-#define M_E2l                  7.389056098930650227230L
-#define M_E3l                  20.08553692318766774093L
-#define M_2_SQRT_PIl           3.5449077018110320545963L       /* 2 sqrt (M_PIl)  */
-#define M_SQRT_PIl             1.77245385090551602729817L      /* sqrt (M_PIl)  */
+#define M_PI_6l                        .52359877559829887307710723054658383L
+#define M_E2l                  7.389056098930650227230427460575008L
+#define M_E3l                  20.085536923187667740928529654581719L
+#define M_2_SQRT_PIl           3.5449077018110320545963349666822903L   /* 2 sqrt (M_PIl)  */
+#define M_SQRT_PIl             1.7724538509055160272981674833411451L   /* sqrt (M_PIl)  */
 #define M_LOG_SQRT_PIl         0.572364942924700087072L        /* log(sqrt(M_PIl))  */
 #define M_LOG_2_SQRT_PIl       1.265512123484645396489L        /* log(2*sqrt(M_PIl))  */
 #define M_PI_34l               (M_PIl - M_PI_4l)               /* 3*pi/4 */
@@ -188,6 +189,7 @@ init_max_error (void)
   max_error = 0;
   real_max_error = 0;
   imag_max_error = 0;
+  feclearexcept (FE_ALL_EXCEPT);
 }
 
 static void
@@ -692,6 +694,12 @@ matherr (struct exception *x __attribute__ ((unused)))
 static void
 acos_test (void)
 {
+  errno = 0;
+  FUNC(acos) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (acos);
 
   TEST_f_f (acos, plus_infty, nan_value, INVALID_EXCEPTION);
@@ -716,6 +724,12 @@ acos_test (void)
 static void
 acosh_test (void)
 {
+  errno = 0;
+  FUNC(acosh) (7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (acosh);
 
   TEST_f_f (acosh, plus_infty, plus_infty);
@@ -733,6 +747,11 @@ acosh_test (void)
 static void
 asin_test (void)
 {
+  errno = 0;
+  FUNC(asin) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (asin);
 
@@ -758,6 +777,12 @@ asin_test (void)
 static void
 asinh_test (void)
 {
+  errno = 0;
+  FUNC(asinh) (0.7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (asinh);
 
   TEST_f_f (asinh, 0, 0);
@@ -775,6 +800,11 @@ asinh_test (void)
 static void
 atan_test (void)
 {
+  errno = 0;
+  FUNC(atan) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (atan);
 
@@ -798,6 +828,11 @@ atan_test (void)
 static void
 atanh_test (void)
 {
+  errno = 0;
+  FUNC(atanh) (0.7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (atanh);
 
@@ -821,6 +856,12 @@ atanh_test (void)
 static void
 atan2_test (void)
 {
+  errno = 0;
+  FUNC(atan2) (-0, 1);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (atan2);
 
   /* atan2 (0,x) == 0 for x > 0.  */
@@ -887,6 +928,12 @@ atan2_test (void)
 static void
 cabs_test (void)
 {
+  errno = 0;
+  FUNC(cabs) (BUILD_COMPLEX (0.7, 12.4));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cabs);
 
   /* cabs (x + iy) is specified as hypot (x,y) */
@@ -927,6 +974,12 @@ cabs_test (void)
 static void
 cacos_test (void)
 {
+  errno = 0;
+  FUNC(cacos) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cacos);
 
 
@@ -985,6 +1038,12 @@ cacos_test (void)
 static void
 cacosh_test (void)
 {
+  errno = 0;
+  FUNC(cacosh) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cacosh);
 
 
@@ -1108,6 +1167,12 @@ carg_test (void)
 static void
 casin_test (void)
 {
+  errno = 0;
+  FUNC(casin) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (casin);
 
   TEST_c_c (casin, 0, 0, 0.0, 0.0);
@@ -1166,6 +1231,12 @@ casin_test (void)
 static void
 casinh_test (void)
 {
+  errno = 0;
+  FUNC(casinh) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (casinh);
 
   TEST_c_c (casinh, 0, 0, 0.0, 0.0);
@@ -1224,6 +1295,12 @@ casinh_test (void)
 static void
 catan_test (void)
 {
+  errno = 0;
+  FUNC(catan) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (catan);
 
   TEST_c_c (catan, 0, 0, 0, 0);
@@ -1286,6 +1363,12 @@ catan_test (void)
 static void
 catanh_test (void)
 {
+  errno = 0;
+  FUNC(catanh) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (catanh);
 
   TEST_c_c (catanh, 0, 0, 0.0, 0.0);
@@ -1346,6 +1429,12 @@ catanh_test (void)
 static void
 cbrt_test (void)
 {
+  errno = 0;
+  FUNC(cbrt) (8);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cbrt);
 
   TEST_f_f (cbrt, 0.0, 0.0);
@@ -1355,10 +1444,10 @@ cbrt_test (void)
   TEST_f_f (cbrt, minus_infty, minus_infty);
   TEST_f_f (cbrt, nan_value, nan_value);
 
-  TEST_f_f (cbrt, -0.001, -0.1);
+  TEST_f_f (cbrt, -0.001, -0.1L);
   TEST_f_f (cbrt, 8, 2);
   TEST_f_f (cbrt, -27.0, -3.0);
-  TEST_f_f (cbrt, 0.970299, 0.99);
+  TEST_f_f (cbrt, 0.970299, 0.99L);
   TEST_f_f (cbrt, 0.7, 0.8879040017426007084L);
 
   END (cbrt);
@@ -1367,6 +1456,11 @@ cbrt_test (void)
 static void
 ccos_test (void)
 {
+  errno = 0;
+  FUNC(ccos) (BUILD_COMPLEX (0, 0));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (ccos);
 
@@ -1431,6 +1525,11 @@ ccos_test (void)
 static void
 ccosh_test (void)
 {
+  errno = 0;
+  FUNC(ccosh) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (ccosh);
 
@@ -1513,6 +1612,12 @@ ceil_test (void)
 static void
 cexp_test (void)
 {
+  errno = 0;
+  FUNC(cexp) (BUILD_COMPLEX (0, 0));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cexp);
 
   TEST_c_c (cexp, plus_zero, plus_zero, 1, 0.0);
@@ -1585,6 +1690,12 @@ cimag_test (void)
 static void
 clog_test (void)
 {
+  errno = 0;
+  FUNC(clog) (BUILD_COMPLEX (-2, -3));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (clog);
 
   TEST_c_c (clog, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
@@ -1644,6 +1755,12 @@ clog_test (void)
 static void
 clog10_test (void)
 {
+  errno = 0;
+  FUNC(clog10) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (clog10);
 
   TEST_c_c (clog10, minus_zero, 0, minus_infty, M_PIl, DIVIDE_BY_ZERO_EXCEPTION);
@@ -1748,6 +1865,12 @@ copysign_test (void)
 static void
 cos_test (void)
 {
+  errno = 0;
+  FUNC(cos) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cos);
 
   TEST_f_f (cos, 0, 1);
@@ -1768,6 +1891,12 @@ cos_test (void)
 static void
 cosh_test (void)
 {
+  errno = 0;
+  FUNC(cosh) (0.7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cosh);
   TEST_f_f (cosh, 0, 1);
   TEST_f_f (cosh, minus_zero, 1);
@@ -1786,6 +1915,12 @@ cosh_test (void)
 static void
 cpow_test (void)
 {
+  errno = 0;
+  FUNC(cpow) (BUILD_COMPLEX (1, 0), BUILD_COMPLEX (0, 0));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (cpow);
 
   TEST_cc_c (cpow, 1, 0, 0, 0, 1.0, 0.0);
@@ -1839,6 +1974,11 @@ creal_test (void)
 static void
 csin_test (void)
 {
+  errno = 0;
+  FUNC(csin) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (csin);
 
@@ -1903,6 +2043,11 @@ csin_test (void)
 static void
 csinh_test (void)
 {
+  errno = 0;
+  FUNC(csinh) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (csinh);
 
@@ -1965,6 +2110,12 @@ csinh_test (void)
 static void
 csqrt_test (void)
 {
+  errno = 0;
+  FUNC(csqrt) (BUILD_COMPLEX (-1, 0));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (csqrt);
 
   TEST_c_c (csqrt, 0, 0, 0.0, 0.0);
@@ -2025,6 +2176,12 @@ csqrt_test (void)
 static void
 ctan_test (void)
 {
+  errno = 0;
+  FUNC(ctan) (BUILD_COMPLEX (0.7, 1.2));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (ctan);
 
   TEST_c_c (ctan, 0, 0, 0.0, 0.0);
@@ -2077,6 +2234,12 @@ ctan_test (void)
 static void
 ctanh_test (void)
 {
+  errno = 0;
+  FUNC(ctanh) (BUILD_COMPLEX (0, 0));
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (ctanh);
 
   TEST_c_c (ctanh, 0, 0, 0.0, 0.0);
@@ -2184,6 +2347,12 @@ erfc_test (void)
 static void
 exp_test (void)
 {
+  errno = 0;
+  FUNC(exp) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (exp);
 
   TEST_f_f (exp, 0, 1);
@@ -2222,7 +2391,7 @@ exp10_test (void)
   TEST_f_f (exp10, minus_infty, 0);
   TEST_f_f (exp10, nan_value, nan_value);
   TEST_f_f (exp10, 3, 1000);
-  TEST_f_f (exp10, -1, 0.1);
+  TEST_f_f (exp10, -1, 0.1L);
   TEST_f_f (exp10, 1e6, plus_infty);
   TEST_f_f (exp10, -1e6, 0);
   TEST_f_f (exp10, 0.7, 5.0118723362727228500L);
@@ -2259,6 +2428,12 @@ exp2_test (void)
 static void
 expm1_test (void)
 {
+  errno = 0;
+  FUNC(expm1) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (expm1);
 
   TEST_f_f (expm1, 0, 0);
@@ -2450,6 +2625,11 @@ fmin_test (void)
 static void
 fmod_test (void)
 {
+  errno = 0;
+  FUNC(fmod) (6.5, 2.3);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (fmod);
 
@@ -2475,10 +2655,10 @@ fmod_test (void)
 
   TEST_ff_f (fmod, nan_value, nan_value, nan_value);
 
-  TEST_ff_f (fmod, 6.5, 2.3, 1.9);
-  TEST_ff_f (fmod, -6.5, 2.3, -1.9);
-  TEST_ff_f (fmod, 6.5, -2.3, 1.9);
-  TEST_ff_f (fmod, -6.5, -2.3, -1.9);
+  TEST_ff_f (fmod, 6.5, 2.3, 1.9L);
+  TEST_ff_f (fmod, -6.5, 2.3, -1.9L);
+  TEST_ff_f (fmod, 6.5, -2.3, 1.9L);
+  TEST_ff_f (fmod, -6.5, -2.3, -1.9L);
 
   END (fmod);
 }
@@ -2513,8 +2693,8 @@ frexp_test (void)
   TEST_fI_f1 (frexp, 0.0, 0.0, 0.0);
   TEST_fI_f1 (frexp, minus_zero, minus_zero, 0.0);
 
-  TEST_fI_f1 (frexp, 12.8, 0.8, 4);
-  TEST_fI_f1 (frexp, -27.34, -0.854375, 5);
+  TEST_fI_f1 (frexp, 12.8, 0.8L, 4);
+  TEST_fI_f1 (frexp, -27.34, -0.854375L, 5);
 
   END (frexp);
 }
@@ -2551,6 +2731,12 @@ gamma_test (void)
 static void
 hypot_test (void)
 {
+  errno = 0;
+  FUNC(hypot) (0.7, 12.4);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (hypot);
 
   TEST_ff_f (hypot, plus_infty, 1, plus_infty, IGNORE_ZERO_INF_SIGN);
@@ -2576,9 +2762,9 @@ hypot_test (void)
   TEST_ff_f (hypot, -12.4, -0.7, 12.41974234837422060118L);
 
   /*  hypot (x,0) == fabs (x)  */
-  TEST_ff_f (hypot, 0.7, 0, 0.7);
-  TEST_ff_f (hypot, -0.7, 0, 0.7);
-  TEST_ff_f (hypot, -5.7e7, 0, 5.7e7);
+  TEST_ff_f (hypot, 0.7, 0, 0.7L);
+  TEST_ff_f (hypot, -0.7, 0, 0.7L);
+  TEST_ff_f (hypot, -5.7e7, 0, 5.7e7L);
 
   TEST_ff_f (hypot, 0.7, 1.2, 1.3892443989449804508L);
 
@@ -2879,6 +3065,11 @@ llrint_test (void)
 static void
 log_test (void)
 {
+  errno = 0;
+  FUNC(log) (1);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
   START (log);
 
   TEST_f_f (log, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
@@ -2902,6 +3093,12 @@ log_test (void)
 static void
 log10_test (void)
 {
+  errno = 0;
+  FUNC(log10) (1);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (log10);
 
   TEST_f_f (log10, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
@@ -2929,6 +3126,12 @@ log10_test (void)
 static void
 log1p_test (void)
 {
+  errno = 0;
+  FUNC(log1p) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (log1p);
 
   TEST_f_f (log1p, 0, 0);
@@ -2951,6 +3154,12 @@ log1p_test (void)
 static void
 log2_test (void)
 {
+  errno = 0;
+  FUNC(log2) (1);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (log2);
 
   TEST_f_f (log2, 0, minus_infty, DIVIDE_BY_ZERO_EXCEPTION);
@@ -3101,6 +3310,12 @@ nearbyint_test (void)
   TEST_f_f (nearbyint, minus_infty, minus_infty);
   TEST_f_f (nearbyint, nan_value, nan_value);
 
+  /* Default rounding mode is round to nearest.  */
+  TEST_f_f (nearbyint, 0.5, 0.0);
+  TEST_f_f (nearbyint, 1.5, 2.0);
+  TEST_f_f (nearbyint, -0.5, minus_zero);
+  TEST_f_f (nearbyint, -1.5, -2.0);
+
   END (nearbyint);
 }
 
@@ -3160,6 +3375,12 @@ static void
 pow_test (void)
 {
 
+  errno = 0;
+  FUNC(pow) (0, 0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (pow);
 
   TEST_ff_f (pow, 0, 0, 1);
@@ -3330,6 +3551,11 @@ pow_test (void)
 static void
 remainder_test (void)
 {
+  errno = 0;
+  FUNC(remainder) (1.625, 1.0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
 
   START (remainder);
 
@@ -3355,6 +3581,12 @@ remquo_test (void)
   /* x is needed.  */
   int x;
 
+  errno = 0;
+  FUNC(remquo) (1.625, 1.0, &x);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (remquo);
 
   TEST_ffI_f1 (remquo, 1, 0, nan_value, IGNORE, INVALID_EXCEPTION);
@@ -3384,6 +3616,12 @@ rint_test (void)
   TEST_f_f (rint, plus_infty, plus_infty);
   TEST_f_f (rint, minus_infty, minus_infty);
 
+  /* Default rounding mode is round to nearest.  */
+  TEST_f_f (rint, 0.5, 0.0);
+  TEST_f_f (rint, 1.5, 2.0);
+  TEST_f_f (rint, -0.5, -0.0);
+  TEST_f_f (rint, -1.5, -2.0);
+
   END (rint);
 }
 
@@ -3533,6 +3771,12 @@ signbit_test (void)
 static void
 sin_test (void)
 {
+  errno = 0;
+  FUNC(sin) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (sin);
 
   TEST_f_f (sin, 0, 0);
@@ -3556,6 +3800,12 @@ sincos_test (void)
 {
   FLOAT sin_res, cos_res;
 
+  errno = 0;
+  FUNC(sincos) (0, &sin_res, &cos_res);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (sincos);
 
   /* sincos is treated differently because it returns void.  */
@@ -3577,6 +3827,12 @@ sincos_test (void)
 static void
 sinh_test (void)
 {
+  errno = 0;
+  FUNC(sinh) (0.7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (sinh);
   TEST_f_f (sinh, 0, 0);
   TEST_f_f (sinh, minus_zero, minus_zero);
@@ -3595,6 +3851,12 @@ sinh_test (void)
 static void
 sqrt_test (void)
 {
+  errno = 0;
+  FUNC(sqrt) (1);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (sqrt);
 
   TEST_f_f (sqrt, 0, 0);
@@ -3613,7 +3875,7 @@ sqrt_test (void)
   TEST_f_f (sqrt, 2, M_SQRT2l);
   TEST_f_f (sqrt, 0.25, 0.5);
   TEST_f_f (sqrt, 6642.25, 81.5);
-  TEST_f_f (sqrt, 15239.9025, 123.45);
+  TEST_f_f (sqrt, 15239.9025, 123.45L);
   TEST_f_f (sqrt, 0.7, 0.83666002653407554798L);
 
   END (sqrt);
@@ -3622,6 +3884,12 @@ sqrt_test (void)
 static void
 tan_test (void)
 {
+  errno = 0;
+  FUNC(tan) (0);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (tan);
 
   TEST_f_f (tan, 0, 0);
@@ -3639,6 +3907,12 @@ tan_test (void)
 static void
 tanh_test (void)
 {
+  errno = 0;
+  FUNC(tanh) (0.7);
+  if (errno == ENOSYS)
+    /* Function not implemented.  */
+    return;
+
   START (tanh);
 
   TEST_f_f (tanh, 0, 0);
@@ -3965,7 +4239,7 @@ main (int argc, char **argv)
   output_points = 1;
   /* XXX set to 0 for releases.  */
   ignore_max_ulp = 0;
-  
+
   /* Parse and process arguments.  */
   argp_parse (&argp, argc, argv, 0, &remaining, NULL);
 
@@ -3994,7 +4268,7 @@ main (int argc, char **argv)
   check_ulp ();
 #endif
 
-  /* Keep the tests a wee bit ordered (according to ISO 9X).  */
+  /* Keep the tests a wee bit ordered (according to ISO C99).  */
   /* Classification macros:  */
   fpclassify_test ();
   isfinite_test ();