/* Test file for mpfr_set_z. Copyright 1999, 2001, 2002 Free Software Foundation, Inc. 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 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 Lesser General Public License for more details. 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. */ #include #include #include #include "gmp.h" #include "mpfr.h" #include "mpfr-test.h" void check _PROTO((long, unsigned char)); void check_large _PROTO((void)); void check(long i, unsigned char rnd) { mpfr_t f; mpz_t z; mpfr_init2(f, 53); mpz_init(z); mpz_set_ui(z, i); mpfr_set_z(f, z, rnd); if ((long)mpfr_get_d1 (f) != i) { printf("Error in mpfr_set_z for i=%ld rnd_mode=%d\n",i,rnd); exit(1); } mpfr_clear(f); mpz_clear(z); } void check_large (void) { mpz_t z; mpfr_t x,y; mpz_init(z); mpfr_init2(x, 160); mpfr_init2(y, 160); mpz_set_str(z, "77031627725494291259359895954016675357279104942148788042", 10); mpfr_set_z(x, z, GMP_RNDN); mpfr_set_str_raw(y, "0.1100100100001111110110101010001000100001011010001100001000110100110001001100011001100010100010111000000011011100000111001101000100101001000000100100111000001001E186"); if (mpfr_cmp(x,y)) { fprintf(stderr, "Error in mpfr_set_z on large input\n"); exit(1); } mpz_clear(z); mpfr_clear(x); mpfr_clear(y); } /* tset_z z rnd prec */ int main (int argc, char *argv[]) { long j; check_large(); SEED_RAND (time(NULL)); check(0, 0); for (j=0; j<1000000; j++) check(LONG_RAND(), LONG_RAND()%4); return 0; }