static void
coverage_01032011 (void)
{
mpfr_t val, cval, sval, svalf;
int status_f, status;
mpfr_init2 (val, MPFR_PREC_MIN);
mpfr_init2 (cval, MPFR_PREC_MIN);
mpfr_init2 (sval, MPFR_PREC_MIN);
mpfr_init2 (svalf, MPFR_PREC_MIN);
mpfr_set_str1 (val, "-0.7");
status_f = mpfr_sincos_fast (svalf, NULL, val, MPFR_RNDN);
status = mpfr_sin_cos (sval, cval, val, MPFR_RNDN);
if (! mpfr_equal_p (svalf, sval) || SIGN (status_f) != SIGN (status))
{
printf ("mpfr_sincos_fast differ from mpfr_sin_cos result:\n"
" sin fast is ");
mpfr_dump (svalf);
printf (" sin is ");
mpfr_dump (sval);
printf ("status_f = %d, status = %d\n", status_f, status);
exit (1);
}
mpfr_clears(val, cval, sval, svalf, (mpfr_ptr) 0);
}
/* check that mpfr_sin_cos and test_mpfr_sincos_fast agree */
static void
test_mpfr_sincos_fast (void)
{
mpfr_t x, y, z, yref, zref, h;
mpfr_prec_t p = 1000;
int i, inex, inexref;
mpfr_rnd_t r;
mpfr_init2 (x, p);
mpfr_init2 (y, p);
mpfr_init2 (z, p);
mpfr_init2 (yref, p);
mpfr_init2 (zref, p);
mpfr_init2 (h, p);
mpfr_set_ui (x, 0, MPFR_RNDN);
/* we generate a random value x, compute sin(x) and cos(x) with both
mpfr_sin_cos and mpfr_sincos_fast, and check the values and the flags
agree */
for (i = 0; i < 100; i++)
{
mpfr_urandomb (h, RANDS);
mpfr_add (x, x, h, MPFR_RNDN);
r = RND_RAND ();
inexref = mpfr_sin_cos (yref, zref, x, r);
inex = mpfr_sincos_fast (y, z, x, r);
if (mpfr_cmp (y, yref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree\n");
printf ("x="); mpfr_dump (x);
printf ("rnd=%s\n", mpfr_print_rnd_mode (r));
printf ("yref="); mpfr_dump (yref);
printf ("y="); mpfr_dump (y);
exit (1);
}
if (mpfr_cmp (z, zref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree\n");
printf ("x="); mpfr_dump (x);
printf ("rnd=%s\n", mpfr_print_rnd_mode (r));
printf ("zref="); mpfr_dump (zref);
printf ("z="); mpfr_dump (z);
exit (1);
}
if (inex != inexref)
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree\n");
printf ("x="); mpfr_dump (x);
printf ("rnd=%s\n", mpfr_print_rnd_mode (r));
printf ("inexref=%d inex=%d\n", inexref, inex);
exit (1);
}
}
mpfr_clear (x);
mpfr_clear (y);
mpfr_clear (z);
mpfr_clear (yref);
mpfr_clear (zref);
mpfr_clear (h);
}
static int
mpfr_cos_fast (mpfr_ptr y, mpfr_srcptr x, mpfr_rnd_t rnd_mode)
{
int inex;
inex = mpfr_sincos_fast (NULL, y, x, rnd_mode);
inex = inex >> 2; /* 0: exact, 1: rounded up, 2: rounded down */
return (inex == 2) ? -1 : inex;
}
static void
bug20091013 (void)
{
mpfr_t x, y, z, yref, zref;
mpfr_prec_t p = 1000;
int inex, inexref;
mpfr_rnd_t r = MPFR_RNDN;
mpfr_init2 (x, p);
mpfr_init2 (y, p);
mpfr_init2 (z, p);
mpfr_init2 (yref, p);
mpfr_init2 (zref, p);
mpfr_set_str (x, "3.240ff3fdcb1ee7cd667b96287593ae24e20fb63ed7c2d5bf4bd0f2cc5509283b04e7628e66382605f14ed5967cef15296041539a1bdaa626c777c7fbb6f2068414759b78cee14f37848689b3a170f583656be4e0837f464d8210556a3a822d4ecfdd59f4e0d5fdb76bf7e15b8a57234e2160b98e14c17bbdf27c4643b8@1", 16, MPFR_RNDN);
inexref = mpfr_sin_cos (yref, zref, x, r);
inex = mpfr_sincos_fast (y, z, x, r);
if (mpfr_cmp (y, yref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091013)\n");
printf ("yref="); mpfr_dump (yref);
printf ("y="); mpfr_dump (y);
exit (1);
}
if (mpfr_cmp (z, zref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091013)\n");
printf ("zref="); mpfr_dump (zref);
printf ("z="); mpfr_dump (z);
exit (1);
}
/* sin(x) is rounded down and cos(x) is rounded down, thus we should get
2+4*2 = 10 as return value */
if (inex != inexref)
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091013)\n");
printf ("inexref=%d inex=%d (10 expected)\n", inexref, inex);
exit (1);
}
mpfr_clear (x);
mpfr_clear (y);
mpfr_clear (z);
mpfr_clear (yref);
mpfr_clear (zref);
}
/* Note: with the sin_cos.c code before r6507, the disagreement occurs
only on the return ("inexact") value, which is new in r6444. */
static void
bug20091008 (void)
{
mpfr_t x, y, z, yref, zref;
mpfr_prec_t p = 1000;
int inex, inexref;
mpfr_rnd_t r = MPFR_RNDN;
mpfr_init2 (x, p);
mpfr_init2 (y, p);
mpfr_init2 (z, p);
mpfr_init2 (yref, p);
mpfr_init2 (zref, p);
mpfr_set_str (x, "c.91813724e28ef6a711d33e6505984699daef7fe93636c1ed5d0168bc96989cc6802f7f9e405c902ec62fb90cd39c9d21084c8ad8b5af4c4aa87bf402e2e4a78e6fe1ffeb6dbbbdbbc2983c196c518966ccc1e094ed39ee77984ef2428069d65de37928e75247edbe7007245e682616b5ebbf05f2fdefc74ad192024f10", 16, MPFR_RNDN);
inexref = mpfr_sin_cos (yref, zref, x, r);
inex = mpfr_sincos_fast (y, z, x, r);
if (mpfr_cmp (y, yref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091008)\n");
printf ("yref="); mpfr_dump (yref);
printf ("y="); mpfr_dump (y);
exit (1);
}
if (mpfr_cmp (z, zref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091008)\n");
printf ("zref="); mpfr_dump (zref);
printf ("z="); mpfr_dump (z);
exit (1);
}
/* sin(x) is rounded up, cos(x) is rounded up too, thus we should get 5
for the return value */
if (inex != inexref)
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091008)\n");
printf ("inexref=%d inex=%d (5 expected)\n", inexref, inex);
exit (1);
}
mpfr_clear (x);
mpfr_clear (y);
mpfr_clear (z);
mpfr_clear (yref);
mpfr_clear (zref);
}
static void
bug20091007 (void)
{
mpfr_t x, y, z, yref, zref;
mpfr_prec_t p = 1000;
int inex, inexref;
mpfr_rnd_t r = MPFR_RNDZ;
mpfr_init2 (x, p);
mpfr_init2 (y, p);
mpfr_init2 (z, p);
mpfr_init2 (yref, p);
mpfr_init2 (zref, p);
mpfr_set_str (x, "1.9ecdc22ba77a5ab2560f7e84289e2a328906f47377ea3fd4c82d1bb2f13ee05c032cffc1933eadab7b0a5498e03e3bd0508968e59c25829d97a0b54f20cd4662c8dfffa54e714de41fc8ee3e0e0b244d110a194db05b70022b7d77f88955d415b09f17dd404576098dc51a583a3e49c35839551646e880c7eb790a01a4@1", 16, MPFR_RNDN);
inexref = mpfr_sin_cos (yref, zref, x, r);
inex = mpfr_sincos_fast (y, z, x, r);
if (mpfr_cmp (y, yref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091007)\n");
printf ("yref="); mpfr_dump (yref);
printf ("y="); mpfr_dump (y);
exit (1);
}
if (mpfr_cmp (z, zref))
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091007)\n");
printf ("zref="); mpfr_dump (zref);
printf ("z="); mpfr_dump (z);
exit (1);
}
if (inex != inexref)
{
printf ("mpfr_sin_cos and mpfr_sincos_fast disagree (bug20091007)\n");
printf ("inexref=%d inex=%d\n", inexref, inex);
exit (1);
}
mpfr_clear (x);
mpfr_clear (y);
mpfr_clear (z);
mpfr_clear (yref);
mpfr_clear (zref);
}
