version 1.1.1.1, 2000/09/09 14:12:19 |
version 1.1.1.2, 2003/08/25 16:06:10 |
|
|
/* Test file for mpfr_mul_2exp. |
/* 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. |
This file is part of the MPFR Library. |
|
|
The MPFR Library is free software; you can redistribute it and/or modify |
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 |
it under the terms of the GNU Lesser General Public License as published by |
the Free Software Foundation; either version 2 of the License, or (at your |
the Free Software Foundation; either version 2.1 of the License, or (at your |
option) any later version. |
option) any later version. |
|
|
The MPFR Library is distributed in the hope that it will be useful, but |
The MPFR Library is distributed in the hope that it will be useful, but |
WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY |
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. |
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 |
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, |
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, |
MA 02111-1307, USA. */ |
MA 02111-1307, USA. */ |
Line 24 MA 02111-1307, USA. */ |
|
Line 24 MA 02111-1307, USA. */ |
|
#include <time.h> |
#include <time.h> |
#include "gmp.h" |
#include "gmp.h" |
#include "mpfr.h" |
#include "mpfr.h" |
|
#include "mpfr-impl.h" |
|
#include "mpfr-test.h" |
|
|
/* checks that x*y gives the same results in double |
/* checks that x*y gives the same results in double |
and with mpfr with 53 bits of precision */ |
and with mpfr with 53 bits of precision */ |
|
|
int |
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); |
mpfr_init2(w, 53); |
|
|
srand48(time(NULL)); |
mpfr_set_inf (w, 1); |
x = drand48(); |
mpfr_mul_2exp (w, w, 10, GMP_RNDZ); |
mpfr_set_d(w, x, 0); |
if (!MPFR_IS_INF(w)) { fprintf(stderr, "Inf != Inf"); exit(-1); } |
mpfr_mul_2exp(w, w, 10, GMP_RNDZ); |
|
if (x != (z = mpfr_get_d(w)/1024)) |
mpfr_set_nan (w); |
{ |
mpfr_mul_2exp (w, w, 10, GMP_RNDZ); |
fprintf(stderr, "%f != %f\n", x, z); |
if (!MPFR_IS_NAN(w)) { fprintf(stderr, "NaN != NaN"); exit(-1); } |
exit (1); |
|
} |
|
exit (0); |
|
} |
|
|
|
|
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; |
|
} |