| /* Test file for mpfr_sec. |
| |
| Copyright 2005, 2006, 2007, 2008, 2009 Free Software Foundation, Inc. |
| Contributed by the Arenaire and Cacao projects, INRIA. |
| |
| This file is part of the GNU MPFR Library. |
| |
| The GNU MPFR Library is free software; you can redistribute it and/or modify |
| it under the terms of the GNU Lesser General Public License as published by |
| the Free Software Foundation; either version 2.1 of the License, or (at your |
| option) any later version. |
| |
| The GNU MPFR Library is distributed in the hope that it will be useful, but |
| WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
| or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public |
| License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public License |
| along with the GNU MPFR Library; see the file COPYING.LIB. If not, write to |
| the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, |
| MA 02110-1301, USA. */ |
| |
| #include <stdio.h> |
| #include <stdlib.h> |
| |
| #include "mpfr-test.h" |
| |
| #define TEST_FUNCTION mpfr_sec |
| #define REDUCE_EMAX 262143 /* otherwise arg. reduction is too expensive */ |
| #include "tgeneric.c" |
| |
| 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, GMP_RNDN); |
| if (! mpfr_nan_p (y)) |
| { |
| printf ("Error: sec(NaN) != NaN\n"); |
| exit (1); |
| } |
| |
| mpfr_set_inf (x, 1); |
| mpfr_sec (y, x, GMP_RNDN); |
| if (! mpfr_nan_p (y)) |
| { |
| printf ("Error: sec(Inf) != NaN\n"); |
| exit (1); |
| } |
| |
| mpfr_set_inf (x, -1); |
| mpfr_sec (y, x, GMP_RNDN); |
| if (! mpfr_nan_p (y)) |
| { |
| printf ("Error: sec(-Inf) != NaN\n"); |
| exit (1); |
| } |
| |
| /* sec(+/-0) = 1 */ |
| mpfr_set_ui (x, 0, GMP_RNDN); |
| mpfr_sec (y, x, GMP_RNDN); |
| if (mpfr_cmp_ui (y, 1)) |
| { |
| printf ("Error: sec(+0) != 1\n"); |
| exit (1); |
| } |
| mpfr_neg (x, x, GMP_RNDN); |
| mpfr_sec (y, x, GMP_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; |
| mp_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, GMP_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), GMP_RNDN); |
| mpfr_clear_flags (); |
| inex = mpfr_sec (x, x, (mp_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 ((mp_rnd_t) rnd)); |
| err = 1; |
| } |
| if (rnd == GMP_RNDZ || rnd == GMP_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 ((mp_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 ((mp_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 ((mp_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 ((mp_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); |
| } |
| |
| int |
| main (int argc, char *argv[]) |
| { |
| tests_start_mpfr (); |
| |
| check_specials (); |
| |
| test_generic (2, 200, 10); |
| overflowed_sec0 (); |
| |
| tests_end_mpfr (); |
| return 0; |
| } |