diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2009-03-03 16:51:07 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 2009-03-03 16:51:07 +0000 |
commit | f5f7883fc7cf47d38ed493a882fd7109bdd589bb (patch) | |
tree | 1058cf0d07593fa2cec9f32b4988ac40611b4bff /tests/tfms.c | |
parent | cbd40e1cb659b579e608afbe7a05da5b15a860be (diff) | |
download | mpfr-f5f7883fc7cf47d38ed493a882fd7109bdd589bb.tar.gz |
GMP_RNDX -> MPFR_RNDX
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@6053 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'tests/tfms.c')
-rw-r--r-- | tests/tfms.c | 134 |
1 files changed, 67 insertions, 67 deletions
diff --git a/tests/tfms.c b/tests/tfms.c index b7d9d4827..0cfb3688a 100644 --- a/tests/tfms.c +++ b/tests/tfms.c @@ -43,9 +43,9 @@ test_exact (void) for (k = 0; k < sv; k++) RND_LOOP (rnd) { - if (mpfr_set_str (a, val[i], 10, GMP_RNDN) || - mpfr_set_str (b, val[j], 10, GMP_RNDN) || - mpfr_set_str (c, val[k], 10, GMP_RNDN) || + if (mpfr_set_str (a, val[i], 10, MPFR_RNDN) || + mpfr_set_str (b, val[j], 10, MPFR_RNDN) || + mpfr_set_str (c, val[k], 10, MPFR_RNDN) || mpfr_mul (r1, a, b, (mp_rnd_t) rnd) || mpfr_sub (r1, r1, c, (mp_rnd_t) rnd)) { @@ -70,9 +70,9 @@ test_exact (void) if (mpfr_cmp (r1, r2) || MPFR_SIGN (r1) != MPFR_SIGN (r2)) { printf ("test_exact(%d,%d,%d,%d):\nexpected ", i, j, k, rnd); - mpfr_out_str (stdout, 10, 0, r1, GMP_RNDN); + mpfr_out_str (stdout, 10, 0, r1, MPFR_RNDN); printf ("\n got "); - mpfr_out_str (stdout, 10, 0, r2, GMP_RNDN); + mpfr_out_str (stdout, 10, 0, r2, MPFR_RNDN); printf ("\n"); exit (1); } @@ -90,18 +90,18 @@ test_overflow1 (void) mpfr_inits2 (8, x, y, z, r, (mpfr_ptr) 0); MPFR_SET_POS (x); mpfr_setmax (x, mpfr_get_emax ()); /* x = 2^emax - ulp */ - mpfr_set_ui (y, 2, GMP_RNDN); /* y = 2 */ - mpfr_set (z, x, GMP_RNDN); /* z = x = 2^emax - ulp */ + mpfr_set_ui (y, 2, MPFR_RNDN); /* y = 2 */ + mpfr_set (z, x, MPFR_RNDN); /* z = x = 2^emax - ulp */ mpfr_clear_flags (); /* The intermediate multiplication x * y overflows, but x * y - z = x is representable. */ - inex = mpfr_fms (r, x, y, z, GMP_RNDN); + inex = mpfr_fms (r, x, y, z, MPFR_RNDN); if (inex || ! mpfr_equal_p (r, x)) { printf ("Error in test_overflow1\nexpected "); - mpfr_out_str (stdout, 2, 0, x, GMP_RNDN); + mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); printf (" with inex = 0\n got "); - mpfr_out_str (stdout, 2, 0, r, GMP_RNDN); + mpfr_out_str (stdout, 2, 0, r, MPFR_RNDN); printf (" with inex = %d\n", inex); exit (1); } @@ -123,7 +123,7 @@ test_overflow2 (void) MPFR_SET_POS (x); mpfr_setmin (x, mpfr_get_emax ()); /* x = 0.1@emax */ - mpfr_set_si (y, -2, GMP_RNDN); /* y = -2 */ + mpfr_set_si (y, -2, MPFR_RNDN); /* y = -2 */ /* The intermediate multiplication x * y will overflow. */ for (i = -9; i <= 9; i++) @@ -131,10 +131,10 @@ test_overflow2 (void) { int inf, overflow; - inf = rnd == GMP_RNDN || rnd == GMP_RNDD || rnd == GMP_RNDA; + inf = rnd == MPFR_RNDN || rnd == MPFR_RNDD || rnd == MPFR_RNDA; overflow = inf || i <= 0; - inex = mpfr_set_si_2exp (z, -i, mpfr_get_emin (), GMP_RNDN); + inex = mpfr_set_si_2exp (z, -i, mpfr_get_emin (), MPFR_RNDN); MPFR_ASSERTN (inex == 0); mpfr_clear_flags (); @@ -209,12 +209,12 @@ test_underflow1 (void) for (signy = -1; signy <= 1; signy += 2) { - mpfr_set_si_2exp (y, signy, -1, GMP_RNDN); /* |y| = 1/2 */ + mpfr_set_si_2exp (y, signy, -1, MPFR_RNDN); /* |y| = 1/2 */ for (signz = -3; signz <= 3; signz += 2) { RND_LOOP (rnd) { - mpfr_set_si (z, signz, GMP_RNDN); + mpfr_set_si (z, signz, MPFR_RNDN); if (ABS (signz) != 1) mpfr_setmax (z, mpfr_get_emax ()); /* |z| = 1 or 2^emax - ulp */ @@ -227,7 +227,7 @@ test_underflow1 (void) mpfr_print_rnd_mode ((mp_rnd_t) rnd)); err = 1; } - mpfr_neg (z, z, GMP_RNDN); + mpfr_neg (z, z, MPFR_RNDN); if (signy < 0 && MPFR_IS_LIKE_RNDD(rnd, -signz)) mpfr_nextbelow (z); if (signy > 0 && MPFR_IS_LIKE_RNDU(rnd, -signz)) @@ -254,18 +254,18 @@ test_underflow1 (void) printf ("\n"); err = 1; } - if (inex >= 0 && (rnd == GMP_RNDD || - (rnd == GMP_RNDZ && signz < 0) || - (rnd == GMP_RNDN && signy > 0))) + if (inex >= 0 && (rnd == MPFR_RNDD || + (rnd == MPFR_RNDZ && signz < 0) || + (rnd == MPFR_RNDN && signy > 0))) { printf (ERRTU1 "ternary value = %d instead of < 0\n", signy, signz, mpfr_print_rnd_mode ((mp_rnd_t) rnd), inex); err = 1; } - if (inex <= 0 && (rnd == GMP_RNDU || - (rnd == GMP_RNDZ && signz > 0) || - (rnd == GMP_RNDN && signy < 0))) + if (inex <= 0 && (rnd == MPFR_RNDU || + (rnd == MPFR_RNDZ && signz > 0) || + (rnd == MPFR_RNDN && signy < 0))) { printf (ERRTU1 "ternary value = %d instead of > 0\n", signy, signz, mpfr_print_rnd_mode ((mp_rnd_t) rnd), @@ -289,21 +289,21 @@ test_underflow2 (void) mpfr_inits2 (32, x, y, z, r, (mpfr_ptr) 0); - mpfr_set_si_2exp (z, -1, mpfr_get_emin (), GMP_RNDN); /* z = -2^emin */ - mpfr_set_si_2exp (x, 1, mpfr_get_emin (), GMP_RNDN); /* x = 2^emin */ + mpfr_set_si_2exp (z, -1, mpfr_get_emin (), MPFR_RNDN); /* z = -2^emin */ + mpfr_set_si_2exp (x, 1, mpfr_get_emin (), MPFR_RNDN); /* x = 2^emin */ for (b = 0; b <= 1; b++) { for (i = 15; i <= 17; i++) { - mpfr_set_si_2exp (y, i, -4 - MPFR_PREC (z), GMP_RNDN); + mpfr_set_si_2exp (y, i, -4 - MPFR_PREC (z), MPFR_RNDN); /* z = -1.000...00b * xy = 01111 * or 10000 * or 10001 */ mpfr_clear_flags (); - inex = mpfr_fms (r, x, y, z, GMP_RNDN); + inex = mpfr_fms (r, x, y, z, MPFR_RNDN); #define ERRTU2 "Error in test_underflow2 (b = %d, i = %d)\n " if (__gmpfr_flags != MPFR_FLAGS_INEXACT) { @@ -318,7 +318,7 @@ test_underflow2 (void) b, i, inex, same ? '<' : '>'); err = 1; } - mpfr_neg (y, z, GMP_RNDN); + mpfr_neg (y, z, MPFR_RNDN); if (!same) mpfr_nextabove (y); if (! mpfr_equal_p (r, y)) @@ -356,10 +356,10 @@ main (int argc, char *argv[]) mpfr_set_prec (y, 2); mpfr_set_prec (z, 2); mpfr_set_prec (s, 2); - mpfr_set_str (x, "-0.75", 10, GMP_RNDN); - mpfr_set_str (y, "0.5", 10, GMP_RNDN); - mpfr_set_str (z, "-0.375", 10, GMP_RNDN); - mpfr_fms (s, x, y, z, GMP_RNDU); /* result is 0 */ + mpfr_set_str (x, "-0.75", 10, MPFR_RNDN); + mpfr_set_str (y, "0.5", 10, MPFR_RNDN); + mpfr_set_str (z, "-0.375", 10, MPFR_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDU); /* result is 0 */ if (mpfr_cmp_ui(s, 0)) { printf("Error: -0.75 * 0.5 - -0.375 should be equal to 0 for prec=2\n"); @@ -371,9 +371,9 @@ main (int argc, char *argv[]) mpfr_set_prec (z, 27); mpfr_set_prec (s, 27); mpfr_set_str_binary (x, "1.11111111111111111111111111e-1"); - mpfr_set (y, x, GMP_RNDN); + mpfr_set (y, x, MPFR_RNDN); mpfr_set_str_binary (z, "1.00011110100011001011001001e-1"); - if (mpfr_fms (s, x, y, z, GMP_RNDN) >= 0) + if (mpfr_fms (s, x, y, z, MPFR_RNDN) >= 0) { printf ("Wrong inexact flag for x=y=1-2^(-27)\n"); exit (1); @@ -382,7 +382,7 @@ main (int argc, char *argv[]) mpfr_set_nan (x); mpfr_urandomb (y, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in x=NAN does not return NAN"); @@ -392,7 +392,7 @@ main (int argc, char *argv[]) mpfr_set_nan (y); mpfr_urandomb (x, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p(s)) { printf ("evaluation of function in y=NAN does not return NAN"); @@ -402,7 +402,7 @@ main (int argc, char *argv[]) mpfr_set_nan (z); mpfr_urandomb (y, RANDS); mpfr_urandomb (x, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in z=NAN does not return NAN"); @@ -412,7 +412,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, 1); mpfr_set_inf (y, 1); mpfr_set_inf (z, -1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) < 0) { printf ("Error for (+inf) * (+inf) - (-inf)\n"); @@ -422,7 +422,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, -1); mpfr_set_inf (y, -1); mpfr_set_inf (z, -1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) < 0) { printf ("Error for (-inf) * (-inf) - (-inf)\n"); @@ -432,7 +432,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, 1); mpfr_set_inf (y, -1); mpfr_set_inf (z, 1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) > 0) { printf ("Error for (+inf) * (-inf) - (+inf)\n"); @@ -442,7 +442,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, -1); mpfr_set_inf (y, 1); mpfr_set_inf (z, 1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) > 0) { printf ("Error for (-inf) * (+inf) - (+inf)\n"); @@ -450,9 +450,9 @@ main (int argc, char *argv[]) } mpfr_set_inf (x, 1); - mpfr_set_ui (y, 0, GMP_RNDN); + mpfr_set_ui (y, 0, MPFR_RNDN); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in x=INF y=0 does not return NAN"); @@ -460,9 +460,9 @@ main (int argc, char *argv[]) } mpfr_set_inf (y, 1); - mpfr_set_ui (x, 0, GMP_RNDN); + mpfr_set_ui (x, 0, MPFR_RNDN); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in x=0 y=INF does not return NAN"); @@ -472,7 +472,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, 1); mpfr_urandomb (y, RANDS); /* always positive */ mpfr_set_inf (z, 1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in x=INF y>0 z=INF does not return NAN"); @@ -482,7 +482,7 @@ main (int argc, char *argv[]) mpfr_set_inf (y, 1); mpfr_urandomb (x, RANDS); mpfr_set_inf (z, 1); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_nan_p (s)) { printf ("evaluation of function in x>0 y=INF z=INF does not return NAN"); @@ -492,7 +492,7 @@ main (int argc, char *argv[]) mpfr_set_inf (x, 1); mpfr_urandomb (y, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) < 0) { printf ("evaluation of function in x=INF does not return INF"); @@ -502,7 +502,7 @@ main (int argc, char *argv[]) mpfr_set_inf (y, 1); mpfr_urandomb (x, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) < 0) { printf ("evaluation of function in y=INF does not return INF"); @@ -512,29 +512,29 @@ main (int argc, char *argv[]) mpfr_set_inf (z, -1); mpfr_urandomb (x, RANDS); mpfr_urandomb (y, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); if (!mpfr_inf_p (s) || mpfr_sgn (s) < 0) { printf ("evaluation of function in z=-INF does not return INF"); exit (1); } - mpfr_set_ui (x, 0, GMP_RNDN); + mpfr_set_ui (x, 0, MPFR_RNDN); mpfr_urandomb (y, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); - mpfr_neg (z, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); + mpfr_neg (z, z, MPFR_RNDN); if (mpfr_cmp (s, z)) { printf ("evaluation of function in x=0 does not return -z\n"); exit (1); } - mpfr_set_ui (y, 0, GMP_RNDN); + mpfr_set_ui (y, 0, MPFR_RNDN); mpfr_urandomb (x, RANDS); mpfr_urandomb (z, RANDS); - mpfr_fms (s, x, y, z, GMP_RNDN); - mpfr_neg (z, z, GMP_RNDN); + mpfr_fms (s, x, y, z, MPFR_RNDN); + mpfr_neg (z, z, MPFR_RNDN); if (mpfr_cmp (s, z)) { printf ("evaluation of function in y=0 does not return -z\n"); @@ -570,11 +570,11 @@ main (int argc, char *argv[]) mpfr_urandomb (z, RANDS); if (randlimb () % 2) - mpfr_neg (x, x, GMP_RNDN); + mpfr_neg (x, x, MPFR_RNDN); if (randlimb () % 2) - mpfr_neg (y, y, GMP_RNDN); + mpfr_neg (y, y, MPFR_RNDN); if (randlimb () % 2) - mpfr_neg (z, z, GMP_RNDN); + mpfr_neg (z, z, MPFR_RNDN); rnd = RND_RAND (); mpfr_set_prec (slong, 2 * prec); @@ -588,18 +588,18 @@ main (int argc, char *argv[]) if (mpfr_cmp (s, t)) { printf ("results differ for x="); - mpfr_out_str (stdout, 2, prec, x, GMP_RNDN); + mpfr_out_str (stdout, 2, prec, x, MPFR_RNDN); printf (" y="); - mpfr_out_str (stdout, 2, prec, y, GMP_RNDN); + mpfr_out_str (stdout, 2, prec, y, MPFR_RNDN); printf (" z="); - mpfr_out_str (stdout, 2, prec, z, GMP_RNDN); + mpfr_out_str (stdout, 2, prec, z, MPFR_RNDN); printf (" prec=%u rnd_mode=%s\n", (unsigned int) prec, mpfr_print_rnd_mode (rnd)); printf ("got "); - mpfr_out_str (stdout, 2, prec, s, GMP_RNDN); + mpfr_out_str (stdout, 2, prec, s, MPFR_RNDN); puts (""); printf ("expected "); - mpfr_out_str (stdout, 2, prec, t, GMP_RNDN); + mpfr_out_str (stdout, 2, prec, t, MPFR_RNDN); puts (""); printf ("approx "); mpfr_print_binary (slong); @@ -612,10 +612,10 @@ main (int argc, char *argv[]) { printf ("Wrong inexact flag for rnd=%s: expected %d, got %d\n", mpfr_print_rnd_mode (rnd), compare, inexact); - printf (" x="); mpfr_out_str (stdout, 2, 0, x, GMP_RNDN); - printf (" y="); mpfr_out_str (stdout, 2, 0, y, GMP_RNDN); - printf (" z="); mpfr_out_str (stdout, 2, 0, z, GMP_RNDN); - printf (" s="); mpfr_out_str (stdout, 2, 0, s, GMP_RNDN); + printf (" x="); mpfr_out_str (stdout, 2, 0, x, MPFR_RNDN); + printf (" y="); mpfr_out_str (stdout, 2, 0, y, MPFR_RNDN); + printf (" z="); mpfr_out_str (stdout, 2, 0, z, MPFR_RNDN); + printf (" s="); mpfr_out_str (stdout, 2, 0, s, MPFR_RNDN); printf ("\n"); exit (1); } |