summaryrefslogtreecommitdiff
path: root/tests/tadd.c
diff options
context:
space:
mode:
authorzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2009-03-03 16:51:07 +0000
committerzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2009-03-03 16:51:07 +0000
commitf5f7883fc7cf47d38ed493a882fd7109bdd589bb (patch)
tree1058cf0d07593fa2cec9f32b4988ac40611b4bff /tests/tadd.c
parentcbd40e1cb659b579e608afbe7a05da5b15a860be (diff)
downloadmpfr-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/tadd.c')
-rw-r--r--tests/tadd.c352
1 files changed, 176 insertions, 176 deletions
diff --git a/tests/tadd.c b/tests/tadd.c
index fee6eebe4..f1efb7a9a 100644
--- a/tests/tadd.c
+++ b/tests/tadd.c
@@ -38,7 +38,7 @@ test_add (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mp_rnd_t rnd_mode)
{
int res;
#ifdef CHECK_EXTERNAL
- int ok = rnd_mode == GMP_RNDN && mpfr_number_p (b) && mpfr_number_p (c);
+ int ok = rnd_mode == MPFR_RNDN && mpfr_number_p (b) && mpfr_number_p (c);
if (ok)
{
mpfr_print_raw (b);
@@ -94,7 +94,7 @@ pcheck (const char *xs, const char *ys, const char *zs, mp_rnd_t rnd_mode,
if (mpfr_cmp_str1 (zz, zs) )
{
printf ("expected sum is %s, got ", zs);
- mpfr_out_str(stdout, 10, 0, zz, GMP_RNDN);
+ mpfr_out_str(stdout, 10, 0, zz, MPFR_RNDN);
printf ("mpfr_add failed for x=%s y=%s with rnd_mode=%s\n",
xs, ys, mpfr_print_rnd_mode (rnd_mode));
exit (1);
@@ -114,12 +114,12 @@ check2b (const char *xs, int px,
mpfr_set_str_binary (xx, xs);
mpfr_set_str_binary (yy, ys);
test_add (zz, xx, yy, rnd_mode);
- if (mpfr_cmp_str (zz, rs, 2, GMP_RNDN))
+ if (mpfr_cmp_str (zz, rs, 2, MPFR_RNDN))
{
printf ("(2) x=%s,%d y=%s,%d pz=%d,rnd=%s\n",
xs, px, ys, py, pz, mpfr_print_rnd_mode (rnd_mode));
printf ("got "); mpfr_print_binary(zz); puts ("");
- mpfr_set_str(zz, rs, 2, GMP_RNDN);
+ mpfr_set_str(zz, rs, 2, MPFR_RNDN);
printf ("instead of "); mpfr_print_binary(zz); puts ("");
exit (1);
}
@@ -140,14 +140,14 @@ check64 (void)
mpfr_set_prec (t, 58);
mpfr_set_str_binary (t, "0.11100010011111001001100110010111110110011000000100101E-1");
mpfr_set_prec (u, 29);
- test_add (u, x, t, GMP_RNDD);
+ test_add (u, x, t, MPFR_RNDD);
mpfr_set_str_binary (t, "1.0101011100001000011100111110e-1");
if (mpfr_cmp (u, t))
{
printf ("mpfr_add(u, x, t) failed for prec(x)=29, prec(t)=58\n");
- printf ("expected "); mpfr_out_str (stdout, 2, 29, t, GMP_RNDN);
+ printf ("expected "); mpfr_out_str (stdout, 2, 29, t, MPFR_RNDN);
puts ("");
- printf ("got "); mpfr_out_str (stdout, 2, 29, u, GMP_RNDN);
+ printf ("got "); mpfr_out_str (stdout, 2, 29, u, MPFR_RNDN);
puts ("");
exit(1);
}
@@ -157,7 +157,7 @@ check64 (void)
mpfr_set_prec (t, 2);
mpfr_set_str_binary (t, "-1.1e-2");
mpfr_set_prec (u, 2);
- test_add (u, x, t, GMP_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
if (MPFR_MANT(u)[0] << 2)
{
printf ("result not normalized for prec=2\n");
@@ -169,7 +169,7 @@ check64 (void)
{
printf ("mpfr_add(u, x, t) failed for prec(x)=4, prec(t)=2\n");
printf ("expected -1.0e-1\n");
- printf ("got "); mpfr_out_str (stdout, 2, 4, u, GMP_RNDN);
+ printf ("got "); mpfr_out_str (stdout, 2, 4, u, MPFR_RNDN);
puts ("");
exit (1);
}
@@ -179,12 +179,12 @@ check64 (void)
mpfr_set_prec (t, 4);
mpfr_set_str_binary (t, "-1.110e-5"); /* -7/128 */
mpfr_set_prec (u, 4);
- test_add (u, x, t, GMP_RNDN); /* should give -5/8 */
+ test_add (u, x, t, MPFR_RNDN); /* should give -5/8 */
mpfr_set_str_binary (t, "-1.010e-1");
if (mpfr_cmp (u, t)) {
printf ("mpfr_add(u, x, t) failed for prec(x)=8, prec(t)=4\n");
printf ("expected -1.010e-1\n");
- printf ("got "); mpfr_out_str (stdout, 2, 4, u, GMP_RNDN);
+ printf ("got "); mpfr_out_str (stdout, 2, 4, u, MPFR_RNDN);
puts ("");
exit (1);
}
@@ -192,7 +192,7 @@ check64 (void)
mpfr_set_prec (x, 112); mpfr_set_prec (t, 98); mpfr_set_prec (u, 54);
mpfr_set_str_binary (x, "-0.11111100100000000011000011100000101101010001000111E-401");
mpfr_set_str_binary (t, "0.10110000100100000101101100011111111011101000111000101E-464");
- test_add (u, x, t, GMP_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
if (mpfr_cmp (u, x))
{
printf ("mpfr_add(u, x, t) failed for prec(x)=112, prec(t)=98\n");
@@ -200,9 +200,9 @@ check64 (void)
}
mpfr_set_prec (x, 92); mpfr_set_prec (t, 86); mpfr_set_prec (u, 53);
- mpfr_set_str (x, "-5.03525136761487735093e-74", 10, GMP_RNDN);
- mpfr_set_str (t, "8.51539046314262304109e-91", 10, GMP_RNDN);
- test_add (u, x, t, GMP_RNDN);
+ mpfr_set_str (x, "-5.03525136761487735093e-74", 10, MPFR_RNDN);
+ mpfr_set_str (t, "8.51539046314262304109e-91", 10, MPFR_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
if (mpfr_cmp_str1 (u, "-5.0352513676148773509283672e-74") )
{
printf ("mpfr_add(u, x, t) failed for prec(x)=92, prec(t)=86\n");
@@ -212,7 +212,7 @@ check64 (void)
mpfr_set_prec(x, 53); mpfr_set_prec(t, 76); mpfr_set_prec(u, 76);
mpfr_set_str_binary(x, "-0.10010010001001011011110000000000001010011011011110001E-32");
mpfr_set_str_binary(t, "-0.1011000101110010000101111111011111010001110011110111100110101011110010011111");
- mpfr_sub(u, x, t, GMP_RNDU);
+ mpfr_sub(u, x, t, MPFR_RNDU);
mpfr_set_str_binary(t, "0.1011000101110010000101111111011100111111101010011011110110101011101000000100");
if (mpfr_cmp(u,t))
{
@@ -223,7 +223,7 @@ check64 (void)
mpfr_set_prec(x, 53); mpfr_set_prec(t, 108); mpfr_set_prec(u, 108);
mpfr_set_str_binary(x, "-0.10010010001001011011110000000000001010011011011110001E-32");
mpfr_set_str_binary(t, "-0.101100010111001000010111111101111101000111001111011110011010101111001001111000111011001110011000000000111111");
- mpfr_sub(u, x, t, GMP_RNDU);
+ mpfr_sub(u, x, t, MPFR_RNDU);
mpfr_set_str_binary(t, "0.101100010111001000010111111101110011111110101001101111011010101110100000001011000010101110011000000000111111");
if (mpfr_cmp(u,t))
{
@@ -233,8 +233,8 @@ check64 (void)
}
mpfr_set_prec(x, 97); mpfr_set_prec(t, 97); mpfr_set_prec(u, 97);
mpfr_set_str_binary(x, "0.1111101100001000000001011000110111101000001011111000100001000101010100011111110010000000000000000E-39");
- mpfr_set_ui(t, 1, GMP_RNDN);
- test_add (u, x, t, GMP_RNDN);
+ mpfr_set_ui(t, 1, MPFR_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
mpfr_set_str_binary(x, "0.1000000000000000000000000000000000000000111110110000100000000101100011011110100000101111100010001E1");
if (mpfr_cmp(u,x))
{
@@ -243,17 +243,17 @@ check64 (void)
}
mpfr_set_prec(x, 128); mpfr_set_prec(t, 128); mpfr_set_prec(u, 128);
mpfr_set_str_binary(x, "0.10101011111001001010111011001000101100111101000000111111111011010100001100011101010001010111111101111010100110111111100101100010E-4");
- mpfr_set(t, x, GMP_RNDN);
- mpfr_sub(u, x, t, GMP_RNDN);
+ mpfr_set(t, x, MPFR_RNDN);
+ mpfr_sub(u, x, t, MPFR_RNDN);
mpfr_set_prec(x, 96); mpfr_set_prec(t, 96); mpfr_set_prec(u, 96);
mpfr_set_str_binary(x, "0.111000000001110100111100110101101001001010010011010011100111100011010100011001010011011011000010E-4");
- mpfr_set(t, x, GMP_RNDN);
- mpfr_sub(u, x, t, GMP_RNDN);
+ mpfr_set(t, x, MPFR_RNDN);
+ mpfr_sub(u, x, t, MPFR_RNDN);
mpfr_set_prec(x, 85); mpfr_set_prec(t, 85); mpfr_set_prec(u, 85);
mpfr_set_str_binary(x, "0.1111101110100110110110100010101011101001100010100011110110110010010011101100101111100E-4");
mpfr_set_str_binary(t, "0.1111101110100110110110100010101001001000011000111000011101100101110100001110101010110E-4");
- mpfr_sub(u, x, t, GMP_RNDU);
- mpfr_sub(x, x, t, GMP_RNDU);
+ mpfr_sub(u, x, t, MPFR_RNDU);
+ mpfr_sub(x, x, t, MPFR_RNDU);
if (mpfr_cmp(x, u) != 0)
{
printf ("Error in mpfr_sub: u=x-t and x=x-t give different results\n");
@@ -268,7 +268,7 @@ check64 (void)
mpfr_set_prec(x, 65); mpfr_set_prec(t, 65); mpfr_set_prec(u, 65);
mpfr_set_str_binary(x, "0.10011010101000110101010000000011001001001110001011101011111011101E623");
mpfr_set_str_binary(t, "0.10011010101000110101010000000011001001001110001011101011111011100E623");
- mpfr_sub(u, x, t, GMP_RNDU);
+ mpfr_sub(u, x, t, MPFR_RNDU);
if (mpfr_cmp_ui_2exp(u, 1, 558))
{ /* 2^558 */
printf ("Error (1) in mpfr_sub\n");
@@ -278,10 +278,10 @@ check64 (void)
mpfr_set_prec(x, 64); mpfr_set_prec(t, 64); mpfr_set_prec(u, 64);
mpfr_set_str_binary(x, "0.1000011110101111011110111111000011101011101111101101101100000100E-220");
mpfr_set_str_binary(t, "0.1000011110101111011110111111000011101011101111101101010011111101E-220");
- test_add (u, x, t, GMP_RNDU);
+ test_add (u, x, t, MPFR_RNDU);
if ((MPFR_MANT(u)[0] & 1) != 1)
{
- printf ("error in mpfr_add with rnd_mode=GMP_RNDU\n");
+ printf ("error in mpfr_add with rnd_mode=MPFR_RNDU\n");
printf ("b= "); mpfr_print_binary(x); puts ("");
printf ("c= "); mpfr_print_binary(t); puts ("");
printf ("b+c="); mpfr_print_binary(u); puts ("");
@@ -292,31 +292,31 @@ check64 (void)
mpfr_set_prec(x, 56); mpfr_set_prec(t, 83); mpfr_set_prec(u, 10);
mpfr_set_str_binary(x, "0.10001001011011001111101100110100000101111010010111010111E-7");
mpfr_set_str_binary(t, "0.10001001011011001111101100110100000101111010010111010111000000000111110110110000100E-7");
- mpfr_sub(u, x, t, GMP_RNDU);
+ mpfr_sub(u, x, t, MPFR_RNDU);
/* array bound write found by Norbert Mueller, 26 Sep 2000 */
mpfr_set_prec(x, 109); mpfr_set_prec(t, 153); mpfr_set_prec(u, 95);
mpfr_set_str_binary(x,"0.1001010000101011101100111000110001111111111111111111111111111111111111111111111111111111111111100000000000000E33");
mpfr_set_str_binary(t,"-0.100101000010101110110011100011000000000000000000000000000000000000000000000000000000000000000000000000000000000000000011100101101000000100100001100110111E33");
- test_add (u, x, t, GMP_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
/* array bound writes found by Norbert Mueller, 27 Sep 2000 */
mpfr_set_prec(x, 106); mpfr_set_prec(t, 53); mpfr_set_prec(u, 23);
mpfr_set_str_binary(x, "-0.1000011110101111111001010001000100001011000000000000000000000000000000000000000000000000000000000000000000E-59");
mpfr_set_str_binary(t, "-0.10000111101011111110010100010001101100011100110100000E-59");
- mpfr_sub(u, x, t, GMP_RNDN);
+ mpfr_sub(u, x, t, MPFR_RNDN);
mpfr_set_prec(x, 177); mpfr_set_prec(t, 217); mpfr_set_prec(u, 160);
mpfr_set_str_binary(x, "-0.111010001011010000111001001010010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000E35");
mpfr_set_str_binary(t, "0.1110100010110100001110010010100100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000111011010011100001111001E35");
- test_add (u, x, t, GMP_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
mpfr_set_prec(x, 214); mpfr_set_prec(t, 278); mpfr_set_prec(u, 207);
mpfr_set_str_binary(x, "0.1000100110100110101101101101000000010000100111000001001110001000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000E66");
mpfr_set_str_binary(t, "-0.10001001101001101011011011010000000100001001110000010011100010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001111011111001001100011E66");
- test_add (u, x, t, GMP_RNDN);
+ test_add (u, x, t, MPFR_RNDN);
mpfr_set_prec(x, 32); mpfr_set_prec(t, 247); mpfr_set_prec(u, 223);
mpfr_set_str_binary(x, "0.10000000000000000000000000000000E1");
mpfr_set_str_binary(t, "0.1111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111100000110001110100000100011110000101110110011101110100110110111111011010111100100000000000000000000000000E0");
- mpfr_sub(u, x, t, GMP_RNDN);
+ mpfr_sub(u, x, t, MPFR_RNDN);
if ((MPFR_MANT(u)[(MPFR_PREC(u)-1)/mp_bits_per_limb] &
((mp_limb_t)1<<(mp_bits_per_limb-1)))==0)
{
@@ -330,7 +330,7 @@ check64 (void)
mpfr_set_prec (u, 65);
mpfr_set_str_binary (x, "0.11100100101101001100111011111111110001101001000011101001001010010E-35");
mpfr_set_str_binary (t, "0.10000000000000000000000000000000000001110010010110100110011110000E1");
- mpfr_sub (u, t, x, GMP_RNDU);
+ mpfr_sub (u, t, x, MPFR_RNDU);
if ((MPFR_MANT(u)[(MPFR_PREC(u)-1)/mp_bits_per_limb] &
((mp_limb_t)1<<(mp_bits_per_limb-1)))==0)
{
@@ -344,7 +344,7 @@ check64 (void)
mpfr_set_prec (u, 107);
mpfr_set_str_binary (x, "0.10111001001111010010001000000010111111011011011101000001001000101000000000000000000000000000000000000000000E315");
mpfr_set_str_binary (t, "0.10000000000000000000000000000000000101110100100101110110000001100101011111001000011101111100100100111011000E350");
- mpfr_sub (u, x, t, GMP_RNDU);
+ mpfr_sub (u, x, t, MPFR_RNDU);
if ((MPFR_MANT(u)[(MPFR_PREC(u)-1)/mp_bits_per_limb] &
((mp_limb_t)1<<(mp_bits_per_limb-1)))==0)
{
@@ -353,14 +353,14 @@ check64 (void)
}
/* checks that NaN flag is correctly reset */
- mpfr_set_ui (t, 1, GMP_RNDN);
- mpfr_set_ui (u, 1, GMP_RNDN);
+ mpfr_set_ui (t, 1, MPFR_RNDN);
+ mpfr_set_ui (u, 1, MPFR_RNDN);
mpfr_set_nan (x);
- test_add (x, t, u, GMP_RNDN);
+ test_add (x, t, u, MPFR_RNDN);
if (mpfr_cmp_ui (x, 2))
{
printf ("Error in mpfr_add: 1+1 gives ");
- mpfr_out_str(stdout, 10, 0, x, GMP_RNDN);
+ mpfr_out_str(stdout, 10, 0, x, MPFR_RNDN);
exit (1);
}
@@ -387,20 +387,20 @@ check_case_1b (void)
dif = prec_b - prec_a;
mpfr_set_prec (b, prec_b);
/* b = 1 - 2^(-prec_a) + 2^(-prec_b) */
- mpfr_set_ui (b, 1, GMP_RNDN);
- mpfr_div_2exp (b, b, dif, GMP_RNDN);
- mpfr_sub_ui (b, b, 1, GMP_RNDN);
- mpfr_div_2exp (b, b, prec_a, GMP_RNDN);
- mpfr_add_ui (b, b, 1, GMP_RNDN);
+ mpfr_set_ui (b, 1, MPFR_RNDN);
+ mpfr_div_2exp (b, b, dif, MPFR_RNDN);
+ mpfr_sub_ui (b, b, 1, MPFR_RNDN);
+ mpfr_div_2exp (b, b, prec_a, MPFR_RNDN);
+ mpfr_add_ui (b, b, 1, MPFR_RNDN);
for (prec_c = dif; prec_c <= 64; prec_c++)
{
/* c = 2^(-prec_a) - 2^(-prec_b) */
mpfr_set_prec (c, prec_c);
- mpfr_set_si (c, -1, GMP_RNDN);
- mpfr_div_2exp (c, c, dif, GMP_RNDN);
- mpfr_add_ui (c, c, 1, GMP_RNDN);
- mpfr_div_2exp (c, c, prec_a, GMP_RNDN);
- test_add (a, b, c, GMP_RNDN);
+ mpfr_set_si (c, -1, MPFR_RNDN);
+ mpfr_div_2exp (c, c, dif, MPFR_RNDN);
+ mpfr_add_ui (c, c, 1, MPFR_RNDN);
+ mpfr_div_2exp (c, c, prec_a, MPFR_RNDN);
+ test_add (a, b, c, MPFR_RNDN);
if (mpfr_cmp_ui (a, 1) != 0)
{
printf ("case (1b) failed for prec_a=%u, prec_b=%u,"
@@ -433,15 +433,15 @@ check_case_2 (void)
mpfr_set_str_binary(a, "1E110"); /* a = 2^110 */
mpfr_set_str_binary(b, "1E900"); /* b = 2^900 */
mpfr_set_str_binary(c, "1E500"); /* c = 2^500 */
- test_add (c, c, a, GMP_RNDZ); /* c = 2^500 + 2^110 */
- mpfr_sub (d, b, c, GMP_RNDZ); /* d = 2^900 - 2^500 - 2^110 */
- test_add (b, b, c, GMP_RNDZ); /* b = 2^900 + 2^500 + 2^110 */
- test_add (a, b, d, GMP_RNDZ); /* a = 2^901 */
+ test_add (c, c, a, MPFR_RNDZ); /* c = 2^500 + 2^110 */
+ mpfr_sub (d, b, c, MPFR_RNDZ); /* d = 2^900 - 2^500 - 2^110 */
+ test_add (b, b, c, MPFR_RNDZ); /* b = 2^900 + 2^500 + 2^110 */
+ test_add (a, b, d, MPFR_RNDZ); /* a = 2^901 */
if (mpfr_cmp_ui_2exp (a, 1, 901))
{
printf ("b + d fails for b=2^900+2^500+2^110, d=2^900-2^500-2^110\n");
printf ("expected 1.0e901, got ");
- mpfr_out_str (stdout, 2, 0, a, GMP_RNDN);
+ mpfr_out_str (stdout, 2, 0, a, MPFR_RNDN);
printf ("\n");
exit (1);
}
@@ -458,8 +458,8 @@ check_same (void)
{
mpfr_t x;
- mpfr_init(x); mpfr_set_ui(x, 1, GMP_RNDZ);
- test_add (x, x, x, GMP_RNDZ);
+ mpfr_init(x); mpfr_set_ui(x, 1, MPFR_RNDZ);
+ test_add (x, x, x, MPFR_RNDZ);
if (mpfr_cmp_ui (x, 2))
{
printf ("Error when all 3 operands are equal\n");
@@ -490,7 +490,7 @@ check_inexact (void)
mpfr_set_prec (u, 33);
mpfr_set_str_binary (u, "0.101110100101101100000000111100000E-1");
mpfr_set_prec (y, 31);
- if ((inexact = test_add (y, x, u, GMP_RNDN)))
+ if ((inexact = test_add (y, x, u, MPFR_RNDN)))
{
printf ("Wrong inexact flag (2): expected 0, got %d\n", inexact);
exit (1);
@@ -501,7 +501,7 @@ check_inexact (void)
mpfr_set_prec (u, 33);
mpfr_set_str_binary (u, "0.101110100101101100000000111100000E-1");
mpfr_set_prec (y, 28);
- if ((inexact = test_add (y, x, u, GMP_RNDN)))
+ if ((inexact = test_add (y, x, u, MPFR_RNDN)))
{
printf ("Wrong inexact flag (1): expected 0, got %d\n", inexact);
exit (1);
@@ -581,34 +581,34 @@ check_nans (void)
/* +inf + -inf == nan */
mpfr_set_inf (x, 1);
mpfr_set_inf (y, -1);
- test_add (s, x, y, GMP_RNDN);
+ test_add (s, x, y, MPFR_RNDN);
MPFR_ASSERTN (mpfr_nan_p (s));
/* +inf + 1 == +inf */
mpfr_set_inf (x, 1);
- mpfr_set_ui (y, 1L, GMP_RNDN);
- test_add (s, x, y, GMP_RNDN);
+ mpfr_set_ui (y, 1L, MPFR_RNDN);
+ test_add (s, x, y, MPFR_RNDN);
MPFR_ASSERTN (mpfr_inf_p (s));
MPFR_ASSERTN (mpfr_sgn (s) > 0);
/* -inf + 1 == -inf */
mpfr_set_inf (x, -1);
- mpfr_set_ui (y, 1L, GMP_RNDN);
- test_add (s, x, y, GMP_RNDN);
+ mpfr_set_ui (y, 1L, MPFR_RNDN);
+ test_add (s, x, y, MPFR_RNDN);
MPFR_ASSERTN (mpfr_inf_p (s));
MPFR_ASSERTN (mpfr_sgn (s) < 0);
/* 1 + +inf == +inf */
- mpfr_set_ui (x, 1L, GMP_RNDN);
+ mpfr_set_ui (x, 1L, MPFR_RNDN);
mpfr_set_inf (y, 1);
- test_add (s, x, y, GMP_RNDN);
+ test_add (s, x, y, MPFR_RNDN);
MPFR_ASSERTN (mpfr_inf_p (s));
MPFR_ASSERTN (mpfr_sgn (s) > 0);
/* 1 + -inf == -inf */
- mpfr_set_ui (x, 1L, GMP_RNDN);
+ mpfr_set_ui (x, 1L, MPFR_RNDN);
mpfr_set_inf (y, -1);
- test_add (s, x, y, GMP_RNDN);
+ test_add (s, x, y, MPFR_RNDN);
MPFR_ASSERTN (mpfr_inf_p (s));
MPFR_ASSERTN (mpfr_sgn (s) < 0);
@@ -624,11 +624,11 @@ check_alloc (void)
mpfr_init2 (a, 10000);
mpfr_set_prec (a, 53);
- mpfr_set_ui (a, 15236, GMP_RNDN);
- test_add (a, a, a, GMP_RNDN);
- mpfr_mul (a, a, a, GMP_RNDN);
- mpfr_div (a, a, a, GMP_RNDN);
- mpfr_sub (a, a, a, GMP_RNDN);
+ mpfr_set_ui (a, 15236, MPFR_RNDN);
+ test_add (a, a, a, MPFR_RNDN);
+ mpfr_mul (a, a, a, MPFR_RNDN);
+ mpfr_div (a, a, a, MPFR_RNDN);
+ mpfr_sub (a, a, a, MPFR_RNDN);
mpfr_clear (a);
}
@@ -643,9 +643,9 @@ check_overflow (void)
mpfr_init2 (b, 256);
mpfr_init2 (c, 256);
- mpfr_set_ui (b, 1, GMP_RNDN);
+ mpfr_set_ui (b, 1, MPFR_RNDN);
mpfr_setmax (b, mpfr_get_emax ());
- mpfr_set_ui (c, 1, GMP_RNDN);
+ mpfr_set_ui (c, 1, MPFR_RNDN);
mpfr_set_exp (c, mpfr_get_emax () - 192);
RND_LOOP(r)
for (prec_a = 128; prec_a < 512; prec_a += 64)
@@ -663,7 +663,7 @@ check_overflow (void)
mpfr_set_exp (c, mpfr_get_emax () - 512);
mpfr_set_prec (a, 256);
mpfr_clear_overflow ();
- test_add (a, b, c, GMP_RNDU);
+ test_add (a, b, c, MPFR_RNDU);
if (!mpfr_overflow_p ())
{
printf ("No overflow in check_overflow\n");
@@ -682,7 +682,7 @@ check_1111 (void)
long n;
mpfr_init2 (one, MPFR_PREC_MIN);
- mpfr_set_ui (one, 1, GMP_RNDN);
+ mpfr_set_ui (one, 1, MPFR_RNDN);
for (n = 0; n < NUM; n++)
{
mp_prec_t prec_a, prec_b, prec_c;
@@ -703,25 +703,25 @@ check_1111 (void)
if (sb != 0)
{
tb = 1 + (randlimb () % (prec_b - (sb != 2)));
- mpfr_div_2ui (b, one, tb, GMP_RNDN);
+ mpfr_div_2ui (b, one, tb, MPFR_RNDN);
if (sb == 2)
- mpfr_neg (b, b, GMP_RNDN);
- test_add (b, b, one, GMP_RNDN);
+ mpfr_neg (b, b, MPFR_RNDN);
+ test_add (b, b, one, MPFR_RNDN);
}
else
- mpfr_set (b, one, GMP_RNDN);
+ mpfr_set (b, one, MPFR_RNDN);
tc = 1 + (randlimb () % (prec_c - 1));
- mpfr_div_2ui (c, one, tc, GMP_RNDN);
+ mpfr_div_2ui (c, one, tc, MPFR_RNDN);
sc = randlimb () % 2;
if (sc)
- mpfr_neg (c, c, GMP_RNDN);
- test_add (c, c, one, GMP_RNDN);
+ mpfr_neg (c, c, MPFR_RNDN);
+ test_add (c, c, one, MPFR_RNDN);
diff = (randlimb () % (2*m)) - m;
- mpfr_mul_2si (c, c, diff, GMP_RNDN);
+ mpfr_mul_2si (c, c, diff, MPFR_RNDN);
rnd_mode = RND_RAND ();
inex_a = test_add (a, b, c, rnd_mode);
mpfr_init2 (s, MPFR_PREC_MIN + 2*m);
- inex_s = test_add (s, b, c, GMP_RNDN); /* exact */
+ inex_s = test_add (s, b, c, MPFR_RNDN); /* exact */
if (inex_s)
{
printf ("check_1111: result should have been exact.\n");
@@ -778,21 +778,21 @@ check_1minuseps (void)
mpfr_init2 (a, prec_a[ia]);
mpfr_init2 (b, prec_b);
- mpfr_set_ui (c, 1, GMP_RNDN);
- mpfr_div_ui (b, c, prec_a[ia], GMP_RNDN);
- mpfr_sub (b, c, b, GMP_RNDN); /* b = 1 - 2^(-prec_a) */
+ mpfr_set_ui (c, 1, MPFR_RNDN);
+ mpfr_div_ui (b, c, prec_a[ia], MPFR_RNDN);
+ mpfr_sub (b, c, b, MPFR_RNDN); /* b = 1 - 2^(-prec_a) */
for (ic = 0; ic < numberof(supp_b); ic++)
- for (rnd_mode = 0; rnd_mode < GMP_RND_MAX; rnd_mode++)
+ for (rnd_mode = 0; rnd_mode < MPFR_RND_MAX; rnd_mode++)
{
mpfr_t s;
int inex_a, inex_s;
- mpfr_set_ui (c, 1, GMP_RNDN);
- mpfr_div_ui (c, c, prec_a[ia] + supp_b[ic], GMP_RNDN);
+ mpfr_set_ui (c, 1, MPFR_RNDN);
+ mpfr_div_ui (c, c, prec_a[ia] + supp_b[ic], MPFR_RNDN);
inex_a = test_add (a, b, c, (mp_rnd_t) rnd_mode);
mpfr_init2 (s, 256);
- inex_s = test_add (s, b, c, GMP_RNDN); /* exact */
+ inex_s = test_add (s, b, c, MPFR_RNDN); /* exact */
if (inex_s)
{
printf ("check_1minuseps: result should have been exact "
@@ -829,239 +829,239 @@ tests (void)
check_case_2 ();
check64();
- check("293607738.0", "1.9967571564050541e-5", GMP_RNDU, 64, 53, 53,
+ check("293607738.0", "1.9967571564050541e-5", MPFR_RNDU, 64, 53, 53,
"2.9360773800002003e8");
- check("880524.0", "-2.0769715792901673e-5", GMP_RNDN, 64, 53, 53,
+ check("880524.0", "-2.0769715792901673e-5", MPFR_RNDN, 64, 53, 53,
"8.8052399997923023e5");
- check("1196426492.0", "-1.4218093058435347e-3", GMP_RNDN, 64, 53, 53,
+ check("1196426492.0", "-1.4218093058435347e-3", MPFR_RNDN, 64, 53, 53,
"1.1964264919985781e9");
- check("982013018.0", "-8.941829477291838e-7", GMP_RNDN, 64, 53, 53,
+ check("982013018.0", "-8.941829477291838e-7", MPFR_RNDN, 64, 53, 53,
"9.8201301799999905e8");
- check("1092583421.0", "1.0880649218158844e9", GMP_RNDN, 64, 53, 53,
+ check("1092583421.0", "1.0880649218158844e9", MPFR_RNDN, 64, 53, 53,
"2.1806483428158846e9");
- check("1.8476886419022969e-6", "961494401.0", GMP_RNDN, 53, 64, 53,
+ check("1.8476886419022969e-6", "961494401.0", MPFR_RNDN, 53, 64, 53,
"9.6149440100000179e8");
- check("-2.3222118418069868e5", "1229318102.0", GMP_RNDN, 53, 64, 53,
+ check("-2.3222118418069868e5", "1229318102.0", MPFR_RNDN, 53, 64, 53,
"1.2290858808158193e9");
- check("-3.0399171300395734e-6", "874924868.0", GMP_RNDN, 53, 64, 53,
+ check("-3.0399171300395734e-6", "874924868.0", MPFR_RNDN, 53, 64, 53,
"8.749248679999969e8");
- check("9.064246624706179e1", "663787413.0", GMP_RNDN, 53, 64, 53,
+ check("9.064246624706179e1", "663787413.0", MPFR_RNDN, 53, 64, 53,
"6.6378750364246619e8");
- check("-1.0954322421551264e2", "281806592.0", GMP_RNDD, 53, 64, 53,
+ check("-1.0954322421551264e2", "281806592.0", MPFR_RNDD, 53, 64, 53,
"2.8180648245677572e8");
- check("5.9836930386056659e-8", "1016217213.0", GMP_RNDN, 53, 64, 53,
+ check("5.9836930386056659e-8", "1016217213.0", MPFR_RNDN, 53, 64, 53,
"1.0162172130000001e9");
- check("-1.2772161928500301e-7", "1237734238.0", GMP_RNDN, 53, 64, 53,
+ check("-1.2772161928500301e-7", "1237734238.0", MPFR_RNDN, 53, 64, 53,
"1.2377342379999998e9");
- check("-4.567291988483277e8", "1262857194.0", GMP_RNDN, 53, 64, 53,
+ check("-4.567291988483277e8", "1262857194.0", MPFR_RNDN, 53, 64, 53,
"8.0612799515167236e8");
- check("4.7719471752925262e7", "196089880.0", GMP_RNDN, 53, 53, 53,
+ check("4.7719471752925262e7", "196089880.0", MPFR_RNDN, 53, 53, 53,
"2.4380935175292528e8");
- check("4.7719471752925262e7", "196089880.0", GMP_RNDN, 53, 64, 53,
+ check("4.7719471752925262e7", "196089880.0", MPFR_RNDN, 53, 64, 53,
"2.4380935175292528e8");
- check("-1.716113812768534e-140", "1271212614.0", GMP_RNDZ, 53, 64, 53,
+ check("-1.716113812768534e-140", "1271212614.0", MPFR_RNDZ, 53, 64, 53,
"1.2712126139999998e9");
- check("-1.2927455200185474e-50", "1675676122.0", GMP_RNDD, 53, 64, 53,
+ check("-1.2927455200185474e-50", "1675676122.0", MPFR_RNDD, 53, 64, 53,
"1.6756761219999998e9");
- check53("1.22191250737771397120e+20", "948002822.0", GMP_RNDN,
+ check53("1.22191250737771397120e+20", "948002822.0", MPFR_RNDN,
"122191250738719408128.0");
- check53("9966027674114492.0", "1780341389094537.0", GMP_RNDN,
+ check53("9966027674114492.0", "1780341389094537.0", MPFR_RNDN,
"11746369063209028.0");
check53("2.99280481918991653800e+272", "5.34637717585790933424e+271",
- GMP_RNDN, "3.5274425367757071711e272");
+ MPFR_RNDN, "3.5274425367757071711e272");
check_same();
check53("6.14384195492641560499e-02", "-6.14384195401037683237e-02",
- GMP_RNDU, "9.1603877261370314499e-12");
+ MPFR_RNDU, "9.1603877261370314499e-12");
check53("1.16809465359248765399e+196", "7.92883212101990665259e+196",
- GMP_RNDU, "9.0969267746123943065e196");
- check53("3.14553393112021279444e-67", "3.14553401015952024126e-67", GMP_RNDU,
+ MPFR_RNDU, "9.0969267746123943065e196");
+ check53("3.14553393112021279444e-67", "3.14553401015952024126e-67", MPFR_RNDU,
"6.2910679412797336946e-67");
- check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",GMP_RNDN,
+ check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDN,
"5.4388530464436950905e185");
- check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",GMP_RNDZ,
+ check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDZ,
"5.4388530464436944867e185");
- check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",GMP_RNDU,
+ check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDU,
"5.4388530464436950905e185");
- check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",GMP_RNDD,
+ check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDD,
"5.4388530464436944867e185");
check2b("1.001010101110011000000010100101110010111001010000001000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e358",187,
"-1.11100111001101100010001111111110101101110001000000000000000000000000000000000000000000e160",87,
"1.001010101110011000000010100101110010111001010000000111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111e358",178,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.111100100011100111010101010101001010100100000111001000000000000000000e481",70,
"1.1111000110100011110101111110110010010000000110101000000000000000e481",65,
"-1.001010111111101011010000001100011101100101000000000000000000e472",61,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("1.0100010111010000100101000000111110011100011001011010000000000000000000000000000000e516",83,
"-1.1001111000100001011100000001001100110011110010111111000000e541",59,
"-1.1001111000100001011011110111000001001011100000011110100000110001110011010011000000000000000000000000000000000000000000000000e541",125,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("-1.0010111100000100110001011011010000000011000111101000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e261",155,
"-1.00111110100011e239",15,
"-1.00101111000001001100101010101110001100110001111010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e261",159,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.110111000011111011000000001001111101101001010100111000000000000000000000000e880",76,
"-1.1010010e-634",8,
"-1.11011100001111101100000000100111110110100101010011100000000000000000000000e880",75,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("1.00100100110110101001010010101111000001011100100101010000000000000000000000000000e-530",81,
"-1.101101111100000111000011001010110011001011101001110100000e-908",58,
"1.00100100110110101001010010101111000001011100100101010e-530",54,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("1.0101100010010111101000000001000010010010011000111011000000000000000000000000000000000000000000000000000000000000000000e374",119,
"1.11100101100101e358",15,
"1.01011000100110011000010110100100100100100110001110110000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e374",150,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("-1.10011001000010100000010100100110110010011111101111000000000000000000000000000000000000000000000000000000000000000000e-172",117,
"1.111011100000101010110000100100110100100001001000011100000000e-173",61,
"-1.0100010000001001010110011011101001001011101011110001000000000000000e-173",68,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("-1.011110000111101011100001100110100011100101000011011000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e-189",175,
"1.1e631",2,
"1.011111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111e631",115,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("-1.101011101001011101100011001000001100010100001101011000000000000000000000000000000000000000000e-449",94,
"-1.01111101010111000011000110011101000111001100110111100000000000000e-429",66,
"-1.01111101010111000100110010000110100100101111111111101100010100001101011000000000000000000000000000000000000000e-429",111,
- GMP_RNDU);
+ MPFR_RNDU);
check2b("-1.101011101001011101100011001000001100010100001101011000000000000000000000000000000000000000000e-449",94,
"-1.01111101010111000011000110011101000111001100110111100000000000000e-429",66,
"-1.01111101010111000100110010000110100100101111111111101100010100001101011000000000000000000000000000000000000000e-429",111,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.1001000011101000110000111110010100100101110101111100000000000000000000000000000000000000000000000000000000e-72",107,
"-1.001100011101100100010101101010101011010010111111010000000000000000000000000000e521",79,
"-1.00110001110110010001010110101010101101001011111101000000000000000000000000000000000000000000000001e521",99,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.01010001111000000101010100100100110101011011100001110000000000e498",63,
"1.010000011010101111000100111100011100010101011110010100000000000e243",64,
"-1.010100011110000001010101001001001101010110111000011100000000000e498",64,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("1.00101100010101000011010000011000111111011110010111000000000000000000000000000000000000000000000000000000000e178",108,
"-1.10101101010101000110011011111001001101111111110000100000000e160",60,
"1.00101100010100111100100011000011111001000010011101110010000000001111100000000000000000000000000000000000e178",105,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("1.00110011010100111110011010110100111101110101100100110000000000000000000000000000000000000000000000e559",99,
"-1.011010110100111011100110100110011100000000111010011000000000000000e559",67,
"-1.101111111101011111111111001001100100011100001001100000000000000000000000000000000000000000000e556",94,
- GMP_RNDU);
+ MPFR_RNDU);
check2b("-1.100000111100101001100111011100011011000001101001111100000000000000000000000000e843",79,
"-1.1101101010110000001001000100001100110011000110110111000000000000000000000000000000000000000000e414",95,
"-1.1000001111001010011001110111000110110000011010100000e843",53,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.110110010110100010100011000110111001010000010111110000000000e-415",61,
"-1.0000100101100001111100110011111111110100011101101011000000000000000000e751",71,
"-1.00001001011000011111001100111111111101000111011010110e751",54,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("-1.1011011011110001001101010101001000010100010110111101000000000000000000000e258",74,
"-1.00011100010110110101001011000100100000100010101000010000000000000000000000000000000000000000000000e268",99,
"-1.0001110011001001000011110001000111010110101011110010011011110100000000000000000000000000000000000000e268",101,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.1011101010011101011000000100100110101101101110000001000000000e629",62,
"1.111111100000011100100011100000011101100110111110111000000000000000000000000000000000000000000e525",94,
"-1.101110101001110101100000010010011010110110111000000011111111111111111111111111111111111111111111111111101e629",106,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("1.111001000010001100010000001100000110001011110111011000000000000000000000000000000000000e152",88,
"1.111110111001100100000100111111010111000100111111001000000000000000e152",67,
"1.1110111111011110000010101001011011101010000110110100e153",53,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("1.000001100011110010110000110100001010101101111011110100e696",55,
"-1.1011001111011100100001011110100101010101110111010101000000000000000000000000000000000000000000000000000000000000e730",113,
"-1.1011001111011100100001011110100100010100010011100010e730",53,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("-1.11010111100001001111000001110101010010001111111001100000000000000000000000000000000000000000000000000000000000e530",111,
"1.01110100010010000000010110111101011101000001111101100000000000000000000000000000000000000000000000e530",99,
"-1.1000110011110011101010101101111101010011011111000000000000000e528",62,
- GMP_RNDD);
+ MPFR_RNDD);
check2b("-1.0001100010010100111101101011101000100100010011100011000000000000000000000000000000000000000000000000000000000e733",110,
"-1.001000000111110010100101010100110111001111011011001000000000000000000000000000000000000000000000000000000000e710",109,
"-1.000110001001010011111000111110110001110110011000110110e733",55,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("-1.1101011110000100111100000111010101001000111111100110000000000000000000000e530",74,
"1.01110100010010000000010110111101011101000001111101100000000000000000000000000000000000000000000000000000000000e530",111,
"-1.10001100111100111010101011011111010100110111110000000000000000000000000000e528",75,
- GMP_RNDU);
+ MPFR_RNDU);
check2b("1.00110011010100111110011010110100111101110101100100110000000000000000000000000000000000000000000000e559",99,
"-1.011010110100111011100110100110011100000000111010011000000000000000e559",67,
"-1.101111111101011111111111001001100100011100001001100000000000000000000000000000000000000000000e556",94,
- GMP_RNDU);
+ MPFR_RNDU);
check2b("-1.100101111110110000000110111111011010011101101111100100000000000000e-624",67,
"1.10111010101110100000010110101000000000010011100000100000000e-587",60,
"1.1011101010111010000001011010011111110100011110001011111111001000000100101100010010000011100000000000000000000e-587",110,
- GMP_RNDU);
+ MPFR_RNDU);
check2b("-1.10011001000010100000010100100110110010011111101111000000000000000000000000000000000000000000000000000000000000000000e-172",117,
"1.111011100000101010110000100100110100100001001000011100000000e-173",61,
"-1.0100010000001001010110011011101001001011101011110001000000000000000e-173",68,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("1.1000111000110010101001010011010011101100010110001001000000000000000000000000000000000000000000000000e167",101,
"1.0011110010000110000000101100100111000001110110110000000000000000000000000e167",74,
"1.01100101010111000101001111111111010101110001100111001000000000000000000000000000000000000000000000000000e168",105,
- GMP_RNDZ);
+ MPFR_RNDZ);
check2b("1.100101111111110010100101110111100001110000100001010000000000000000000000000000000000000000000000e808",97,
"-1.1110011001100000100000111111110000110010100111001011000000000000000000000000000000e807",83,
"1.01001001100110001100011111000000000001011010010111010000000000000000000000000000000000000000000e807",96,
- GMP_RNDN);
+ MPFR_RNDN);
check2b("1e128",128,
"1e0",128,
"100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001e0",256,
- GMP_RNDN);
+ MPFR_RNDN);
/* Checking double precision (53 bits) */
- check53("-8.22183238641455905806e-19", "7.42227178769761587878e-19",GMP_RNDD,
+ check53("-8.22183238641455905806e-19", "7.42227178769761587878e-19",MPFR_RNDD,
"-7.9956059871694317927e-20");
- check53("5.82106394662028628236e+234","-5.21514064202368477230e+89",GMP_RNDD,
+ check53("5.82106394662028628236e+234","-5.21514064202368477230e+89",MPFR_RNDD,
"5.8210639466202855763e234");
check53("5.72931679569871602371e+122","-5.72886070363264321230e+122",
- GMP_RNDN, "4.5609206607281141508e118");
+ MPFR_RNDN, "4.5609206607281141508e118");
check53("-5.09937369394650450820e+238", "2.70203299854862982387e+250",
- GMP_RNDD, "2.7020329985435301323e250");
- check53("-2.96695924472363684394e+27", "1.22842938251111500000e+16",GMP_RNDD,
+ MPFR_RNDD, "2.7020329985435301323e250");
+ check53("-2.96695924472363684394e+27", "1.22842938251111500000e+16",MPFR_RNDD,
"-2.96695924471135255027e27");
check53("1.74693641655743793422e-227", "-7.71776956366861843469e-229",
- GMP_RNDN, "1.669758720920751867e-227");
+ MPFR_RNDN, "1.669758720920751867e-227");
/* x = -7883040437021647.0; for (i=0; i<468; i++) x = x / 2.0;*/
check53("-1.03432206392780011159e-125", "1.30127034799251347548e-133",
- GMP_RNDN,
+ MPFR_RNDN,
"-1.0343220509150965661100887242027378881805094180354e-125");
- check53("1.05824655795525779205e+71", "-1.06022698059744327881e+71",GMP_RNDZ,
+ check53("1.05824655795525779205e+71", "-1.06022698059744327881e+71",MPFR_RNDZ,
"-1.9804226421854867632e68");
check53("-5.84204911040921732219e+240", "7.26658169050749590763e+240",
- GMP_RNDD, "1.4245325800982785854e240");
- check53("1.00944884131046636376e+221","2.33809162651471520268e+215",GMP_RNDN,
+ MPFR_RNDD, "1.4245325800982785854e240");
+ check53("1.00944884131046636376e+221","2.33809162651471520268e+215",MPFR_RNDN,
"1.0094511794020929787e221");
/*x = 7045852550057985.0; for (i=0; i<986; i++) x = x / 2.0;*/
check53("4.29232078932667367325e-278",
"1.0773525047389793833221116707010783793203080117586e-281"
- , GMP_RNDU, "4.2933981418314132787e-278");
- check53("5.27584773801377058681e-80", "8.91207657803547196421e-91", GMP_RNDN,
+ , MPFR_RNDU, "4.2933981418314132787e-278");
+ check53("5.27584773801377058681e-80", "8.91207657803547196421e-91", MPFR_RNDN,
"5.2758477381028917269e-80");
check53("2.99280481918991653800e+272", "5.34637717585790933424e+271",
- GMP_RNDN, "3.5274425367757071711e272");
+ MPFR_RNDN, "3.5274425367757071711e272");
check53("4.67302514390488041733e-184", "2.18321376145645689945e-190",
- GMP_RNDN, "4.6730273271186420541e-184");
- check53("5.57294120336300389254e+71", "2.60596167942024924040e+65", GMP_RNDZ,
+ MPFR_RNDN, "4.6730273271186420541e-184");
+ check53("5.57294120336300389254e+71", "2.60596167942024924040e+65", MPFR_RNDZ,
"5.5729438093246831053e71");
- check53("6.6052588496951015469e24", "4938448004894539.0", GMP_RNDU,
+ check53("6.6052588496951015469e24", "4938448004894539.0", MPFR_RNDU,
"6.6052588546335505068e24");
check53("1.23056185051606761523e-190", "1.64589756643433857138e-181",
- GMP_RNDU, "1.6458975676649006598e-181");
+ MPFR_RNDU, "1.6458975676649006598e-181");
check53("2.93231171510175981584e-280", "3.26266919161341483877e-273",
- GMP_RNDU, "3.2626694848445867288e-273");
- check53("5.76707395945001907217e-58", "4.74752971449827687074e-51", GMP_RNDD,
+ MPFR_RNDU, "3.2626694848445867288e-273");
+ check53("5.76707395945001907217e-58", "4.74752971449827687074e-51", MPFR_RNDD,
"4.747530291205672325e-51");
- check53("277363943109.0", "11.0", GMP_RNDN, "277363943120.0");
+ check53("277363943109.0", "11.0", MPFR_RNDN, "277363943120.0");
check53("1.44791789689198883921e-140", "-1.90982880222349071284e-121",
- GMP_RNDN, "-1.90982880222349071e-121");
+ MPFR_RNDN, "-1.90982880222349071e-121");
/* tests for particular cases (Vincent Lefevre, 22 Aug 2001) */
- check53("9007199254740992.0", "1.0", GMP_RNDN, "9007199254740992.0");
- check53("9007199254740994.0", "1.0", GMP_RNDN, "9007199254740996.0");
- check53("9007199254740992.0", "-1.0", GMP_RNDN, "9007199254740991.0");
- check53("9007199254740994.0", "-1.0", GMP_RNDN, "9007199254740992.0");
- check53("9007199254740996.0", "-1.0", GMP_RNDN, "9007199254740996.0");
+ check53("9007199254740992.0", "1.0", MPFR_RNDN, "9007199254740992.0");
+ check53("9007199254740994.0", "1.0", MPFR_RNDN, "9007199254740996.0");
+ check53("9007199254740992.0", "-1.0", MPFR_RNDN, "9007199254740991.0");
+ check53("9007199254740994.0", "-1.0", MPFR_RNDN, "9007199254740992.0");
+ check53("9007199254740996.0", "-1.0", MPFR_RNDN, "9007199254740996.0");
check_overflow ();
check_1111 ();