int
main (int argc, char *argv[])
{
tests_start_mpfr ();
check_nans ();
/* integer part is exact, frac. part is exact: return value should be 0 */
check ("61680","3.52935791015625e-1", "61680.352935791015625",
53, 53, 53, 0, MPFR_RNDZ);
/* integer part is rounded up, fractional part is rounded up: return value
should be 1+4*1=5 */
check ("-53968","-3.529052734375e-1", "-53970.352935791015625",
13, 13, 53, 5, MPFR_RNDZ);
/* integer part is rounded down, fractional part is rounded down:
return value should be 2+4*2=10 */
check ("61632","3.525390625e-1", "61648.352935791015625",
10, 10, 53, 10, MPFR_RNDZ);
check ("61680", "0", "61680", 53, 53, 53, 0, MPFR_RNDZ);
/* integer part is rounded up, fractional part is exact: 1 */
check ("-53968","0", "-53970", 13, 13, 53, 1, MPFR_RNDZ);
/* integer part is rounded up, fractional part is exact: 1 */
check ("-43392","0", "-43399", 13, 13, 53, 1, MPFR_RNDU);
/* integer part is rounded down, fractional part is exact: 2 */
check ("-52720","0", "-52719", 13, 13, 53, 2, MPFR_RNDD);
/* integer part is rounded down, fractional part is exact: 2 */
check ("61632", "0", "61648", 10, 10, 53, 2, MPFR_RNDZ);
check_special_exprange ();
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
int i;
mpfr_t x;
tests_start_mpfr ();
check_nans ();
/* check against values given by sinh(x), cosh(x) */
mpfr_init2 (x, 53);
mpfr_set_str (x, "FEDCBA987654321p-48", 16, MPFR_RNDN);
for (i = 0; i < 10; ++i)
{
/* x = i - x / 2 : boggle sign and bits */
mpfr_ui_sub (x, i, x, MPFR_RNDD);
mpfr_div_2ui (x, x, 2, MPFR_RNDD);
check (x, MPFR_RNDN);
check (x, MPFR_RNDU);
check (x, MPFR_RNDD);
}
mpfr_clear (x);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
mp_prec_t p;
unsigned k;
tests_start_mpfr ();
check_nans ();
special ();
for (p=2; p<100; p++)
for (k=0; k<100; k++)
check_two_sum (p);
check(1196426492, "1.4218093058435347e-3", GMP_RNDN,
"1.1964264919985781e9");
check(1092583421, "-1.0880649218158844e9", GMP_RNDN,
"2.1806483428158845901e9");
check(948002822, "1.22191250737771397120e+20", GMP_RNDN,
"-1.2219125073682338611e20");
check(832100416, "4.68311314939691330000e-215", GMP_RNDD,
"8.3210041599999988079e8");
check(1976245324, "1.25296395864546893357e+232", GMP_RNDZ,
"-1.2529639586454686577e232");
check(2128997392, "-1.08496826129284207724e+187", GMP_RNDU,
"1.0849682612928422704e187");
check(293607738, "-1.9967571564050541e-5", GMP_RNDU,
"2.9360773800002003e8");
check(354270183, "2.9469161763489528e3", GMP_RNDN,
"3.5426723608382362e8");
tests_end_mpfr ();
return 0;
}
/* tsin_cos prec [N] performs N tests with prec bits */
int
main (int argc, char *argv[])
{
tests_start_mpfr ();
if (argc > 1)
{
if (argc != 3 && argc != 4)
{
fprintf (stderr, "Usage: tsin_cos x prec [n]\n");
exit (1);
}
large_test (argv[1], atoi (argv[2]), (argc > 3) ? atoi (argv[3]) : 1);
goto end;
}
bug20091013 ();
bug20091008 ();
bug20091007 ();
bug20091122 ();
consistency ();
test_mpfr_sincos_fast ();
check_nans ();
/* worst case from PhD thesis of Vincent Lefe`vre: x=8980155785351021/2^54 */
check53 ("4.984987858808754279e-1", "4.781075595393330379e-1",
"8.783012931285841817e-1", MPFR_RNDN);
check53 ("4.984987858808754279e-1", "4.781075595393329824e-1",
"8.783012931285840707e-1", MPFR_RNDD);
check53 ("4.984987858808754279e-1", "4.781075595393329824e-1",
"8.783012931285840707e-1", MPFR_RNDZ);
check53 ("4.984987858808754279e-1", "4.781075595393330379e-1",
"8.783012931285841817e-1", MPFR_RNDU);
check53 ("1.00031274099908640274", "8.416399183372403892e-1",
"0.540039116973283217504", MPFR_RNDN);
check53 ("1.00229256850978698523", "8.427074524447979442e-1",
"0.538371757797526551137", MPFR_RNDZ);
check53 ("1.00288304857059840103", "8.430252033025980029e-1",
"0.537874062022526966409", MPFR_RNDZ);
check53 ("1.00591265847407274059", "8.446508805292128885e-1",
"0.53531755997839769456", MPFR_RNDN);
/* check one argument only */
check53sin ("1.00591265847407274059", "8.446508805292128885e-1", MPFR_RNDN);
check53cos ("1.00591265847407274059", "0.53531755997839769456", MPFR_RNDN);
overflowed_sin_cos0 ();
tiny ();
test20071214 ();
coverage_01032011 ();
end:
tests_end_mpfr ();
return 0;
}
int
main (void)
{
tests_start_mpfr ();
check_nans ();
check_regulars ();
test_generic (2, 1000, 100);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
tests_start_mpfr ();
check_nans ();
check_exact ();
check_float ();
check53("6.9314718055994530941514e-1", "0.0", MPFR_RNDZ, "0.0");
check53("0.0", "6.9314718055994530941514e-1", MPFR_RNDZ, "0.0");
check_sign();
check53("-4.165000000e4", "-0.00004801920768307322868063274915", MPFR_RNDN,
"2.0");
check53("2.71331408349172961467e-08", "-6.72658901114033715233e-165",
MPFR_RNDZ, "-1.8251348697787782844e-172");
check53("2.71331408349172961467e-08", "-6.72658901114033715233e-165",
MPFR_RNDA, "-1.8251348697787786e-172");
check53("0.31869277231188065", "0.88642843322303122", MPFR_RNDZ,
"2.8249833483992453642e-1");
check("8.47622108205396074254e-01", "3.24039313247872939883e-01", MPFR_RNDU,
28, 45, 2, "0.375");
check("8.47622108205396074254e-01", "3.24039313247872939883e-01", MPFR_RNDA,
28, 45, 2, "0.375");
check("2.63978122803639081440e-01", "6.8378615379333496093e-1", MPFR_RNDN,
34, 23, 31, "0.180504585267044603");
check("1.0", "0.11835170935876249132", MPFR_RNDU, 6, 41, 36,
"0.1183517093595583");
check53("67108865.0", "134217729.0", MPFR_RNDN, "9.007199456067584e15");
check("1.37399642157394197284e-01", "2.28877275604219221350e-01", MPFR_RNDN,
49, 15, 32, "0.0314472340833162888");
check("4.03160720978664954828e-01", "5.854828e-1"
/*"5.85483042917246621073e-01"*/, MPFR_RNDZ,
51, 22, 32, "0.2360436821472831");
check("3.90798504668055102229e-14", "9.85394674650308388664e-04", MPFR_RNDN,
46, 22, 12, "0.385027296503914762e-16");
check("4.58687081072827851358e-01", "2.20543551472118792844e-01", MPFR_RNDN,
49, 3, 2, "0.09375");
check_max();
check_min();
check_regression ();
test_generic (2, 500, 100);
data_check ("data/mulpi", mpfr_mulpi, "mpfr_mulpi");
valgrind20110503 ();
tests_end_mpfr ();
return 0;
}
int
main (void)
{
mpfr_t x, y, z;
double d;
int inexact;
tests_start_mpfr ();
/* check with enough precision */
mpfr_init2 (x, IEEE_DBL_MANT_DIG);
mpfr_init2 (y, IEEE_DBL_MANT_DIG);
mpfr_init2 (z, IEEE_DBL_MANT_DIG);
mpfr_set_str (y, "4096", 10, MPFR_RNDN);
d = 0.125;
mpfr_clear_flags ();
inexact = mpfr_d_sub (x, d, y, MPFR_RNDN);
if (inexact != 0)
{
printf ("Inexact flag error in mpfr_d_sub\n");
exit (1);
}
mpfr_set_str (z, "-4095.875", 10, MPFR_RNDN);
if (mpfr_cmp (z, x))
{
printf ("Error in mpfr_d_sub (");
mpfr_out_str (stdout, 10, 7, y, MPFR_RNDN);
printf (" + %.20g)\nexpected ", d);
mpfr_out_str (stdout, 10, 0, z, MPFR_RNDN);
printf ("\ngot ");
mpfr_out_str (stdout, 10, 0, x, MPFR_RNDN);
printf ("\n");
exit (1);
}
mpfr_clears (x, y, z, (mpfr_ptr) 0);
check_nans ();
test_generic (2, 1000, 100);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
tests_start_mpfr ();
check_nans ();
special ();
check3 ("-1.716113812768534e-140", 1271212614, MPFR_RNDZ,
"1.27121261399999976e9");
check3 ("1.22191250737771397120e+20", 948002822, MPFR_RNDN,
"122191250738719408128.0");
check3 ("-6.72658901114033715233e-165", 2000878121, MPFR_RNDZ,
"2.0008781209999997615e9");
check3 ("-2.0769715792901673e-5", 880524, MPFR_RNDN,
"8.8052399997923023e5");
test_generic_ui (2, 1000, 100);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
mpfr_prec_t p;
int k;
tests_start_mpfr ();
check_nans ();
for (p=2; p<200; p++)
for (k=0; k<200; k++)
check_two_sum (p);
check3 ("0.9999999999", 1, MPFR_RNDN,
"-10000000827403709990903735160827636718750e-50");
test_generic_ui (2, 1000, 100);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
tests_start_mpfr ();
check_nans ();
check ("61680","3.52935791015625e-1", "61680.352935791015625",
53, 53, 53, 0, GMP_RNDZ);
check ("-53968","-3.529052734375e-1", "-53970.352935791015625",
13, 13, 53, 2, GMP_RNDZ);
check ("61632","3.525390625e-1", "61648.352935791015625",
10, 10, 53, -2, GMP_RNDZ);
check ("61680", "0", "61680", 53, 53, 53, 0, GMP_RNDZ);
check ("-53968","0", "-53970", 13, 13, 53, 1, GMP_RNDZ);
check ("-43392","0", "-43399", 13, 13, 53, 1, GMP_RNDU);
check ("-52720","0", "-52719", 13, 13, 53, -1, GMP_RNDD);
check ("61632", "0", "61648", 10, 10, 53, -1, GMP_RNDZ);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char* argv[])
{
tests_start_mpfr ();
check_nans ();
check_large ();
check4 ("2.0", "1.0", MPFR_RNDN, "1.456791031046906869", -1);
check4 ("6.0", "4.0", MPFR_RNDN, "4.949360872472608925", 1);
check4 ("62.0", "61.0", MPFR_RNDN, "61.498983718845075902", -1);
check4 ("0.5", "1.0", MPFR_RNDN, "0.72839551552345343459", -1);
check4 ("1.0", "2.0", MPFR_RNDN, "1.456791031046906869", -1);
check4 ("234375765.0", "234375000.0", MPFR_RNDN, "234375382.49984394025", 1);
check4 ("8.0", "1.0", MPFR_RNDU, "3.61575617759736274873", 1);
check4 ("1.0", "44.0", MPFR_RNDU, "13.3658354512981243907", 1);
check4 ("1.0", "3.7252902984619140625e-9", MPFR_RNDU,
"0.07553933569711989657765", 1);
test_generic (2, 300, 17);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char* argv[])
{
tests_start_mpfr ();
check_nans ();
check_large ();
check4 ("2.0", "1.0", MPFR_RNDN, "1.45679103104690677029");
check4 ("6.0", "4.0", MPFR_RNDN, "4.94936087247260925182");
check4 ("62.0", "61.0", MPFR_RNDN, "6.14989837188450749750e+01");
check4 ("0.5", "1.0", MPFR_RNDN, "7.28395515523453385143e-01");
check4 ("1.0", "2.0", MPFR_RNDN, "1.45679103104690677029");
check4 ("234375765.0", "234375000.0", MPFR_RNDN, "2.3437538249984395504e8");
check4 ("8.0", "1.0", MPFR_RNDU, "3.615756177597362786");
check4 ("1.0", "44.0", MPFR_RNDU, "1.33658354512981247808e1");
check4 ("1.0", "3.7252902984619140625e-9", MPFR_RNDU,
"7.55393356971199025907e-02");
test_generic (2, 300, 17);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
mpfr_t x, y;
int inex;
tests_start_mpfr ();
special_overflow ();
check_nans ();
mpfr_init (x);
mpfr_init (y);
mpfr_set_prec (x, 53);
mpfr_set_prec (y, 2);
mpfr_set_str (x, "9.81333845856942e-1", 10, MPFR_RNDN);
test_cos (y, x, MPFR_RNDN);
mpfr_set_prec (x, 30);
mpfr_set_prec (y, 30);
mpfr_set_str_binary (x, "1.00001010001101110010100010101e-1");
test_cos (y, x, MPFR_RNDU);
mpfr_set_str_binary (x, "1.10111100010101011110101010100e-1");
if (mpfr_cmp (y, x))
{
printf ("Error for prec=30, rnd=MPFR_RNDU\n");
printf ("expected "); mpfr_print_binary (x); puts ("");
printf (" got "); mpfr_print_binary (y); puts ("");
exit (1);
}
mpfr_set_prec (x, 59);
mpfr_set_prec (y, 59);
mpfr_set_str_binary (x, "1.01101011101111010011111110111111111011011101100111100011e-3");
test_cos (y, x, MPFR_RNDU);
mpfr_set_str_binary (x, "1.1111011111110010001001001011100111101110100010000010010011e-1");
if (mpfr_cmp (y, x))
{
printf ("Error for prec=59, rnd=MPFR_RNDU\n");
printf ("expected "); mpfr_print_binary (x); puts ("");
printf (" got "); mpfr_print_binary (y); puts ("");
exit (1);
}
mpfr_set_prec (x, 5);
mpfr_set_prec (y, 5);
mpfr_set_str_binary (x, "1.1100e-2");
test_cos (y, x, MPFR_RNDD);
mpfr_set_str_binary (x, "1.1100e-1");
if (mpfr_cmp (y, x))
{
printf ("Error for x=1.1100e-2, rnd=MPFR_RNDD\n");
printf ("expected 1.1100e-1, got "); mpfr_print_binary (y); puts ("");
exit (1);
}
mpfr_set_prec (x, 32);
mpfr_set_prec (y, 32);
mpfr_set_str_binary (x, "0.10001000001001011000100001E-6");
mpfr_set_str_binary (y, "0.1111111111111101101111001100001");
test_cos (x, x, MPFR_RNDN);
if (mpfr_cmp (x, y))
{
printf ("Error for prec=32 (1)\n");
exit (1);
}
mpfr_set_str_binary (x, "-0.1101011110111100111010011001011E-1");
mpfr_set_str_binary (y, "0.11101001100110111011011010100011");
test_cos (x, x, MPFR_RNDN);
if (mpfr_cmp (x, y))
{
printf ("Error for prec=32 (2)\n");
exit (1);
}
/* huge argument reduction */
mpfr_set_str_binary (x, "0.10000010000001101011101111001011E40");
mpfr_set_str_binary (y, "0.10011000001111010000101011001011E-1");
test_cos (x, x, MPFR_RNDN);
if (mpfr_cmp (x, y))
{
printf ("Error for prec=32 (3)\n");
exit (1);
}
mpfr_set_prec (x, 3);
mpfr_set_prec (y, 3);
mpfr_set_str_binary (x, "0.110E60");
inex = mpfr_cos (y, x, MPFR_RNDD);
MPFR_ASSERTN(inex < 0);
/* worst case from PhD thesis of Vincent Lefe`vre: x=8980155785351021/2^54 */
check53 ("4.984987858808754279e-1", "8.783012931285841817e-1", MPFR_RNDN);
check53 ("4.984987858808754279e-1", "8.783012931285840707e-1", MPFR_RNDD);
check53 ("4.984987858808754279e-1", "8.783012931285840707e-1", MPFR_RNDZ);
check53 ("4.984987858808754279e-1", "8.783012931285841817e-1", MPFR_RNDU);
check53 ("1.00031274099908640274", "0.540039116973283217504", MPFR_RNDN);
check53 ("1.00229256850978698523", "0.538371757797526551137", MPFR_RNDZ);
check53 ("1.00288304857059840103", "0.537874062022526966409", MPFR_RNDZ);
check53 ("1.00591265847407274059", "0.53531755997839769456", MPFR_RNDN);
check53 ("1.00591265847407274059", "0.53531755997839769456", MPFR_RNDN);
overflowed_cos0 ();
test_generic (2, 100, 15);
/* check inexact flag */
mpfr_set_prec (x, 3);
mpfr_set_prec (y, 13);
mpfr_set_str_binary (x, "-0.100E196");
inex = mpfr_cos (y, x, MPFR_RNDU);
mpfr_set_prec (x, 13);
mpfr_set_str_binary (x, "0.1111111100101");
MPFR_ASSERTN (inex > 0 && mpfr_equal_p (x, y));
mpfr_clear (x);
mpfr_clear (y);
bug20091030 ();
data_check ("data/cos", mpfr_cos, "mpfr_cos");
bad_cases (mpfr_cos, mpfr_acos, "mpfr_cos", 256, -40, 0, 4, 128, 800, 50);
tests_end_mpfr ();
return 0;
}
int
main (int argc, char *argv[])
{
mpfr_t x, c, s, c2, s2;
tests_start_mpfr ();
check_regression ();
check_nans ();
/* worst case from PhD thesis of Vincent Lefe`vre: x=8980155785351021/2^54 */
check53 ("4.984987858808754279e-1", "4.781075595393330379e-1", MPFR_RNDN);
check53 ("4.984987858808754279e-1", "4.781075595393329824e-1", MPFR_RNDD);
check53 ("4.984987858808754279e-1", "4.781075595393329824e-1", MPFR_RNDZ);
check53 ("4.984987858808754279e-1", "4.781075595393330379e-1", MPFR_RNDU);
check53 ("1.00031274099908640274", "8.416399183372403892e-1", MPFR_RNDN);
check53 ("1.00229256850978698523", "8.427074524447979442e-1", MPFR_RNDZ);
check53 ("1.00288304857059840103", "8.430252033025980029e-1", MPFR_RNDZ);
check53 ("1.00591265847407274059", "8.446508805292128885e-1", MPFR_RNDN);
/* Other worst cases showing a bug introduced on 2005-01-29 in rev 3248 */
check53b ("1.0111001111010111010111111000010011010001110001111011e-21",
"1.0111001111010111010111111000010011010001101001110001e-21",
MPFR_RNDU);
check53b ("1.1011101111111010000001010111000010000111100100101101",
"1.1111100100101100001111100000110011110011010001010101e-1",
MPFR_RNDU);
mpfr_init2 (x, 2);
mpfr_set_str (x, "0.5", 10, MPFR_RNDN);
test_sin (x, x, MPFR_RNDD);
if (mpfr_cmp_ui_2exp (x, 3, -3)) /* x != 0.375 = 3/8 */
{
printf ("mpfr_sin(0.5, MPFR_RNDD) failed with precision=2\n");
exit (1);
}
/* bug found by Kevin Ryde */
mpfr_const_pi (x, MPFR_RNDN);
mpfr_mul_ui (x, x, 3L, MPFR_RNDN);
mpfr_div_ui (x, x, 2L, MPFR_RNDN);
test_sin (x, x, MPFR_RNDN);
if (mpfr_cmp_ui (x, 0) >= 0)
{
printf ("Error: wrong sign for sin(3*Pi/2)\n");
exit (1);
}
/* Can fail on an assert */
mpfr_set_prec (x, 53);
mpfr_set_str (x, "77291789194529019661184401408", 10, MPFR_RNDN);
mpfr_init2 (c, 4); mpfr_init2 (s, 42);
mpfr_init2 (c2, 4); mpfr_init2 (s2, 42);
test_sin (s, x, MPFR_RNDN);
mpfr_cos (c, x, MPFR_RNDN);
mpfr_sin_cos (s2, c2, x, MPFR_RNDN);
if (mpfr_cmp (c2, c))
{
printf("cos differs for x=77291789194529019661184401408");
exit (1);
}
if (mpfr_cmp (s2, s))
{
printf("sin differs for x=77291789194529019661184401408");
exit (1);
}
mpfr_set_str_binary (x, "1.1001001000011111101101010100010001000010110100010011");
test_sin (x, x, MPFR_RNDZ);
if (mpfr_cmp_str (x, "1.1111111111111111111111111111111111111111111111111111e-1", 2, MPFR_RNDN))
{
printf ("Error for x= 1.1001001000011111101101010100010001000010110100010011\nGot ");
mpfr_dump (x);
exit (1);
}
mpfr_set_prec (s, 9);
mpfr_set_prec (x, 190);
mpfr_const_pi (x, MPFR_RNDN);
mpfr_sin (s, x, MPFR_RNDZ);
if (mpfr_cmp_str (s, "0.100000101e-196", 2, MPFR_RNDN))
{
printf ("Error for x ~= pi\n");
mpfr_dump (s);
exit (1);
}
mpfr_clear (s2);
mpfr_clear (c2);
mpfr_clear (s);
mpfr_clear (c);
mpfr_clear (x);
test_generic (2, 100, 15);
test_generic (MPFR_SINCOS_THRESHOLD-1, MPFR_SINCOS_THRESHOLD+1, 2);
test_sign ();
check_tiny ();
data_check ("data/sin", mpfr_sin, "mpfr_sin");
bad_cases (mpfr_sin, mpfr_asin, "mpfr_sin", 256, -40, 0, 4, 128, 800, 50);
tests_end_mpfr ();
return 0;
}
static void
tests (void)
{
check_alloc ();
check_nans ();
check_inexact ();
check_case_1b ();
check_case_2 ();
check64();
coverage_bk_eq_0 ();
check("293607738.0", "1.9967571564050541e-5", MPFR_RNDU, 64, 53, 53,
"2.9360773800002003e8");
check("880524.0", "-2.0769715792901673e-5", MPFR_RNDN, 64, 53, 53,
"8.8052399997923023e5");
check("1196426492.0", "-1.4218093058435347e-3", MPFR_RNDN, 64, 53, 53,
"1.1964264919985781e9");
check("982013018.0", "-8.941829477291838e-7", MPFR_RNDN, 64, 53, 53,
"9.8201301799999905e8");
check("1092583421.0", "1.0880649218158844e9", MPFR_RNDN, 64, 53, 53,
"2.1806483428158846e9");
check("1.8476886419022969e-6", "961494401.0", MPFR_RNDN, 53, 64, 53,
"9.6149440100000179e8");
check("-2.3222118418069868e5", "1229318102.0", MPFR_RNDN, 53, 64, 53,
"1.2290858808158193e9");
check("-3.0399171300395734e-6", "874924868.0", MPFR_RNDN, 53, 64, 53,
"8.749248679999969e8");
check("9.064246624706179e1", "663787413.0", MPFR_RNDN, 53, 64, 53,
"6.6378750364246619e8");
check("-1.0954322421551264e2", "281806592.0", MPFR_RNDD, 53, 64, 53,
"2.8180648245677572e8");
check("5.9836930386056659e-8", "1016217213.0", MPFR_RNDN, 53, 64, 53,
"1.0162172130000001e9");
check("-1.2772161928500301e-7", "1237734238.0", MPFR_RNDN, 53, 64, 53,
"1.2377342379999998e9");
check("-4.567291988483277e8", "1262857194.0", MPFR_RNDN, 53, 64, 53,
"8.0612799515167236e8");
check("4.7719471752925262e7", "196089880.0", MPFR_RNDN, 53, 53, 53,
"2.4380935175292528e8");
check("4.7719471752925262e7", "196089880.0", MPFR_RNDN, 53, 64, 53,
"2.4380935175292528e8");
check("-1.716113812768534e-140", "1271212614.0", MPFR_RNDZ, 53, 64, 53,
"1.2712126139999998e9");
check("-1.2927455200185474e-50", "1675676122.0", MPFR_RNDD, 53, 64, 53,
"1.6756761219999998e9");
check53("1.22191250737771397120e+20", "948002822.0", MPFR_RNDN,
"122191250738719408128.0");
check53("9966027674114492.0", "1780341389094537.0", MPFR_RNDN,
"11746369063209028.0");
check53("2.99280481918991653800e+272", "5.34637717585790933424e+271",
MPFR_RNDN, "3.5274425367757071711e272");
check_same();
check53("6.14384195492641560499e-02", "-6.14384195401037683237e-02",
MPFR_RNDU, "9.1603877261370314499e-12");
check53("1.16809465359248765399e+196", "7.92883212101990665259e+196",
MPFR_RNDU, "9.0969267746123943065e196");
check53("3.14553393112021279444e-67", "3.14553401015952024126e-67", MPFR_RNDU,
"6.2910679412797336946e-67");
check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDN,
"5.4388530464436950905e185");
check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDZ,
"5.4388530464436944867e185");
check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDU,
"5.4388530464436950905e185");
check53("5.43885304644369509058e+185","-1.87427265794105342763e-57",MPFR_RNDD,
"5.4388530464436944867e185");
check2b("1.001010101110011000000010100101110010111001010000001000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e358",187,
"-1.11100111001101100010001111111110101101110001000000000000000000000000000000000000000000e160",87,
"1.001010101110011000000010100101110010111001010000000111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111e358",178,
MPFR_RNDD);
check2b("-1.111100100011100111010101010101001010100100000111001000000000000000000e481",70,
"1.1111000110100011110101111110110010010000000110101000000000000000e481",65,
"-1.001010111111101011010000001100011101100101000000000000000000e472",61,
MPFR_RNDD);
check2b("1.0100010111010000100101000000111110011100011001011010000000000000000000000000000000e516",83,
"-1.1001111000100001011100000001001100110011110010111111000000e541",59,
"-1.1001111000100001011011110111000001001011100000011110100000110001110011010011000000000000000000000000000000000000000000000000e541",125,
MPFR_RNDZ);
check2b("-1.0010111100000100110001011011010000000011000111101000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e261",155,
"-1.00111110100011e239",15,
"-1.00101111000001001100101010101110001100110001111010000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e261",159,
MPFR_RNDD);
check2b("-1.110111000011111011000000001001111101101001010100111000000000000000000000000e880",76,
"-1.1010010e-634",8,
"-1.11011100001111101100000000100111110110100101010011100000000000000000000000e880",75,
MPFR_RNDZ);
check2b("1.00100100110110101001010010101111000001011100100101010000000000000000000000000000e-530",81,
"-1.101101111100000111000011001010110011001011101001110100000e-908",58,
"1.00100100110110101001010010101111000001011100100101010e-530",54,
MPFR_RNDN);
check2b("1.0101100010010111101000000001000010010010011000111011000000000000000000000000000000000000000000000000000000000000000000e374",119,
"1.11100101100101e358",15,
"1.01011000100110011000010110100100100100100110001110110000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e374",150,
MPFR_RNDZ);
check2b("-1.10011001000010100000010100100110110010011111101111000000000000000000000000000000000000000000000000000000000000000000e-172",117,
"1.111011100000101010110000100100110100100001001000011100000000e-173",61,
"-1.0100010000001001010110011011101001001011101011110001000000000000000e-173",68,
MPFR_RNDZ);
check2b("-1.011110000111101011100001100110100011100101000011011000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000e-189",175,
"1.1e631",2,
"1.011111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111e631",115,
MPFR_RNDZ);
check2b("-1.101011101001011101100011001000001100010100001101011000000000000000000000000000000000000000000e-449",94,
"-1.01111101010111000011000110011101000111001100110111100000000000000e-429",66,
"-1.01111101010111000100110010000110100100101111111111101100010100001101011000000000000000000000000000000000000000e-429",111,
MPFR_RNDU);
check2b("-1.101011101001011101100011001000001100010100001101011000000000000000000000000000000000000000000e-449",94,
"-1.01111101010111000011000110011101000111001100110111100000000000000e-429",66,
"-1.01111101010111000100110010000110100100101111111111101100010100001101011000000000000000000000000000000000000000e-429",111,
MPFR_RNDD);
check2b("-1.1001000011101000110000111110010100100101110101111100000000000000000000000000000000000000000000000000000000e-72",107,
"-1.001100011101100100010101101010101011010010111111010000000000000000000000000000e521",79,
"-1.00110001110110010001010110101010101101001011111101000000000000000000000000000000000000000000000001e521",99,
MPFR_RNDD);
check2b("-1.01010001111000000101010100100100110101011011100001110000000000e498",63,
"1.010000011010101111000100111100011100010101011110010100000000000e243",64,
"-1.010100011110000001010101001001001101010110111000011100000000000e498",64,
MPFR_RNDN);
check2b("1.00101100010101000011010000011000111111011110010111000000000000000000000000000000000000000000000000000000000e178",108,
"-1.10101101010101000110011011111001001101111111110000100000000e160",60,
"1.00101100010100111100100011000011111001000010011101110010000000001111100000000000000000000000000000000000e178",105,
MPFR_RNDN);
check2b("1.00110011010100111110011010110100111101110101100100110000000000000000000000000000000000000000000000e559",99,
"-1.011010110100111011100110100110011100000000111010011000000000000000e559",67,
"-1.101111111101011111111111001001100100011100001001100000000000000000000000000000000000000000000e556",94,
MPFR_RNDU);
check2b("-1.100000111100101001100111011100011011000001101001111100000000000000000000000000e843",79,
"-1.1101101010110000001001000100001100110011000110110111000000000000000000000000000000000000000000e414",95,
"-1.1000001111001010011001110111000110110000011010100000e843",53,
MPFR_RNDD);
check2b("-1.110110010110100010100011000110111001010000010111110000000000e-415",61,
"-1.0000100101100001111100110011111111110100011101101011000000000000000000e751",71,
"-1.00001001011000011111001100111111111101000111011010110e751",54,
MPFR_RNDN);
check2b("-1.1011011011110001001101010101001000010100010110111101000000000000000000000e258",74,
"-1.00011100010110110101001011000100100000100010101000010000000000000000000000000000000000000000000000e268",99,
"-1.0001110011001001000011110001000111010110101011110010011011110100000000000000000000000000000000000000e268",101,
MPFR_RNDD);
check2b("-1.1011101010011101011000000100100110101101101110000001000000000e629",62,
"1.111111100000011100100011100000011101100110111110111000000000000000000000000000000000000000000e525",94,
"-1.101110101001110101100000010010011010110110111000000011111111111111111111111111111111111111111111111111101e629",106,
MPFR_RNDD);
check2b("1.111001000010001100010000001100000110001011110111011000000000000000000000000000000000000e152",88,
"1.111110111001100100000100111111010111000100111111001000000000000000e152",67,
"1.1110111111011110000010101001011011101010000110110100e153",53,
MPFR_RNDN);
check2b("1.000001100011110010110000110100001010101101111011110100e696",55,
"-1.1011001111011100100001011110100101010101110111010101000000000000000000000000000000000000000000000000000000000000e730",113,
"-1.1011001111011100100001011110100100010100010011100010e730",53,
MPFR_RNDN);
check2b("-1.11010111100001001111000001110101010010001111111001100000000000000000000000000000000000000000000000000000000000e530",111,
"1.01110100010010000000010110111101011101000001111101100000000000000000000000000000000000000000000000e530",99,
"-1.1000110011110011101010101101111101010011011111000000000000000e528",62,
MPFR_RNDD);
check2b("-1.0001100010010100111101101011101000100100010011100011000000000000000000000000000000000000000000000000000000000e733",110,
"-1.001000000111110010100101010100110111001111011011001000000000000000000000000000000000000000000000000000000000e710",109,
"-1.000110001001010011111000111110110001110110011000110110e733",55,
MPFR_RNDN);
check2b("-1.1101011110000100111100000111010101001000111111100110000000000000000000000e530",74,
"1.01110100010010000000010110111101011101000001111101100000000000000000000000000000000000000000000000000000000000e530",111,
"-1.10001100111100111010101011011111010100110111110000000000000000000000000000e528",75,
MPFR_RNDU);
check2b("1.00110011010100111110011010110100111101110101100100110000000000000000000000000000000000000000000000e559",99,
"-1.011010110100111011100110100110011100000000111010011000000000000000e559",67,
"-1.101111111101011111111111001001100100011100001001100000000000000000000000000000000000000000000e556",94,
MPFR_RNDU);
check2b("-1.100101111110110000000110111111011010011101101111100100000000000000e-624",67,
"1.10111010101110100000010110101000000000010011100000100000000e-587",60,
"1.1011101010111010000001011010011111110100011110001011111111001000000100101100010010000011100000000000000000000e-587",110,
MPFR_RNDU);
check2b("-1.10011001000010100000010100100110110010011111101111000000000000000000000000000000000000000000000000000000000000000000e-172",117,
"1.111011100000101010110000100100110100100001001000011100000000e-173",61,
"-1.0100010000001001010110011011101001001011101011110001000000000000000e-173",68,
MPFR_RNDZ);
check2b("1.1000111000110010101001010011010011101100010110001001000000000000000000000000000000000000000000000000e167",101,
"1.0011110010000110000000101100100111000001110110110000000000000000000000000e167",74,
"1.01100101010111000101001111111111010101110001100111001000000000000000000000000000000000000000000000000000e168",105,
MPFR_RNDZ);
check2b("1.100101111111110010100101110111100001110000100001010000000000000000000000000000000000000000000000e808",97,
"-1.1110011001100000100000111111110000110010100111001011000000000000000000000000000000e807",83,
"1.01001001100110001100011111000000000001011010010111010000000000000000000000000000000000000000000e807",96,
MPFR_RNDN);
check2b("1e128",128,
"1e0",128,
"100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001e0",256,
MPFR_RNDN);
/* Checking double precision (53 bits) */
check53("-8.22183238641455905806e-19", "7.42227178769761587878e-19",MPFR_RNDD,
"-7.9956059871694317927e-20");
check53("5.82106394662028628236e+234","-5.21514064202368477230e+89",MPFR_RNDD,
"5.8210639466202855763e234");
check53("5.72931679569871602371e+122","-5.72886070363264321230e+122",
MPFR_RNDN, "4.5609206607281141508e118");
check53("-5.09937369394650450820e+238", "2.70203299854862982387e+250",
MPFR_RNDD, "2.7020329985435301323e250");
check53("-2.96695924472363684394e+27", "1.22842938251111500000e+16",MPFR_RNDD,
"-2.96695924471135255027e27");
check53("1.74693641655743793422e-227", "-7.71776956366861843469e-229",
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",
MPFR_RNDN,
"-1.0343220509150965661100887242027378881805094180354e-125");
check53("1.05824655795525779205e+71", "-1.06022698059744327881e+71",MPFR_RNDZ,
"-1.9804226421854867632e68");
check53("-5.84204911040921732219e+240", "7.26658169050749590763e+240",
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"
, MPFR_RNDU, "4.2933981418314132787e-278");
check53("5.27584773801377058681e-80", "8.91207657803547196421e-91", MPFR_RNDN,
"5.2758477381028917269e-80");
check53("2.99280481918991653800e+272", "5.34637717585790933424e+271",
MPFR_RNDN, "3.5274425367757071711e272");
check53("4.67302514390488041733e-184", "2.18321376145645689945e-190",
MPFR_RNDN, "4.6730273271186420541e-184");
check53("5.57294120336300389254e+71", "2.60596167942024924040e+65", MPFR_RNDZ,
"5.5729438093246831053e71");
check53("6.6052588496951015469e24", "4938448004894539.0", MPFR_RNDU,
"6.6052588546335505068e24");
check53("1.23056185051606761523e-190", "1.64589756643433857138e-181",
MPFR_RNDU, "1.6458975676649006598e-181");
check53("2.93231171510175981584e-280", "3.26266919161341483877e-273",
MPFR_RNDU, "3.2626694848445867288e-273");
check53("5.76707395945001907217e-58", "4.74752971449827687074e-51", MPFR_RNDD,
"4.747530291205672325e-51");
check53("277363943109.0", "11.0", MPFR_RNDN, "277363943120.0");
check53("1.44791789689198883921e-140", "-1.90982880222349071284e-121",
MPFR_RNDN, "-1.90982880222349071e-121");
/* tests for particular cases (Vincent Lefevre, 22 Aug 2001) */
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 ();
check_1minuseps ();
}
int
main (int argc, char *argv[])
{
mpfr_t x;
unsigned int i;
unsigned int prec[10] = {14, 15, 19, 22, 23, 24, 25, 40, 41, 52};
unsigned int prec2[10] = {4, 5, 6, 19, 70, 95, 100, 106, 107, 108};
tests_start_mpfr ();
check_nans ();
mpfr_init (x);
mpfr_set_prec (x, 2);
mpfr_set_str (x, "0.5", 10, MPFR_RNDN);
mpfr_tan (x, x, MPFR_RNDD);
if (mpfr_cmp_ui_2exp(x, 1, -1))
{
printf ("mpfr_tan(0.5, MPFR_RNDD) failed\n"
"expected 0.5, got");
mpfr_print_binary(x);
putchar('\n');
exit (1);
}
/* check that tan(3*Pi/4) ~ -1 */
for (i=0; i<10; i++)
{
mpfr_set_prec (x, prec[i]);
mpfr_const_pi (x, MPFR_RNDN);
mpfr_mul_ui (x, x, 3, MPFR_RNDN);
mpfr_div_ui (x, x, 4, MPFR_RNDN);
mpfr_tan (x, x, MPFR_RNDN);
if (mpfr_cmp_si (x, -1))
{
printf ("tan(3*Pi/4) fails for prec=%u\n", prec[i]);
exit (1);
}
}
/* check that tan(7*Pi/4) ~ -1 */
for (i=0; i<10; i++)
{
mpfr_set_prec (x, prec2[i]);
mpfr_const_pi (x, MPFR_RNDN);
mpfr_mul_ui (x, x, 7, MPFR_RNDN);
mpfr_div_ui (x, x, 4, MPFR_RNDN);
mpfr_tan (x, x, MPFR_RNDN);
if (mpfr_cmp_si (x, -1))
{
printf ("tan(3*Pi/4) fails for prec=%u\n", prec2[i]);
exit (1);
}
}
mpfr_clear (x);
test_generic (2, 100, 10);
data_check ("data/tan", mpfr_tan, "mpfr_tan");
bad_cases (mpfr_tan, mpfr_atan, "mpfr_tan", 256, -256, 255, 4, 128, 800, 40);
tests_end_mpfr ();
return 0;
}
