math/test-misc.c

   1 /* Miscellaneous tests which don't fit anywhere else.
   2    Copyright (C) 2000, 2001 Free Software Foundation, Inc.
   3    This file is part of the GNU C Library.
   4
   5    The GNU C Library is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU Library General Public License as
   7    published by the Free Software Foundation; either version 2 of the
   8    License, or (at your option) any later version.
   9
  10    The GNU C Library is distributed in the hope that it will be useful,
  11    but WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13    Library General Public License for more details.
  14
  15    You should have received a copy of the GNU Library General Public
  16    License along with the GNU C Library; see the file COPYING.LIB.  If not,
  17    write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  18    Boston, MA 02111-1307, USA.  */
  19
  20 #include <fenv.h>
  21 #include <float.h>
  22 #include <ieee754.h>
  23 #include <math.h>
  24 #include <stdio.h>
  25 #include <string.h>
  26
  27
  28 int
  29 main (void)
  30 {
  31   int result = 0;
  32
  33 #ifndef NO_LONG_DOUBLE
  34   {
  35     long double x = 0x100000001ll + (long double) 0.5;
  36     long double q;
  37     long double r;
  38
  39     r = modfl (x, &q);
  40     if (q != (long double) 0x100000001ll || r != 0.5)
  41       {
  42         printf ("modfl (%Lg, ...) failed\n", x);
  43         result = 1;
  44       }
  45   }
  46
  47 # if __GNUC__ >= 3 || __GNUC_MINOR__ >= 96
  48   {
  49     long double x;
  50     long double m;
  51     long double r;
  52     int e;
  53     int i;
  54
  55 #  if LDBL_MANT_DIG == 64
  56     m = 0xf.fffffffffffffffp-4L;
  57 #  elif LDBL_MANT_DIG == 113
  58     m = 0x1.ffffffffffffffffffffffffffffp-1L;
  59 #  else
  60 #   error "Please adjust"
  61 #  endif
  62
  63     for (i = LDBL_MAX_EXP, x = LDBL_MAX; i >= LDBL_MIN_EXP; --i, x /= 2.0L)
  64       {
  65         printf ("2^%d: ", i);
  66
  67         r = frexpl (x, &e);
  68         if (r != m)
  69           {
  70             printf ("mantissa incorrect: %.20La\n", r);
  71             result = 1;
  72             continue;
  73           }
  74         if (e != i)
  75           {
  76             printf ("exponent wrong %d (%.20Lg)\n", e, x);
  77             result = 1;
  78             continue;
  79           }
  80         puts ("ok");
  81       }
  82   }
  83 # endif
  84
  85 #if 0
  86   {
  87     int e;
  88     long double r = frexpl (LDBL_MIN * LDBL_EPSILON, &e);
  89
  90     if (r != 0.5)
  91       {
  92         printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): mantissa wrong: %Lg\n",
  93                 r);
  94         result = 1;
  95       }
  96     else if (e != -16444)
  97       {
  98         printf ("frexpl (LDBL_MIN * LDBL_EPSILON, ...): exponent wrong: %d\n",
  99                 e);
 100         result = 1;
 101       }
 102   }
 103 #endif
 104 #endif
 105
 106   {
 107     double x = 0x100000001ll + (double) 0.5;
 108     double q;
 109     double r;
 110
 111     r = modf (x, &q);
 112     if (q != (double) 0x100000001ll || r != 0.5)
 113       {
 114         printf ("modf (%g, ...) failed\n", x);
 115         result = 1;
 116       }
 117   }
 118
 119   {
 120     union ieee754_float v1;
 121     union ieee754_float v2;
 122     float f;
 123
 124     v1.f = f = FLT_MIN;
 125     if (fpclassify (f) != FP_NORMAL)
 126       {
 127         printf ("fpclassify (FLT_MIN) failed: %d\n", fpclassify (f));
 128         result = 1;
 129       }
 130     f = nextafterf (f, FLT_MIN / 2.0f);
 131     if (fpclassify (f) != FP_SUBNORMAL)
 132       {
 133         printf ("fpclassify (FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
 134         result = 1;
 135       }
 136     v2.f = f = nextafterf (f, FLT_MIN);
 137     if (fpclassify (f) != FP_NORMAL)
 138       {
 139         printf ("fpclassify (FLT_MIN-epsilon+epsilon) failed: %d\n",
 140                 fpclassify (f));
 141         result = 1;
 142       }
 143
 144     if (v1.ieee.mantissa != v2.ieee.mantissa)
 145       {
 146         printf ("FLT_MIN: mantissa differs: %8x vs %8x\n",
 147                 v1.ieee.mantissa, v2.ieee.mantissa);
 148         result = 1;
 149       }
 150     if (v1.ieee.exponent != v2.ieee.exponent)
 151       {
 152         printf ("FLT_MIN: exponent differs: %4x vs %4x\n",
 153                 v1.ieee.exponent, v2.ieee.exponent);
 154         result = 1;
 155       }
 156     if (v1.ieee.negative != v2.ieee.negative)
 157       {
 158         printf ("FLT_MIN: negative differs: %d vs %d\n",
 159                 v1.ieee.negative, v2.ieee.negative);
 160         result = 1;
 161       }
 162
 163     v1.f = f = -FLT_MIN;
 164     if (fpclassify (f) != FP_NORMAL)
 165       {
 166         printf ("fpclassify (-FLT_MIN) failed: %d\n", fpclassify (f));
 167         result = 1;
 168       }
 169     f = nextafterf (f, -FLT_MIN / 2.0f);
 170     if (fpclassify (f) != FP_SUBNORMAL)
 171       {
 172         printf ("fpclassify (-FLT_MIN-epsilon) failed: %d\n", fpclassify (f));
 173         result = 1;
 174       }
 175     v2.f = f = nextafterf (f, -FLT_MIN);
 176     if (fpclassify (f) != FP_NORMAL)
 177       {
 178         printf ("fpclassify (-FLT_MIN-epsilon+epsilon) failed: %d\n",
 179                 fpclassify (f));
 180         result = 1;
 181       }
 182
 183     if (v1.ieee.mantissa != v2.ieee.mantissa)
 184       {
 185         printf ("-FLT_MIN: mantissa differs: %8x vs %8x\n",
 186                 v1.ieee.mantissa, v2.ieee.mantissa);
 187         result = 1;
 188       }
 189     if (v1.ieee.exponent != v2.ieee.exponent)
 190       {
 191         printf ("-FLT_MIN: exponent differs: %4x vs %4x\n",
 192                 v1.ieee.exponent, v2.ieee.exponent);
 193         result = 1;
 194       }
 195     if (v1.ieee.negative != v2.ieee.negative)
 196       {
 197         printf ("-FLT_MIN: negative differs: %d vs %d\n",
 198                 v1.ieee.negative, v2.ieee.negative);
 199         result = 1;
 200       }
 201
 202     f = FLT_MAX;
 203     if (fpclassify (f) != FP_NORMAL)
 204       {
 205         printf ("fpclassify (FLT_MAX) failed: %d\n", fpclassify (f));
 206         result = 1;
 207       }
 208     f = nextafterf (f, INFINITY);
 209     if (fpclassify (f) != FP_INFINITE)
 210       {
 211         printf ("fpclassify (FLT_MAX+epsilon) failed: %d\n", fpclassify (f));
 212         result = 1;
 213       }
 214
 215     f = -FLT_MAX;
 216     if (fpclassify (f) != FP_NORMAL)
 217       {
 218         printf ("fpclassify (-FLT_MAX) failed: %d\n", fpclassify (f));
 219         result = 1;
 220       }
 221     f = nextafterf (f, -INFINITY);
 222     if (fpclassify (f) != FP_INFINITE)
 223       {
 224         printf ("fpclassify (-FLT_MAX-epsilon) failed: %d\n", fpclassify (f));
 225         result = 1;
 226       }
 227
 228     v1.f = f = 0.0625;
 229     f = nextafterf (f, 0.0);
 230     v2.f = f = nextafterf (f, 1.0);
 231
 232     if (v1.ieee.mantissa != v2.ieee.mantissa)
 233       {
 234         printf ("0.0625f down: mantissa differs: %8x vs %8x\n",
 235                 v1.ieee.mantissa, v2.ieee.mantissa);
 236         result = 1;
 237       }
 238     if (v1.ieee.exponent != v2.ieee.exponent)
 239       {
 240         printf ("0.0625f down: exponent differs: %4x vs %4x\n",
 241                 v1.ieee.exponent, v2.ieee.exponent);
 242         result = 1;
 243       }
 244     if (v1.ieee.negative != v2.ieee.negative)
 245       {
 246         printf ("0.0625f down: negative differs: %d vs %d\n",
 247                 v1.ieee.negative, v2.ieee.negative);
 248         result = 1;
 249       }
 250
 251     v1.f = f = 0.0625;
 252     f = nextafterf (f, 1.0);
 253     v2.f = f = nextafterf (f, 0.0);
 254
 255     if (v1.ieee.mantissa != v2.ieee.mantissa)
 256       {
 257         printf ("0.0625f up: mantissa differs: %8x vs %8x\n",
 258                 v1.ieee.mantissa, v2.ieee.mantissa);
 259         result = 1;
 260       }
 261     if (v1.ieee.exponent != v2.ieee.exponent)
 262       {
 263         printf ("0.0625f up: exponent differs: %4x vs %4x\n",
 264                 v1.ieee.exponent, v2.ieee.exponent);
 265         result = 1;
 266       }
 267     if (v1.ieee.negative != v2.ieee.negative)
 268       {
 269         printf ("0.0625f up: negative differs: %d vs %d\n",
 270                 v1.ieee.negative, v2.ieee.negative);
 271         result = 1;
 272       }
 273
 274     v1.f = f = -0.0625;
 275     f = nextafterf (f, 0.0);
 276     v2.f = f = nextafterf (f, -1.0);
 277
 278     if (v1.ieee.mantissa != v2.ieee.mantissa)
 279       {
 280         printf ("-0.0625f up: mantissa differs: %8x vs %8x\n",
 281                 v1.ieee.mantissa, v2.ieee.mantissa);
 282         result = 1;
 283       }
 284     if (v1.ieee.exponent != v2.ieee.exponent)
 285       {
 286         printf ("-0.0625f up: exponent differs: %4x vs %4x\n",
 287                 v1.ieee.exponent, v2.ieee.exponent);
 288         result = 1;
 289       }
 290     if (v1.ieee.negative != v2.ieee.negative)
 291       {
 292         printf ("-0.0625f up: negative differs: %d vs %d\n",
 293                 v1.ieee.negative, v2.ieee.negative);
 294         result = 1;
 295       }
 296
 297     v1.f = f = -0.0625;
 298     f = nextafterf (f, -1.0);
 299     v2.f = f = nextafterf (f, 0.0);
 300
 301     if (v1.ieee.mantissa != v2.ieee.mantissa)
 302       {
 303         printf ("-0.0625f down: mantissa differs: %8x vs %8x\n",
 304                 v1.ieee.mantissa, v2.ieee.mantissa);
 305         result = 1;
 306       }
 307     if (v1.ieee.exponent != v2.ieee.exponent)
 308       {
 309         printf ("-0.0625f down: exponent differs: %4x vs %4x\n",
 310                 v1.ieee.exponent, v2.ieee.exponent);
 311         result = 1;
 312       }
 313     if (v1.ieee.negative != v2.ieee.negative)
 314       {
 315         printf ("-0.0625f down: negative differs: %d vs %d\n",
 316                 v1.ieee.negative, v2.ieee.negative);
 317         result = 1;
 318       }
 319
 320     v1.f = f = 0.0f;
 321     f = nextafterf (f, 1.0);
 322     v2.f = nextafterf (f, -1.0);
 323
 324     if (v1.ieee.mantissa != v2.ieee.mantissa)
 325       {
 326         printf ("0.0f up: mantissa differs: %8x vs %8x\n",
 327                 v1.ieee.mantissa, v2.ieee.mantissa);
 328         result = 1;
 329       }
 330     if (v1.ieee.exponent != v2.ieee.exponent)
 331       {
 332         printf ("0.0f up: exponent differs: %4x vs %4x\n",
 333                 v1.ieee.exponent, v2.ieee.exponent);
 334         result = 1;
 335       }
 336     if (0 != v2.ieee.negative)
 337       {
 338         printf ("0.0f up: negative differs: 0 vs %d\n",
 339                 v2.ieee.negative);
 340         result = 1;
 341       }
 342
 343     v1.f = f = 0.0f;
 344     f = nextafterf (f, -1.0);
 345     v2.f = nextafterf (f, 1.0);
 346
 347     if (v1.ieee.mantissa != v2.ieee.mantissa)
 348       {
 349         printf ("0.0f down: mantissa differs: %8x vs %8x\n",
 350                 v1.ieee.mantissa, v2.ieee.mantissa);
 351         result = 1;
 352       }
 353     if (v1.ieee.exponent != v2.ieee.exponent)
 354       {
 355         printf ("0.0f down: exponent differs: %4x vs %4x\n",
 356                 v1.ieee.exponent, v2.ieee.exponent);
 357         result = 1;
 358       }
 359     if (1 != v2.ieee.negative)
 360       {
 361         printf ("0.0f down: negative differs: 1 vs %d\n",
 362                 v2.ieee.negative);
 363         result = 1;
 364       }
 365   }
 366
 367   {
 368     union ieee754_double v1;
 369     union ieee754_double v2;
 370     double d;
 371
 372     v1.d = d = DBL_MIN;
 373     if (fpclassify (d) != FP_NORMAL)
 374       {
 375         printf ("fpclassify (DBL_MIN) failed: %d\n", fpclassify (d));
 376         result = 1;
 377       }
 378     d = nextafter (d, DBL_MIN / 2.0);
 379     if (fpclassify (d) != FP_SUBNORMAL)
 380       {
 381         printf ("fpclassify (DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
 382         result = 1;
 383       }
 384     v2.d = d = nextafter (d, DBL_MIN);
 385     if (fpclassify (d) != FP_NORMAL)
 386       {
 387         printf ("fpclassify (DBL_MIN-epsilon+epsilon) failed: %d\n",
 388                 fpclassify (d));
 389         result = 1;
 390       }
 391
 392     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 393       {
 394         printf ("DBL_MIN: mantissa0 differs: %8x vs %8x\n",
 395                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 396         result = 1;
 397       }
 398     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 399       {
 400         printf ("DBL_MIN: mantissa1 differs: %8x vs %8x\n",
 401                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 402         result = 1;
 403       }
 404     if (v1.ieee.exponent != v2.ieee.exponent)
 405       {
 406         printf ("DBL_MIN: exponent differs: %4x vs %4x\n",
 407                 v1.ieee.exponent, v2.ieee.exponent);
 408         result = 1;
 409       }
 410     if (v1.ieee.negative != v2.ieee.negative)
 411       {
 412         printf ("DBL_MIN: negative differs: %d vs %d\n",
 413                 v1.ieee.negative, v2.ieee.negative);
 414         result = 1;
 415       }
 416
 417     v1.d = d = -DBL_MIN;
 418     if (fpclassify (d) != FP_NORMAL)
 419       {
 420         printf ("fpclassify (-DBL_MIN) failed: %d\n", fpclassify (d));
 421         result = 1;
 422       }
 423     d = nextafter (d, -DBL_MIN / 2.0);
 424     if (fpclassify (d) != FP_SUBNORMAL)
 425       {
 426         printf ("fpclassify (-DBL_MIN-epsilon) failed: %d\n", fpclassify (d));
 427         result = 1;
 428       }
 429     v2.d = d = nextafter (d, -DBL_MIN);
 430     if (fpclassify (d) != FP_NORMAL)
 431       {
 432         printf ("fpclassify (-DBL_MIN-epsilon+epsilon) failed: %d\n",
 433                 fpclassify (d));
 434         result = 1;
 435       }
 436
 437     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 438       {
 439         printf ("-DBL_MIN: mantissa0 differs: %8x vs %8x\n",
 440                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 441         result = 1;
 442       }
 443     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 444       {
 445         printf ("-DBL_MIN: mantissa1 differs: %8x vs %8x\n",
 446                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 447         result = 1;
 448       }
 449     if (v1.ieee.exponent != v2.ieee.exponent)
 450       {
 451         printf ("-DBL_MIN: exponent differs: %4x vs %4x\n",
 452                 v1.ieee.exponent, v2.ieee.exponent);
 453         result = 1;
 454       }
 455     if (v1.ieee.negative != v2.ieee.negative)
 456       {
 457         printf ("-DBL_MIN: negative differs: %d vs %d\n",
 458                 v1.ieee.negative, v2.ieee.negative);
 459         result = 1;
 460       }
 461
 462     d = DBL_MAX;
 463     if (fpclassify (d) != FP_NORMAL)
 464       {
 465         printf ("fpclassify (DBL_MAX) failed: %d\n", fpclassify (d));
 466         result = 1;
 467       }
 468     d = nextafter (d, INFINITY);
 469     if (fpclassify (d) != FP_INFINITE)
 470       {
 471         printf ("fpclassify (DBL_MAX+epsilon) failed: %d\n", fpclassify (d));
 472         result = 1;
 473       }
 474
 475     d = -DBL_MAX;
 476     if (fpclassify (d) != FP_NORMAL)
 477       {
 478         printf ("fpclassify (-DBL_MAX) failed: %d\n", fpclassify (d));
 479         result = 1;
 480       }
 481     d = nextafter (d, -INFINITY);
 482     if (fpclassify (d) != FP_INFINITE)
 483       {
 484         printf ("fpclassify (-DBL_MAX-epsilon) failed: %d\n", fpclassify (d));
 485         result = 1;
 486       }
 487
 488     v1.d = d = 0.0625;
 489     d = nextafter (d, 0.0);
 490     v2.d = d = nextafter (d, 1.0);
 491
 492     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 493       {
 494         printf ("0.0625 down: mantissa0 differs: %8x vs %8x\n",
 495                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 496         result = 1;
 497       }
 498     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 499       {
 500         printf ("0.0625 down: mantissa1 differs: %8x vs %8x\n",
 501                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 502         result = 1;
 503       }
 504     if (v1.ieee.exponent != v2.ieee.exponent)
 505       {
 506         printf ("0.0625 down: exponent differs: %4x vs %4x\n",
 507                 v1.ieee.exponent, v2.ieee.exponent);
 508         result = 1;
 509       }
 510     if (v1.ieee.negative != v2.ieee.negative)
 511       {
 512         printf ("0.0625 down: negative differs: %d vs %d\n",
 513                 v1.ieee.negative, v2.ieee.negative);
 514         result = 1;
 515       }
 516
 517     v1.d = d = 0.0625;
 518     d = nextafter (d, 1.0);
 519     v2.d = d = nextafter (d, 0.0);
 520
 521     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 522       {
 523         printf ("0.0625 up: mantissa0 differs: %8x vs %8x\n",
 524                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 525         result = 1;
 526       }
 527     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 528       {
 529         printf ("0.0625 up: mantissa1 differs: %8x vs %8x\n",
 530                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 531         result = 1;
 532       }
 533     if (v1.ieee.exponent != v2.ieee.exponent)
 534       {
 535         printf ("0.0625 up: exponent differs: %4x vs %4x\n",
 536                 v1.ieee.exponent, v2.ieee.exponent);
 537         result = 1;
 538       }
 539     if (v1.ieee.negative != v2.ieee.negative)
 540       {
 541         printf ("0.0625 up: negative differs: %d vs %d\n",
 542                 v1.ieee.negative, v2.ieee.negative);
 543         result = 1;
 544       }
 545
 546     v1.d = d = -0.0625;
 547     d = nextafter (d, 0.0);
 548     v2.d = d = nextafter (d, -1.0);
 549
 550     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 551       {
 552         printf ("-0.0625 up: mantissa0 differs: %8x vs %8x\n",
 553                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 554         result = 1;
 555       }
 556     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 557       {
 558         printf ("-0.0625 up: mantissa1 differs: %8x vs %8x\n",
 559                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 560         result = 1;
 561       }
 562     if (v1.ieee.exponent != v2.ieee.exponent)
 563       {
 564         printf ("-0.0625 up: exponent differs: %4x vs %4x\n",
 565                 v1.ieee.exponent, v2.ieee.exponent);
 566         result = 1;
 567       }
 568     if (v1.ieee.negative != v2.ieee.negative)
 569       {
 570         printf ("-0.0625 up: negative differs: %d vs %d\n",
 571                 v1.ieee.negative, v2.ieee.negative);
 572         result = 1;
 573       }
 574
 575     v1.d = d = -0.0625;
 576     d = nextafter (d, -1.0);
 577     v2.d = d = nextafter (d, 0.0);
 578
 579     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 580       {
 581         printf ("-0.0625 down: mantissa0 differs: %8x vs %8x\n",
 582                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 583         result = 1;
 584       }
 585     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 586       {
 587         printf ("-0.0625 down: mantissa1 differs: %8x vs %8x\n",
 588                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 589         result = 1;
 590       }
 591     if (v1.ieee.exponent != v2.ieee.exponent)
 592       {
 593         printf ("-0.0625 down: exponent differs: %4x vs %4x\n",
 594                 v1.ieee.exponent, v2.ieee.exponent);
 595         result = 1;
 596       }
 597     if (v1.ieee.negative != v2.ieee.negative)
 598       {
 599         printf ("-0.0625 down: negative differs: %d vs %d\n",
 600                 v1.ieee.negative, v2.ieee.negative);
 601         result = 1;
 602       }
 603
 604     v1.d = d = 0.0;
 605     d = nextafter (d, 1.0);
 606     v2.d = nextafter (d, -1.0);
 607
 608     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 609       {
 610         printf ("0.0 up: mantissa0 differs: %8x vs %8x\n",
 611                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 612         result = 1;
 613       }
 614     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 615       {
 616         printf ("0.0 up: mantissa1 differs: %8x vs %8x\n",
 617                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 618         result = 1;
 619       }
 620     if (v1.ieee.exponent != v2.ieee.exponent)
 621       {
 622         printf ("0.0 up: exponent differs: %4x vs %4x\n",
 623                 v1.ieee.exponent, v2.ieee.exponent);
 624         result = 1;
 625       }
 626     if (0 != v2.ieee.negative)
 627       {
 628         printf ("0.0 up: negative differs: 0 vs %d\n",
 629                 v2.ieee.negative);
 630         result = 1;
 631       }
 632
 633     v1.d = d = 0.0;
 634     d = nextafter (d, -1.0);
 635     v2.d = nextafter (d, 1.0);
 636
 637     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 638       {
 639         printf ("0.0 down: mantissa0 differs: %8x vs %8x\n",
 640                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 641         result = 1;
 642       }
 643     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 644       {
 645         printf ("0.0 down: mantissa1 differs: %8x vs %8x\n",
 646                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 647         result = 1;
 648       }
 649     if (v1.ieee.exponent != v2.ieee.exponent)
 650       {
 651         printf ("0.0 down: exponent differs: %4x vs %4x\n",
 652                 v1.ieee.exponent, v2.ieee.exponent);
 653         result = 1;
 654       }
 655     if (1 != v2.ieee.negative)
 656       {
 657         printf ("0.0 down: negative differs: 1 vs %d\n",
 658                 v2.ieee.negative);
 659         result = 1;
 660       }
 661   }
 662
 663 #ifndef NO_LONG_DOUBLE
 664   {
 665     union ieee854_long_double v1;
 666     union ieee854_long_double v2;
 667     long double ld;
 668
 669     v1.d = ld = LDBL_MIN;
 670     if (fpclassify (ld) != FP_NORMAL)
 671       {
 672         printf ("fpclassify (LDBL_MIN) failed: %d\n", fpclassify (ld));
 673         result = 1;
 674       }
 675     ld = nextafterl (ld, LDBL_MIN / 2.0);
 676     if (fpclassify (ld) != FP_SUBNORMAL)
 677       {
 678         printf ("fpclassify (LDBL_MIN-epsilon) failed: %d (%La)\n",
 679                 fpclassify (ld), ld);
 680         result = 1;
 681       }
 682     v2.d = ld = nextafterl (ld, LDBL_MIN);
 683     if (fpclassify (ld) != FP_NORMAL)
 684       {
 685         printf ("fpclassify (LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
 686                 fpclassify (ld), ld);
 687         result = 1;
 688       }
 689
 690     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 691       {
 692         printf ("LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
 693                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 694         result = 1;
 695       }
 696     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 697       {
 698         printf ("LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
 699                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 700         result = 1;
 701       }
 702     if (v1.ieee.exponent != v2.ieee.exponent)
 703       {
 704         printf ("LDBL_MIN: exponent differs: %4x vs %4x\n",
 705                 v1.ieee.exponent, v2.ieee.exponent);
 706         result = 1;
 707       }
 708     if (v1.ieee.negative != v2.ieee.negative)
 709       {
 710         printf ("LDBL_MIN: negative differs: %d vs %d\n",
 711                 v1.ieee.negative, v2.ieee.negative);
 712         result = 1;
 713       }
 714
 715     v1.d = ld = -LDBL_MIN;
 716     if (fpclassify (ld) != FP_NORMAL)
 717       {
 718         printf ("fpclassify (-LDBL_MIN) failed: %d\n", fpclassify (ld));
 719         result = 1;
 720       }
 721     ld = nextafterl (ld, -LDBL_MIN / 2.0);
 722     if (fpclassify (ld) != FP_SUBNORMAL)
 723       {
 724         printf ("fpclassify (-LDBL_MIN-epsilon) failed: %d (%La)\n",
 725                 fpclassify (ld), ld);
 726         result = 1;
 727       }
 728     v2.d = ld = nextafterl (ld, -LDBL_MIN);
 729     if (fpclassify (ld) != FP_NORMAL)
 730       {
 731         printf ("fpclassify (-LDBL_MIN-epsilon+epsilon) failed: %d (%La)\n",
 732                 fpclassify (ld), ld);
 733         result = 1;
 734       }
 735
 736     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 737       {
 738         printf ("-LDBL_MIN: mantissa0 differs: %8x vs %8x\n",
 739                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 740         result = 1;
 741       }
 742     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 743       {
 744         printf ("-LDBL_MIN: mantissa1 differs: %8x vs %8x\n",
 745                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 746         result = 1;
 747       }
 748     if (v1.ieee.exponent != v2.ieee.exponent)
 749       {
 750         printf ("-LDBL_MIN: exponent differs: %4x vs %4x\n",
 751                 v1.ieee.exponent, v2.ieee.exponent);
 752         result = 1;
 753       }
 754     if (v1.ieee.negative != v2.ieee.negative)
 755       {
 756         printf ("-LDBL_MIN: negative differs: %d vs %d\n",
 757                 v1.ieee.negative, v2.ieee.negative);
 758         result = 1;
 759       }
 760
 761     ld = LDBL_MAX;
 762     if (fpclassify (ld) != FP_NORMAL)
 763       {
 764         printf ("fpclassify (LDBL_MAX) failed: %d\n", fpclassify (ld));
 765         result = 1;
 766       }
 767     ld = nextafterl (ld, INFINITY);
 768     if (fpclassify (ld) != FP_INFINITE)
 769       {
 770         printf ("fpclassify (LDBL_MAX+epsilon) failed: %d\n", fpclassify (ld));
 771         result = 1;
 772       }
 773
 774     ld = -LDBL_MAX;
 775     if (fpclassify (ld) != FP_NORMAL)
 776       {
 777         printf ("fpclassify (-LDBL_MAX) failed: %d\n", fpclassify (ld));
 778         result = 1;
 779       }
 780     ld = nextafterl (ld, -INFINITY);
 781     if (fpclassify (ld) != FP_INFINITE)
 782       {
 783         printf ("fpclassify (-LDBL_MAX-epsilon) failed: %d\n",
 784                 fpclassify (ld));
 785         result = 1;
 786       }
 787
 788     v1.d = ld = 0.0625;
 789     ld = nextafterl (ld, 0.0);
 790     v2.d = ld = nextafterl (ld, 1.0);
 791
 792     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 793       {
 794         printf ("0.0625L down: mantissa0 differs: %8x vs %8x\n",
 795                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 796         result = 1;
 797       }
 798     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 799       {
 800         printf ("0.0625L down: mantissa1 differs: %8x vs %8x\n",
 801                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 802         result = 1;
 803       }
 804     if (v1.ieee.exponent != v2.ieee.exponent)
 805       {
 806         printf ("0.0625L down: exponent differs: %4x vs %4x\n",
 807                 v1.ieee.exponent, v2.ieee.exponent);
 808         result = 1;
 809       }
 810     if (v1.ieee.negative != v2.ieee.negative)
 811       {
 812         printf ("0.0625L down: negative differs: %d vs %d\n",
 813                 v1.ieee.negative, v2.ieee.negative);
 814         result = 1;
 815       }
 816
 817     v1.d = ld = 0.0625;
 818     ld = nextafterl (ld, 1.0);
 819     v2.d = ld = nextafterl (ld, 0.0);
 820
 821     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 822       {
 823         printf ("0.0625L up: mantissa0 differs: %8x vs %8x\n",
 824                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 825         result = 1;
 826       }
 827     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 828       {
 829         printf ("0.0625L up: mantissa1 differs: %8x vs %8x\n",
 830                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 831         result = 1;
 832       }
 833     if (v1.ieee.exponent != v2.ieee.exponent)
 834       {
 835         printf ("0.0625L up: exponent differs: %4x vs %4x\n",
 836                 v1.ieee.exponent, v2.ieee.exponent);
 837         result = 1;
 838       }
 839     if (v1.ieee.negative != v2.ieee.negative)
 840       {
 841         printf ("0.0625L up: negative differs: %d vs %d\n",
 842                 v1.ieee.negative, v2.ieee.negative);
 843         result = 1;
 844       }
 845
 846     v1.d = ld = -0.0625;
 847     ld = nextafterl (ld, 0.0);
 848     v2.d = ld = nextafterl (ld, -1.0);
 849
 850     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 851       {
 852         printf ("-0.0625L up: mantissa0 differs: %8x vs %8x\n",
 853                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 854         result = 1;
 855       }
 856     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 857       {
 858         printf ("-0.0625L up: mantissa1 differs: %8x vs %8x\n",
 859                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 860         result = 1;
 861       }
 862     if (v1.ieee.exponent != v2.ieee.exponent)
 863       {
 864         printf ("-0.0625L up: exponent differs: %4x vs %4x\n",
 865                 v1.ieee.exponent, v2.ieee.exponent);
 866         result = 1;
 867       }
 868     if (v1.ieee.negative != v2.ieee.negative)
 869       {
 870         printf ("-0.0625L up: negative differs: %d vs %d\n",
 871                 v1.ieee.negative, v2.ieee.negative);
 872         result = 1;
 873       }
 874
 875     v1.d = ld = -0.0625;
 876     ld = nextafterl (ld, -1.0);
 877     v2.d = ld = nextafterl (ld, 0.0);
 878
 879     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 880       {
 881         printf ("-0.0625L down: mantissa0 differs: %8x vs %8x\n",
 882                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 883         result = 1;
 884       }
 885     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 886       {
 887         printf ("-0.0625L down: mantissa1 differs: %8x vs %8x\n",
 888                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 889         result = 1;
 890       }
 891     if (v1.ieee.exponent != v2.ieee.exponent)
 892       {
 893         printf ("-0.0625L down: exponent differs: %4x vs %4x\n",
 894                 v1.ieee.exponent, v2.ieee.exponent);
 895         result = 1;
 896       }
 897     if (v1.ieee.negative != v2.ieee.negative)
 898       {
 899         printf ("-0.0625L down: negative differs: %d vs %d\n",
 900                 v1.ieee.negative, v2.ieee.negative);
 901         result = 1;
 902       }
 903
 904     v1.d = ld = 0.0;
 905     ld = nextafterl (ld, 1.0);
 906     v2.d = nextafterl (ld, -1.0);
 907
 908     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 909       {
 910         printf ("0.0L up: mantissa0 differs: %8x vs %8x\n",
 911                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 912         result = 1;
 913       }
 914     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 915       {
 916         printf ("0.0L up: mantissa1 differs: %8x vs %8x\n",
 917                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 918         result = 1;
 919       }
 920     if (v1.ieee.exponent != v2.ieee.exponent)
 921       {
 922         printf ("0.0L up: exponent differs: %4x vs %4x\n",
 923                 v1.ieee.exponent, v2.ieee.exponent);
 924         result = 1;
 925       }
 926     if (0 != v2.ieee.negative)
 927       {
 928         printf ("0.0L up: negative differs: 0 vs %d\n",
 929                 v2.ieee.negative);
 930         result = 1;
 931       }
 932
 933     v1.d = ld = 0.0;
 934     ld = nextafterl (ld, -1.0);
 935     v2.d = nextafterl (ld, 1.0);
 936
 937     if (v1.ieee.mantissa0 != v2.ieee.mantissa0)
 938       {
 939         printf ("0.0L down: mantissa0 differs: %8x vs %8x\n",
 940                 v1.ieee.mantissa0, v2.ieee.mantissa0);
 941         result = 1;
 942       }
 943     if (v1.ieee.mantissa1 != v2.ieee.mantissa1)
 944       {
 945         printf ("0.0L down: mantissa1 differs: %8x vs %8x\n",
 946                 v1.ieee.mantissa1, v2.ieee.mantissa1);
 947         result = 1;
 948       }
 949     if (v1.ieee.exponent != v2.ieee.exponent)
 950       {
 951         printf ("0.0L down: exponent differs: %4x vs %4x\n",
 952                 v1.ieee.exponent, v2.ieee.exponent);
 953         result = 1;
 954       }
 955     if (1 != v2.ieee.negative)
 956       {
 957         printf ("0.0L down: negative differs: 1 vs %d\n",
 958                 v2.ieee.negative);
 959         result = 1;
 960       }
 961   }
 962 #endif
 963
 964   if (! isnormal (FLT_MIN))
 965     {
 966       puts ("isnormal (FLT_MIN) failed");
 967       result = 1;
 968     }
 969   if (! isnormal (DBL_MIN))
 970     {
 971       puts ("isnormal (DBL_MIN) failed");
 972       result = 1;
 973     }
 974 #ifndef NO_LONG_DOUBLE
 975   if (! isnormal (LDBL_MIN))
 976     {
 977       puts ("isnormal (LDBL_MIN) failed");
 978       result = 1;
 979     }
 980 #endif
 981
 982 #ifdef __i386__
 983   /* This is a test for the strange long doubles in x86 FPUs.  */
 984   {
 985     union
 986     {
 987       char b[10];
 988       long double d;
 989     } u =
 990       { .b = { 0, 0, 0, 0, 0, 0, 0, 0x80, 0, 0 } };
 991
 992     if (fpclassify (u.d) != FP_NORMAL)
 993       {
 994         printf ("fpclassify (0x00008000000000000000) failed: %d (%Lg)\n",
 995                 fpclassify (u.d), u.d);
 996         result = 1;
 997       }
 998   }
 999
1000   /* Special NaNs in x86 long double.  Test for scalbl.  */
1001   {
1002     union
1003     {
1004       char b[10];
1005       long double d;
1006     } u =
1007       { .b = { 0, 1, 0, 0, 0, 0, 0, 0xc0, 0xff, 0x7f } };
1008     long double r;
1009
1010     r = scalbl (u.d, 0.0);
1011     if (!isnan (r))
1012       {
1013         puts ("scalbl(NaN, 0) does not return NaN");
1014         result = 1;
1015       }
1016     else if (memcmp (&r, &u.d, sizeof (double)) != 0)
1017       {
1018         puts ("scalbl(NaN, 0) does not return the same NaN");
1019         result = 1;
1020       }
1021   }
1022 #endif
1023
1024 #ifndef NO_LONG_DOUBLE
1025   {
1026     long double r;
1027
1028     feclearexcept (FE_ALL_EXCEPT);
1029     r = scalbl (LDBL_MIN, 2147483647);
1030     if (! isinf (r))
1031       {
1032         puts ("scalbl (LDBL_MIN, 2147483647) does not return Inf");
1033         result = 1;
1034       }
1035     else if (signbit (r) != 0)
1036       {
1037         puts ("scalbl (LDBL_MIN, 2147483647) returns -Inf");
1038         result = 1;
1039       }
1040     else if (fetestexcept (FE_UNDERFLOW))
1041       {
1042         puts ("scalbl(NaN, 0) raises underflow exception");
1043         result = 1;
1044       }
1045
1046     feclearexcept (FE_ALL_EXCEPT);
1047     r = scalbl (LDBL_MAX, -2147483647);
1048     if (r != 0.0)
1049       {
1050         puts ("scalbl (LDBL_MAX, -2147483647) does not return 0");
1051         result = 1;
1052       }
1053     else if (signbit (r) != 0)
1054       {
1055         puts ("scalbl (LDBL_MAX, -2147483647) returns -Inf");
1056         result = 1;
1057       }
1058     else if (fetestexcept (FE_OVERFLOW))
1059       {
1060         puts ("scalbl(NaN, 0) raises overflow exception");
1061         result = 1;
1062       }
1063   }
1064 #endif
1065
1066   return result;
1067 }