Return to tset_q.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpfr / tests

Annotation of OpenXM_contrib/gmp/mpfr/tests/tset_q.c, Revision 1.1

1.1     ! ohara       1: /* Test file for mpfr_set_q.
        !             2:
        !             3: Copyright 2000, 2001, 2002 Free Software Foundation, Inc.
        !             4:
        !             5: This file is part of the MPFR Library.
        !             6:
        !             7: The MPFR Library is free software; you can redistribute it and/or modify
        !             8: it under the terms of the GNU Lesser General Public License as published by
        !             9: the Free Software Foundation; either version 2.1 of the License, or (at your
        !            10: option) any later version.
        !            11:
        !            12: The MPFR Library is distributed in the hope that it will be useful, but
        !            13: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
        !            14: or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU Lesser General Public
        !            15: License for more details.
        !            16:
        !            17: You should have received a copy of the GNU Lesser General Public License
        !            18: along with the MPFR Library; see the file COPYING.LIB.  If not, write to
        !            19: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
        !            20: MA 02111-1307, USA. */
        !            21:
        !            22: #include <stdio.h>
        !            23: #include <stdlib.h>
        !            24: #include <time.h>
        !            25: #include "gmp.h"
        !            26: #include "mpfr.h"
        !            27: #include "mpfr-test.h"
        !            28:
        !            29: void check _PROTO((long int, long int, mp_rnd_t, double));
        !            30:
        !            31: void
        !            32: check (long int n, long int d, mp_rnd_t rnd, double y)
        !            33: {
        !            34:   mpq_t q;
        !            35:   mpfr_t x, t;
        !            36:   double z;
        !            37:   int inexact, compare;
        !            38:
        !            39:   mpfr_init2 (x, 53);
        !            40:   mpfr_init2 (t, mpfr_get_prec (x) + mp_bits_per_limb);
        !            41:   mpq_init (q);
        !            42:   mpq_set_si (q, n, d);
        !            43: #ifdef MPFR_HAVE_FESETROUND
        !            44:   mpfr_set_machine_rnd_mode (rnd);
        !            45:   y = (double) n / d;
        !            46: #endif
        !            47:   inexact = mpfr_set_q (x, q, rnd);
        !            48:   z = mpfr_get_d1 (x);
        !            49:
        !            50:   /* check values */
        !            51:   if (y != z)
        !            52:     {
        !            53:     fprintf (stderr, "Error for q=%ld/%lu and rnd=%s\n", n, d,
        !            54:             mpfr_print_rnd_mode (rnd));
        !            55:     fprintf (stderr, "libm.a gives %1.20e, mpfr_set_q gives %1.20e\n", y, z);
        !            56:     exit (1);
        !            57:     }
        !            58:
        !            59:   /* check inexact flag */
        !            60:   if (mpfr_mul_ui (t, x, (d < 0) ? (-d) : d, rnd))
        !            61:     {
        !            62:       fprintf (stderr, "t <- x * d should be exact\n");
        !            63:       exit (1);
        !            64:     }
        !            65:   compare = mpfr_cmp_si (t, n);
        !            66:   if (((inexact == 0) && (compare != 0)) ||
        !            67:       ((inexact < 0) && (compare >= 0)) ||
        !            68:       ((inexact > 0) && (compare <= 0)))
        !            69:     {
        !            70:       fprintf (stderr, "wrong inexact flag: expected %d, got %d\n", compare,
        !            71:               inexact);
        !            72:       exit (1);
        !            73:     }
        !            74:
        !            75:   mpfr_clear (x);
        !            76:   mpfr_clear (t);
        !            77:   mpq_clear (q);
        !            78: }
        !            79:
        !            80: int
        !            81: main (void)
        !            82: {
        !            83: #ifdef MPFR_HAVE_FESETROUND
        !            84:   long int i, n;
        !            85:   unsigned long int d;
        !            86:   double y;
        !            87:   unsigned char rnd;
        !            88:
        !            89:   mpfr_test_init ();
        !            90:
        !            91:   SEED_RAND(time(NULL));
        !            92:   for (i=0;i<1000000;i++) {
        !            93:     n = LONG_RAND();
        !            94:     d = LONG_RAND();
        !            95:     if (LONG_RAND()%2) n = -n;
        !            96:     rnd = LONG_RAND() % 4;
        !            97:     y = (double) n / d;
        !            98:     check(n, d, rnd, y);
        !            99:   }
        !           100: #endif
        !           101:   check(-1647229822, 40619231, GMP_RNDZ, -4.055295438754120596e1);
        !           102:   check(-148939696, 1673285490, GMP_RNDZ, -8.9010331404953485501e-2);
        !           103:   check(-441322590, 273662545, GMP_RNDZ, -1.6126525096812205362);
        !           104:   check(-1631156895, 1677687197, GMP_RNDU, -9.722652100563177191e-1);
        !           105:   check(2141332571, 3117601, GMP_RNDZ, 6.8685267004982347316e2);
        !           106:   check(75504803, 400207282, GMP_RNDU, 1.8866424074712365155e-1);
        !           107:   check(643562308, 23100894, GMP_RNDD, 2.7858762002890447462e1);
        !           108:   check(632549085, 1831935802, GMP_RNDN, 3.4528998467600230393e-1);
        !           109:
        !           110:   return 0;
        !           111: }