about summary refs log tree commit diff
path: root/math
diff options
context:
space:
mode:
Diffstat (limited to 'math')
-rw-r--r--math/libm-test.c464
-rw-r--r--math/math.h66
2 files changed, 268 insertions, 262 deletions
diff --git a/math/libm-test.c b/math/libm-test.c
index 3dffcad6dc..8fa473eaf5 100644
--- a/math/libm-test.c
+++ b/math/libm-test.c
@@ -137,9 +137,9 @@
 #define NO_PRINT 0
 
 /* Various constants (we must supply them precalculated for accuracy).  */
-#define M_PI_6  .52359877559829887308L
-#define M_E2	7.389056098930650227230L
-#define M_E3	20.08553692318766774093L
+#define M_PI_6l  .52359877559829887308L
+#define M_E2l	7.389056098930650227230L
+#define M_E3l	20.08553692318766774093L
 
 static int noErrors;   /* number of errors */
 static int noTests;    /* number of tests (without testing exceptions) */
@@ -840,15 +840,15 @@ acos_test (void)
 		   FUNC(acos) (x),
 		   INVALID_EXCEPTION);
 #endif
-  check ("acos (0) == pi/2", FUNC(acos) (0), M_PI_2);
-  check ("acos (-0) == pi/2", FUNC(acos) (minus_zero), M_PI_2);
+  check ("acos (0) == pi/2", FUNC(acos) (0), M_PI_2l);
+  check ("acos (-0) == pi/2", FUNC(acos) (minus_zero), M_PI_2l);
 
   check ("acos (1) == 0", FUNC(acos) (1), 0);
-  check ("acos (-1) == pi", FUNC(acos) (-1), M_PI);
+  check ("acos (-1) == pi", FUNC(acos) (-1), M_PIl);
 
-  check_eps ("acos (0.5) == pi/3", FUNC(acos) (0.5), M_PI_6 * 2.0,
+  check_eps ("acos (0.5) == pi/3", FUNC(acos) (0.5), M_PI_6l * 2.0,
 	     CHOOSE (1e-18, 0, 0));
-  check_eps ("acos (-0.5) == 2*pi/3", FUNC(acos) (-0.5), M_PI_6 * 4.0,
+  check_eps ("acos (-0.5) == 2*pi/3", FUNC(acos) (-0.5), M_PI_6l * 4.0,
 	     CHOOSE (1e-17, 0, 0));
 
   check_eps ("acos (0.7) == 0.795398830...", FUNC(acos) (0.7),
@@ -895,12 +895,12 @@ asin_test (void)
 
   check ("asin (0) == 0", FUNC(asin) (0), 0);
   check ("asin (-0) == -0", FUNC(asin) (minus_zero), minus_zero);
-  check_eps ("asin (0.5) ==  pi/6", FUNC(asin) (0.5), M_PI_6,
+  check_eps ("asin (0.5) ==  pi/6", FUNC(asin) (0.5), M_PI_6l,
 	     CHOOSE(3.5e-18, 0, 2e-7));
-  check_eps ("asin (-0.5) ==  -pi/6", FUNC(asin) (-0.5), -M_PI_6,
+  check_eps ("asin (-0.5) ==  -pi/6", FUNC(asin) (-0.5), -M_PI_6l,
 	     CHOOSE(3.5e-18, 0, 2e-7));
-  check ("asin (1.0) ==  pi/2", FUNC(asin) (1.0), M_PI_2);
-  check ("asin (-1.0) ==  -pi/2", FUNC(asin) (-1.0), -M_PI_2);
+  check ("asin (1.0) ==  pi/2", FUNC(asin) (1.0), M_PI_2l);
+  check ("asin (-1.0) ==  -pi/2", FUNC(asin) (-1.0), -M_PI_2l);
   check_eps ("asin (0.7) ==  0.775397496...", FUNC(asin) (0.7),
 	     0.7753974966107530637L, CHOOSE(7e-17L, 2e-16, 0));
 }
@@ -928,12 +928,12 @@ atan_test (void)
   check ("atan (0) == 0", FUNC(atan) (0), 0);
   check ("atan (-0) == -0", FUNC(atan) (minus_zero), minus_zero);
 
-  check ("atan (+inf) == pi/2", FUNC(atan) (plus_infty), M_PI_2);
-  check ("atan (-inf) == -pi/2", FUNC(atan) (minus_infty), -M_PI_2);
+  check ("atan (+inf) == pi/2", FUNC(atan) (plus_infty), M_PI_2l);
+  check ("atan (-inf) == -pi/2", FUNC(atan) (minus_infty), -M_PI_2l);
 
-  check_eps ("atan (1) == pi/4", FUNC(atan) (1), M_PI_4,
+  check_eps ("atan (1) == pi/4", FUNC(atan) (1), M_PI_4l,
 	     CHOOSE (1e-18, 0, 0));
-  check_eps ("atan (-1) == -pi/4", FUNC(atan) (1), M_PI_4,
+  check_eps ("atan (-1) == -pi/4", FUNC(atan) (1), M_PI_4l,
 	     CHOOSE (1e-18, 0, 0));
 
   check_eps ("atan (0.7) == 0.610725964...", FUNC(atan) (0.7),
@@ -957,25 +957,27 @@ atan2_test (void)
   check ("atan2 (-0,+0) == -0", FUNC(atan2) (minus_zero, 0), minus_zero);
 
   x = -random_greater (0);
-  check ("atan2 (+0,x) == +pi for x < 0", FUNC(atan2) (0, x), M_PI);
+  check ("atan2 (+0,x) == +pi for x < 0", FUNC(atan2) (0, x), M_PIl);
 
   x = -random_greater (0);
-  check ("atan2 (-0,x) == -pi for x < 0", FUNC(atan2) (minus_zero, x), -M_PI);
+  check ("atan2 (-0,x) == -pi for x < 0", FUNC(atan2) (minus_zero, x), -M_PIl);
 
-  check ("atan2 (+0,-0) == +pi", FUNC(atan2) (0, minus_zero), M_PI);
-  check ("atan2 (-0,-0) == -pi", FUNC(atan2) (minus_zero, minus_zero), -M_PI);
+  check ("atan2 (+0,-0) == +pi", FUNC(atan2) (0, minus_zero), M_PIl);
+  check ("atan2 (-0,-0) == -pi", FUNC(atan2) (minus_zero, minus_zero), -M_PIl);
 
   x = random_greater (0);
-  check ("atan2 (y,+0) == pi/2 for y > 0", FUNC(atan2) (x, 0), M_PI_2);
+  check ("atan2 (y,+0) == pi/2 for y > 0", FUNC(atan2) (x, 0), M_PI_2l);
 
   x = random_greater (0);
-  check ("atan2 (y,-0) == pi/2 for y > 0", FUNC(atan2) (x, minus_zero), M_PI_2);
+  check ("atan2 (y,-0) == pi/2 for y > 0", FUNC(atan2) (x, minus_zero),
+	 M_PI_2l);
 
   x = random_less (0);
-  check ("atan2 (y,+0) == -pi/2 for y < 0", FUNC(atan2) (x, 0), -M_PI_2);
+  check ("atan2 (y,+0) == -pi/2 for y < 0", FUNC(atan2) (x, 0), -M_PI_2l);
 
   x = random_less (0);
-  check ("atan2 (y,-0) == -pi/2 for y < 0", FUNC(atan2) (x, minus_zero), -M_PI_2);
+  check ("atan2 (y,-0) == -pi/2 for y < 0", FUNC(atan2) (x, minus_zero),
+	 -M_PI_2l);
 
   x = random_greater (0);
   check ("atan2 (y,inf) == +0 for finite y > 0",
@@ -987,31 +989,31 @@ atan2_test (void)
 
   x = random_value (-1e4, 1e4);
   check ("atan2(+inf, x) == pi/2 for finite x",
-         FUNC(atan2) (plus_infty, x), M_PI_2);
+         FUNC(atan2) (plus_infty, x), M_PI_2l);
 
   x = random_value (-1e4, 1e4);
   check ("atan2(-inf, x) == -pi/2 for finite x",
-         FUNC(atan2) (minus_infty, x), -M_PI_2);
+         FUNC(atan2) (minus_infty, x), -M_PI_2l);
 
   x = random_greater (0);
   check ("atan2 (y,-inf) == +pi for finite y > 0",
-	 FUNC(atan2) (x, minus_infty), M_PI);
+	 FUNC(atan2) (x, minus_infty), M_PIl);
 
   x = -random_greater (0);
   check ("atan2 (y,-inf) == -pi for finite y < 0",
-	 FUNC(atan2) (x, minus_infty), -M_PI);
+	 FUNC(atan2) (x, minus_infty), -M_PIl);
 
   check ("atan2 (+inf,+inf) == +pi/4",
-	 FUNC(atan2) (plus_infty, plus_infty), M_PI_4);
+	 FUNC(atan2) (plus_infty, plus_infty), M_PI_4l);
 
   check ("atan2 (-inf,+inf) == -pi/4",
-	 FUNC(atan2) (minus_infty, plus_infty), -M_PI_4);
+	 FUNC(atan2) (minus_infty, plus_infty), -M_PI_4l);
 
   check ("atan2 (+inf,-inf) == +3*pi/4",
-	 FUNC(atan2) (plus_infty, minus_infty), 3 * M_PI_4);
+	 FUNC(atan2) (plus_infty, minus_infty), 3 * M_PI_4l);
 
   check ("atan2 (-inf,-inf) == -3*pi/4",
-	 FUNC(atan2) (minus_infty, minus_infty), -3 * M_PI_4);
+	 FUNC(atan2) (minus_infty, minus_infty), -3 * M_PI_4l);
 
   /* FIXME: Add some specific tests */
   check_eps ("atan2 (0.7,1) == 0.61072...", FUNC(atan2) (0.7,1),
@@ -1084,8 +1086,8 @@ ceil_test (void)
   check_isinfp ("ceil (+inf) == +inf", FUNC(ceil) (plus_infty));
   check_isinfn ("ceil (-inf) == -inf", FUNC(ceil) (minus_infty));
 
-  check ("ceil (pi) == 4", FUNC(ceil) (M_PI), 4.0);
-  check ("ceil (-pi) == -3", FUNC(ceil) (-M_PI), -3.0);
+  check ("ceil (pi) == 4", FUNC(ceil) (M_PIl), 4.0);
+  check ("ceil (-pi) == -3", FUNC(ceil) (-M_PIl), -3.0);
 }
 
 
@@ -1102,11 +1104,11 @@ cos_test (void)
 		   FUNC(cos) (minus_infty),
 		   INVALID_EXCEPTION);
 
-  check_eps ("cos (pi/3) == 0.5", FUNC(cos) (M_PI_6 * 2.0),
+  check_eps ("cos (pi/3) == 0.5", FUNC(cos) (M_PI_6l * 2.0),
 	     0.5, CHOOSE (4e-18L, 1e-15L, 1e-7L));
-  check_eps ("cos (2*pi/3) == -0.5", FUNC(cos) (M_PI_6 * 4.0),
+  check_eps ("cos (2*pi/3) == -0.5", FUNC(cos) (M_PI_6l * 4.0),
 	     -0.5, CHOOSE (4e-18L, 1e-15L, 1e-7L));
-  check_eps ("cos (pi/2) == 0", FUNC(cos) (M_PI_2),
+  check_eps ("cos (pi/2) == 0", FUNC(cos) (M_PI_2l),
 	     0, CHOOSE (1e-19L, 1e-16L, 1e-7L));
 
   check_eps ("cos (0.7) == 0.7648421872...", FUNC(cos) (0.7),
@@ -1177,11 +1179,11 @@ exp_test (void)
   check_isinfp ("exp (+inf) == +inf", FUNC(exp) (plus_infty));
   check ("exp (-inf) == 0", FUNC(exp) (minus_infty), 0);
 #endif
-  check_eps ("exp (1) == e", FUNC(exp) (1), M_E, CHOOSE (4e-18L, 0, 0));
+  check_eps ("exp (1) == e", FUNC(exp) (1), M_El, CHOOSE (4e-18L, 0, 0));
 
-  check_eps ("exp (2) == e^2", FUNC(exp) (2), M_E2,
+  check_eps ("exp (2) == e^2", FUNC(exp) (2), M_E2l,
 	     CHOOSE (1e-18, 0, 0));
-  check_eps ("exp (3) == e^3", FUNC(exp) (3), M_E3,
+  check_eps ("exp (3) == e^3", FUNC(exp) (3), M_E3l,
 	     CHOOSE (1.5e-17, 0, 0));
   check_eps ("exp (0.7) == 2.0137527074...", FUNC(exp) (0.7),
 	     2.0137527074704765216L, CHOOSE(9e-17L, 0, 0));
@@ -1222,7 +1224,7 @@ expm1_test (void)
   check ("expm1 (-inf) == -1", FUNC(expm1) (minus_infty), -1);
 #endif
 
-  check_eps ("expm1 (1) == e-1", FUNC(expm1) (1), M_E - 1.0,
+  check_eps ("expm1 (1) == e-1", FUNC(expm1) (1), M_El - 1.0,
 	     CHOOSE (4e-18L, 0, 2e-7));
 
   check_eps ("expm1 (0.7) == 1.01375...", FUNC(expm1) (0.7),
@@ -1418,17 +1420,17 @@ gamma_test (void)
   check_int ("gamma (0) sets signgam to 1", signgam, 1);
 
   signgam = 0;
-  check ("gamma (3) == M_LN2", FUNC(gamma) (3), M_LN2);
+  check ("gamma (3) == M_LN2", FUNC(gamma) (3), M_LN2l);
   check_int ("gamma (3) sets signgam to 1", signgam, 1);
 
   signgam = 0;
   check_eps ("gamma (0.5) == log(sqrt(pi))", FUNC(gamma) (0.5),
-             FUNC(log) (FUNC(sqrt) (M_PI)), CHOOSE (0, 1e-15, 1e-7));
+             FUNC(log) (FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 1e-7));
   check_int ("gamma (0.5) sets signgam to 1", signgam, 1);
 
   signgam = 0;
   check_eps ("gamma (-0.5) == log(2*sqrt(pi))", FUNC(gamma) (-0.5),
-             FUNC(log) (2*FUNC(sqrt) (M_PI)), CHOOSE (0, 1e-15, 0));
+             FUNC(log) (2*FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 0));
 
   check_int ("gamma (-0.5) sets signgam to -1", signgam, -1);
 
@@ -1445,11 +1447,11 @@ gamma_test (void)
                    FUNC(gamma) (minus_infty), INVALID_EXCEPTION);
 
 #ifdef TODO
-  check_eps ("gamma (0.5) == sqrt(pi)", FUNC(gamma) (0.5), FUNC(sqrt) (M_PI),
+  check_eps ("gamma (0.5) == sqrt(pi)", FUNC(gamma) (0.5), FUNC(sqrt) (M_PIl),
 	     CHOOSE (0, 5e-16, 2e-7));
 #endif
   check_eps ("gamma (-0.5) == -2*sqrt(pi)", FUNC(gamma) (-0.5),
-	     -2*FUNC(sqrt) (M_PI), CHOOSE (0, 5e-16, 3e-7));
+	     -2*FUNC(sqrt) (M_PIl), CHOOSE (0, 5e-16, 3e-7));
 
   check ("gamma (1) == 1", FUNC(gamma) (1), 1);
   check ("gamma (4) == 6", FUNC(gamma) (4), 6);
@@ -1486,17 +1488,17 @@ lgamma_test (void)
   check_int ("lgamma (0) sets signgam to 1", signgam, 1);
 
   signgam = 0;
-  check ("lgamma (3) == M_LN2", FUNC(lgamma) (3), M_LN2);
+  check ("lgamma (3) == M_LN2", FUNC(lgamma) (3), M_LN2l);
   check_int ("lgamma (3) sets signgam to 1", signgam, 1);
 
   signgam = 0;
   check_eps ("lgamma (0.5) == log(sqrt(pi))", FUNC(lgamma) (0.5),
-             FUNC(log) (FUNC(sqrt) (M_PI)), CHOOSE (0, 1e-15, 1e-7));
+             FUNC(log) (FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 1e-7));
   check_int ("lgamma (0.5) sets signgam to 1", signgam, 1);
 
   signgam = 0;
   check_eps ("lgamma (-0.5) == log(2*sqrt(pi))", FUNC(lgamma) (-0.5),
-             FUNC(log) (2*FUNC(sqrt) (M_PI)), CHOOSE (0, 1e-15, 0));
+             FUNC(log) (2*FUNC(sqrt) (M_PIl)), CHOOSE (0, 1e-15, 0));
 
   check_int ("lgamma (-0.5) sets signgam to -1", signgam, -1);
 
@@ -1519,7 +1521,7 @@ ilogb_test (void)
   int i;
 
   check_int ("ilogb (1) == 0", FUNC(ilogb) (1), 0);
-  check_int ("ilogb (e) == 1", FUNC(ilogb) (M_E), 1);
+  check_int ("ilogb (e) == 1", FUNC(ilogb) (M_El), 1);
   check_int ("ilogb (1024) == 10", FUNC(ilogb) (1024), 10);
   check_int ("ilogb (-2000) == 10", FUNC(ilogb) (-2000), 10);
 
@@ -1569,11 +1571,11 @@ log_test (void)
 		   FUNC(log) (-1), INVALID_EXCEPTION);
   check_isinfp ("log (+inf) == +inf", FUNC(log) (plus_infty));
 
-  check_eps ("log (e) == 1", FUNC(log) (M_E), 1, CHOOSE (1e-18L, 0, 9e-8L));
-  check_eps ("log (1/e) == -1", FUNC(log) (1.0 / M_E), -1,
+  check_eps ("log (e) == 1", FUNC(log) (M_El), 1, CHOOSE (1e-18L, 0, 9e-8L));
+  check_eps ("log (1/e) == -1", FUNC(log) (1.0 / M_El), -1,
 	     CHOOSE (2e-18L, 0, 0));
-  check ("log (2) == M_LN2", FUNC(log) (2), M_LN2);
-  check_eps ("log (10) == M_LN10", FUNC(log) (10), M_LN10,
+  check ("log (2) == M_LN2", FUNC(log) (2), M_LN2l);
+  check_eps ("log (10) == M_LN10", FUNC(log) (10), M_LN10l,
 	     CHOOSE (1e-18L, 0, 0));
   check_eps ("log (0.7) == -0.3566749439...", FUNC(log) (0.7),
 	     -0.35667494393873237891L, CHOOSE(7e-17L, 6e-17, 3e-8));
@@ -1602,7 +1604,7 @@ log10_test (void)
   check_eps ("log10 (100) == 2", FUNC(log10) (100.0), 2,
 	     CHOOSE (1e-18L, 0, 0));
   check ("log10 (10000) == 4", FUNC(log10) (10000.0), 4);
-  check_eps ("log10 (e) == M_LOG10E", FUNC(log10) (M_E), M_LOG10E,
+  check_eps ("log10 (e) == M_LOG10E", FUNC(log10) (M_El), M_LOG10El,
 	     CHOOSE (1e-18, 0, 9e-8));
   check_eps ("log10 (0.7) == -0.1549019599...", FUNC(log10) (0.7),
 	     -0.15490195998574316929L, CHOOSE(3e-17L, 3e-17, 0));
@@ -1622,7 +1624,7 @@ log1p_test (void)
 
   check_isinfp ("log1p (+inf) == +inf", FUNC(log1p) (plus_infty));
 
-  check_eps ("log1p (e-1) == 1", FUNC(log1p) (M_E - 1.0), 1,
+  check_eps ("log1p (e-1) == 1", FUNC(log1p) (M_El - 1.0), 1,
 	     CHOOSE (1e-18L, 0, 6e-8));
 
   check_eps ("log1p (-0.3) == -0.35667...", FUNC(log1p) (-0.3),
@@ -1645,7 +1647,7 @@ log2_test (void)
 
   check_isinfp ("log2 (+inf) == +inf", FUNC(log2) (plus_infty));
 
-  check_eps ("log2 (e) == M_LOG2E", FUNC(log2) (M_E), M_LOG2E,
+  check_eps ("log2 (e) == M_LOG2E", FUNC(log2) (M_El), M_LOG2El,
 	     CHOOSE (1e-18L, 0, 0));
   check ("log2 (2) == 1", FUNC(log2) (2.0), 1);
   check_eps ("log2 (16) == 4", FUNC(log2) (16.0), 4, CHOOSE (1e-18L, 0, 0));
@@ -1669,7 +1671,7 @@ logb_test (void)
 		    FUNC(logb) (minus_zero), DIVIDE_BY_ZERO_EXCEPTION);
 
   check ("logb (1) == 0", FUNC(logb) (1), 0);
-  check ("logb (e) == 1", FUNC(logb) (M_E), 1);
+  check ("logb (e) == 1", FUNC(logb) (M_El), 1);
   check ("logb (1024) == 10", FUNC(logb) (1024), 10);
   check ("logb (-2000) == 10", FUNC(logb) (-2000), 10);
 
@@ -1826,12 +1828,12 @@ sin_test (void)
 		   FUNC(sin) (minus_infty),
 		   INVALID_EXCEPTION);
 
-  check_eps ("sin (pi/6) == 0.5", FUNC(sin) (M_PI_6),
+  check_eps ("sin (pi/6) == 0.5", FUNC(sin) (M_PI_6l),
 	     0.5, CHOOSE (4e-18L, 0, 0));
-  check_eps ("sin (-pi/6) == -0.5", FUNC(sin) (-M_PI_6),
+  check_eps ("sin (-pi/6) == -0.5", FUNC(sin) (-M_PI_6l),
 	     -0.5, CHOOSE (4e-18L, 0, 0));
-  check ("sin (pi/2) == 1", FUNC(sin) (M_PI_2), 1);
-  check ("sin (-pi/2) == -1", FUNC(sin) (-M_PI_2), -1);
+  check ("sin (pi/2) == 1", FUNC(sin) (M_PI_2l), 1);
+  check ("sin (-pi/2) == -1", FUNC(sin) (-M_PI_2l), -1);
   check_eps ("sin (0.7) == 0.6442176872...", FUNC(sin) (0.7),
 	     0.64421768723769105367L, CHOOSE(4e-17L, 0, 0));
 }
@@ -1888,18 +1890,18 @@ sincos_test (void)
   check_isnan_exc ("sincos (-inf,&sin, &cos) puts NaN in cos plus invalid exception",
 		   cos_res, INVALID_EXCEPTION);
 
-  FUNC(sincos) (M_PI_2, &sin_res, &cos_res);
+  FUNC(sincos) (M_PI_2l, &sin_res, &cos_res);
   fegetenv (&fenv);
   check ("sincos (pi/2, &sin, &cos) puts 1 in sin", sin_res, 1);
   fesetenv (&fenv);
   check_eps ("sincos (pi/2, &sin, &cos) puts 0 in cos", cos_res, 0,
 	     CHOOSE (1e-18L, 1e-16, 1e-7));
 
-  FUNC(sincos) (M_PI_6, &sin_res, &cos_res);
+  FUNC(sincos) (M_PI_6l, &sin_res, &cos_res);
   check_eps ("sincos (pi/6, &sin, &cos) puts 0.5 in sin", sin_res, 0.5,
 	     CHOOSE (5e-18L, 0, 0));
 
-  FUNC(sincos) (M_PI_6*2.0, &sin_res, &cos_res);
+  FUNC(sincos) (M_PI_6l*2.0, &sin_res, &cos_res);
   check_eps ("sincos (pi/3, &sin, &cos) puts 0.5 in cos", cos_res, 0.5,
 	     CHOOSE (5e-18L, 1e-15, 1e-7));
 
@@ -1921,7 +1923,7 @@ tan_test (void)
   check_isnan_exc ("tan (-inf) == NaN plus invalid exception",
 		   FUNC(tan) (minus_infty), INVALID_EXCEPTION);
 
-  check_eps ("tan (pi/4) == 1", FUNC(tan) (M_PI_4), 1,
+  check_eps ("tan (pi/4) == 1", FUNC(tan) (M_PI_4l), 1,
 	     CHOOSE (2e-18L, 1e-15L, 2e-7));
   check_eps ("tan (0.7) == 0.8422883804...", FUNC(tan) (0.7),
 	     0.84228838046307944813L, CHOOSE(8e-17L, 0, 0));
@@ -1953,7 +1955,7 @@ fabs_test (void)
   check_isinfp ("fabs (-inf) == +inf", FUNC(fabs) (minus_infty));
 
   check ("fabs (+38) == 38", FUNC(fabs) (38.0), 38.0);
-  check ("fabs (-e) == e", FUNC(fabs) (-M_E), M_E);
+  check ("fabs (-e) == e", FUNC(fabs) (-M_El), M_El);
 }
 
 
@@ -1965,8 +1967,8 @@ floor_test (void)
   check_isinfp ("floor (+inf) == +inf", FUNC(floor) (plus_infty));
   check_isinfn ("floor (-inf) == -inf", FUNC(floor) (minus_infty));
 
-  check ("floor (pi) == 3", FUNC(floor) (M_PI), 3.0);
-  check ("floor (-pi) == -4", FUNC(floor) (-M_PI), -4.0);
+  check ("floor (pi) == 3", FUNC(floor) (M_PIl), 3.0);
+  check ("floor (-pi) == -4", FUNC(floor) (-M_PIl), -4.0);
 }
 
 
@@ -3440,62 +3442,64 @@ cacos_test (void)
   __complex__ MATHTYPE result;
 
   result = FUNC(cacos) (BUILD_COMPLEX (0, 0));
-  check ("real(cacos(0 + i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0 + i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(cacos(0 + i0)) = -0", __imag__ result, minus_zero);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_zero, 0));
-  check ("real(cacos(-0 + i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(-0 + i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(cacos(-0 + i0)) = -0", __imag__ result, minus_zero);
   result = FUNC(cacos) (BUILD_COMPLEX (0, minus_zero));
-  check ("real(cacos(0 - i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0 - i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(cacos(0 - i0)) = 0", __imag__ result, 0);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_zero, minus_zero));
-  check ("real(cacos(-0 - i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(-0 - i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(cacos(-0 - i0)) = 0", __imag__ result, 0);
 
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, plus_infty));
-  check ("real(cacos(-Inf + i Inf)) = 3*pi/4", __real__ result, M_PI - M_PI_4);
+  check ("real(cacos(-Inf + i Inf)) = 3*pi/4", __real__ result,
+	 M_PIl - M_PI_4l);
   check_isinfn ("imag(cacos(-Inf + i Inf)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, minus_infty));
-  check ("real(cacos(-Inf - i Inf)) = 3*pi/4", __real__ result, M_PI - M_PI_4);
+  check ("real(cacos(-Inf - i Inf)) = 3*pi/4", __real__ result,
+ M_PIl - M_PI_4l);
   check_isinfp ("imag(cacos(-Inf - i Inf)) = +Inf", __imag__ result);
 
   result = FUNC(cacos) (BUILD_COMPLEX (plus_infty, plus_infty));
-  check ("real(cacos(+Inf + i Inf)) = pi/4", __real__ result, M_PI_4);
+  check ("real(cacos(+Inf + i Inf)) = pi/4", __real__ result, M_PI_4l);
   check_isinfn ("imag(cacos(+Inf + i Inf)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (plus_infty, minus_infty));
-  check ("real(cacos(+Inf - i Inf)) = pi/4", __real__ result, M_PI_4);
+  check ("real(cacos(+Inf - i Inf)) = pi/4", __real__ result, M_PI_4l);
   check_isinfp ("imag(cacos(+Inf - i Inf)) = +Inf", __imag__ result);
 
   result = FUNC(cacos) (BUILD_COMPLEX (-10.0, plus_infty));
-  check ("real(cacos(-10.0 + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(-10.0 + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfn ("imag(cacos(-10.0 + i Inf)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (-10.0, minus_infty));
-  check ("real(cacos(-10.0 - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(-10.0 - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfp ("imag(cacos(-10.0 - i Inf)) = +Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (0, plus_infty));
-  check ("real(cacos(0 + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0 + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfn ("imag(cacos(0 + i Inf)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (0, minus_infty));
-  check ("real(cacos(0 - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0 - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfp ("imag(cacos(0 - i Inf)) = +Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (0.1, plus_infty));
-  check ("real(cacos(0.1 + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0.1 + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfn ("imag(cacos(0.1 + i Inf)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (0.1, minus_infty));
-  check ("real(cacos(0.1 - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0.1 - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check_isinfp ("imag(cacos(0.1 - i Inf)) = +Inf", __imag__ result);
 
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, 0));
-  check ("real(cacos(-Inf + i0)) = pi", __real__ result, M_PI);
+  check ("real(cacos(-Inf + i0)) = pi", __real__ result, M_PIl);
   check_isinfn ("imag(cacos(-Inf + i0)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, minus_zero));
-  check ("real(cacos(-Inf - i0)) = pi", __real__ result, M_PI);
+  check ("real(cacos(-Inf - i0)) = pi", __real__ result, M_PIl);
   check_isinfp ("imag(cacos(-Inf - i0)) = +Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, 100));
-  check ("real(cacos(-Inf + i100)) = pi", __real__ result, M_PI);
+  check ("real(cacos(-Inf + i100)) = pi", __real__ result, M_PIl);
   check_isinfn ("imag(cacos(-Inf + i100)) = -Inf", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_infty, -100));
-  check ("real(cacos(-Inf - i100)) = pi", __real__ result, M_PI);
+  check ("real(cacos(-Inf - i100)) = pi", __real__ result, M_PIl);
   check_isinfp ("imag(cacos(-Inf - i100)) = +Inf", __imag__ result);
 
   result = FUNC(cacos) (BUILD_COMPLEX (plus_infty, 0));
@@ -3521,10 +3525,10 @@ cacos_test (void)
 		FUNC(fabs) (__imag__ result));
 
   result = FUNC(cacos) (BUILD_COMPLEX (0, nan_value));
-  check ("real(cacos(0 + i NaN)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(0 + i NaN)) = pi/2", __real__ result, M_PI_2l);
   check_isnan ("imag(cacos(0 + i NaN)) = NaN", __imag__ result);
   result = FUNC(cacos) (BUILD_COMPLEX (minus_zero, nan_value));
-  check ("real(cacos(-0 + i NaN)) = pi/2", __real__ result, M_PI_2);
+  check ("real(cacos(-0 + i NaN)) = pi/2", __real__ result, M_PI_2l);
   check_isnan ("imag(cacos(-0 + i NaN)) = NaN", __imag__ result);
 
   result = FUNC(cacos) (BUILD_COMPLEX (nan_value, plus_infty));
@@ -3581,64 +3585,64 @@ cacosh_test (void)
 
   result = FUNC(cacosh) (BUILD_COMPLEX (0, 0));
   check ("real(cacosh(0 + i0)) = 0", __real__ result, 0);
-  check ("imag(cacosh(0 + i0)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(cacosh(0 + i0)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_zero, 0));
   check ("real(cacosh(-0 + i0)) = 0", __real__ result, 0);
-  check ("imag(cacosh(-0 + i0)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(cacosh(-0 + i0)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (0, minus_zero));
   check ("real(cacosh(0 - i0)) = 0", __real__ result, 0);
-  check ("imag(cacosh(0 - i0)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(cacosh(0 - i0)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_zero, minus_zero));
   check ("real(cacosh(-0 - i0)) = 0", __real__ result, 0);
-  check ("imag(cacosh(-0 - i0)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(cacosh(-0 - i0)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, plus_infty));
   check_isinfp ("real(cacosh(-Inf + i Inf)) = +Inf", __real__ result);
   check ("imag(cacosh(-Inf + i Inf)) = 3*pi/4", __imag__ result,
-	 M_PI - M_PI_4);
+	 M_PIl - M_PI_4l);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, minus_infty));
   check_isinfp ("real(cacosh(-Inf - i Inf)) = +Inf", __real__ result);
   check ("imag(cacosh(-Inf - i Inf)) = -3*pi/4", __imag__ result,
-	 M_PI_4 - M_PI);
+	 M_PI_4l - M_PIl);
 
   result = FUNC(cacosh) (BUILD_COMPLEX (plus_infty, plus_infty));
   check_isinfp ("real(cacosh(+Inf + i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4);
+  check ("imag(cacosh(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4l);
   result = FUNC(cacosh) (BUILD_COMPLEX (plus_infty, minus_infty));
   check_isinfp ("real(cacosh(+Inf - i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4);
+  check ("imag(cacosh(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4l);
 
   result = FUNC(cacosh) (BUILD_COMPLEX (-10.0, plus_infty));
   check_isinfp ("real(cacosh(-10.0 + i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(cacosh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (-10.0, minus_infty));
   check_isinfp ("real(cacosh(-10.0 - i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(cacosh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (0, plus_infty));
   check_isinfp ("real(cacosh(0 + i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(cacosh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (0, minus_infty));
   check_isinfp ("real(cacosh(0 - i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(cacosh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (0.1, plus_infty));
   check_isinfp ("real(cacosh(0.1 + i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(cacosh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(cacosh) (BUILD_COMPLEX (0.1, minus_infty));
   check_isinfp ("real(cacosh(0.1 - i Inf)) = +Inf", __real__ result);
-  check ("imag(cacosh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(cacosh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, 0));
   check_isinfp ("real(cacosh(-Inf + i0)) = +Inf", __real__ result);
-  check ("imag(cacosh(-Inf + i0)) = pi", __imag__ result, M_PI);
+  check ("imag(cacosh(-Inf + i0)) = pi", __imag__ result, M_PIl);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, minus_zero));
   check_isinfp ("real(cacosh(-Inf - i0)) = +Inf", __real__ result);
-  check ("imag(cacosh(-Inf - i0)) = -pi", __imag__ result, -M_PI);
+  check ("imag(cacosh(-Inf - i0)) = -pi", __imag__ result, -M_PIl);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, 100));
   check_isinfp ("real(cacosh(-Inf + i100)) = +Inf", __real__ result);
-  check ("imag(cacosh(-Inf + i100)) = pi", __imag__ result, M_PI);
+  check ("imag(cacosh(-Inf + i100)) = pi", __imag__ result, M_PIl);
   result = FUNC(cacosh) (BUILD_COMPLEX (minus_infty, -100));
   check_isinfp ("real(cacosh(-Inf - i100)) = +Inf", __real__ result);
-  check ("imag(cacosh(-Inf - i100)) = -pi", __imag__ result, -M_PI);
+  check ("imag(cacosh(-Inf - i100)) = -pi", __imag__ result, -M_PIl);
 
   result = FUNC(cacosh) (BUILD_COMPLEX (plus_infty, 0));
   check_isinfp ("real(cacosh(+Inf + i0)) = +Inf", __real__ result);
@@ -3733,16 +3737,16 @@ casin_test (void)
   check ("imag(casin(-0 - i0)) = -0", __imag__ result, minus_zero);
 
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, plus_infty));
-  check ("real(casin(+Inf + i Inf)) = pi/4", __real__ result, M_PI_4);
+  check ("real(casin(+Inf + i Inf)) = pi/4", __real__ result, M_PI_4l);
   check_isinfp ("imag(casin(+Inf + i Inf)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, minus_infty));
-  check ("real(casin(+Inf - i Inf)) = pi/4", __real__ result, M_PI_4);
+  check ("real(casin(+Inf - i Inf)) = pi/4", __real__ result, M_PI_4l);
   check_isinfn ("imag(casin(+Inf - i Inf)) = -Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, plus_infty));
-  check ("real(casin(-Inf + i Inf)) = -pi/4", __real__ result, -M_PI_4);
+  check ("real(casin(-Inf + i Inf)) = -pi/4", __real__ result, -M_PI_4l);
   check_isinfp ("imag(casin(-Inf + i Inf)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, minus_infty));
-  check ("real(casin(-Inf - i Inf)) = -pi/4", __real__ result, -M_PI_4);
+  check ("real(casin(-Inf - i Inf)) = -pi/4", __real__ result, -M_PI_4l);
   check_isinfn ("imag(casin(-Inf - i Inf)) = -Inf", __imag__ result);
 
   result = FUNC(casin) (BUILD_COMPLEX (-10.0, plus_infty));
@@ -3771,29 +3775,29 @@ casin_test (void)
   check_isinfn ("imag(casin(0.1 - i Inf)) = -Inf", __imag__ result);
 
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, 0));
-  check ("real(casin(-Inf + i0)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(casin(-Inf + i0)) = -pi/2", __real__ result, -M_PI_2l);
   check_isinfp ("imag(casin(-Inf + i0)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, minus_zero));
-  check ("real(casin(-Inf - i0)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(casin(-Inf - i0)) = -pi/2", __real__ result, -M_PI_2l);
   check_isinfn ("imag(casin(-Inf - i0)) = -Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, 100));
-  check ("real(casin(-Inf + i100)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(casin(-Inf + i100)) = -pi/2", __real__ result, -M_PI_2l);
   check_isinfp ("imag(casin(-Inf + i100)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (minus_infty, -100));
-  check ("real(casin(-Inf - i100)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(casin(-Inf - i100)) = -pi/2", __real__ result, -M_PI_2l);
   check_isinfn ("imag(casin(-Inf - i100)) = -Inf", __imag__ result);
 
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, 0));
-  check ("real(casin(+Inf + i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(casin(+Inf + i0)) = pi/2", __real__ result, M_PI_2l);
   check_isinfp ("imag(casin(+Inf + i0)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, minus_zero));
-  check ("real(casin(+Inf - i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(casin(+Inf - i0)) = pi/2", __real__ result, M_PI_2l);
   check_isinfn ("imag(casin(+Inf - i0)) = -Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, 0.5));
-  check ("real(casin(+Inf + i0.5)) = pi/2", __real__ result, M_PI_2);
+  check ("real(casin(+Inf + i0.5)) = pi/2", __real__ result, M_PI_2l);
   check_isinfp ("imag(casin(+Inf + i0.5)) = +Inf", __imag__ result);
   result = FUNC(casin) (BUILD_COMPLEX (plus_infty, -0.5));
-  check ("real(casin(+Inf - i0.5)) = pi/2", __real__ result, M_PI_2);
+  check ("real(casin(+Inf - i0.5)) = pi/2", __real__ result, M_PI_2l);
   check_isinfn ("imag(casin(+Inf - i0.5)) = -Inf", __imag__ result);
 
   result = FUNC(casin) (BUILD_COMPLEX (nan_value, plus_infty));
@@ -3879,41 +3883,41 @@ casinh_test (void)
 
   result = FUNC(casinh) (BUILD_COMPLEX (plus_infty, plus_infty));
   check_isinfp ("real(casinh(+Inf + i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4);
+  check ("imag(casinh(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4l);
   result = FUNC(casinh) (BUILD_COMPLEX (plus_infty, minus_infty));
   check_isinfp ("real(casinh(+Inf - i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4);
+  check ("imag(casinh(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4l);
   result = FUNC(casinh) (BUILD_COMPLEX (minus_infty, plus_infty));
   check_isinfn ("real(casinh(-Inf + i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-Inf + i Inf)) = pi/4", __imag__ result, M_PI_4);
+  check ("imag(casinh(-Inf + i Inf)) = pi/4", __imag__ result, M_PI_4l);
   result = FUNC(casinh) (BUILD_COMPLEX (minus_infty, minus_infty));
   check_isinfn ("real(casinh(-Inf - i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4);
+  check ("imag(casinh(-Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4l);
 
   result = FUNC(casinh) (BUILD_COMPLEX (-10.0, plus_infty));
   check_isinfn ("real(casinh(-10.0 + i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(casinh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (-10.0, minus_infty));
   check_isinfn ("real(casinh(-10.0 - i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(casinh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (0, plus_infty));
   check_isinfp ("real(casinh(0 + i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(casinh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (0, minus_infty));
   check_isinfp ("real(casinh(0 - i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(casinh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (minus_zero, plus_infty));
   check_isinfn ("real(casinh(-0 + i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(casinh(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (minus_zero, minus_infty));
   check_isinfn ("real(casinh(-0 - i Inf)) = -Inf", __real__ result);
-  check ("imag(casinh(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(casinh(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (0.1, plus_infty));
   check_isinfp ("real(casinh(0.1 + i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(casinh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(casinh) (BUILD_COMPLEX (0.1, minus_infty));
   check_isinfp ("real(casinh(0.1 - i Inf)) = +Inf", __real__ result);
-  check ("imag(casinh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(casinh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(casinh) (BUILD_COMPLEX (minus_infty, 0));
   check_isinfn ("real(casinh(-Inf + i0)) = -Inf", __real__ result);
@@ -4023,67 +4027,67 @@ catan_test (void)
   check ("imag(catan(-0 - i0)) = -0", __imag__ result, minus_zero);
 
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, plus_infty));
-  check ("real(catan(+Inf + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(+Inf + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(+Inf + i Inf)) = 0", __imag__ result, 0);
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, minus_infty));
-  check ("real(catan(+Inf - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(+Inf - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(+Inf - i Inf)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, plus_infty));
-  check ("real(catan(-Inf + i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf + i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf + i Inf)) = 0", __imag__ result, 0.0);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, minus_infty));
-  check ("real(catan(-Inf - i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf - i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf - i Inf)) = -0", __imag__ result, minus_zero);
 
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, -10.0));
-  check ("real(catan(+Inf - i10.0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(+Inf - i10.0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(+Inf - i10.0)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, -10.0));
-  check ("real(catan(-Inf - i10.0)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf - i10.0)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf - i10.0)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, minus_zero));
-  check ("real(catan(Inf - i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(Inf - i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(Inf - i0)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, minus_zero));
-  check ("real(catan(-Inf - i0)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf - i0)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf - i0)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, 0.0));
-  check ("real(catan(Inf + i0)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(Inf + i0)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(Inf + i0)) = 0", __imag__ result, 0.0);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, 0.0));
-  check ("real(catan(-Inf + i0)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf + i0)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf + i0)) = 0", __imag__ result, 0.0);
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, 0.1));
-  check ("real(catan(+Inf + i0.1)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(+Inf + i0.1)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(+Inf + i0.1)) = 0", __imag__ result, 0);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, 0.1));
-  check ("real(catan(-Inf + i0.1)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf + i0.1)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf + i0.1)) = 0", __imag__ result, 0);
 
   result = FUNC(catan) (BUILD_COMPLEX (0.0, minus_infty));
-  check ("real(catan(0 - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(0 - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(0 - i Inf)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (minus_zero, minus_infty));
-  check ("real(catan(-0 - i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-0 - i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-0 - i Inf)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (100.0, minus_infty));
-  check ("real(catan(100 - i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(100 - i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(100 - i Inf)) = -0", __imag__ result, minus_zero);
   result = FUNC(catan) (BUILD_COMPLEX (-100.0, minus_infty));
-  check ("real(catan(-100 - i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-100 - i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-100 - i Inf)) = -0", __imag__ result, minus_zero);
 
   result = FUNC(catan) (BUILD_COMPLEX (0.0, plus_infty));
-  check ("real(catan(0 + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(0 + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(0 + i Inf)) = 0", __imag__ result, 0);
   result = FUNC(catan) (BUILD_COMPLEX (minus_zero, plus_infty));
-  check ("real(catan(-0 + i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-0 + i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-0 + i Inf)) = 0", __imag__ result, 0);
   result = FUNC(catan) (BUILD_COMPLEX (0.5, plus_infty));
-  check ("real(catan(0.5 + i Inf)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(0.5 + i Inf)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(0.5 + i Inf)) = 0", __imag__ result, 0);
   result = FUNC(catan) (BUILD_COMPLEX (-0.5, plus_infty));
-  check ("real(catan(-0.5 + i Inf)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-0.5 + i Inf)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-0.5 + i Inf)) = 0", __imag__ result, 0);
 
   result = FUNC(catan) (BUILD_COMPLEX (nan_value, 0.0));
@@ -4108,10 +4112,10 @@ catan_test (void)
   check_isnan ("imag(catan(-0 + i NaN)) = NaN", __imag__ result);
 
   result = FUNC(catan) (BUILD_COMPLEX (plus_infty, nan_value));
-  check ("real(catan(+Inf + i NaN)) = pi/2", __real__ result, M_PI_2);
+  check ("real(catan(+Inf + i NaN)) = pi/2", __real__ result, M_PI_2l);
   check ("imag(catan(+Inf + i NaN)) = +-0", FUNC(fabs) (__imag__ result), 0);
   result = FUNC(catan) (BUILD_COMPLEX (minus_infty, nan_value));
-  check ("real(catan(-Inf + i NaN)) = -pi/2", __real__ result, -M_PI_2);
+  check ("real(catan(-Inf + i NaN)) = -pi/2", __real__ result, -M_PI_2l);
   check ("imag(catan(-Inf + i NaN)) = +-0", FUNC(fabs) (__imag__ result), 0);
 
   result = FUNC(catan) (BUILD_COMPLEX (nan_value, 10.5));
@@ -4174,67 +4178,67 @@ catanh_test (void)
 
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, plus_infty));
   check ("real(catanh(+Inf + i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(+Inf + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, minus_infty));
   check ("real(catanh(+Inf - i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(+Inf - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, plus_infty));
   check ("real(catanh(-Inf + i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(-Inf + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, minus_infty));
   check ("real(catanh(-Inf - i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(-Inf - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(catanh) (BUILD_COMPLEX (-10.0, plus_infty));
   check ("real(catanh(-10.0 + i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(-10.0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (-10.0, minus_infty));
   check ("real(catanh(-10.0 - i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(-10.0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_zero, plus_infty));
   check ("real(catanh(-0 + i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_zero, minus_infty));
   check ("real(catanh(-0 - i Inf)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (0, plus_infty));
   check ("real(catanh(0 + i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (0, minus_infty));
   check ("real(catanh(0 - i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (0.1, plus_infty));
   check ("real(catanh(0.1 + i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(0.1 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (0.1, minus_infty));
   check ("real(catanh(0.1 - i Inf)) = 0", __real__ result, 0);
-  check ("imag(catanh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(0.1 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, 0));
   check ("real(catanh(-Inf + i0)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf + i0)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(-Inf + i0)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, minus_zero));
   check ("real(catanh(-Inf - i0)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf - i0)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(-Inf - i0)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, 100));
   check ("real(catanh(-Inf + i100)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf + i100)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(-Inf + i100)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (minus_infty, -100));
   check ("real(catanh(-Inf - i100)) = -0", __real__ result, minus_zero);
-  check ("imag(catanh(-Inf - i100)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(-Inf - i100)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, 0));
   check ("real(catanh(+Inf + i0)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf + i0)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(+Inf + i0)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, minus_zero));
   check ("real(catanh(+Inf - i0)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf - i0)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(+Inf - i0)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, 0.5));
   check ("real(catanh(+Inf + i0.5)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf + i0.5)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(+Inf + i0.5)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (plus_infty, -0.5));
   check ("real(catanh(+Inf - i0.5)) = 0", __real__ result, 0);
-  check ("imag(catanh(+Inf - i0.5)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(+Inf - i0.5)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(catanh) (BUILD_COMPLEX (0, nan_value));
   check ("real(catanh(0 + i NaN)) = 0", __real__ result, 0);
@@ -4259,10 +4263,10 @@ catanh_test (void)
 
   result = FUNC(catanh) (BUILD_COMPLEX (nan_value, plus_infty));
   check ("real(catanh(NaN + i Inf)) = +-0", FUNC(fabs) (__real__ result), 0);
-  check ("imag(catanh(NaN + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(catanh(NaN + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(catanh) (BUILD_COMPLEX (nan_value, minus_infty));
   check ("real(catanh(NaN - i Inf)) = +-0", FUNC(fabs) (__real__ result), 0);
-  check ("imag(catanh(NaN - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(catanh(NaN - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(catanh) (BUILD_COMPLEX (10.5, nan_value));
   check_isnan_maybe_exc ("real(catanh(10.5 + i NaN)) = NaN plus maybe invalid exception",
@@ -4588,7 +4592,7 @@ ctanh_test (void)
   check_isnan ("real(ctanh(NaN + i NaN)) = NaN", __real__ result);
   check_isnan ("imag(ctanh(NaN + i NaN)) = NaN", __imag__ result);
 
-  result = FUNC(ctanh) (BUILD_COMPLEX (0, M_PI_4));
+  result = FUNC(ctanh) (BUILD_COMPLEX (0, M_PI_4l));
   check ("real(ctanh (0 + i pi/4)) == 0", __real__ result, 0);
   check_eps ("imag(ctanh (0 + i pi/4)) == 1", __imag__ result, 1,
 	     CHOOSE (0, 2e-16, 2e-7));
@@ -4616,12 +4620,12 @@ clog_test (void)
   check_isinfn_exc ("real(clog(-0 + i0)) = -Inf plus divide-by-zero exception",
 		    __real__ result, DIVIDE_BY_ZERO_EXCEPTION);
   check ("imag(clog(-0 + i0)) = pi plus divide-by-zero exception",
-	 __imag__ result, M_PI);
+	 __imag__ result, M_PIl);
   result = FUNC(clog) (BUILD_COMPLEX (minus_zero, minus_zero));
   check_isinfn_exc ("real(clog(-0 - i0)) = -Inf plus divide-by-zero exception",
 		    __real__ result, DIVIDE_BY_ZERO_EXCEPTION);
   check ("imag(clog(-0 - i0)) = -pi plus divide-by-zero exception",
-	 __imag__ result, -M_PI);
+	 __imag__ result, -M_PIl);
 
   result = FUNC(clog) (BUILD_COMPLEX (0, 0));
   check_isinfn_exc ("real(clog(0 + i0)) = -Inf plus divide-by-zero exception",
@@ -4636,55 +4640,57 @@ clog_test (void)
 
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, plus_infty));
   check_isinfp ("real(clog(-Inf + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf + i Inf)) = 3*pi/4", __imag__ result, M_PI - M_PI_4);
+  check ("imag(clog(-Inf + i Inf)) = 3*pi/4", __imag__ result,
+	 M_PIl - M_PI_4l);
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, minus_infty));
   check_isinfp ("real(clog(-Inf - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf - i Inf)) = -3*pi/4", __imag__ result, M_PI_4 - M_PI);
+  check ("imag(clog(-Inf - i Inf)) = -3*pi/4", __imag__ result,
+	 M_PI_4l - M_PIl);
 
   result = FUNC(clog) (BUILD_COMPLEX (plus_infty, plus_infty));
   check_isinfp ("real(clog(+Inf + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4);
+  check ("imag(clog(+Inf + i Inf)) = pi/4", __imag__ result, M_PI_4l);
   result = FUNC(clog) (BUILD_COMPLEX (plus_infty, minus_infty));
   check_isinfp ("real(clog(+Inf - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4);
+  check ("imag(clog(+Inf - i Inf)) = -pi/4", __imag__ result, -M_PI_4l);
 
   result = FUNC(clog) (BUILD_COMPLEX (0, plus_infty));
   check_isinfp ("real(clog(0 + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(clog(0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (3, plus_infty));
   check_isinfp ("real(clog(3 + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(3 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(clog(3 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (minus_zero, plus_infty));
   check_isinfp ("real(clog(-0 + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(clog(-0 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (-3, plus_infty));
   check_isinfp ("real(clog(-3 + i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-3 + i Inf)) = pi/2", __imag__ result, M_PI_2);
+  check ("imag(clog(-3 + i Inf)) = pi/2", __imag__ result, M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (0, minus_infty));
   check_isinfp ("real(clog(0 - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(clog(0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (3, minus_infty));
   check_isinfp ("real(clog(3 - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(3 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(clog(3 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (minus_zero, minus_infty));
   check_isinfp ("real(clog(-0 - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(clog(-0 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
   result = FUNC(clog) (BUILD_COMPLEX (-3, minus_infty));
   check_isinfp ("real(clog(-3 - i Inf)) = +Inf", __real__ result);
-  check ("imag(clog(-3 - i Inf)) = -pi/2", __imag__ result, -M_PI_2);
+  check ("imag(clog(-3 - i Inf)) = -pi/2", __imag__ result, -M_PI_2l);
 
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, 0));
   check_isinfp ("real(clog(-Inf + i0)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf + i0)) = pi", __imag__ result, M_PI);
+  check ("imag(clog(-Inf + i0)) = pi", __imag__ result, M_PIl);
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, 1));
   check_isinfp ("real(clog(-Inf + i1)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf + i1)) = pi", __imag__ result, M_PI);
+  check ("imag(clog(-Inf + i1)) = pi", __imag__ result, M_PIl);
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, minus_zero));
   check_isinfp ("real(clog(-Inf - i0)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf - i0)) = -pi", __imag__ result, -M_PI);
+  check ("imag(clog(-Inf - i0)) = -pi", __imag__ result, -M_PIl);
   result = FUNC(clog) (BUILD_COMPLEX (minus_infty, -1));
   check_isinfp ("real(clog(-Inf - i1)) = +Inf", __real__ result);
-  check ("imag(clog(-Inf - i1)) = -pi", __imag__ result, -M_PI);
+  check ("imag(clog(-Inf - i1)) = -pi", __imag__ result, -M_PIl);
 
   result = FUNC(clog) (BUILD_COMPLEX (plus_infty, 0));
   check_isinfp ("real(clog(+Inf + i0)) = +Inf", __real__ result);
@@ -4782,12 +4788,12 @@ clog10_test (void)
   check_isinfn_exc ("real(clog10(-0 + i0)) = -Inf plus divide-by-zero exception",
 		    __real__ result, DIVIDE_BY_ZERO_EXCEPTION);
   check ("imag(clog10(-0 + i0)) = pi plus divide-by-zero exception",
-	 __imag__ result, M_PI);
+	 __imag__ result, M_PIl);
   result = FUNC(clog10) (BUILD_COMPLEX (minus_zero, minus_zero));
   check_isinfn_exc ("real(clog10(-0 - i0)) = -Inf plus divide-by-zero exception",
 		    __real__ result, DIVIDE_BY_ZERO_EXCEPTION);
   check ("imag(clog10(-0 - i0)) = -pi plus divide-by-zero exception",
-	 __imag__ result, -M_PI);
+	 __imag__ result, -M_PIl);
 
   result = FUNC(clog10) (BUILD_COMPLEX (0, 0));
   check_isinfn_exc ("real(clog10(0 + i0)) = -Inf plus divide-by-zero exception",
@@ -4803,70 +4809,70 @@ clog10_test (void)
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, plus_infty));
   check_isinfp ("real(clog10(-Inf + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf + i Inf)) = 3*pi/4*M_LOG10E", __imag__ result,
-	     (M_PI - M_PI_4) * M_LOG10E, CHOOSE (0, 3e-16, 0));
+	     (M_PIl - M_PI_4l) * M_LOG10El, CHOOSE (0, 3e-16, 0));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, minus_infty));
   check_isinfp ("real(clog10(-Inf - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf - i Inf)) = -3*pi/4*M_LOG10E", __imag__ result,
-	     (M_PI_4 - M_PI) * M_LOG10E, CHOOSE (0, 3e-16, 0));
+	     (M_PI_4l - M_PIl) * M_LOG10El, CHOOSE (0, 3e-16, 0));
 
   result = FUNC(clog10) (BUILD_COMPLEX (plus_infty, plus_infty));
   check_isinfp ("real(clog10(+Inf + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(+Inf + i Inf)) = pi/4*M_LOG10E", __imag__ result,
-	     M_PI_4 * M_LOG10E, CHOOSE (0, 6e-17, 3e-8));
+	     M_PI_4l * M_LOG10El, CHOOSE (0, 6e-17, 3e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (plus_infty, minus_infty));
   check_isinfp ("real(clog10(+Inf - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(+Inf - i Inf)) = -pi/4*M_LOG10E", __imag__ result,
-	     -M_PI_4 * M_LOG10E, CHOOSE (0, 6e-17, 3e-8));
+	     -M_PI_4l * M_LOG10El, CHOOSE (0, 6e-17, 3e-8));
 
   result = FUNC(clog10) (BUILD_COMPLEX (0, plus_infty));
   check_isinfp ("real(clog10(0 + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(0 + i Inf)) = pi/2*M_LOG10E", __imag__ result,
-	     M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (3, plus_infty));
   check_isinfp ("real(clog10(3 + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(3 + i Inf)) = pi/2*M_LOG10E", __imag__ result,
-	     M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_zero, plus_infty));
   check_isinfp ("real(clog10(-0 + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-0 + i Inf)) = pi/2*M_LOG10E", __imag__ result,
-	     M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (-3, plus_infty));
   check_isinfp ("real(clog10(-3 + i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-3 + i Inf)) = pi/2*M_LOG10E", __imag__ result,
-	     M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (0, minus_infty));
   check_isinfp ("real(clog10(0 - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(0 - i Inf)) = -pi/2*M_LOG10E", __imag__ result,
-	     -M_PI_2*M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     -M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (3, minus_infty));
   check_isinfp ("real(clog10(3 - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(3 - i Inf)) = -pi/2*M_LOG10E", __imag__ result,
-	     -M_PI_2*M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     -M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_zero, minus_infty));
   check_isinfp ("real(clog10(-0 - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-0 - i Inf)) = -pi/2*M_LOG10E", __imag__ result,
-	     -M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     -M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
   result = FUNC(clog10) (BUILD_COMPLEX (-3, minus_infty));
   check_isinfp ("real(clog10(-3 - i Inf)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-3 - i Inf)) = -pi/2*M_LOG10E", __imag__ result,
-	     -M_PI_2 * M_LOG10E, CHOOSE (0, 2e-16, 6e-8));
+	     -M_PI_2l * M_LOG10El, CHOOSE (0, 2e-16, 6e-8));
 
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, 0));
   check_isinfp ("real(clog10(-Inf + i0)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf + i0)) = pi*M_LOG10E", __imag__ result,
-	     M_PI * M_LOG10E, CHOOSE (0, 3e-16, 2e-7));
+	     M_PIl * M_LOG10El, CHOOSE (0, 3e-16, 2e-7));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, 1));
   check_isinfp ("real(clog10(-Inf + i1)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf + i1)) = pi*M_LOG10E", __imag__ result,
-	     M_PI * M_LOG10E, CHOOSE (0, 3e-16, 2e-7));
+	     M_PIl * M_LOG10El, CHOOSE (0, 3e-16, 2e-7));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, minus_zero));
   check_isinfp ("real(clog10(-Inf - i0)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf - i0)) = -pi*M_LOG10E", __imag__ result,
-	     -M_PI * M_LOG10E, CHOOSE (0, 3e-16, 2e-7));
+	     -M_PIl * M_LOG10El, CHOOSE (0, 3e-16, 2e-7));
   result = FUNC(clog10) (BUILD_COMPLEX (minus_infty, -1));
   check_isinfp ("real(clog10(-Inf - i1)) = +Inf", __real__ result);
   check_eps ("imag(clog10(-Inf - i1)) = -pi*M_LOG10E", __imag__ result,
-	     -M_PI * M_LOG10E, CHOOSE (0, 3e-16, 2e-7));
+	     -M_PIl * M_LOG10El, CHOOSE (0, 3e-16, 2e-7));
 
   result = FUNC(clog10) (BUILD_COMPLEX (plus_infty, 0));
   check_isinfp ("real(clog10(+Inf + i0)) = +Inf", __real__ result);
@@ -5135,7 +5141,7 @@ cpow_test (void)
 	     CHOOSE (2e-16L, 0, 0));
   check ("imag(cpow (2 + i0), (10 + i0)) == 0", __imag__ result, 0);
 
-  result = FUNC (cpow) (BUILD_COMPLEX (M_E, 0), BUILD_COMPLEX (0, 2*M_PI));
+  result = FUNC (cpow) (BUILD_COMPLEX (M_El, 0), BUILD_COMPLEX (0, 2 * M_PIl));
   check_eps ("real(cpow (e + i0), (0 + i 2*PI)) == 1", __real__ result, 1,
 	     CHOOSE (0, 0, 6e-8));
   check_eps ("imag(cpow (e + i0), (0 + i 2*PI)) == 0", __imag__ result, 0,
@@ -5212,32 +5218,32 @@ carg_test (void)
 
   x = -random_greater (0);
   check ("carg (x + i 0) == +pi for x < 0", FUNC(carg) (BUILD_COMPLEX(x, 0)),
-	 M_PI);
+	 M_PIl);
 
   x = -random_greater (0);
   check ("carg (x - i 0) == -pi for x < 0",
-	 FUNC(carg) (BUILD_COMPLEX(x, minus_zero)), -M_PI);
+	 FUNC(carg) (BUILD_COMPLEX(x, minus_zero)), -M_PIl);
 
   check ("carg (-0 + i 0) == +pi", FUNC(carg) (BUILD_COMPLEX(minus_zero, 0)),
-	 M_PI);
+	 M_PIl);
   check ("carg (-0 - i 0) == -pi",
-	 FUNC(carg) (BUILD_COMPLEX(minus_zero, minus_zero)), -M_PI);
+	 FUNC(carg) (BUILD_COMPLEX(minus_zero, minus_zero)), -M_PIl);
 
   x = random_greater (0);
   check ("carg (+0 + i y) == pi/2 for y > 0", FUNC(carg) (BUILD_COMPLEX(0, x)),
-	 M_PI_2);
+	 M_PI_2l);
 
   x = random_greater (0);
   check ("carg (-0 + i y) == pi/2 for y > 0",
-	 FUNC(carg) (BUILD_COMPLEX(minus_zero, x)), M_PI_2);
+	 FUNC(carg) (BUILD_COMPLEX(minus_zero, x)), M_PI_2l);
 
   x = random_less (0);
   check ("carg (+0 + i y) == -pi/2 for y < 0", FUNC(carg) (BUILD_COMPLEX(0, x)),
-	 -M_PI_2);
+	 -M_PI_2l);
 
   x = random_less (0);
   check ("carg (-0 + i y) == -pi/2 for y < 0",
-	 FUNC(carg) (BUILD_COMPLEX(minus_zero, x)), -M_PI_2);
+	 FUNC(carg) (BUILD_COMPLEX(minus_zero, x)), -M_PI_2l);
 
   x = random_greater (0);
   check ("carg (inf + i y) == +0 for finite y > 0",
@@ -5249,31 +5255,31 @@ carg_test (void)
 
   x = random_value (-1e4, 1e4);
   check ("carg(x + i inf) == pi/2 for finite x",
-         FUNC(carg) (BUILD_COMPLEX(x, plus_infty)), M_PI_2);
+         FUNC(carg) (BUILD_COMPLEX(x, plus_infty)), M_PI_2l);
 
   x = random_value (-1e4, 1e4);
   check ("carg(x - i inf) == -pi/2 for finite x",
-         FUNC(carg) (BUILD_COMPLEX(x, minus_infty)), -M_PI_2);
+         FUNC(carg) (BUILD_COMPLEX(x, minus_infty)), -M_PI_2l);
 
   x = random_greater (0);
   check ("carg (-inf + i y) == +pi for finite y > 0",
-	 FUNC(carg) (BUILD_COMPLEX(minus_infty, x)), M_PI);
+	 FUNC(carg) (BUILD_COMPLEX(minus_infty, x)), M_PIl);
 
   x = -random_greater (0);
   check ("carg (-inf + i y) == -pi for finite y < 0",
-	 FUNC(carg) (BUILD_COMPLEX(minus_infty, x)), -M_PI);
+	 FUNC(carg) (BUILD_COMPLEX(minus_infty, x)), -M_PIl);
 
   check ("carg (+inf + i inf) == +pi/4",
-	 FUNC(carg) (BUILD_COMPLEX(plus_infty, plus_infty)), M_PI_4);
+	 FUNC(carg) (BUILD_COMPLEX(plus_infty, plus_infty)), M_PI_4l);
 
   check ("carg (+inf -i inf) == -pi/4",
-	 FUNC(carg) (BUILD_COMPLEX(plus_infty, minus_infty)), -M_PI_4);
+	 FUNC(carg) (BUILD_COMPLEX(plus_infty, minus_infty)), -M_PI_4l);
 
   check ("carg (-inf +i inf) == +3*pi/4",
-	 FUNC(carg) (BUILD_COMPLEX(minus_infty, plus_infty)), 3 * M_PI_4);
+	 FUNC(carg) (BUILD_COMPLEX(minus_infty, plus_infty)), 3 * M_PI_4l);
 
   check ("carg (-inf -i inf) == -3*pi/4",
-	 FUNC(carg) (BUILD_COMPLEX(minus_infty, minus_infty)), -3 * M_PI_4);
+	 FUNC(carg) (BUILD_COMPLEX(minus_infty, minus_infty)), -3 * M_PI_4l);
 
 }
 
diff --git a/math/math.h b/math/math.h
index ae83262338..5103f2ef3c 100644
--- a/math/math.h
+++ b/math/math.h
@@ -1,5 +1,5 @@
 /* Declarations for math functions.
-   Copyright (C) 1991, 92, 93, 95, 96, 97 Free Software Foundation, Inc.
+   Copyright (C) 1991, 92, 93, 95, 96, 97, 98 Free Software Foundation, Inc.
    This file is part of the GNU C Library.
 
    The GNU C Library is free software; you can redistribute it and/or
@@ -262,41 +262,41 @@ extern int matherr __P ((struct exception *__exc));
 #endif	/* SVID */
 
 
-#if defined __USE_BSD || defined __USE_UNIX98
-
 /* Some useful constants.  */
-# define M_E		_Mldbl(2.7182818284590452354)	/* e */
-# define M_LOG2E	_Mldbl(1.4426950408889634074)	/* log_2 e */
-# define M_LOG10E	_Mldbl(0.43429448190325182765)	/* log_10 e */
-# define M_LN2		_Mldbl(0.69314718055994530942)	/* log_e 2 */
-# define M_LN10		_Mldbl(2.30258509299404568402)	/* log_e 10 */
-# define M_PI		_Mldbl(3.14159265358979323846)	/* pi */
-# define M_PI_2		_Mldbl(1.57079632679489661923)	/* pi/2 */
-# define M_PI_4		_Mldbl(0.78539816339744830962)	/* pi/4 */
-# define M_1_PI		_Mldbl(0.31830988618379067154)	/* 1/pi */
-# define M_2_PI		_Mldbl(0.63661977236758134308)	/* 2/pi */
-# define M_2_SQRTPI	_Mldbl(1.12837916709551257390)	/* 2/sqrt(pi) */
-# define M_SQRT2	_Mldbl(1.41421356237309504880)	/* sqrt(2) */
-# define M_SQRT1_2	_Mldbl(0.70710678118654752440)	/* 1/sqrt(2) */
-
+#if defined __USE_BSD || defined __USE_UNIX98
+# define M_E		2.7182818284590452354	/* e */
+# define M_LOG2E	1.4426950408889634074	/* log_2 e */
+# define M_LOG10E	0.43429448190325182765	/* log_10 e */
+# define M_LN2		0.69314718055994530942	/* log_e 2 */
+# define M_LN10		2.30258509299404568402	/* log_e 10 */
+# define M_PI		3.14159265358979323846	/* pi */
+# define M_PI_2		1.57079632679489661923	/* pi/2 */
+# define M_PI_4		0.78539816339744830962	/* pi/4 */
+# define M_1_PI		0.31830988618379067154	/* 1/pi */
+# define M_2_PI		0.63661977236758134308	/* 2/pi */
+# define M_2_SQRTPI	1.12837916709551257390	/* 2/sqrt(pi) */
+# define M_SQRT2	1.41421356237309504880	/* sqrt(2) */
+# define M_SQRT1_2	0.70710678118654752440	/* 1/sqrt(2) */
 #endif
 
-/* Our constants might specify more precision than `double' can represent.
-   Use `long double' constants in standard and GNU C, where they are
-   supported and the cast to `double'.
-
-   If the constants are use in code which does not use prototypes, one
-   might get problems if a function takes a `double' argument and any
-   of the constants are provided as the argument.  In this case, cast
-   the argument to `double'.
-
-   Please note we define the macro even if the constants are not defined.
-   This helps us to use the macros in other places.  */
-#if (__STDC__ - 0 || __GNUC__ - 0) && defined __USE_GNU
-# define _Mldbl(x) x##L
-#else	/* Traditional C.  */
-# define _Mldbl(x) x
-#endif	/* Standard or GNU C.  */
+/* The above constants are not adequate for computation using `long double's.
+   Therefore we provide as an extension constants with similar names as a
+   GNU extension.  */
+#ifdef __USE_GNU
+# define M_El		2.7182818284590452354L	/* e */
+# define M_LOG2El	1.4426950408889634074L	/* log_2 e */
+# define M_LOG10El	0.43429448190325182765L	/* log_10 e */
+# define M_LN2l		0.69314718055994530942L	/* log_e 2 */
+# define M_LN10l	2.30258509299404568402L	/* log_e 10 */
+# define M_PIl		3.14159265358979323846L	/* pi */
+# define M_PI_2l	1.57079632679489661923L	/* pi/2 */
+# define M_PI_4l	0.78539816339744830962L	/* pi/4 */
+# define M_1_PIl	0.31830988618379067154L	/* 1/pi */
+# define M_2_PIl	0.63661977236758134308L	/* 2/pi */
+# define M_2_SQRTPIl	1.12837916709551257390L	/* 2/sqrt(pi) */
+# define M_SQRT2l	1.41421356237309504880L	/* sqrt(2) */
+# define M_SQRT1_2l	0.70710678118654752440L	/* 1/sqrt(2) */
+#endif
 
 
 /* Get machine-dependent inline versions (if there are any).  */