static void check_specials (void) { mpfr_t x, y; mpfr_init2 (x, 123L); mpfr_init2 (y, 123L); mpfr_set_nan (x); mpfr_sec (y, x, MPFR_RNDN); if (! mpfr_nan_p (y)) { printf ("Error: sec(NaN) != NaN\n"); exit (1); } mpfr_set_inf (x, 1); mpfr_sec (y, x, MPFR_RNDN); if (! mpfr_nan_p (y)) { printf ("Error: sec(Inf) != NaN\n"); exit (1); } mpfr_set_inf (x, -1); mpfr_sec (y, x, MPFR_RNDN); if (! mpfr_nan_p (y)) { printf ("Error: sec(-Inf) != NaN\n"); exit (1); } /* sec(+/-0) = 1 */ mpfr_set_ui (x, 0, MPFR_RNDN); mpfr_sec (y, x, MPFR_RNDN); if (mpfr_cmp_ui (y, 1)) { printf ("Error: sec(+0) != 1\n"); exit (1); } mpfr_neg (x, x, MPFR_RNDN); mpfr_sec (y, x, MPFR_RNDN); if (mpfr_cmp_ui (y, 1)) { printf ("Error: sec(-0) != 1\n"); exit (1); } mpfr_clear (x); mpfr_clear (y); }
static void overflowed_sec0 (void) { mpfr_t x, y; int emax, i, inex, rnd, err = 0; mpfr_exp_t old_emax; old_emax = mpfr_get_emax (); mpfr_init2 (x, 8); mpfr_init2 (y, 8); for (emax = -1; emax <= 0; emax++) { mpfr_set_ui_2exp (y, 1, emax, MPFR_RNDN); mpfr_nextbelow (y); set_emax (emax); /* 1 is not representable. */ for (i = -1; i <= 1; i++) RND_LOOP (rnd) { mpfr_set_si_2exp (x, i, -512 * ABS (i), MPFR_RNDN); mpfr_clear_flags (); inex = mpfr_sec (x, x, (mpfr_rnd_t) rnd); if (! mpfr_overflow_p ()) { printf ("Error in overflowed_sec0 (i = %d, rnd = %s):\n" " The overflow flag is not set.\n", i, mpfr_print_rnd_mode ((mpfr_rnd_t) rnd)); err = 1; } if (rnd == MPFR_RNDZ || rnd == MPFR_RNDD) { if (inex >= 0) { printf ("Error in overflowed_sec0 (i = %d, rnd = %s):\n" " The inexact value must be negative.\n", i, mpfr_print_rnd_mode ((mpfr_rnd_t) rnd)); err = 1; } if (! mpfr_equal_p (x, y)) { printf ("Error in overflowed_sec0 (i = %d, rnd = %s):\n" " Got ", i, mpfr_print_rnd_mode ((mpfr_rnd_t) rnd)); mpfr_print_binary (x); printf (" instead of 0.11111111E%d.\n", emax); err = 1; } } else { if (inex <= 0) { printf ("Error in overflowed_sec0 (i = %d, rnd = %s):\n" " The inexact value must be positive.\n", i, mpfr_print_rnd_mode ((mpfr_rnd_t) rnd)); err = 1; } if (! (mpfr_inf_p (x) && MPFR_SIGN (x) > 0)) { printf ("Error in overflowed_sec0 (i = %d, rnd = %s):\n" " Got ", i, mpfr_print_rnd_mode ((mpfr_rnd_t) rnd)); mpfr_print_binary (x); printf (" instead of +Inf.\n"); err = 1; } } } set_emax (old_emax); } if (err) exit (1); mpfr_clear (x); mpfr_clear (y); }
MpfrFloat MpfrFloat::sec(const MpfrFloat& value) { MpfrFloat retval(MpfrFloat::kNoInitialization); mpfr_sec(retval.mData->mFloat, value.mData->mFloat, GMP_RNDN); return retval; }
void bvisit(const Sec &x) { apply(result_, *(x.get_arg())); mpfr_sec(result_, result_, rnd_); }