Annotation of OpenXM_contrib/gmp/mpfr/tests/tfma.c, Revision 1.1
1.1 ! ohara 1: /* Test file for mpfr_fma.
! 2:
! 3: Copyright 2001 Free Software Foundation.
! 4: Adapted from tarctan.c.
! 5:
! 6: This file is part of the MPFR Library.
! 7:
! 8: The MPFR Library is free software; you can redistribute it and/or modify
! 9: it under the terms of the GNU Lesser General Public License as published by
! 10: the Free Software Foundation; either version 2.1 of the License, or (at your
! 11: option) any later version.
! 12:
! 13: The MPFR Library is distributed in the hope that it will be useful, but
! 14: WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
! 15: or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
! 16: License for more details.
! 17:
! 18: You should have received a copy of the GNU Lesser General Public License
! 19: along with the MPFR Library; see the file COPYING.LIB. If not, write to
! 20: the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
! 21: MA 02111-1307, USA. */
! 22:
! 23: #include <stdio.h>
! 24: #include <limits.h>
! 25: #include <stdlib.h>
! 26: #include "gmp.h"
! 27: #include "gmp-impl.h"
! 28: #include "mpfr.h"
! 29: #include "mpfr-impl.h"
! 30: #include "mpfr-test.h"
! 31:
! 32:
! 33: int
! 34: main (int argc, char *argv[])
! 35: {
! 36: mpfr_t x, y,z,s;
! 37: mpfr_init (x);
! 38: mpfr_init (s);
! 39: mpfr_init (y);
! 40: mpfr_init (z);
! 41:
! 42: /* check special cases */
! 43: mpfr_set_prec (x, 2);
! 44: mpfr_set_prec (y, 2);
! 45: mpfr_set_prec (z, 2);
! 46: mpfr_set_prec (s, 2);
! 47: mpfr_set_d (x, -0.75, GMP_RNDN);
! 48: mpfr_set_d (y, 0.5, GMP_RNDN);
! 49: mpfr_set_d (z, 0.375, GMP_RNDN);
! 50: mpfr_fma (s, x, y, z, GMP_RNDU); /* result is 0 */
! 51:
! 52: mpfr_set_prec (x, 27);
! 53: mpfr_set_prec (y, 27);
! 54: mpfr_set_prec (z, 27);
! 55: mpfr_set_prec (s, 27);
! 56: mpfr_set_str_raw (x, "1.11111111111111111111111111e-1");
! 57: mpfr_set (y, x, GMP_RNDN);
! 58: mpfr_set_str_raw (z, "-1.00011110100011001011001001e-1");
! 59: if (mpfr_fma (s, x, y, z, GMP_RNDN) >= 0)
! 60: {
! 61: fprintf (stderr, "Wrong inexact flag for x=y=1-2^(-27)\n");
! 62: exit (1);
! 63: }
! 64:
! 65: MPFR_SET_NAN(x);
! 66: mpfr_random(y);
! 67: mpfr_random(z);
! 68: mpfr_fma (s,x, y,z, GMP_RNDN);
! 69: if(!MPFR_IS_NAN(s))
! 70: {
! 71: fprintf (stderr, "evaluation of function in x=NAN does not return NAN");
! 72: exit (1);
! 73: }
! 74:
! 75: MPFR_CLEAR_FLAGS(x);
! 76: MPFR_CLEAR_FLAGS(y);
! 77: MPFR_CLEAR_FLAGS(z);
! 78:
! 79: MPFR_SET_NAN(y);
! 80: mpfr_random(x);
! 81: mpfr_random(z);
! 82: mpfr_fma (s,x, y,z, GMP_RNDN);
! 83: if(!MPFR_IS_NAN(s))
! 84: {
! 85: printf ("evaluation of function in y=NAN does not return NAN");
! 86: exit (1);
! 87: }
! 88:
! 89: MPFR_CLEAR_FLAGS(x);
! 90: MPFR_CLEAR_FLAGS(y);
! 91: MPFR_CLEAR_FLAGS(z);
! 92:
! 93: MPFR_SET_NAN(z);
! 94: mpfr_random(y);
! 95: mpfr_random(x);
! 96: mpfr_fma (s,x, y,z, GMP_RNDN);
! 97: if(!MPFR_IS_NAN(s))
! 98: {
! 99: printf ("evaluation of function in z=NAN does not return NAN");
! 100: exit (1);
! 101: }
! 102:
! 103: MPFR_CLEAR_FLAGS(x);
! 104: MPFR_CLEAR_FLAGS(y);
! 105: MPFR_CLEAR_FLAGS(z);
! 106:
! 107: MPFR_SET_INF(x);
! 108: MPFR_SET_ZERO(y);
! 109: mpfr_random(z);
! 110: mpfr_fma (s,x, y,z, GMP_RNDN);
! 111: if(!MPFR_IS_NAN(s))
! 112: {
! 113: printf ("evaluation of function in x=INF y=0 does not return NAN");
! 114: exit (1);
! 115: }
! 116:
! 117: MPFR_CLEAR_FLAGS(z);
! 118: MPFR_CLEAR_FLAGS(y);
! 119: MPFR_CLEAR_FLAGS(x);
! 120:
! 121: MPFR_SET_INF(y);
! 122: MPFR_SET_ZERO(x);
! 123: mpfr_random(z);
! 124: mpfr_fma (s,x, y,z, GMP_RNDN);
! 125: if(!MPFR_IS_NAN(s))
! 126: {
! 127: printf ("evaluation of function in x=0 y=INF does not return NAN");
! 128: exit (1);
! 129: }
! 130:
! 131: MPFR_CLEAR_FLAGS(z);
! 132: MPFR_CLEAR_FLAGS(y);
! 133: MPFR_CLEAR_FLAGS(x);
! 134:
! 135: MPFR_SET_INF(x);
! 136: mpfr_random(y);
! 137: MPFR_SET_INF(z);
! 138: if((MPFR_SIGN(x) * MPFR_SIGN(y)) == MPFR_SIGN(z))
! 139: MPFR_CHANGE_SIGN(z);
! 140: mpfr_fma (s,x, y,z, GMP_RNDN);
! 141: if(!MPFR_IS_NAN(s))
! 142: {
! 143: printf ("evaluation of function in x=INF z=(-sign(x)*sign(y))INF does not return NAN");
! 144: exit (1);
! 145: }
! 146:
! 147: MPFR_CLEAR_FLAGS(z);
! 148: MPFR_CLEAR_FLAGS(y);
! 149: MPFR_CLEAR_FLAGS(x);
! 150:
! 151: MPFR_SET_INF(y);
! 152: mpfr_random(x);
! 153: MPFR_SET_INF(z);
! 154: if((MPFR_SIGN(x) * MPFR_SIGN(y)) == MPFR_SIGN(z))
! 155: MPFR_CHANGE_SIGN(z);
! 156:
! 157: mpfr_fma (s,x, y,z, GMP_RNDN);
! 158: if(!MPFR_IS_NAN(s))
! 159: {
! 160: printf ("evaluation of function in y=INF z=(-sign(x)*sign(y))INF does not return NAN");
! 161: exit (1);
! 162: }
! 163:
! 164: MPFR_CLEAR_FLAGS(z);
! 165: MPFR_CLEAR_FLAGS(y);
! 166: MPFR_CLEAR_FLAGS(x);
! 167:
! 168: MPFR_SET_INF(x);
! 169: mpfr_random(y);
! 170: mpfr_random(z);
! 171: mpfr_fma (s,x, y,z, GMP_RNDN);
! 172: if(!MPFR_IS_INF(s))
! 173: {
! 174: printf ("evaluation of function in x=INF does not return INF");
! 175: exit (1);
! 176: }
! 177:
! 178: MPFR_CLEAR_FLAGS(z);
! 179: MPFR_CLEAR_FLAGS(y);
! 180: MPFR_CLEAR_FLAGS(x);
! 181:
! 182: MPFR_SET_INF(y);
! 183: mpfr_random(x);
! 184: mpfr_random(z);
! 185: mpfr_fma (s,x, y,z, GMP_RNDN);
! 186: if(!MPFR_IS_INF(s))
! 187: {
! 188: printf ("evaluation of function in y=INF does not return INF");
! 189: exit (1);
! 190: }
! 191:
! 192: MPFR_CLEAR_FLAGS(z);
! 193: MPFR_CLEAR_FLAGS(y);
! 194: MPFR_CLEAR_FLAGS(x);
! 195:
! 196: MPFR_SET_INF(z);
! 197: mpfr_random(x);
! 198: mpfr_random(y);
! 199: mpfr_fma (s,x, y,z, GMP_RNDN);
! 200: if(!MPFR_IS_INF(s))
! 201: {
! 202: printf ("evaluation of function in z=INF does not return INF");
! 203: exit (1);
! 204: }
! 205:
! 206: MPFR_CLEAR_FLAGS(z);
! 207: MPFR_CLEAR_FLAGS(y);
! 208: MPFR_CLEAR_FLAGS(x);
! 209:
! 210: MPFR_SET_ZERO(x);
! 211: mpfr_random(y);
! 212: mpfr_random(z);
! 213: mpfr_fma (s,x, y,z, GMP_RNDN);
! 214: if(mpfr_cmp(s,z)!=0)
! 215: {
! 216: printf ("evaluation of function in x=0 does not return z\n");
! 217: exit (1);
! 218: }
! 219:
! 220: MPFR_CLEAR_FLAGS(z);
! 221: MPFR_CLEAR_FLAGS(y);
! 222: MPFR_CLEAR_FLAGS(x);
! 223:
! 224: MPFR_SET_ZERO(y);
! 225: mpfr_random(x);
! 226: mpfr_random(z);
! 227: mpfr_fma (s,x, y,z, GMP_RNDN);
! 228: if(mpfr_cmp(s,z)!=0)
! 229: {
! 230: printf ("evaluation of function in y=0 does not return z\n");
! 231: exit (1);
! 232: }
! 233:
! 234: MPFR_CLEAR_FLAGS(z);
! 235: MPFR_CLEAR_FLAGS(y);
! 236: MPFR_CLEAR_FLAGS(x);
! 237:
! 238: {
! 239: mp_prec_t prec;
! 240: mpfr_t t, slong;
! 241: mp_rnd_t rnd;
! 242: int inexact, compare;
! 243: unsigned int n;
! 244:
! 245: int p0=2;
! 246: int p1=200;
! 247: int N=200;
! 248:
! 249: mpfr_init (t);
! 250: mpfr_init (slong);
! 251:
! 252: /* generic test */
! 253: for (prec = p0; prec <= p1; prec++)
! 254: {
! 255: mpfr_set_prec (x, prec);
! 256: mpfr_set_prec (y, prec);
! 257: mpfr_set_prec (z, prec);
! 258: mpfr_set_prec (s, prec);
! 259: mpfr_set_prec (t, prec);
! 260:
! 261: for (n=0; n<N; n++)
! 262: {
! 263: mpfr_random (x);
! 264: mpfr_random (y);
! 265: mpfr_random (z);
! 266:
! 267: if (random() % 2)
! 268: mpfr_neg (x, x, GMP_RNDN);
! 269: if (random() % 2)
! 270: mpfr_neg (y, y, GMP_RNDN);
! 271: if (random() % 2)
! 272: mpfr_neg (z, z, GMP_RNDN);
! 273:
! 274: rnd = random () % 4;
! 275: mpfr_set_prec (slong, 2 * prec);
! 276: if (mpfr_mul (slong, x, y, rnd))
! 277: {
! 278: fprintf (stderr, "x*y should be exact\n");
! 279: exit (1);
! 280: }
! 281: compare = mpfr_add (t, slong, z, rnd);
! 282: inexact = mpfr_fma (s, x, y, z, rnd);
! 283: if (mpfr_cmp (s, t))
! 284: {
! 285: printf ("results differ for x=");
! 286: mpfr_out_str (stdout, 2, prec, x, GMP_RNDN);
! 287: printf (" y=");
! 288: mpfr_out_str (stdout, 2, prec, y, GMP_RNDN);
! 289: printf (" z=");
! 290: mpfr_out_str (stdout, 2, prec, z, GMP_RNDN);
! 291: printf (" prec=%u rnd_mode=%s\n", (unsigned) prec,
! 292: mpfr_print_rnd_mode (rnd));
! 293: printf ("got ");
! 294: mpfr_out_str (stdout, 2, prec, s, GMP_RNDN);
! 295: putchar ('\n');
! 296: printf ("expected ");
! 297: mpfr_out_str (stdout, 2, prec, t, GMP_RNDN);
! 298: putchar ('\n');
! 299: printf ("approx ");
! 300: mpfr_print_binary (slong);
! 301: putchar ('\n');
! 302: exit (1);
! 303: }
! 304: if (((inexact == 0) && (compare != 0)) ||
! 305: ((inexact < 0) && (compare >= 0)) ||
! 306: ((inexact > 0) && (compare <= 0)))
! 307: {
! 308: fprintf (stderr, "Wrong inexact flag for rnd=%s: expected %d, got %d\n",
! 309: mpfr_print_rnd_mode (rnd), compare, inexact);
! 310: fprintf (stderr, "x="); mpfr_out_str (stderr, 2, 0, x, GMP_RNDN);
! 311: fprintf (stderr, " y="); mpfr_out_str (stderr, 2, 0, y, GMP_RNDN);
! 312: fprintf (stderr, " z="); mpfr_out_str (stderr, 2, 0, z, GMP_RNDN);
! 313: fprintf (stderr, " s="); mpfr_out_str (stderr, 2, 0, s, GMP_RNDN);
! 314: fprintf (stderr, "\n");
! 315: fprintf (stderr, "z=%1.20e s=%1.20e\n", mpfr_get_d1 (z),
! 316: mpfr_get_d1 (s));
! 317: exit (1);
! 318: }
! 319: }
! 320: }
! 321: mpfr_clear (t);
! 322: mpfr_clear (slong);
! 323:
! 324: }
! 325: mpfr_clear (x);
! 326: mpfr_clear (y);
! 327: mpfr_clear (z);
! 328: mpfr_clear (s);
! 329:
! 330: return 0;
! 331: }
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to tfma.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib / gmp / mpfr / tests