=================================================================== RCS file: /home/cvs/OpenXM_contrib/gmp/mpfr/tests/Attic/tmul_2exp.c,v retrieving revision 1.1.1.1 retrieving revision 1.1.1.2 diff -u -p -r1.1.1.1 -r1.1.1.2 --- OpenXM_contrib/gmp/mpfr/tests/Attic/tmul_2exp.c 2000/09/09 14:12:19 1.1.1.1 +++ OpenXM_contrib/gmp/mpfr/tests/Attic/tmul_2exp.c 2003/08/25 16:06:10 1.1.1.2 @@ -1,20 +1,20 @@ /* Test file for mpfr_mul_2exp. -Copyright (C) 1999 PolKA project, Inria Lorraine and Loria +Copyright 1999, 2001, 2002 Free Software Foundation. This file is part of the MPFR Library. The MPFR Library is free software; you can redistribute it and/or modify -it under the terms of the GNU Library General Public License as published by -the Free Software Foundation; either version 2 of the License, or (at your +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 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 Library General Public +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 Library General Public License +You should have received a copy of the GNU Lesser General Public License along with the MPFR Library; see the file COPYING.LIB. If not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. */ @@ -24,26 +24,49 @@ MA 02111-1307, USA. */ #include #include "gmp.h" #include "mpfr.h" +#include "mpfr-impl.h" +#include "mpfr-test.h" /* checks that x*y gives the same results in double and with mpfr with 53 bits of precision */ int -main(argc,argv) int argc; char *argv[]; +main (int argc, char *argv[]) { - double x, z; mpfr_t w; + double x, z; mpfr_t w; unsigned long k; mpfr_init2(w, 53); - srand48(time(NULL)); - x = drand48(); - mpfr_set_d(w, x, 0); - mpfr_mul_2exp(w, w, 10, GMP_RNDZ); - if (x != (z = mpfr_get_d(w)/1024)) - { - fprintf(stderr, "%f != %f\n", x, z); - exit (1); - } - exit (0); -} + mpfr_set_inf (w, 1); + mpfr_mul_2exp (w, w, 10, GMP_RNDZ); + if (!MPFR_IS_INF(w)) { fprintf(stderr, "Inf != Inf"); exit(-1); } + + mpfr_set_nan (w); + mpfr_mul_2exp (w, w, 10, GMP_RNDZ); + if (!MPFR_IS_NAN(w)) { fprintf(stderr, "NaN != NaN"); exit(-1); } + SEED_RAND (time(NULL)); + for (k = 0; k < 100000; k++) { + x = DBL_RAND (); + mpfr_set_d (w, x, 0); + mpfr_mul_2exp (w, w, 10, GMP_RNDZ); + if (x != (z = mpfr_get_d1 (w)/1024)) + { + fprintf(stderr, "%f != %f\n", x, z); + return -1; + } + + mpfr_set_d(w, x, 0); + mpfr_div_2exp(w, w, 10, GMP_RNDZ); + if (x != (z = mpfr_get_d1 (w)*1024)) + { + fprintf(stderr, "%f != %f\n", x, z); + mpfr_clear(w); + return -1; + } + } + + mpfr_clear(w); + + return 0; +}